ALBERT

Description: Ice nuclei in marine air: bioparticles or dust? Atmospheric Chemistry and Physics Discussions, 12, 4373-4416, 2012 Author(s): S. M. Burrows, C. Hoose, U. Pöschl, and M. G. Lawrence Ice nuclei can influence the properties of clouds and precipitation, but their sources and distribution in the atmosphere are still not well known. Particularly little attention has been paid to IN sources in marine environments, although anecdotal evidence suggests that IN populations in remote marine regions may be dominated by biological particles associated with sea spray. In this exploratory model study, we aim to bring attention to this long-neglected topic and identify promising target regions for future field campaigns. We assess the likely global distribution of marine biological ice nuclei using a combination of historical observations, satellite data and model output. By comparing simulated marine biological IN distributions and dust IN distributions, we predict strong regional differences in the importance of marine biological IN relative to dust IN. Our analysis suggests that marine biological IN are most likely to play a dominant role in determining IN concentrations over the Southern Ocean, so future field campaigns aimed at investigating marine biological IN should target that region. Climate-related changes in the abundance and emission of biological marine IN could affect marine cloud properties, thereby introducing previously unconsidered feedbacks that influence the hydrological cycle and the Earth's energy balance. Furthermore, marine biological IN may be an important aspect to consider in proposals for marine cloud brightening by artificial sea spray production.

Description: An observation-based approach to identify local natural dust events from routine aerosol ground monitoring Atmospheric Chemistry and Physics Discussions, 12, 4279-4310, 2012 Author(s): D. Q. Tong, M. Dan, T. Wang, and P. Lee Dust is a major component of atmospheric aerosols in many parts of the world. Although there exist many routine aerosol monitoring networks, it is often difficult to obtain dust records from these networks, because these monitors are either deployed far away from dust active regions (most likely collocated with dense population) or contaminated by anthropogenic sources and other natural sources, such as wildfires and vegetation detritus. Here we propose a new approach to identify local dust events relying solely on aerosol mass and composition from general-purpose aerosol measurements. Through analyzing the chemical and physical characteristics of aerosol observations during satellite-detected dust episodes, we select five indicators to be used to identify local dust records: (1) high PM 10 concentrations; (2) low PM 2.5 /PM 10 ratio; (3) higher concentrations and percentage of crustal elements; (4) lower percentage of anthropogenic pollutants; and (5) low enrichment factors of anthropogenic elements. After establishing these identification criteria, we conduct hierarchical cluster analysis for all validated aerosol measurement data over 68 IMPROVE sites in the Western United States. A total of 182 local dust events were identified over 30 of the 68 locations from 2000 to 2007. These locations are either close to the four US Deserts, namely the Great Basin Desert, the Mojave Desert, the Sonoran Desert, and the Chihuahuan Desert, or in the high wind power region (Colorado). During the eight-year study period, the total number of dust events displays an interesting four-year activity cycle (one in 2000–2003 and the other in 2004–2007). The years of 2003, 2002 and 2007 are the three most active dust periods, with 46, 31 and 24 recorded dust events, respectively, while the years of 2000, 2004 and 2005 are the calmest periods, all with single digit dust records. Among these deserts, the Chihuahua Desert (59 cases) and the Sonoran Desert (62 cases) are by far the most active source regions. In general, the Chihuahua Desert dominates dust activities in the first half of the eight-year period while the Sonoran Desert in the second half. The monthly frequency of dust events shows a peak from March to July and a second peak in autumn from September to November. The large quantity of dust events occurring in summertime also suggests the prevailing impact of windblown dust across the year. This seasonal variation is consistent with previous model simulations over the United States.

Description: Quantifying retrieval uncertainties in the CM-SAF cloud physical property algorithm with simulated SEVIRI observations Atmospheric Chemistry and Physics Discussions, 12, 4311-4340, 2012 Author(s): B. J. Jonkheid, R. A. Roebeling, and E. van Meijgaard The uncertainties in the cloud physical properties derived from satellite observations make it difficult to interpret model evaluation studies. In this paper, the uncertainties in the cloud water path (CWP) retrievals derived with the cloud physical properties retrieval algorithm (CPP) of the climate monitoring satellite application facility (CM-SAF) are investigated. To this end, a numerical simulator of MSG-SEVIRI observations was developed that calculates the reflectances at 0.64 and 1.63 μm for a wide range of cloud parameters, satellite viewing geometries and surface albedos. These reflectances are used as input to CPP, and the retrieved values of CWP are compared to the original input of the simulator. It is shown that the CWP retrievals are very sensitive to the assumptions made in the CPP code. The CWP retrieval errors are generally small for unbroken single-phase clouds with COT 〉10, with retrieval errors of ~3% for liquid water clouds to ~10% for ice clouds. When both liquid water and ice clouds are present in a pixel, the CWP retrieval errors increase dramatically; depending on the cloud, this can lead to uncertainties of 40–80%. CWP retrievals also become more uncertain when the cloud does not cover the entire pixel, leading to errors of ~50% for cloud fractions of 0.75 and even larger errors for smaller cloud fractions. Thus, the satellite retrieval of cloud physical properties of broken clouds and multi-phase clouds is complicated by inherent difficulties, and the proper interpretation of such retrievals requires extra care.

Description: Physically-based modeling of topographic effects on spatial evapotranspiration and soil moisture patterns through radiation and wind Hydrology and Earth System Sciences, 16, 357-373, 2012 Author(s): M. Liu, A. Bárdossy, J. Li, and Y. Jiang In this paper, simulations with the Soil Water Atmosphere Plant (SWAP) model are performed to quantify the spatial variability of both potential and actual evapotranspiration (ET), and soil moisture content (SMC) caused by topography-induced spatial wind and radiation differences. To obtain the spatially distributed ET/SMC patterns, the field scale SWAP model is applied in a distributed way for both pointwise and catchment wide simulations. An adapted radiation model from r.sun and the physically-based meso-scale wind model METRAS PC are applied to obtain the spatial radiation and wind patterns respectively, which show significant spatial variation and correlation with aspect and elevation respectively. Such topographic dependences and spatial variations further propagate to ET/SMC. A strong spatial, seasonal-dependent, scale-relevant intra-catchment variability in daily/annual ET and less variability in SMC can be observed from the numerical experiments. The study concludes that topography has a significant effect on ET/SMC in the humid region where ET is a energy limited rather than water availability limited process. It affects the spatial runoff generation through spatial radiation and wind, therefore should be applied to inform hydrological model development. In addition, the methodology used in the study can serve as a general method for physically-based ET estimation for data sparse regions.

Description: Enhanced cold-season warming in semi-arid regions Atmospheric Chemistry and Physics Discussions, 12, 4627-4653, 2012 Author(s): J. Huang, X. Guan, and F. Ji This study examined surface air temperature trends over global land from 1901–2009. It is found that the warming trend was particularly enhanced, in the boreal cold season (November to March) over semi-arid regions (with precipitation of 200–600 mm yr −1 ), showing a temperature increase of 1.53 °C as compared to the global annual mean temperature increase of 1.13 °C over land. In mid-latitude semi-arid areas of Europe, Asia, and North America, temperatures in the cold season increased by 1.41, 2.42, and 1.5 °C, respectively. The semi-arid regions contribute 44.46% to global annual-mean land-surface temperature trend. The mid-latitude semi-arid regions in the Northern Hemisphere accounting contribute by 27.0% of the total, with the mid-latitude semi-arid areas in Europe, Asia, and North America accounting for 6.29%, 13.81%, and 6.85%, respectively. Such enhanced semi-arid warming (ESAW) may cause these regions to become drier and warmer.

Description: Humidity-dependent phase state of SOA particles from biogenic and anthropogenic precursors Atmospheric Chemistry and Physics Discussions, 12, 4447-4476, 2012 Author(s): E. Saukko, A. T. Lambe, P. Massoli, Koop, T., J. P. Wright, D. R. Croasdale, D. A. Pedernera, T. B. Onasch, A. Laaksonen, P. Davidovits, D. R. Worsnop, and A. Virtanen The physical phase state (solid, semi-solid, or liquid) of secondary organic aerosol (SOA) particles has important implications for a number of atmospheric processes. We report the phase state of SOA particles spanning a wide range of oxygen to carbon ratios (O/C), used here as a surrogate for SOA oxidation level, produced in a flow tube reactor by photo-oxidation of various atmospherically relevant surrogate anthropogenic and biogenic volatile organic compounds (VOCs). The phase state of laboratory-generated SOA was determined by the particle bounce behavior after inertial impaction on a polished steel substrate. The measured bounce fraction was evaluated as a function of relative humidity and SOA oxidation level (O/C) measured by an Aerodyne high resolution time of flight aerosol mass spectrometer (HR-ToF AMS). The main findings of the study are: (1) Biogenic and anthropogenic SOA particles are found to be solid or semi-solid until a relative humidity of at least 50 % RH at impaction is reached. (2) Long-chain alkanes produce liquid SOA particles when generated at low oxidation level of O/C

Description: Carbonaceous aerosol AAE inferred from in-situ aerosol measurements at the Gosan ABC super site, and the implications for brown carbon aerosol Atmospheric Chemistry and Physics Discussions, 12, 4507-4539, 2012 Author(s): C. E. Chung, S.-W. Kim, M. Lee, S.-C. Yoon, and S. Lee Carbon mass of aerosols and its division between EC and OC sources were continuously measured at hourly intervals from October 2009 to June 2010. During this 9-month period, we also measured the aerosol absorption coefficient at 7 wavelengths and obtained PM mass density data at 1-h resolution. The measurement was made at the Gosan ABC super site, which is an ideal location for monitoring long-range transported aerosols from China. We remove the absorption data with significant dust influence using the mass ratio of PM 10 to PM 2.5 . The remaining data shows an Absorption Ångström Exponent (AAE) of about 1.27, which we suggest represent the average carbonaceous aerosol (CA) AAE at Gosan. CA AAE is highest in winter, in which the monthly value is near 1.4. We find a positive correlation between the mass ratio of OC to EC and CA AAE, and successfully increase the correlation by filtering out data associated with weak absorption signal. After the filtering, absorption coefficient is regressed on OC and EC mass densities. Black carbon (BC) and organic aerosol (OA) absorption cross sections per carbon mass are found to be 5.1 (4.2–6.0) and 1.4 (1.1–1.8) m 2 g −1 at 520 nm respectively. From the estimated BC & OA MAC, we find that OA contributes about 45% to CA absorption at 520 nm. OA AAE is found to be 1.7 (1.4–2.1). Compared with a previous estimate of OA MAC and AAE, our estimates at Gosan strongly suggest that the strongly-absorbing so-called brown carbon spheres are either unrelated to biomass burning or absent near the emission source.

Description: The importance of parameter resampling for soil moisture data assimilation into hydrologic models using the particle filter Hydrology and Earth System Sciences, 16, 375-390, 2012 Author(s): D. A. Plaza, R. De Keyser, G. J. M. De Lannoy, L. Giustarini, P. Matgen, and V. R. N. Pauwels The Sequential Importance Sampling with Resampling (SISR) particle filter and the SISR with parameter resampling particle filter (SISR-PR) are evaluated for their performance in soil moisture assimilation and the consequent effect on baseflow generation. With respect to the resulting soil moisture time series, both filters perform appropriately. However, the SISR filter has a negative effect on the baseflow due to inconsistency between the parameter values and the states after the assimilation. In order to overcome this inconsistency, parameter resampling is applied along with the SISR filter, to obtain consistent parameter values with the analyzed soil moisture state. Extreme parameter replication, which could lead to a particle collapse, is avoided by the perturbation of the parameters with white noise. Both the modeled soil moisture and baseflow are improved if the complementary parameter resampling is applied. The SISR filter with parameter resampling offers an efficient way to deal with biased observations. The robustness of the methodology is evaluated for 3 model parameter sets and 3 assimilation frequencies. Overall, the results in this paper indicate that the particle filter is a promising tool for hydrologic modeling purposes, but that an additional parameter resampling may be necessary to consistently update all state variables and fluxes within the model.

Description: Organic molecular markers and signature from wood combustion particles in winter ambient aerosols: aerosol mass spectrometer (AMS) and high time-resolved GC-MS measurements in Augsburg, Germany Atmospheric Chemistry and Physics Discussions, 12, 4831-4866, 2012 Author(s): M. Elsasser, M. Crippa, J. Orasche, P. F. DeCarlo, M. Oster, M. Pitz, T. L. Gustafson, J. B. C. Pettersson, J. Schnelle-Kreis, A. S. H. Prévôt, and R. Zimmermann The impact of wood combustion on ambient aerosols was investigated in Augsburg, Germany during a winter measurement campaign of a six-week period. Special attention was paid to the high time resolution observations of wood combustion with different mass spectrometric methods. Here we present and compare the results from an Aerodyne aerosol mass spectrometer (AMS) and gas chromatographic – mass spectrometric (GC-MS) analysed PM 1 filters on an hourly basis. This includes source apportionment of the AMS derived organic mass by using positive matrix factorisation (PMF) and analysis of levoglucosan as wood combustion marker, respectively. In the measurement period nitrate and organics are the main contributors to the defined submicron particle mass with 28% and 35%, respectively. To the latter wood combustion organic aerosol (WCOA) contributes 23% on average and 27% in the evening and night-time. Conclusively, wood combustion has a strong influence on the organics and overall aerosol composition. Levoglucosan accounts for 14% of WCOA mass with a higher percentage in comparison to other studies. The ratio between the mass of levoglucosan and organic carbon amounts to 0.06. This study is unique in the one-hour time resolution comparison between the wood combustion results of the AMS and the GC-MS analysed filter method at a PM 1 particle size range. This comparison of the concentration courses of the PMF WCOA factor, levoglucosan estimated by the AMS data and the levoglucosan measured by GC-MS is highly correlated, and a detailed discussion on the contributors to the wood combustion marker ion at mass-to-charge ratio 60 will be given. This offers a suitable application possibility for the description of the wood combustion course by the WCOA factor and the levoglucosan concentration estimated by AMS data. However, quantitative description of the levoglucosan concentration estimated by the AMS data is difficult due to the offset of latter compared to measured levoglucosan by the GC-MS.

Description: Mean winds, temperatures and the 16- and 5-day planetary waves in the mesosphere and lower thermosphere over Bear Lake Observatory (42° N, 111° W) Atmospheric Chemistry and Physics, 12, 1571-1585, 2012 Author(s): K. A. Day, M. J. Taylor, and N. J. Mitchell Atmospheric temperatures and winds in the mesosphere and lower thermosphere have been measured simultaneously using the Aura satellite and a meteor radar at Bear Lake Observatory (42° N, 111° W), respectively. The data presented in this study is from the interval March 2008 to July 2011. The mean winds observed in the summer-time over Bear Lake Observatory show the meridional winds to be equatorward at meteor heights during April−August and to reach monthly-mean velocities of −12 m s −1 . The mean winds are closely related to temperatures in this region of the atmosphere and in the summer the coldest mesospheric temperatures occur about the same time as the strongest equatorward meridional winds. The zonal winds are eastward through most of the year and in the summer strong eastward zonal wind shears of up to ~4.5 m s −1 km −1 are present. However, westward winds are observed at the upper heights in winter and sometimes during the equinoxes. Considerable inter-annual variability is observed in the mean winds and temperatures. Comparisons of the observed winds with URAP and HWM-07 reveal some large differences. Our radar zonal wind observations are generally more eastward than predicted by the URAP model zonal winds. Considering the radar meridional winds, in comparison to HWM-07 our observations reveal equatorward flow at all meteor heights in the summer whereas HWM-07 suggests that only weakly equatorward, or even poleward flows occur at the lower heights. However, the zonal winds observed by the radar and modelled by HWM-07 are generally similar in structure and strength. Signatures of the 16- and 5-day planetary waves are clearly evident in both the radar-wind data and Aura-temperature data. Short-lived wave events can reach large amplitudes of up to ~15 m s −1 and 8 K and 20 m s −1 and 10 K for the 16- and 5-day waves, respectively. A clear seasonal and short-term variability are observed in the 16- and 5-day planetary wave amplitudes. The 16-day wave reaches largest amplitude in winter and is also present in summer, but with smaller amplitudes. The 5-day wave reaches largest amplitude in winter and in late summer. An inter-annual variability in the amplitude of the planetary waves is evident in the four years of observations. Some 41 episodes of large-amplitude wave occurrence are identified. Temperature and wind amplitudes for these episodes, A T and A W , that passed the Student T-test were found to be related by, A T = 0.34 A W and A T = 0.62 A W for the 16- and 5-day wave, respectively.

Description: Remote sensing-based estimation of gross primary production in a subalpine grassland Biogeosciences Discussions, 9, 1711-1758, 2012 Author(s): M. Rossini, S. Cogliati, M. Meroni, M. Migliavacca, M. Galvagno, L. Busetto, E. Cremonese, T. Julitta, C. Siniscalco, U. Morra di Cella, and R. Colombo This study investigates the performances in a terrestrial ecosystem of gross primary production (GPP) estimation of a suite of spectral vegetation indexes (VIs) that can be computed from currently orbiting platforms. Vegetation indexes were computed from near-surface field spectroscopy measurements collected using an automatic system designed for high temporal frequency acquisition of spectral measurements in the visible near-infrared region. Spectral observations were collected for two consecutive years in Italy in a subalpine grassland equipped with an Eddy Covariance (EC) flux tower which provides continuous measurements of net ecosystem carbon dioxide (CO 2 ) exchange (NEE) and the derived GPP. Different VIs were calculated based on ESA-MERIS and NASA-MODIS spectral bands and correlated with biophysical (Leaf Area Index, LAI; fraction of photosynthetically active radiation intercepted by green vegetation, f IPAR g ), biochemical (chlorophyll concentration) and ecophysiological (green light-use efficiency, LUE g ) canopy variables. In this study, the normalized difference vegetation index (NDVI) showed better correlations with LAI and f PAR g ( r = 0.90 and 0.95, respectively), the MERIS terrestrial chlorophyll index (MTCI) with leaf chlorophyll content ( r = 0.91) and the Photochemical Reflectance Index (PRI 551 ), computed as ( R 531 − R 551 )/( R 531 + R 551 ) with LUE g ( r = 0.64). Subsequently, these VIs were used to estimate GPP using different modelling solutions based on the light-use efficiency model describing the GPP as driven by the photosynthetically active radiation absorbed by green vegetation (APAR g ) and by the efficiency (ε) with which plants use the absorbed radiation to fix carbon via photosynthesis. Results show that GPP can be successfully modelled with a combination of VIs and meteorological data or VIs only. Vegetation indexes designed to be more sensitive to chlorophyll content explained most of the variability in GPP in the ecosystem investigated, characterized by a strong seasonal dynamic of GPP. Accuracy in GPP estimation slightly improves when taking into account high frequency modulations of GPP driven by incident PAR or modelling LUE g with the PRI in model formulation. Similar results were obtained for both measured daily VIs and VIs obtained as 16-day composite time series and then downscaled from the compositing period to daily scale (resampled data). However, the use of resampled data rather than measured daily input data decreases the accuracy of the total GPP estimation on an annual basis.

Description: Soil organic carbon storage changes in coastal wetlands of the modern Yellow River Delta from 2000 to 2009 Biogeosciences Discussions, 9, 1759-1779, 2012 Author(s): J. Yu, Y. Wang, Y. Li, H. Dong, D. Zhou, G. Han, H. Wu, G. Wang, P. Mao, and Y. Gao Soil carbon sequestration plays an essential role in mitigating CO 2 increases and the subsequently global greenhouse effect. The storages and dynamics of soil organic carbon (SOC) of 0–30 cm soil depth in different landscape types including beaches, reservoir and pond, reed wetland, forest wetland, bush wetland, farmland, building land, bare land (severe saline land) and salt field in the modern Yellow River Delta (YRD), were studied based on the data of the regional survey and laboratory analysis. The landscape types were classified by the interpretation of remote sensing images of 2000 and 2009, which was calibrated by field survey results. The results revealed an increase of 10.59 km 2 in the modem YRD area from 2000 to 2009. The SOC density varied ranging from 0.73 kg m −2 to 21.60 kg m −2 at depth of 30 cm. There were ~3.97 × 10 6 t and 3.98 × 10 6 t SOC stored in the YRD in 2000 and 2009, respectively. The SOC storages changed greatly in beaches, bush wetland, farm land and salt field which were affected dominantly by anthropogenic activities. The area of the YRD increased greatly within 10 yr, however, the small increase of SOC storage in the region was observed due to landscape changes, indicating that the modern YRD was a potential carbon sink and anthropogenic activity was a key factor for SOC change.

Description: Evolution of marine storminess in the Belgian part of the North Sea Natural Hazards and Earth System Science, 12, 305-312, 2012 Author(s): D. Van den Eynde, R. De Sutter, and P. Haerens Severe storms have affected European coast lines in the past but knowledge on changes in storminess for the last decades is still sparse. Climate change is assumed to be a main driving factor with the potential to induce changes on the intensity, duration and frequency of powerful marine storms, including a long-term influence on peak wind speeds, surges and waves. It is, therefore, important to investigate whether in the last decades changes in the magnitude of storms, their duration and frequency could be observed. Understanding trends in storminess in the last decades will help to better prepare coastal managers for future events, taking into account potential changes on storm occurrence and magnitude to improve planning of mitigation and adaptation strategies. The purpose of this study was to focus on the evolution of extreme wind conditions, wave height and storm surge levels in the North Sea Region, especially in the Belgian part of the North Sea (BPNS). Based on the analysis performed it is concluded that no clear trend can be observed for the occurrence of significant increasing extreme wind speeds over the BPNS. Furthermore, one can conclude that not enough scientific evidence is available to support scenarios with increased wave height or storminess.

Description: Inferring snow pack ripening and melt out from distributed ground surface temperature measurements The Cryosphere Discussions, 6, 563-591, 2012 Author(s): M.-O. Schmid, S. Gubler, J. Fiddes, and S. Gruber The seasonal snow cover and its melting are heterogeneous both in space and time. Describing and modelling this variability are important because it affects divers phenomena such as runoff, ground temperatures or slope movements. This study investigates the derivation of melting characteristics based on spatial clusters of temperature measurements. Results are based on data from Switzerland where ground surface temperatures were measured with miniature loggers (iButtons) at 40 locations, referred to as footprints. At each footprint, ten iButtons have been distributed randomly few cm below the ground surface over an area of 10 m × 10 m. Footprints span elevations of 2100–3300 m a.s.l. and slope angles of 0–55°, as well as diverse slope expositions and types of surface cover and ground material. Based on two years of temperature data, the basal ripening date and the melt-out date are determined for each iButton, aggregated to the footprint level and further analysed. The date of melt out could be derived for nearly all iButtons, the ripening date could be extracted for only approximately half of them because it requires ground freezing below the snow pack. The variability within a footprint is often considerable and one to three weeks difference between melting or ripening of the points in one footprint is not uncommon. The correlation of mean annual ground surface temperatures, ripening date and melt-out date is moderate, making them useful intuitive complementary measured for model evaluation.

Description: An isotope view on ionising radiation as a source of sulphuric acid Atmospheric Chemistry and Physics Discussions, 12, 5039-5064, 2012 Author(s): M. B. Enghoff, N. Bork, S. Hattori, C. Meusinger, M. Nakagawa, J. O. P. Pedersen, S. Danielache, Y. Ueno, M. S. Johnson, N. Yoshida, and H. Svensmark Sulphuric acid is an important factor in aerosol nucleation and growth. It has been shown that ions enhance the formation of sulphuric acid aerosols, but the exact mechanism has remained undetermined. Furthermore some studies have found a deficiency in the sulphuric acid budget, suggesting a missing source. In this study the production of sulphuric acid from SO 2 through a number of different pathways is investigated. The production methods are standard gas phase oxidation by OH radicals produced by ozone photolysis with UV light, liquid phase oxidation by ozone, and gas phase oxidation initiated by gamma rays. The distributions of stable sulphur isotopes in the products and substrate were measured using isotope ratio mass spectrometry. All methods produced sulphate enriched in 34 S and we find a δ 34 S value of 8.7 ± 0.4‰ (1 standard deviation) for the UV-initiated OH reaction. Only UV light (Hg emission at 253.65 nm) produced a clear non-mass-dependent excess of 33 S. The pattern of isotopic enrichment produced by gamma rays is similar, but not equal, to that produced by aqueous oxidation of SO 2 by ozone. This, combined with the relative yields of the experiments, suggests a mechanism in which ionising radiation may lead to hydrated ion clusters that serve as nanoreactors for S(IV) to S(VI) conversion.

Description: First intercalibration of column-averaged methane from the Total Carbon Column Observing Network and the Network for the Detection of Atmospheric Composition Change Atmospheric Measurement Techniques Discussions, 5, 1355-1379, 2012 Author(s): F. Forster, R. Sussmann, M. Rettinger, N. M. Deutscher, D. W. T. Griffith, N. Jones, and P. K. Patra We present the intercalibration of dry-air column-averaged mole fractions of methane (XCH 4 ) retrieved from solar FTIR measurements of the Network for the Detection of Atmospheric Composition Change (NDACC) in the mid-infrared (MIR) versus near-infrared (NIR) soundings from the Total Carbon Column Observing Network (TCCON). The study uses multi-annual quasi-coincident MIR and NIR measurements from the stations Garmisch, Germany (47.48° N, 11.06° E, 743 m a.s.l.) and Wollongong, Australia (34.41° S, 150.88° E, 30 m a.s.l.). Direct comparison of the retrieved MIR and NIR time series shows a phase shift in XCH 4 seasonality, i.e. a significant time-dependent bias leading to a standard deviation (stdv) of the difference time series (NIR-MIR) of 8.4 ppb. After eliminating differences in a prioris by using ACTM-simulated profiles as a common prior, the seasonalities of the (corrected) MIR and NIR time series agree within the noise (stdv = 5.2 ppb for the difference time series). The difference time series (NIR-MIR) do not show a significant trend. Therefore it is possible to use a simple scaling factor for the intercalibration without a time-dependent linear or seasonal component. Using the Garmisch and Wollongong data together, we obtain an overall calibration factor MIR/NIR = 0.9926(18). The individual calibration factors per station are 0.9940(14) for Garmisch and 0.9893(40) for Wollongong. They agree within their error bars with the overall calibration factor which can therefore be used for both stations. Our results suggest that after applying the proposed intercalibration concept to all stations performing both NIR and MIR measurements, it should be possible to obtain one refined overall intercalibration factor for the two networks. This would allow to set up a harmonized NDACC and TCCON XCH 4 data set which can be exploited for joint trend studies, satellite validation, or the inverse modeling of sources and sinks.

Description: Revised identification of tropical oceanic cumulus congestus as viewed by CloudSat Atmospheric Chemistry and Physics, 12, 1587-1595, 2012 Author(s): S. P. F. Casey, E. J. Fetzer, and B. H. Kahn Congestus cloud convective features are examined in one year of tropical oceanic cloud observations from the CloudSat/CALIPSO instruments. Two types of convective clouds (cumulus and deep convective, based on classification profiles from radar), and associated differences in radar reflectivity and radar/lidar cloud-top height are considered. Congestus convective features are defined as contiguous convective clouds with heights between 3 and 9 km. Three criteria were used in previous studies to identify congestus: (1) CloudSat and CALIPSO cloud-top heights less than 1 km apart; (2) CloudSat 0 dBZ echo-top height less than 1 km from CloudSat cloud-top height, and (3) CloudSat 10 dBZ echo-top height less than 2 km from CloudSat cloud-top height. A majority of congestus convective features satisfy the second and third requirements. However, over 40% of convective features identified had no associated CALIPSO cloud-top height, predominantly due to the extinguishment of the lidar beam above the CloudSat-reported convective cloud. For the remaining cells, approximately 56% of these satisfy all three requirements; when considering the lidar beam-extinction issue, only 31% of congestus convective features are identified using these criteria. This implies that while previous methods used to identify congestus clouds may be accurate in finding vigorous convection (such as transient congestus rising toward the tropopause), these criteria may miss almost 70% of the total observed congestus convective features, suggesting a more general approach should be used to describe congestus and its surrounding environment.

Description: Wavelet-based spatial comparison technique for analysing and evaluating two-dimensional geophysical model fields Geoscientific Model Development, 5, 223-230, 2012 Author(s): S. Saux Picart, M. Butenschön, and J. D. Shutler Complex numerical models of the Earth's environment, based around 3-D or 4-D time and space domains are routinely used for applications including climate predictions, weather forecasts, fishery management and environmental impact assessments. Quantitatively assessing the ability of these models to accurately reproduce geographical patterns at a range of spatial and temporal scales has always been a difficult problem to address. However, this is crucial if we are to rely on these models for decision making. Satellite data are potentially the only observational dataset able to cover the large spatial domains analysed by many types of geophysical models. Consequently optical wavelength satellite data is beginning to be used to evaluate model hindcast fields of terrestrial and marine environments. However, these satellite data invariably contain regions of occluded or missing data due to clouds, further complicating or impacting on any comparisons with the model. This work builds on a published methodology, that evaluates precipitation forecast using radar observations based on predefined absolute thresholds. It allows model skill to be evaluated at a range of spatial scales and rain intensities. Here we extend the original method to allow its generic application to a range of continuous and discontinuous geophysical data fields, and therefore allowing its use with optical satellite data. This is achieved through two major improvements to the original method: (i) all thresholds are determined based on the statistical distribution of the input data, so no a priori knowledge about the model fields being analysed is required and (ii) occluded data can be analysed without impacting on the metric results. The method can be used to assess a model's ability to simulate geographical patterns over a range of spatial scales. We illustrate how the method provides a compact and concise way of visualising the degree of agreement between spatial features in two datasets. The application of the new method, its handling of bias and occlusion and the advantages of the novel method are demonstrated through the analysis of model fields from a marine ecosystem model.

Description: Moving beyond traditional model calibration or how to better identify realistic model parameters: sub-period calibration Hydrology and Earth System Sciences Discussions, 9, 1885-1918, 2012 Author(s): S. Gharari, M. Hrachowitz, F. Fenicia, and H. H. G. Savenije Conceptual hydrological models often rely on calibration for the identification of their parameters. As these models are typically designed to reflect real catchment processes, a key objective of an appropriate calibration strategy is the determination of parameter sets that reflect a "realistic" model behavior. Previous studies have shown that parameter estimates for different calibration periods can be significantly different. This questions model transposability in time, which is one of the key conditions for the set-up of a "realistic" model. This paper presents a new approach that selects parameter sets that provide a consistent model performance in time. The approach consists of confronting model performance in different periods, and selecting parameter sets that are as close as possible to the optimum of each individual sub-period. While aiding model calibration, the approach is also useful as a diagnostic tool, illustrating tradeoffs in the identification of time consistent parameter sets. The approach is demonstrated in a case study where we illustrate the multi-objective calibration of the HyMod hydrological model to a Luxembourgish catchment.

Description: Technical Note: Analytical sensitivity analysis of a two parameter recursive digital baseflow separation filter Hydrology and Earth System Sciences, 16, 451-455, 2012 Author(s): K. Eckhardt A sensitivity analysis for a well-established baseflow separation technique, a two parameter recursive digital filter, is presented. The sensitivity of the calculated baseflow index to errors or uncertainties of the two filter parameters and of the initial baseflow value is analytically ascertained. It is found that the influence of the initial baseflow value is negligible for long time series. The propagation of errors or uncertainties of the two filter parameters into the baseflow index is expressed by a dimensionless sensitivity index, the ratio between the relative error of the baseflow index and the relative error of the respective parameter. Representative index values are derived by application of the resulting equations to 65 North American catchments. In the mean the parameter a , the recession constant, has a stronger influence on the calculated baseflow index than the second filter parameter BFI max . This is favourable in that a can be determined by a recession analysis and therefore should be less uncertain. Whether this finding also applies for a specific catchment can easily be checked by means of the derived equations.

Description: The 2007 flood in the Sahel: causes, characteristics and its presentation in the media and FEWS NET Natural Hazards and Earth System Science, 12, 313-325, 2012 Author(s): C. Samimi, A. H. Fink, and H. Paeth During the rainy season in 2007, reports about exceptional rains and floodings in the Sahel were published in the media, especially in August and September. Institutions and organizations like the World Food Programme (WFP) and FEWS NET put the events on the agenda and released alerts and requested help. The partly controversial picture was that most of the Sahel faced a crisis caused by widespread floodings. Our study shows that the rainy season in 2007 was exceptional with regard to rainfall amount and return periods. In many areas the event had a return period between 1 and 50 yr with high spatial heterogeneity, with the exception of the Upper Volta basin, which yielded return periods of up to 1200 yr. Despite the strong rainfall, the interpretation of satellite images show that the floods were mainly confined to lakes and river beds. However, the study also proves the difficulties in assessing the meteorological processes and the demarcation of flooded areas in satellite images without ground truthing. These facts and the somewhat vague and controversial reports in the media and FEWS NET demonstrate that it is crucial to thoroughly analyze such events at a regional and local scale involving the local population.

Description: Factors selection in landslide susceptibility modelling on large scale following the gis matrix method: application to the river Beiro basin (Spain) Natural Hazards and Earth System Science, 12, 327-340, 2012 Author(s): D. Costanzo, E. Rotigliano, C. Irigaray, J. D. Jiménez-Perálvarez, and J. Chacón A procedure to select the controlling factors connected to the slope instability has been defined. It allowed us to assess the landslide susceptibility in the Rio Beiro basin (about 10 km 2 ) over the northeastern area of the city of Granada (Spain). Field and remote (Google EarthTM) recognition techniques allowed us to generate a landslide inventory consisting in 127 phenomena. To discriminate between stable and unstable conditions, a diagnostic area had been chosen as the one limited to the crown and the toe of the scarp of the landslide. 15 controlling or determining factors have been defined considering topographic, geologic, geomorphologic and pedologic available data. Univariate tests, using both association coefficients and validation results of single-variable susceptibility models, allowed us to select the best predictors, which were combined for the unique conditions analysis. For each of the five recognised landslide typologies, susceptibility maps for the best models were prepared. In order to verify both the goodness of fit and the prediction skill of the susceptibility models, two different validation procedures were applied and compared. Both procedures are based on a random partition of the landslide archive for producing a test and a training subset. The first method is based on the analysis of the shape of the success and prediction rate curves, which are quantitatively analysed exploiting two morphometric indexes. The second method is based on the analysis of the degree of fit, by considering the relative error between the intersected target landslides by each of the different susceptibility classes in which the study area was partitioned. Both the validation procedures confirmed a very good predictive performance of the susceptibility models and of the actual procedure followed to select the controlling factors.

Description: Spatial patterns of linear and nonparametric long-term trends in Baltic sea-level variability Nonlinear Processes in Geophysics, 19, 95-111, 2012 Author(s): R. V. Donner, R. Ehrcke, S. M. Barbosa, J. Wagner, J. F. Donges, and J. Kurths The study of long-term trends in tide gauge data is important for understanding the present and future risk of changes in sea-level variability for coastal zones, particularly with respect to the ongoing debate on climate change impacts. Traditionally, most corresponding analyses have exclusively focused on trends in mean sea-level. However, such studies are not able to provide sufficient information about changes in the full probability distribution (especially in the more extreme quantiles). As an alternative, in this paper we apply quantile regression (QR) for studying changes in arbitrary quantiles of sea-level variability. For this purpose, we chose two different QR approaches and discuss the advantages and disadvantages of different settings. In particular, traditional linear QR poses very restrictive assumptions that are often not met in reality. For monthly data from 47 tide gauges from along the Baltic Sea coast, the spatial patterns of quantile trends obtained in linear and nonparametric (spline-based) frameworks display marked differences, which need to be understood in order to fully assess the impact of future changes in sea-level variability on coastal areas. In general, QR demonstrates that the general variability of Baltic sea-level has increased over the last decades. Linear quantile trends estimated for sliding windows in time reveal a wide-spread acceleration of trends in the median, but only localised changes in the rates of changes in the lower and upper quantiles.

Description: Impact of spatial resolution on the modelling of the Greenland ice sheet surface mass balance between 1990–2010, using the regional climate model MAR The Cryosphere Discussions, 6, 635-672, 2012 Author(s): B. Franco, X. Fettweis, C. Lang, and M. Erpicum With the aim to force an ice dynamical model, the Greenland ice sheet (GrIS) surface mass balance (SMB) was modelled at different spatial resolutions (15–50 km) for the period 1990–2010, using the regional climate model MAR (Modèle Atmosphérique Régional) forced by the ERA-INTERIM reanalysis. This comparison revealed that (i) the inter-annual variability of the SMB components is consistent within the different spatial resolutions investigated, (ii) the MAR model simulates heavier precipitation on average over the GrIS with diminishing spatial resolution, and (iii) the SMB components (except precipitation) can be derived from a simulation at lower resolution with an ''intelligent'' interpolation. This interpolation can also be used to approximate the SMB components over another topography/ice sheet mask of the GrIS. These results are important for the forcing of an ice dynamical model, needed to enable future projections of the GrIS contribution to sea level rise over the coming centuries.

Description: Large surface meltwater discharge from the Kangerlussuaq sector of the Greenland ice sheet during the record-warm year 2010 explained by detailed energy balance observations The Cryosphere, 6, 199-209, 2012 Author(s): D. van As, A. L. Hubbard, B. Hasholt, A. B. Mikkelsen, M. R. van den Broeke, and R. S. Fausto This study uses data from six on-ice weather stations, calibrated MODIS-derived albedo and proglacial river gauging measurements to drive and validate an energy balance model. We aim to quantify the record-setting positive temperature anomaly in 2010 and its effect on mass balance and runoff from the Kangerlussuaq sector of the Greenland ice sheet. In 2010, the average temperature was 4.9 °C (2.7 standard deviations) above the 1974–2010 average in Kangerlussuaq. High temperatures were also observed over the ice sheet, with the magnitude of the positive anomaly increasing with altitude, particularly in August. Simultaneously, surface albedo was anomalously low in 2010, predominantly in the upper ablation zone. The low albedo was caused by high ablation, which in turn profited from high temperatures and low winter snowfall. Surface energy balance calculations show that the largest melt excess (∼170%) occurred in the upper ablation zone (above 1000 m), where higher temperatures and lower albedo contributed equally to the melt anomaly. At lower elevations the melt excess can be attributed to high atmospheric temperatures alone. In total, we calculate that 6.6 ± 1.0 km 3 of surface meltwater ran off the ice sheet in the Kangerlussuaq catchment in 2010, exceeding the reference year 2009 (based on atmospheric temperature measurements) by ∼150%. During future warm episodes we can expect a melt response of at least the same magnitude, unless a larger wintertime snow accumulation delays and moderates the melt-albedo feedback. Due to the hypsometry of the ice sheet, yielding an increasing surface area with elevation, meltwater runoff will be further amplified by increases in melt forcings such as atmospheric heat.

Description: Characterization and airborne deployment of a new counterflow virtual impactor inlet Atmospheric Measurement Techniques Discussions, 5, 1515-1541, 2012 Author(s): T. Shingler, S. Dey, A. Sorooshian, F. J. Brechtel, Z. Wang, A. Metcalf, M. Coggon, J. Mülmenstädt, L. M. Russell, H. H. Jonsson, and J. H. Seinfeld A new counterflow virtual impactor (CVI) inlet is introduced with details of its design, laboratory characterization tests, and deployment on an aircraft during the 2011 Eastern Pacific Emitted Aerosol Cloud Experiment (E-PEACE). The CVI inlet addresses three key issues in previous designs; in particular, the inlet operates with: (i) negligible organic contamination; (ii) a significant sample flow rate to downstream instruments (~15 l min −1 ) that reduces the need for dilution; and (iii) a high level of accessibility to the probe interior for cleaning. Wind tunnel experiments characterized the cut size of sampled droplets and the particle size-dependent transmission efficiency in various parts of the probe. For a range of counter-flow rates and air velocities, the measured cut size was between 8.7–13.1 μm. The percentage error between cut size measurements and predictions from aerodynamic drag theory are less than 13%. The CVI was deployed on the Center for Interdisciplinary Remotely-Piloted Aircraft Studies (CIRPAS) Twin Otter for thirty flights during E-PEACE to study aerosol-cloud-radiation interactions off the central coast of California between July and August 2011. Results are reported to assess the performance of the inlet including comparisons of particle number concentration downstream of the CVI and cloud drop number concentration measured by two independent aircraft probes. Measurements downstream the CVI are also examined from one representative case flight coordinated with shipboard-emitted smoke that was intercepted in cloud by the Twin Otter.

Description: Implications of observed inconsistencies in carbonate chemistry measurements for ocean acidification studies Biogeosciences Discussions, 9, 1781-1792, 2012 Author(s): C. J. M. Hoppe, G. Langer, S. D. Rokitta, D. A. Wolf-Gladrow, and B. Rost The growing field of ocean acidification research is concerned with the investigation of organisms' responses to increasing p CO 2 values. One important approach in this context is culture work using seawater with adjusted CO 2 levels. As aqueous p CO 2 is difficult to measure directly in small scale experiments, it is generally calculated from two other measured parameters of the carbonate system (often A T , C T or pH). Unfortunately, the overall uncertainties of measured and subsequently calculated values are often unknown. Especially under high p CO 2 , this can become a severe problem with respect to the interpretation of physiological and ecological data. In the few datasets from ocean acidification research where all three of these parameters were measured, p CO 2 values calculated from A T and C T are typically about 30 % lower (i.e. ~300 μatm at a target p CO 2 of 1000 μatm) than those calculated from A T and pH or C T and pH. This study presents and discusses these discrepancies as well as likely consequences for the ocean acidification community. Until this problem is solved, one has to consider that calculated parameters of the carbonate system (e.g. p CO 2 , calcite saturation state) may not be comparable between studies, and that this may have important implications for the interpretation of CO 2 perturbation experiments.

Description: The isotopic composition of precipitation from a winter storm – a case study with the limited-area model COSMO iso Atmospheric Chemistry and Physics, 12, 1629-1648, 2012 Author(s): S. Pfahl, H. Wernli, and K. Yoshimura Stable water isotopes are valuable tracers of the atmospheric water cycle, and potentially provide useful information also on weather-related processes. In order to further explore this potential, the water isotopes H 2 18 O and HDO are incorporated into the limited-area model COSMO. In a first case study, the new COSMO iso model is used for simulating a winter storm event in January 1986 over the eastern United States associated with intense frontal precipitation. The modelled isotope ratios in precipitation and water vapour are compared to spatially distributed δ 18 O observations. COSMO iso very accurately reproduces the statistical distribution of δ 18 O in precipitation, and also the synoptic-scale spatial pattern and temporal evolution agree well with the measurements. Perpendicular to the front that triggers most of the rainfall during the event, the model simulates a gradient in the isotopic composition of the precipitation, with high δ 18 O values in the warm air and lower values in the cold sector behind the front. This spatial pattern is created through an interplay of large scale air mass advection, removal of heavy isotopes by precipitation at the front and microphysical interactions between rain drops and water vapour beneath the cloud base. This investigation illustrates the usefulness of high resolution, event-based model simulations for understanding the complex processes that cause synoptic-scale variability of the isotopic composition of atmospheric waters. In future research, this will be particularly beneficial in combination with laser spectrometric isotope observations with high temporal resolution.

Description: Improved cloud screening in MAIAC aerosol retrievals using spectral and spatial analysis Atmospheric Measurement Techniques Discussions, 5, 1575-1595, 2012 Author(s): A. Lyapustin, Y. Wang, I. Laszlo, and S. Korkin An improved cloud/snow screening technique in the Multi-Angle Implementation of Atmospheric Correction (MAIAC) algorithm is described. It is implemented as part of MAIAC aerosol retrievals based on analysis of spectral residuals and spatial variability. Comparisons with AERONET aerosol measurements and a large-scale MODIS data analysis show strong suppression of aerosol optical depth outliers due to unresolved clouds and snow. At the same time, the developed filter does not reduce the aerosol retrieval capability at high 1 km resolution in strongly inhomogeneous environments, such as near centers of the active fires. Despite significant improvement, the optical depth outliers in high spatial resolution data are and will remain the problem to be addressed by the application-dependent specialized filtering techniques.

Description: Investigating riparian groundwater flow close to a losing river using diurnal temperature oscillations at high vertical resolution Hydrology and Earth System Sciences, 16, 473-487, 2012 Author(s): T. Vogt, M. Schirmer, and O. A. Cirpka River-water infiltration is of high relevance for hyporheic and riparian groundwater ecology as well as for drinking water supply by river-bank filtration. Heat has become a popular natural tracer to estimate exchange rates between rivers and groundwater. However, quantifying flow patterns and velocities is impeded by spatial and temporal variations of exchange fluxes, insufficient sensors spacing during field investigations, or simplifying assumptions for analysis or modeling such as uniform flow. The objective of this study is to investigate lateral shallow groundwater flow upon river-water infiltration at the shoreline of the riverbed and in the adjacent riparian zone of the River Thur in northeast Switzerland. Here we have applied distributed temperature sensing (DTS) along optical fibers wrapped around tubes to measure high-resolution vertical temperature profiles of the unsaturated zone and shallow riparian groundwater. Diurnal temperature oscillations were tracked in the subsurface and analyzed by means of dynamic harmonic regression to extract amplitudes and phase angles. Subsequent calculations of amplitude attenuation and time shift relative to the river signal show in detail vertical and temporal variations of heat transport in shallow riparian groundwater. In addition, we apply a numerical two-dimensional heat transport model for the unsaturated zone and shallow groundwater to obtain a better understanding of the observed heat transport processes in shallow riparian groundwater and to estimate the groundwater flow velocity. Our results show that the observed riparian groundwater temperature distribution cannot be described by uniform flow, but rather by horizontal groundwater flow velocities varying over depth. In addition, heat transfer of diurnal temperature oscillations from the losing river through shallow groundwater is influenced by thermal exchange with the unsaturated zone. Neglecting the influence of the unsaturated zone would cause biased interpretation and underestimation of groundwater flow velocities. The combination of high resolution field data and modeling shows the complex hydraulic and thermal processes occurring in shallow riparian groundwater close to losing river sections as well as potential errors sources for interpreting diurnal temperature oscillations in such environments.

Description: Impacts of inhomogeneous landscapes in oasis interior on the oasis self-maintaining mechanism by integrating numerical model with satellite data Hydrology and Earth System Sciences Discussions, 9, 1979-2004, 2012 Author(s): X. Meng, S. Lu, T. Zhang, Y. Ao, S. Li, Y. Bao, L. Wen, and S. Luo Mesoscale meteorological modeling is an important tool to help understand the energy budget of the oasis. While basic dynamic and thermodynamic processes for oasis self-maintaining in the desert environment is well investigated, influence of heterogeneous landscapes of oasis interior on the processes are still important and remain to be investigated. In this study, two simulations are designed for investigating the influence of inhomogeneity. In the first case, land surface parameters including land-use types, vegetation cover fraction, and surface layer soil moisture are derived by satellite remote sensing data from EOS/MODIS, and then be used specify the respective options in the MM5 model, to describe a real inhomogeneity for the oasis interior. In the other run, land use types are set to MM5 default, in which landscapes in the oasis interior is relative uniform, and then surface layer soil moisture and vegetation fraction is set to be averages of the first case for the respective oasis and desert surface lying, to represent a relative homogeneity. Results show that the inhomogeneity leads to a weaker oasis "cold-wet island" effect and a stronger turbulence over the oasis interior, both of which will reduce the oasis-desert secondary circulation and increase the evaporation over the oasis, resulting in a negative impact on the oasis self-protecting mechanism. The simulation of homogeneity indicates that the oasis may be more stable even with relative lower soil moisture if landscapes in the oasis interior are comparatively uniform.

Description: Measuring variations of δ 18 O and δ 2 H in atmospheric water vapour using laser spectroscopy: an instrument characterisation study Atmospheric Measurement Techniques Discussions, 5, 1597-1655, 2012 Author(s): F. Aemisegger, P. Sturm, P. Graf, H. Sodemann, S. Pfahl, A. Knohl, and H. Wernli Variations of stable water isotopes in water vapour have become measurable at a measurement frequency of about 1 Hz in recent years using novel laser spectroscopic techniques. This enables us to perform continuous measurements for process-based investigations of the atmospheric water cycle at the time scales relevant for synoptic meteorology. An important prerequisite for the interpretation of data from automated field measurements lasting for several weeks or months is a detailed knowledge about instrument properties and the sources of measurement uncertainty. We present here a comprehensive characterisation and comparison study of two commercial laser spectroscopic systems based on cavity ring-down spectroscopy (Picarro) and off-axis integrated cavity output spectroscopy (Los Gatos Research). The uncertainty components of the measurements were first assessed in laboratory experiments, focussing on the effects of (i) water vapour mixing ratio, (ii) measurement stability, (iii) uncertainties due to calibration and (iv) response times of the isotope measurements due to adsorption-desorption processes on the tubing and measurement cavity walls. Based on the experience from our laboratory experiments we set up a one-week field campaign for comparing measurements of the ambient isotope signals of the two laser spectroscopic systems. The optimal calibration strategy determined for both instruments was applied as well as the correction functions for water vapour mixing ratio effects. The root mean square difference between the isotope signals from the two instruments during the field deployment was 2.3‰ for δ 2 H, 0.5‰ for δ 18 O and 3.1‰ for deuterium excess. These uncertainty estimates from field measurements compare well to those found in the laboratory experiments. The present quality of measurements from laser spectroscopic instruments combined with a calibration system opens new possibilities for investigating the atmospheric water cycle and the land-atmosphere moisture fluxes.

Description: In plain sight: the Chesapeake Bay crater ejecta blanket Solid Earth Discussions, 4, 363-428, 2012 Author(s): D. L. Griscom The discovery nearly two decades ago of a 90 km-diameter impact crater below the lower Chesapeake Bay has gone unnoted by the general public because to date all published literature on the subject has described it as "buried". To the contrary, evidence is presented here that the so-called "upland deposits" that blanket ∼5000 km 2 of the U.S. Middle-Atlantic Coastal Plain (M-ACP) display morphologic, lithologic, and stratigraphic features consistent with their being ejecta from the 35.4 Ma Chesapeake Bay Impact Structure (CBIS) and absolutely inconsistent with the prevailing belief that they are of fluvial origin. Specifically supporting impact origin are the facts that (i) a 95 %-pure iron ore endemic to the upland deposits of southern Maryland, eastern Virginia, and the District of Columbia has previously been proven to be impactoclastic in origin, (ii) this iron ore welds together a small percentage of well-rounded quartzite pebbles and cobbles of the upland deposits into brittle sheets interpretable as "spall plates" created in the interference-zone of the CBIS impact, (iii) the predominantly non-welded upland gravels have long ago been shown to be size sorted with an extreme crater-centric gradient far too large to have been the work of rivers, but well explained as atmospheric size-sorted interference-zone ejecta, (iv) new evidence is provided here that ~60 % of the non-welded quartzite pebbles and cobbles of the (lower lying) gravel member of the upland deposits display planar fractures attributable to interference-zone tensile waves, (v) the (overlying) loam member of the upland deposits is attributable to base-surge-type deposition, (vi) several exotic clasts found in a debris flow topographically below the upland deposits can only be explained as jetting-phase crater ejecta, and (vii) an allogenic granite boulder found among the upland deposits is deduced to have been launched into space and sculpted by hypervelocity air friction during reentry. An idealized calculation of the CBIS ejecta-blanket elevation profile minutes after the impact was carried out founded on well established rules for explosion and impact-generated craters. This profile is shown here to match the volume of the upland deposits ≥170 km from the crater center. Closer to the crater, much of the "postdicted" ejecta blanket has clearly been removed by erosion. Nevertheless the Shirley and fossil-free Bacons Castle Formations, located between the upland deposits and the CBIS interior and veneering the present day surface with units ∼10–20 m deep, are respectively identified as curtain- and excavation-phase ejecta. The neritic-fossil-bearing Calvert Formation external to the crater is deduced to be of Eocene age (as opposed to early Miocene as currently believed), preserved by the armoring effects of the overlying CBIS ejecta composed of the (distal) upland deposits and the (proximal) Bacons Castle Formation. The lithofacies of the in-crater Calvert Formation can only have resulted from inward mass wasting of the postdicted ejecta blanket, vestiges of which (i.e. the Bacons Castle and Shirley Formations) still overlap the crater rim and sag into its interior, consistent with this expectation. Because there appear to be a total of ∼10 000 km 2 of CBIS ejecta lying on the present-day surface, future research should center the stratigraphic, lithologic, and petrologic properties of these ejecta versus both radial distance from the crater center (to identify ejecta from different ejection stages) and circumferentially at fixed radial distances (to detect possible anisotropies relating the impact angle and direction of approach of the impactor). The geological units described here may comprise the best preserved, and certainly the most accessible, ejecta blanket of a major crater on the Earth's surface and therefore promise to be a boon to the field of impact geology. As a corollary, a major revision of the current stratigraphic column of the M-ACP will be necessary.

Description: Uncertainty in future solid ice discharge from Antarctica The Cryosphere Discussions, 6, 673-714, 2012 Author(s): R. Winkelmann, A. Levermann, K. Frieler, and M. A. Martin Future solid ice discharge from Antarctica under climate scenarios based on the Extended Concentration Pathways is investigated with the Potsdam Parallel Ice Sheet Model (PISM-PIK), a shallow model with a consistent representation of the ice flow in sheet, shelves and the transition zone. Both the uncertainty in the climate forcing as well as the intra-model uncertainty are combined into a probability distribution for solid ice discharge from Antarctica until the year 2500 under the ECP scenarios: All simulations are performed for a 81-member perturbed-physics ensemble and the likely ranges of surface and ocean warming under the emission pathways derived from the results of 20 CMIP3-AOGCMS. The effects of surface warming, ocean warming and increased precipitation on solid ice discharge are separately considered. We find that solid ice discharge caused by enhanced sub-shelf melting exceeds that caused by surface warming. Increasing precipitation leads to a change from net sea-level rise to sea-level drop. Our results suggest that the history of the ice-sheet plays an important role with respect to projections of solid ice discharge. Although all climate-change-forced simulations begin with the year 1850, the ice discharge around 2000 is significantly smaller than observed. Observed changes in ice discharge are reached around 2077 under the ECP-8.5 scenario. During the subsequent century, ice discharge reaches up to 0.24 m.

Description: Floodwater utilisation values of wetland services – a case study in Northeastern China Natural Hazards and Earth System Science, 12, 341-349, 2012 Author(s): S. B. Lü, S. G. Xu, and F. Feng Water plays a significant role in wetlands. Floodwater utilisation in wetlands brings a wide range of wetland services, from goods production and water regulation to animal protection and aesthetics related to water supply in wetlands. In this study, the floodwater utilisation values of wetland services were estimated within the Momoge wetland and Xianghai wetland in western Jilin province of northeastern China. From 2003 to 2008, the floodwater diverted from the Nenjiang and Tao'er River is 381 million m 3 , which translates into a monetary value of approximately 1.35 billion RMB in 2008 (RMB: Chinese Currency, RMB 6.80 = US$ 1), and the ratio of economic value, eco-environmental value, and social value is 1:12:2. Besides the monetary value of the water itself, excessive floodwater utilisation may bring losses to wetlands; the threshold floodwater utilisation volumes in wetlands are discussed. Floodwater utilisation can alleviate water shortages in wetlands, and the evaluation of floodwater utilisation in wetland services in monetary terms is a guide for the effective use of the floodwater resources and for the conservation of wetlands.

Description: On closure parameter estimation in chaotic systems Nonlinear Processes in Geophysics, 19, 127-143, 2012 Author(s): J. Hakkarainen, A. Ilin, A. Solonen, M. Laine, H. Haario, J. Tamminen, E. Oja, and H. Järvinen Many dynamical models, such as numerical weather prediction and climate models, contain so called closure parameters. These parameters usually appear in physical parameterizations of sub-grid scale processes, and they act as "tuning handles" of the models. Currently, the values of these parameters are specified mostly manually, but the increasing complexity of the models calls for more algorithmic ways to perform the tuning. Traditionally, parameters of dynamical systems are estimated by directly comparing the model simulations to observed data using, for instance, a least squares approach. However, if the models are chaotic, the classical approach can be ineffective, since small errors in the initial conditions can lead to large, unpredictable deviations from the observations. In this paper, we study numerical methods available for estimating closure parameters in chaotic models. We discuss three techniques: off-line likelihood calculations using filtering methods, the state augmentation method, and the approach that utilizes summary statistics from long model simulations. The properties of the methods are studied using a modified version of the Lorenz 95 system, where the effect of fast variables are described using a simple parameterization.

Description: Sea surface temperature anomalies, seasonal cycle and trend regimes in the Eastern Pacific coast Ocean Science, 8, 81-90, 2012 Author(s): A. Ramos-Rodríguez, D. B. Lluch-Cota, S. E. Lluch-Cota, and A. Trasviña-Castro We used the extended reconstruction of sea surface temperature (ERSST) to analyze the variation of surface temperature and the seasonal cycle along the coast of the eastern Pacific (60° N–60° S, 61 pixels alongshore) from 1950 to 2010 (732 months). First, we analyzed the monthly anomalies and looked for a relationship of such anomalies with total solar irradiance (TSI) and then the Regime Shift Detector (RSD) was applied to detect possible temperature regimes in the series. Afterwards, we calculated a yearly temperature range per pixel (amplitude of seasonal cycle) and through the subtraction of a latitudinal theoretical curve of temperature based on solar irradiance, the residuals of the seasonal cycle were obtained. The results showed an almost complete spatial synchrony and dominance of negative anomalies from 1950 to mid-late 1970's, with a switch to near-zero and positive anomalies that lasted up to late 1990's when a new shift to negative values was detected. Such a shift lasted until the early 2000's when positive anomalies appeared again but there was a change to negative anomalies in the late 2000's. These results were supported by the RSD. The TSI variability shows a clear relationship with that of sea surface temperature anomalies and with the regime changes. This is probably due to a difference in the amount of energy received from the sun. Comparing the "cool regime" versus the "warm regime", the second one received 0.39% more energy (approximately 3 × 10 8 J m −2 ) from the sun. Seasonal cycles show larger ranges at northern latitudes (〉40° N), northern tropical-temperate transition zone (20°–26° N) and in the tropical-equatorial band (0°–30° S). The smallest ranges occur at 0°–16° N and 50°–60° S. The residuals (seasonal minus the theoretical curve) indicated a clear modulation due to advection by ocean currents.

Description: The Lumen Gini Coefficient: a satellite imagery derived human development index Social Geography Discussions, 8, 27-59, 2012 Author(s): C. D. Elvidge, K. E. Baugh, S. J. Anderson, P. C. Sutton, and T. Ghosh The "Lumen Gini Coefficient" is a simple, objective, spatially explicit and globally available empirical measurement of human development derived solely from nighttime satellite imagery and population density. There is increasing recognition that the distribution of wealth and income amongst the population in a nation or region correlates strongly with both the overall happiness of that population and the environmental quality of that nation or region. Measuring the distribution of wealth and income at national and regional scales is an interesting and challenging problem. Gini coefficients derived from Lorenz curves are a well-established method of measuring income distribution. Nonetheless, there are many shortcomings of the Gini coefficient as a measure of income or wealth distribution. Gini coefficients are typically calculated using national level data on the distribution of income through the population. Such data are not available for many countries and the results are generally limited to single values representing entire countries. In this paper we develop an alternative measure of the distribution of "human development", called the "Lumen Gini coefficient", that is derived without the use of monetary measures of wealth and is capable of providing a spatial depiction of differences in development within countries.

Description: Constraining fault interpretation through tomographic velocity gradients: application to northern Cascadia Solid Earth, 3, 53-61, 2012 Author(s): K. Ramachandran Spatial gradients of tomographic velocities are seldom used in interpretation of subsurface fault structures. This study shows that spatial velocity gradients can be used effectively in identifying subsurface discontinuities in the horizontal and vertical directions. Three-dimensional velocity models constructed through tomographic inversion of active source and/or earthquake traveltime data are generally built from an initial 1-D velocity model that varies only with depth. Regularized tomographic inversion algorithms impose constraints on the roughness of the model that help to stabilize the inversion process. Final velocity models obtained from regularized tomographic inversions have smooth three-dimensional structures that are required by the data. Final velocity models are usually analyzed and interpreted either as a perturbation velocity model or as an absolute velocity model. Compared to perturbation velocity model, absolute velocity models have an advantage of providing constraints on lithology. Both velocity models lack the ability to provide sharp constraints on subsurface faults. An interpretational approach utilizing spatial velocity gradients applied to northern Cascadia shows that subsurface faults that are not clearly interpretable from velocity model plots can be identified by sharp contrasts in velocity gradient plots. This interpretation resulted in inferring the locations of the Tacoma, Seattle, Southern Whidbey Island, and Darrington Devil's Mountain faults much more clearly. The Coast Range Boundary fault, previously hypothesized on the basis of sedimentological and tectonic observations, is inferred clearly from the gradient plots. Many of the fault locations imaged from gradient data correlate with earthquake hypocenters, indicating their seismogenic nature.

Description: The next generation Nature Reviews Microbiology 10, 157 (2012). doi:10.1038/nrmicro2759 A better understanding of the factors that affect children's desire to follow a career in science is vital if we are to inspire the next generation of scientists.

Description: Bacterial pathogenesis: Controlling Fic proteins Nature Reviews Microbiology 10, 160 (2012). doi:10.1038/nrmicro2757 Author: Sheilagh Molloy A general mechanism for controlling the adenylylation activity of Fic proteins has been identified.

Description: Cellular microbiology: An unconventional exit for Brucella Nature Reviews Microbiology 10, 160 (2012). doi:10.1038/nrmicro2749 Author: Rachel David Brucella abortus hijacks autophagy components to complete its life cycle.

Description: Bacterial genomics: Universal bacterial barcode Nature Reviews Microbiology 10, 161 (2012). doi:10.1038/nrmicro2754 Author: Andrew Jermy Although analysis of 16S ribosomal RNA has revolutionized the ability to classify bacteria and understand their phylogenetic relationships, it is less useful for classifying related strains of bacteria that share similar or identical 16S rRNA gene sequences but exhibit distinct properties. To address this, Jolley

Description: Viral infection: Promiscuous packaging Nature Reviews Microbiology 10, 161 (2012). doi:10.1038/nrmicro2755 Author: Andrew Jermy The high error rate of viral polymerases allows viruses to rapidly adapt to a new host, and each replication cycle results in a heterogeneous population of viruses with distinct genetic content. Routh et al. used an unbiased next-generation sequencing approach to look at the

Description: Biofilms: Biofilms take shape Nature Reviews Microbiology 10, 162 (2012). doi:10.1038/nrmicro2756 Author: Sheilagh Molloy Phenol-soluble modulins have a key role in the structuring and detachment of Staphylococcus aureus biofilms.

Description: Adapting to domesticity Nature Reviews Microbiology 10, 163 (2012). doi:10.1038/nrmicro2752 Author: Adam J. Reid This month's Genome Watch highlights new insights into the impact of food production on the evolution of microorganisms.

Description: In the news Nature Reviews Microbiology 10, 164 (2012). doi:10.1038/nrmicro2760 Our monthly round up of infectious diseases news, which this month includes the spread of a new virus affecting livestock in Europe, the H5N1 controversy and a report that mortality from malaria may be much higher than was previously suspected.

Description: Bacterial pathogenesis: UPEC helps host to exfoliate Nature Reviews Microbiology 10, 159 (2012). doi:10.1038/nrmicro2758 Author: Andrew Jermy UPEC pore-forming toxin α-haemolysin (HlyA) stimulates the proteolysis of host proteins involved in cell adhesion, resulting in uroepithelial cell exfoliation.

Description: A bacterial siren song: intimate interactions between Neisseria and neutrophils Nature Reviews Microbiology 10, 178 (2012). doi:10.1038/nrmicro2713 Authors: Alison K. Criss & H. Steven Seifert Neisseria gonorrhoeae and Neisseria meningitidis are Gram-negative bacterial pathogens that are exquisitely adapted for growth at human mucosal surfaces and for efficient transmission between hosts. One factor that is essential to neisserial pathogenesis is the interaction between the bacteria and neutrophils, which are

Description: Phage biology: T4 to TA, checkmate Nature Reviews Microbiology 10, 160 (2012). doi:10.1038/nrmicro2748 Author: Andrew Jermy Phage T4 encodes an antitoxin that is active against multiple Escherichia coli toxins.

Description: ppGpp: magic beyond RNA polymerase Nature Reviews Microbiology 10, 203 (2012). doi:10.1038/nrmicro2720 Authors: Zachary D. Dalebroux & Michele S. Swanson During stress, bacteria undergo extensive physiological transformations, many of which are coordinated by ppGpp. Although ppGpp is best known for enhancing cellular resilience by redirecting the RNA polymerase (RNAP) to certain genes, it also acts as a signal in many other cellular processes in bacteria.

Description: Structural insights into the coupling of virion assembly and rotavirus replication Nature Reviews Microbiology 10, 165 (2012). doi:10.1038/nrmicro2673 Authors: Shane D. Trask, Sarah M. McDonald & John T. Patton Viral replication is rapid and robust, but it is far from a chaotic process. Instead, successful production of infectious progeny requires that events occur in the correct place and at the correct time. Rotaviruses (segmented double-stranded RNA viruses of the Reoviridae family) seem to

Description: Genome analyses highlight the different biological roles of cellulases Nature Reviews Microbiology 10, 227 (2012). doi:10.1038/nrmicro2729 Authors: Felix Mba Medie, Gideon J. Davies, Michel Drancourt & Bernard Henrissat Cellulolytic enzymes have been the subject of renewed interest owing to their potential role in the conversion of plant lignocellulose to sustainable biofuels. An analysis of ∼1,500 complete bacterial genomes, presented here, reveals that ∼40% of the genomes of sequenced bacteria encode at least one

Description: Host response: Phagocytosis runs like clockwork Nature Reviews Microbiology 10, 162 (2012). doi:10.1038/nrmicro2751 Author: Christina Tobin Kåhrström In Drosophila melanogaster, phagocytosis of specific bacteria is regulated by a circadian protein.

Description: From self sufficiency to dependence: mechanisms and factors important for autotransporter biogenesis Nature Reviews Microbiology 10, 213 (2012). doi:10.1038/nrmicro2733 Authors: Denisse L. Leyton, Amanda E. Rossiter & Ian R. Henderson Autotransporters are a superfamily of proteins that use the type V secretion pathway for their delivery to the surface of Gram-negative bacteria. At first glance, autotransporters look to contain all the functional elements required to promote their own secretion: an amino-terminal signal peptide to mediate

Description: Virophages question the existence of satellites Nature Reviews Microbiology 10, 234 (2012). doi:10.1038/nrmicro2676-c3 Authors: Christelle Desnues & Didier Raoult In a recent Comment (Virophages or satellite viruses? Nature Rev. Microbiol.9 762–763(2011)), Mart Krupovic and Virginija Cvirkaite-Krupovic argued that the recently described virophages, Sputnik and Mavirus, should be classified as satellite viruses. In a response, to which Krupovic and

Description: Evolution: Co-evolution promotes innovation Nature Reviews Microbiology 10, 161 (2012). doi:10.1038/nrmicro2753 Author: Christina Tobin Kåhrström The interaction between a viral ligand and its cellular receptor is highly specific and limits host range. A recent paper describes a mutant form of bacteriophage λ that targets an alternative receptor after co-evolution with Escherichia coli. When phage cI26 was cultured with E.

Description: Towards a more comprehensive classification of satellite viruses Nature Reviews Microbiology 10, 234 (2012). doi:10.1038/nrmicro2676-c4 Authors: Mart Krupovic & Virginija Cvirkaite-Krupovic Our recent Comment on the nature and place of the so-called virophages in the viral world (Virophages or satellite viruses? Nature Rev. Microbiol.9 762–763 (2011)) has caught the attention of researchers favouring the virophage concept and initiated a debate on this

Description: Computational tools for the synthetic design of biochemical pathways Nature Reviews Microbiology 10, 191 (2012). doi:10.1038/nrmicro2717 Authors: Marnix H. Medema, Renske van Raaphorst, Eriko Takano & Rainer Breitling As the field of synthetic biology is developing, the prospects for de novo design of biosynthetic pathways are becoming more and more realistic. Hence, there is an increasing need for computational tools that can support these efforts. A range of algorithms has been developed

Description: Hydrogeology and hydrogeochemistry of an alkaline volcanic area: the NE Mt. Meru slope (East African Rift – Northern Tanzania) Hydrology and Earth System Sciences, 16, 529-541, 2012 Author(s): G. Ghiglieri, D. Pittalis, G. Cerri, and G. Oggiano The objective of this study is to analyze the geochemical conditions associated with the presence of fluoride (F − ) in the groundwater of an area of Northern Tanzania. The studied aquifers are composed of volcanic rocks such as phonolitic and nephelinitic lavas, basalts, lahars of various ages and mantling ash. Sedimentary rocks consisting of fine-grained alluvial and lacustrine deposits occur as well. Samples collected from springs, borehole and surface water, during two monitoring surveys, were analyzed for the various physico-chemical and isotopic parameters. The geochemical composition of water is typically sodium bicarbonate. High values of F − (up to 68 mg l −1 ) were recorded. The highest values of fluoride agreed with the highest values of pH, sodium and bicarbonate. Dissolution of major ions, exchange processes and precipitation of Ca 2+ from super-saturated solutions joined with the local permeability and hydraulic gradients, control the fluoride mobilization and the contamination of the area.

Description: Calibration and evaluation of a semi-distributed watershed model of sub-Saharan Africa using GRACE data Hydrology and Earth System Sciences Discussions, 9, 2071-2120, 2012 Author(s): H. Xie, L. Longuevergne, C. Ringler, and B. Scanlon Irrigation development is rapidly expanding in mostly rainfed Sub-Saharan Africa. This expansion underscores the need for a more comprehensive understanding of water resources beyond surface water. Gravity Recovery and Climate Experiment (GRACE) satellites provide valuable information on spatio-temporal variability of water storage. The objective of this study was to calibrate and evaluate a semi-distributed regional-scale hydrological model, or a large-scale application of the Soil and Water Assessment Tool (SWAT) model, for basins in Sub-Saharan Africa using seven-year (2002–2009) 10-day GRACE data. Multi-site river discharge data were used as well, and the analysis was conducted in a multi-criteria framework. In spite of the uncertainty arising from the tradeoff in optimizing model parameters with respect to two non-commensurable criteria defined for two fluxes, it is concluded that SWAT can perform well in simulating total water storage variability in most areas of Sub-Saharan Africa, which have semi-arid and sub-humid climates, and that among various water storages represented in SWAT, the water storage variations from soil, the vadose zone, and groundwater are dominant. On the other hand, the study also showed that the simulated total water storage variations tend to have less agreement with the GRACE data in arid and equatorial humid regions, and the model-based partition of total water storage variations into different water storage compartments could be highly uncertain. Thus, future work will be needed for model enhancement in these areas with inferior model fit and for uncertainty reduction in component-wise estimation of water storage variations.

Description: Estimation of soil redistribution rates due to snow cover related processes in a mountainous area (Valle d'Aosta, NW Italy) Hydrology and Earth System Sciences, 16, 517-528, 2012 Author(s): E. Ceaglio, K. Meusburger, M. Freppaz, E. Zanini, and C. Alewell Mountain areas are widely affected by soil erosion, which is generally linked to runoff processes occurring in the growing season and snowmelt period. Also processes like snow gliding and full-depth snow avalanches may be important factors that can enhance soil erosion, however the role and importance of snow movements as agents of soil redistribution are not well understood yet. The aim of this study was to provide information on the relative importance of snow related processes in comparison to runoff processes. In the study area, which is an avalanche path characterized by intense snow movements, soil redistribution rates were quantified with two methods: (i) by field measurements of sediment yield in an avalanche deposition area during 2009 and 2010 winter seasons; (ii) by caesium-137 method, which supplies the cumulative net soil loss/gain since 1986, including all the soil erosion processes. The snow related soil accumulation estimated with data from the deposit area (27.5 Mg ha −1 event −1 and 161.0 Mg ha −1 event −1 ) was not only higher than the yearly sediment amounts, reported in literature, due to runoff processes, but it was even more intense than the yearly total deposition rate assessed with 137 Cs (12.6 Mg ha −1 yr −1 ). The snow related soil erosion rates estimated from the sediment yield at the avalanche deposit area (3.7 Mg ha −1 and 20.8 Mg ha −1 ) were greater than the erosion rates reported in literature and related to runoff processes; they were comparable to the yearly total erosion rates assessed with the 137 Cs method (13.4 Mg ha −1 yr −1 and 8.8 Mg ha −1 yr −1 ). The 137 Cs method also showed that, where the ground avalanche does not release, the erosion and deposition of soil particles from the upper part of the basin was considerable and likely related to snow gliding. Even though the comparison of both the approaches is linked to high methodological uncertainties, mainly due to the different spatial and temporal scales considered, we still can deduce, from the similarity of the erosion rates, that soil redistribution in this catchment is driven by snow movement, with a greater impact in comparison to the runoff processes occurring in the snow-free season. Nonetheless, the study highlights that soil erosion processes due to the snow movements should be considered in the assessment of soil vulnerability in mountain areas, as they significantly determine the pattern of soil redistribution.

Description: Development of an operational coastal flooding early warning system Natural Hazards and Earth System Science, 12, 379-390, 2012 Author(s): D.-J. Doong, L. Z.-H. Chuang, L.-C. Wu, Y.-M. Fan, C. C. Kao, and J.-H. Wang Coastal floods are a consistent threat to oceanfront countries, causing major human suffering and substantial economic losses. Climate change is exacerbating the problem. An early warning system is essential to mitigate the loss of life and property from coastal flooding. The purpose of this study is to develop a coastal flooding early warning system (CoFEWs) by integrating existing sea-state monitoring technology, numerical ocean forecasting models, historical database and experiences, as well as computer science. The proposed system has capability of offering data for the past, information for the present and future. The system was developed for the Taiwanese coast due to its frequent threat by typhoons. An operational system without any manual work is the basic requirement of the system. Integration of various data sources is the system kernel. Numerical ocean models play an important role within the system because they provide data for assessment of possible flooding. The regional wave model (SWAN) that nested with the large domain wave model (NWW III) is operationally set up for coastal wave forecasting, in addition to the storm surge predicted by a POM model. Data assimilation technology is incorporated for enhanced accuracy. A warning signal is presented when the storm water level that accumulated from astronomical tide, storm surge, and wave-induced run-up exceeds the alarm sea level. This warning system has been in practical use for coastal flooding damage mitigation in Taiwan for years. An example of the system operation during the Typhoon Haitung which struck Taiwan in 2005 is illustrated in this study.

Description: Extreme summer temperatures in Western Europe Advances in Science and Research, 8, 5-9, 2012 Author(s): C. Simolo, M. Brunetti, M. Maugeri, and T. Nanni We discuss the evolution of summer temperature extremes over Western Europe during 1961–2004 in the context of current climate warming. Using a parametric approach, we investigate the role of properties and changes in probability density functions of daily temperatures in modifying the frequency of severe, isolated events. In this perspective, the recent intensification of extremely warm events over Europe turns out to be well consistent with a pure, nonuniform shift of mean values, with no room for conjectures about increasing temperature variability.

Description: Riverine influence on the tropical Atlantic Ocean biogeochemistry Biogeosciences Discussions, 9, 1945-1969, 2012 Author(s): L. Cotrim da Cunha and E. T. Buitenhuis We assess the role of riverine inputs of N, Si, Fe, organic and inorganic C in the tropical Atlantic Ocean using a global ocean biogeochemistry model. We use two sensitivity tests to investigate the role of the western (South American Rivers) and eastern (African Rivers) riverine nutrient inputs on the tropical Atlantic Ocean biogeochemistry (between 20° S–20° N and 70° W–20°). Increased nutrient availability from river inputs in this area (compared to an extreme scenario with no river nutrients) leads to an increase in 14 % (0.7 Pg C a −1 ) in open ocean primary production (PP), and 21 % (0.2 Pg C a −1 ) in coastal ocean PP. We estimate very modest increases in open and coastal ocean export production and sea-air CO 2 fluxes. Results suggest that in the tropical Atlantic Ocean, the large riverine nutrient inputs on the western side have a larger impact on primary production and sea-air CO 2 exchanges. On the other hand, African river inputs, although smaller than South American inputs, have larger impact on the coastal and open tropical Atlantic Ocean export production. This is probably due to a combination of nutrient trapping in upwelling areas off the Congo River outflow, and differences in delivered nutrient ratios leading to alleviation in limitation conditions mainly for diatoms.

Description: Quantification of diesel exhaust gas phase organics by a thermal desorption proton transfer reaction mass spectrometer Atmospheric Chemistry and Physics Discussions, 12, 5389-5423, 2012 Author(s): M. H. Erickson, H. W. Wallace, and B. T. Jobson A new approach was developed to measure the total abundance of long chain alkanes (C 12 and above) in urban air using thermal desorption with a proton transfer reaction mass spectrometer (PTR-MS). These species are emitted in diesel exhaust and may be important precursors to secondary organic aerosol production in urban areas. Long chain alkanes undergo dissociative proton transfer reactions forming a series of fragment ions with formula C n H 2n+1 . The yield of the fragment ions is a function of drift conditions. At a drift field strength of 80 Townsends, the most abundant ion fragments from C 10 to C 16 n -alkanes were m/z 57, 71 and 85. The PTR-MS is insensitive to n -alkanes less than C 8 but displays an increasing sensitivity for larger alkanes. Higher drift field strengths yield greater normalized sensitivity implying that the proton affinity of the long chain n -alkanes is less than H 2 O. Analysis of diesel fuel shows the mass spectrum was dominated by alkanes (C n H 2 n +1 ), monocyclic aromatics, and an ion group with formula C n H 2 n −1 ( m/z 97, 111, 125, 139). The PTR-MS was deployed in Sacramento, CA during the Carbonaceous Aerosols and Radiative Effects Study field experiment in June 2010. The ratio of the m/z 97 to 85 ion intensities in ambient air matched that found in diesel fuel. Total diesel exhaust alkane concentrations calculated from the measured abundance of m/z 85 ranged from the method detection limit of ~1 μg m −3 to 100 μg m −3 in several air pollution episodes. The total diesel exhaust alkane concentration determined by this method was on average a factor of 10 greater than the sum of alkylbenzenes associated with spark ignition vehicle exhaust.

Description: Arctic climate response to forcing from light-absorbing particles in snow and sea ice in CESM Atmospheric Chemistry and Physics Discussions, 12, 5341-5388, 2012 Author(s): N. Goldenson, S. J. Doherty, C. M. Bitz, M. M. Holland, B. Light, and A. J. Conley The presence of light-absorbing aerosol particles deposited on arctic snow and sea ice influences the surface albedo, causing greater shortwave absorption, warming, and loss of snow and sea ice, lowering the albedo further. The Community Earth System Model version 1 (CESM1) now includes the radiative effects of light-absorbing particles in snow on land and sea ice and in sea ice itself. We investigate the model response to the deposition of black carbon and dust to both snow and sea ice. For these purposes we employ a slab ocean version of CESM1, using the Community Atmosphere Model version 4 (CAM4), run to equilibrium for year 2000 levels of CO 2 and fixed aerosol deposition. We construct experiments with and without aerosol deposition, with dust or black carbon deposition alone, and with varying quantities of black carbon and dust to approximate year 1850 and 2000 deposition fluxes. The year 2000 deposition fluxes of both dust and black carbon cause 1–2 °C of surface warming over large areas of the Arctic Ocean and sub-Arctic seas in autumn and winter and in patches of Northern land in every season. Atmospheric circulation changes are a key component of the surface-warming pattern. Arctic sea ice thins by on average about 30 cm. Simulations with year 1850 aerosol deposition are not substantially different from those with year 2000 deposition, given constant levels of CO 2 . The climatic impact of particulate impurities deposited over land exceeds that of particles deposited over sea ice. Even the surface warming over the sea ice and sea ice thinning depends more upon light-absorbing particles deposited over land. For CO 2 doubled relative to year 2000 levels, the climate impact of particulate impurities in snow and sea ice is substantially lower than for the year 2000 equilibrium simulation.

Description: Consistency between Fourier transform and small-volume few-wave decomposition for spectral and spatial variability of gravity waves above a typhoon Atmospheric Measurement Techniques Discussions, 5, 1763-1793, 2012 Author(s): C. I. Lehmann, Y.-H. Kim, P. Preusse, H.-Y. Chun, M. Ern, and S.-Y. Kim Convective gravity wave (GW) sources are spatially localized and emit at the same time waves with a wide spectrum of phase speeds. Any wave analysis therefore compromises between spectral and spatial resolution. Future satellite borne limb imagers will for a first time provide real 3d volumes of observations. These volumes will be however limited which will impose further constraints on the analysis technique. In this study a three dimensional few-wave appoach fitting sinusoidal waves to limited 3-D volumes is introduced. The method is applied to simulated GWs above typhoon Ewiniar and GW momentum flux is estimated from temperature fluctuations. Phase speed spectra as well as average profiles of positive, negative and net momentum fluxes are compared to momentum flux estimated by Fourier transform as well as spatial averaging of wind fluctuations. The results agree within 10–20%. The few-wave method can also reveal the spatial orientation of the GWs with respect to the source. The relevance of the results for different types of measurements as well as its applicability to model data is discussed.

Description: Modelling of organic aerosols over Europe (2002–2007) using a volatility basis set (VBS) framework with application of different assumptions regarding the formation of secondary organic aerosol Atmospheric Chemistry and Physics Discussions, 12, 5425-5485, 2012 Author(s): R. Bergström, H. A. C. Denier van der Gon, A. S. H. Prévôt, K. E. Yttri, and D. Simpson A new organic aerosol (OA) module has been implemented into the EMEP chemical transport model. Four different volatility basis set (VBS) schemes have been tested in long-term simulations for Europe, covering the six years 2002–2007. Different assumptions regarding partitioning of primary OA (POA) and aging of POA and secondary OA (SOA), have been explored. Model results are compared to filter measurements, AMS-data and source-apportionment studies, as well as to other model studies. The present study indicates that many different sources contribute significantly to OA in Europe. Fossil POA and oxidised POA, biogenic and anthropogenic SOA (BSOA and ASOA), residential burning of biomass fuels and wildfire emissions may all contribute more than 10% each over substantial parts of Europe. Simple VBS based OA models can give reasonably good results for summer OA but more observational studies are needed to constrain the VBS parameterisations and to help improve emission inventories. The volatility distribution of primary emissions is an important issue for further work. This study shows smaller contributions from BSOA to OA in Europe than earlier work, but relatively greater ASOA. BVOC emissions are highly uncertain and need further validation. We can not reproduce winter levels of OA in Europe, and there are many indications that the present emission inventories substantially underestimate emissions from residential wood burning in large parts of Europe.

Description: Profiles of CH 4 , HDO, H 2 O, and N 2 O with improved lower tropospheric vertical resolution from Aura TES radiances Atmospheric Measurement Techniques, 5, 397-411, 2012 Author(s): J. Worden, S. Kulawik, C. Frankenberg, V. Payne, K. Bowman, K. Cady-Peirara, K. Wecht, J.-E. Lee, and D. Noone Thermal infrared (IR) radiances measured near 8 microns contain information about the vertical distribution of water vapor (H 2 O), the water isotopologue HDO, and methane (CH 4 ), key gases in the water and carbon cycles. Previous versions (Version 4 or less) of the TES profile retrieval algorithm used a "spectral-window" approach to minimize uncertainty from interfering species at the expense of reduced vertical resolution and sensitivity. In this manuscript we document changes to the vertical resolution and uncertainties of the TES version 5 retrieval algorithm. In this version (Version 5), joint estimates of H 2 O, HDO, CH 4 and nitrous oxide (N 2 O) are made using radiances from almost the entire spectral region between 1100 cm −1 and 1330 cm −1 . The TES retrieval constraints are also modified in order to better use this information. The new H 2 O estimates show improved vertical resolution in the lower troposphere and boundary layer, while the new HDO/H 2 O estimates can now profile the HDO/H 2 O ratio between 925 hPa and 450 hPa in the tropics and during summertime at high latitudes. The new retrievals are now sensitive to methane in the free troposphere between 800 and 150 mb with peak sensitivity near 500 hPa; whereas in previous versions the sensitivity peaked at 200 hPa. However, the upper troposphere methane concentrations are biased high relative to the lower troposphere by approximately 4% on average. This bias is likely related to temperature, calibration, and/or methane spectroscopy errors. This bias can be mitigated by normalizing the CH 4 estimate by the ratio of the N 2 O estimate relative to the N 2 O prior, under the assumption that the same systematic error affects both the N 2 O and CH 4 estimates. We demonstrate that applying this ratio theoretically reduces the CH 4 estimate for non-retrieved parameters that jointly affect both the N 2 O and CH 4 estimates. The relative upper troposphere to lower troposphere bias is approximately 2.8% after this bias correction. Quality flags based upon the vertical variability of the methane and N 2 O estimates can be used to reduce this bias further. While these new CH 4 , HDO/H 2 O, and H 2 O estimates are consistent with previous TES retrievals in the altitude regions where the sensitivities overlap, future comparisons with independent profile measurement will be required to characterize the biases of these new retrievals and determine if the calculated uncertainties using the new constraints are consistent with actual uncertainties.

Description: Natural agricultural production functions in open system, where energy and matter are exchanged freely between geosphere (especially pedosphere), biosphere, and atmosphere. The self regulation of biophysical processes in these spheres is crucial for global sustainability. However, modern farming practices have stressed the system to the extent that damaged C, N, and P cycles are threatening catastrophic consequences. High food-wastage coupled with high food price, high soil fertility zones but low crop productivity, and highly productive irrigated farming marred with ecological disasters are examples of global paradoxes associated with modern farming. High incidence of pesticide residues in soil and water bodies, low use efficiency of agricultural inputs, vulnerability to climate, and low conversion of energy by crops leaving high amount of residues and their subsequent burden on environment, and hunger and malnutrition in many parts of the world are threatening civilizations to crumble. Nanotechnology promises to break these vicious cycles, because technology is based on applying exact amounts of inputs for use by the crops, and only when they are required. Some of the examples are: nanofabricating nutrient ions to improve nutrient use efficiency, targeted use of nano-pesticides, holding of water and then releasing it at the time of crop need by nps, forecasting and elimination of diseases, packaging of food to improve shelf-life, protecting food during storage, DNA nanotechnology, smart treatment delivery systems, bioanalytical nanosensors, bioselective surfaces, nanobioprocessing, protection of the environment by the reduction and conversion of agricultural materials into valuable products, design and development of new nanocatalysts to convert vegetable oils into biobased fuels and biodegradable industrial solvents, and controlled ecological life support system.

Description: Two instruments based on differential optical absorption spectroscopy (DOAS) to measure accurate ammonia concentrations in the atmosphere Atmospheric Measurement Techniques, 5, 413-427, 2012 Author(s): H. Volten, J. B. Bergwerff, M. Haaima, D. E. Lolkema, A. J. C. Berkhout, G. R. van der Hoff, C. J. M. Potma, R. J. Wichink Kruit, W. A. J. van Pul, and D. P. J. Swart We present two Differential Optical Absorption Spectroscopy (DOAS) instruments built at RIVM: the RIVM DOAS and the miniDOAS. Both instruments provide virtually interference-free measurements of NH 3 concentrations in the atmosphere, since they measure over an open path, without suffering from inlet problems or interference problems by ammonium aerosols dissociating on tubes or filters. They measure concentrations up to at least 200 μg m −3 , have a fast response, low maintenance demands, and a high up-time. The RIVM DOAS has a high accuracy of typically 0.15 μg m −3 for ammonia for 5-min averages and over a total light path of 100 m. The miniDOAS has been developed for application in measurement networks such as the Dutch National Air Quality Monitoring Network (LML). Compared to the RIVM DOAS it has a similar accuracy, but is significantly reduced in size, costs, and handling complexity. The RIVM DOAS and miniDOAS results showed excellent agreement ( R 2 = 0.996) during a field measurement campaign in Vredepeel, the Netherlands. This measurement site is located in an agricultural area and is characterized by highly variable, but on average high ammonia concentrations in the air. The RIVM-DOAS and miniDOAS results were compared to the results of the AMOR instrument, a continuous-flow wet denuder system, which is currently used in the LML. Averaged over longer time spans of typically a day, the (mini)DOAS and AMOR results agree reasonably well, although an offset of the AMOR values compared to the (mini)DOAS results exists. On short time scales, the (mini)DOAS shows a faster response and does not show the memory effects due to inlet tubing and transport of absorption fluids encountered by the AMOR. Due to its high accuracy, high uptime, low maintenance and its open path, the (mini)DOAS shows a good potential for flux measurements by using two (or more) systems in a gradient set-up and applying the aerodynamic gradient technique.

Description: Importance of rare taxa for bacterial diversity in the rhizosphere of Bt- and conventional maize varieties The ISME Journal 7, 37 (January 2013). doi:10.1038/ismej.2012.77 Authors: Anja B Dohrmann, Meike Küting, Sebastian Jünemann, Sebastian Jaenicke, Andreas Schlüter & Christoph C Tebbe

Description: Dispersal of thermophilic Desulfotomaculum endospores into Baltic Sea sediments over thousands of years The ISME Journal 7, 72 (January 2013). doi:10.1038/ismej.2012.83 Authors: Júlia Rosa de Rezende, Kasper Urup Kjeldsen, Casey R J Hubert, Kai Finster, Alexander Loy & Bo Barker Jørgensen

Description: Phylogenetic stratigraphy in the Guerrero Negro hypersaline microbial mat The ISME Journal 7, 50 (January 2013). doi:10.1038/ismej.2012.79 Authors: J Kirk Harris, J Gregory Caporaso, Jeffrey J Walker, John R Spear, Nicholas J Gold, Charles E Robertson, Philip Hugenholtz, Julia Goodrich, Daniel McDonald, Dan Knights, Paul Marshall, Henry Tufo, Rob Knight & Norman R Pace

Description: Metaproteogenomic insights beyond bacterial response to naphthalene exposure and bio-stimulation The ISME Journal 7, 122 (January 2013). doi:10.1038/ismej.2012.82 Authors: María-Eugenia Guazzaroni, Florian-Alexander Herbst, Iván Lores, Javier Tamames, Ana Isabel Peláez, Nieves López-Cortés, María Alcaide, Mercedes V Del Pozo, José María Vieites, Martin von Bergen, José Luis R Gallego, Rafael Bargiela, Arantxa López-López, Dietmar H Pieper, Ramón Rosselló-Móra, Jesús Sánchez, Jana Seifert & Manuel Ferrer

Description: Extracellular Ca2+ is a danger signal activating the NLRP3 inflammasome through G protein-coupled calcium sensing receptors Nature Communications 3, 1329 (2012). doi:10.1038/ncomms2339 Authors: Manuela Rossol, Matthias Pierer, Nora Raulien, Dagmar Quandt, Undine Meusch, Kathrin Rothe, Kristin Schubert, Torsten Schöneberg, Michael Schaefer, Ute Krügel, Sanela Smajilovic, Hans Bräuner-Osborne, Christoph Baerwald & Ulf Wagner

Description: Phospho-dependent ubiquitination and degradation of PAR-1 regulates synaptic morphology and tau-mediated Aβ toxicity in Drosophila Nature Communications 3, 1312 (2012). doi:10.1038/ncomms2278 Authors: Seongsoo Lee, Ji-Wu Wang, Wendou Yu & Bingwei Lu

Description: Atomic-scale engineering of magnetic anisotropy of nanostructures through interfaces and interlines Nature Communications 3, 1313 (2012). doi:10.1038/ncomms2316 Authors: S. Ouazi, S. Vlaic, S. Rusponi, G. Moulas, P. Buluschek, K. Halleux, S. Bornemann, S. Mankovsky, J. Minár, J.B. Staunton, H. Ebert & H. Brune

Description: Organic light-emitting transistors (OLETs) are of great research interest because they combine the advantage of the active channel of a transistor that can control the luminescence of an in-situ light-emitting diode in the same device. Here we report a novel single-crystal OLET (SCLET) that is coupled with single crystal optical feedback resonators. The combination of single-crystal waveguides with native Fabry-Perot cavities, formed by parallel crystal edges, drastically lowers the threshold energy for spectral narrowing and non-linear intensity enhancement. We apply this structure to SCLETs and demonstrate the first fabrication of a SCLET with the optical feedback resonators. Scientific Reports 2 doi: 10.1038/srep00985

Description: Assembly of nanowires into ordered macroscopic structures with new functionalities has been a recent focus. In this Letter, we report a new route for ordering hydrophilic Ag nanowires with high aspect ratio by flowing through a glass capillary. The present glass capillary with well-defined silver nanowire films inside can serve as a portable and reusable substrate for surface-enhanced Raman spectroscopy (SERS), which may provide a versatile and promising platform for detecting mixture pollutions. By controlling the flow parameters of nanowire suspensions, initially random Ag nanowires can be aligned to form nanowire arrays with tunable density, forming cambered nanowire films adhered onto the inner wall of the capillary. Compared with the planar ordered Ag nanowire films by the Langmuir-Blodgett (LB) technique, the cambered nanowire films show better SERS performance. Scientific Reports 2 doi: 10.1038/srep00987

Description: Acervuli are calcified concretions in the pineal gland (PG). Particularly interesting are their incidence and size, which are believed to affect neurological disorders and many physiological functions of PG such as regulating circadian rhythm. Despite long investigations for a century, detailed growth mechanism of acervuli has yet to be studied. Here we study the growth morphology of acervuli in human PGs by a direct visualization in 3-dimension (3-D) using a synchrotron X-ray imaging method. For an entire PG, non-aggregated acervuli show Gaussian distribution in size with 47±28 µm. The 3-D volume rendered images of acervuli reveal that the bumpy surfaces developed by lamination result in the mulberry-like structure. In addition, coalescence of multiple acervuli leads to large-scale lamination on the whole aggregate. We suggest a novel hypothesis on the growth patterns of acervuli by their nucleation density (Nd): i) mulberry-like structure at low Nd, and ii) large-scale lamination on an aggregate at high Nd. Scientific Reports 2 doi: 10.1038/srep00984

Description: The density maximum of water dominates the thermodynamics of the system under ambient conditions, is strongly P-dependent, and disappears at a crossover pressure Pcross ~ 1.8 kbar. We study this variable across a wide area of the T–P phase diagram. We consider old and new data of both the isothermal compressibility KT(T, P) and the coefficient of thermal expansion αP(T, P). We observe that KT(T) shows a minimum at T* ~ 315±5 K for all the studied pressures. We find the behavior of αP to also be surprising: all the αP(T) curves measured at different P cross at T*. The experimental data show a “singular and universal expansivity point” at T* ~ 315 K and αP(T*) ≃ 0.44 10−3 K−1. Unlike other water singularities, we find this temperature to be thermodynamically consistent in the relationship connecting the two response functions. Scientific Reports 2 doi: 10.1038/srep00993

Description: MicroRNAs (miRNAs) are key regulators of sequence-specific gene silencing. However, crucial factors that determine the efficacy of miRNA-mediated target gene silencing are poorly understood. Here we mathematized base-pairing stability and showed that miRNAs with an unstable 5′ terminal duplex and stable seed-target duplex exhibit strong silencing activity. The results are consistent with the previous findings that an RNA strand with unstable 5′ terminal in miRNA duplex easily loads onto the RNA-induced silencing complex (RISC), and miRNA recognizes target mRNAs with seed-complementary sequences to direct posttranscriptional repression. Our results suggested that both the unwinding and target recognition processes of miRNAs could be proficiently controlled by the thermodynamics of base-pairing in protein-free condition. Interestingly, such thermodynamic parameters might be evolutionarily well adapted to the body temperatures of various species. Scientific Reports 2 doi: 10.1038/srep00996

Description: A spin liquid is a novel quantum state of matter with no conventional order parameter where a finite charge gap exists even though the band theory would predict metallic behavior. Finding a stable spin liquid in two or higher spatial dimensions is one of the most challenging and debated issues in condensed matter physics. Very recently, it has been reported that a model of graphene, i.e., the Hubbard model on the honeycomb lattice, can show a spin liquid ground state in a wide region of the phase diagram, between a semi-metal (SM) and an antiferromagnetic insulator (AFMI). Here, by performing numerically exact quantum Monte Carlo simulations, we extend the previous study to much larger clusters (containing up to 2592 sites), and find, if any, a very weak evidence of this spin liquid region. Instead, our calculations strongly indicate a direct and continuous quantum phase transition between SM and AFMI. Scientific Reports 2 doi: 10.1038/srep00992

Description: Snow specific surface area simulation using the one-layer snow model in the Canadian LAnd Surface Scheme (CLASS) The Cryosphere Discussions, 6, 5255-5289, 2012 Author(s): A. Roy, A. Royer, B. Montpetit, P. A. Bartlett, and A. Langlois Snow grain size is a key parameter for modeling microwave snow emission properties and the surface energy balance because of its influence on the snow albedo, thermal conductivity and diffusivity. A model of the specific surface area (SSA) of snow was implemented in the one-layer snow model in the Canadian LAnd Surface Scheme (CLASS) version 3.4. This offline multilayer model (CLASS-SSA) simulates the decrease of SSA based on snow age, snow temperature and the temperature gradient under dry snow conditions, whereas it considers the liquid water content for wet snow metamorphism. We compare the model with ground-based measurements from several sites (alpine, Arctic and sub-Arctic) with different types of snow. The model provides simulated SSA in good agreement with measurements with an overall point-to-point comparison RMSE of 8.1 m 2 kg −1 , and a RMSE of 4.9 m 2 kg −1 for the snowpack average SSA. The model, however, is limited under wet conditions due to the single-layer nature of the CLASS model, leading to a single liquid water content value for the whole snowpack. The SSA simulations are of great interest for satellite passive microwave brightness temperature assimilations, snow mass balance retrievals and surface energy balance calculations with associated climate feedbacks.

Description: Effects of fire on ash thickness in a Lithuanian grassland and short-term spatio-temporal changes Solid Earth Discussions, 4, 1545-1584, 2012 Author(s): P. Pereira, A. Cerdà, X. Úbeda, J. Mataix-Solera, D. Martin, A. Jordán, and M. Burguet Ash thickness is a key variable in the protection of soil against erosion agents after planned and unplanned fires. Thicker ash provides better protection against raindrop impact and reduces the runoff response by retaining water and promoting water infiltration although little is known about the distribution and the evolution of the ash layer after the fires. Ash thickness measurements were conducted along two transects (flat and sloping areas) following a a grid experimental design. Both transects extended from the burned area into an adjacent unburned area. We analysed ash thickness evolution according to time and fire severity. In order to interpolate data with accuracy and identify the techniques with the least bias, several interpolation methods were tested in the grid plot. Overall, the fire had a low severity. The fire significantly reduced the ground cover, especially on sloping areas owing to the higher fire severity and/or less biomass previous to the fire. Ash thickness depends on fire severity and is thin where fire severity was higher and thicker in lower fire severity sites. The ash thickness decreased with time after the fire. Between 4 and 16 days after the fire, ash was transported by wind. The major reduction took place between 16 and 34 days after the fire as a result of rainfall, and was more efficient where fire severity was higher. Between 34 and 45 days after the fire no significant differences in ash thickness were identified among ash colours and only traces of the ash layer remained. The omni-directional experimental variograms shown that variable structure did not change importantly with the time, however, the most accurate interpolation methods were different highlighting the slight different patterns of ash thickness distribution with the time. The ash spatial variability increased with the time, particularly on the slope, as a result of water erosion.

Description: Evaluation of a near-global eddy-resolving ocean model Geoscientific Model Development Discussions, 5, 4305-4354, 2012 Author(s): P. R. Oke, D. A. Griffin, A. Schiller, R. J. Matear, R. Fiedler, J. Mansbridge, A. Lenton, M. Cahill, M. A. Chamberlain, and K. Ridgway Analysis of the variability in an 18-yr run of a near-global, eddy-resolving ocean general circulation model coupled with biogeochemistry is presented. Comparisons between modelled and observed mean sea level (MSL), mixed-layer depth (MLD), sea-level anomaly (SLA), sea-surface temperature (SST), and Chlorophyll a indicate that the model variability is realistic. We find some systematic errors in the modelled MLD, with the model generally deeper than observations, that results in errors in the Chlorophyll a , owing to the strong biophysical coupling. We evaluate several other metrics in the model, including the zonally-averaged seasonal cycle of SST, meridional overturning, volume transports through key Straits and passages, zonal averaged temperature and salinity, and El Nino-related SST indices. We find that the modelled seasonal cycle in SST is 0.5–1.5 °C weaker than observed; volume transports of the Antarctic Circumpolar Current, the East Australian Current, and Indonesian Throughflow are in good agreement with observational estimates; and the correlation between the modelled and observed NINO SST indices exceed 0.91. Most aspects of the model circulation are realistic. We conclude that the model output is suitable for broader analysis to better understand ocean dynamics and ocean variability.

Description: Elusive drought: uncertainty in observed trends and short- and long-term CMIP5 projections Hydrology and Earth System Sciences Discussions, 9, 13773-13803, 2012 Author(s): B. Orlowsky and S. I. Seneviratne Recent years have seen a number of severe droughts in different regions around the world, causing agricultural and economic losses, famines and migration. Despite their devastating consequences, the Standardised Precipitation Index (SPI) of these events lies within the range of internal climate variability, which we estimate from simulations from the 5th phase of the Coupled Model Intercomparison Project (CMIP5). In terms of drought magnitude, regional trends of SPI over the last decades remain mostly inconclusive in observations and CMIP5 simulations, although Soil Moisture Anomalies (SMAs) in CMIP5 simulations hint at increased drought in a few regions (e.g. the Mediterranean, Central America/Mexico, the Amazon, North-East Brazil and South Africa). Also for the future, projections of meteorological (SPI) and agricultural (SMA) drought in CMIP5 display large uncertainties over all time frames, generally impeding trend detection. Analogue analyses of the frequencies rather than magnitudes of future drought display, however, more robust signal-to-noise ratios with detectable trends towards more frequent drought until the end of the 21st century in the Mediterranean, South Africa and Central America/Mexico. Other present-day hot spots are projected to become less drought-prone, or to display unsignificant changes in drought occurrence. A separation of different sources of uncertainty in drought projections reveals that for the near term, internal climate variability is the dominant source, while the formulation of Global Climate Models (GCMs) generally becomes the dominant source of uncertainty by the end of the 21st century, especially for agricultural (soil moisture) drought. In comparison, the uncertainty in Green-House Gas (GHG) concentrations scenarios is negligible for most regions. These findings stand in contrast to respective analyses for a heat wave indicator, for which GHG concentrations scenarios constitute the main source of uncertainty. Our results highlight the inherent difficulty of drought quantification and the uncertainty of drought projections. However, high uncertainty should not be equated with low drought risk, since potential scenarios include large drought increases in key agricultural and ecosystem regions.

Description: Circulation changes in the winter lower atmosphere and long-lasting solar/geomagnetic activity Annales Geophysicae, 30, 1719-1726, 2012 Author(s): J. Bochníček, H. Davídkovová, P. Hejda, and R. Huth The paper describes the association between high long-lasting solar/geomagnetic activity and geopotential height (GPH) changes in the winter lower atmosphere, based on their development in the Northern Hemisphere in the winter periods (December–March) of 1950–1969 and 1970–2002. Solar/geomagnetic activity is characterised by the 60-day mean of the sunspot number R/by the 60-day mean of the daily sum of the Kp index. The GPH distributions in the lower atmosphere are described by 60-day anomalies from their long-term daily average at 20 hPa/850 hPa. The data have been adopted from the NCEP/NCAR reanalysis. The 60-day mean values of solar/geomagnetic activity and GPH anomalies were calculated in five-day steps over the whole winter period. The analysis was carried out using composite maps which represent their distribution of the GPH anomalies during high solar activity ( R ≥ 100) and high geomagnetic activity (ΣKp ≥ 20). Analysis has shown that the distribution of GPH anomalies depends on solar activity, geomagnetic activity and the phase of winter period (early or late winter). The nature of this relationship then depends on the time interval involved, i.e. 1950–1969 or 1970–2002. Positive anomalies in the polar stratosphere (20 hPa) were detected during the whole winter periods of the years 1950–1969. Significant anomalies were detected in the lower troposphere (850 hPa) during the second half of the winter period. The distribution of GPH anomalies on the maps compiled with regard to solar activity was similar to the distribution on maps compiled with regard to geomagnetic activity. In the interval 1970–2002, significant negative GPH anomalies were detected in the stratosphere at high latitudes, and positive anomalies were detected in the region of low latitudes. The distribution of GPH anomalies in the lower troposphere was substantially affected by situations in which, together with high solar activity, also high geomagnetic activity occurred.

Description: Sensitivity studies of dust ice nuclei effect on cirrus clouds with the Community Atmosphere Model CAM5 Atmospheric Chemistry and Physics, 12, 12061-12079, 2012 Author(s): X. Liu, X. Shi, K. Zhang, E. J. Jensen, A. Gettelman, D. Barahona, A. Nenes, and P. Lawson In this study the effect of dust aerosol on upper tropospheric cirrus clouds through heterogeneous ice nucleation is investigated in the Community Atmospheric Model version 5 (CAM5) with two ice nucleation parameterizations. Both parameterizations consider homogeneous and heterogeneous nucleation and the competition between the two mechanisms in cirrus clouds, but differ significantly in the number concentration of heterogeneous ice nuclei (IN) from dust. Heterogeneous nucleation on dust aerosol reduces the occurrence frequency of homogeneous nucleation and thus the ice crystal number concentration in the Northern Hemisphere (NH) cirrus clouds compared to simulations with pure homogeneous nucleation. Global and annual mean shortwave and longwave cloud forcing are reduced by up to 2.0 ± 0.1 W m −2 (1σ uncertainty) and 2.4 ± 0.1 W m −2 , respectively due to the presence of dust IN, with the net cloud forcing change of −0.40 ± 0.20 W m −2 . Comparison of model simulations with in situ aircraft data obtained in NH mid-latitudes suggests that homogeneous ice nucleation may play an important role in the ice nucleation at these regions with temperatures of 205–230 K. However, simulations overestimate observed ice crystal number concentrations in the tropical tropopause regions with temperatures of 190–205 K, and overestimate the frequency of occurrence of high ice crystal number concentration (〉 200 L −1 ) and underestimate the frequency of low ice crystal number concentration ( 〈 30 L −1 ) at NH mid-latitudes. These results highlight the importance of quantifying the number concentrations and properties of heterogeneous IN (including dust aerosol) in the upper troposphere from the global perspective.

Description: Dissolved organic carbon (DOC) and select aldehydes in cloud and fog water: the role of the aqueous phase in impacting trace gas budgets Atmospheric Chemistry and Physics Discussions, 12, 33083-33125, 2012 Author(s): B. Ervens, Y. Wang, J. Eagar, W. R. Leaitch, A. M. Macdonald, K. T. Valsaraj, and P. Herckes Cloud and fog droplets efficiently scavenge and process water-soluble compounds and thus modify the chemical composition of the gas and particle phases. The concentrations of dissolved organic carbon (DOC) in the aqueous phase reach concentrations on the order of ~10 mg C L −1 which is typically on the same order of magnitude as the sum of inorganic anions. Aldehydes and carboxylic acids typically comprise a large fraction of DOC because of their high solubility. The dissolution of species in the aqueous phase can lead to (i) the removal of species from the gas phase preventing their processing by gas phase reactions (e.g. photolysis of aldehydes) and (ii) the formation of unique products that do not have any efficient gas phase sources (e.g. dicarboxylic acids). We present measurements of DOC and select aldehydes in fog water at high elevation and intercepted clouds in a biogenically-impacted location (Whistler, Canada) and in fog water in a more polluted area (Davis, CA). Concentrations of formaldehyde, glyoxal and methylglyoxal were in the micromolar range and comprised ≤2% each individually of the DOC. Comparison of the DOC and aldehyde concentrations to those at other locations shows good agreement and reveals highest levels for both in anthropogenically impacted regions. Based on this overview, we conclude that the fraction of organic carbon (dissolved and insoluble inclusions) in the aqueous phase comprises 1–~40% of total organic carbon. Higher values are observed to be associated with aged air masses where organics are expected to be more highly oxidized and thus more soluble. Accordingly, the aqueous/gas partitioning ratio expressed here as an effective Henry's law constant for DOC ( K H *DOC ) increases by an order of magnitude from 7×10 3 M atm −1 to 7×10 4 M atm −1 during the ageing of air masses. The measurements are accompanied by photochemical box model simulations. They suggest that the scavenging of aldehydes by the aqueous phase can reduce HO 2 gas phase levels by two orders of magnitude due to a weaker net source of HO 2 production from aldehyde photolysis in the gas phase. Despite the high solubility of dialdehydes (glyoxal, methylglyoxal), their impact on the HO 2 budget by scavenging is

Description: Climate versus emission drivers of methane lifetime against loss by tropospheric OH from 1860–2100 Atmospheric Chemistry and Physics, 12, 12021-12036, 2012 Author(s): J. G. John, A. M. Fiore, V. Naik, L. W. Horowitz, and J. P. Dunne With a more-than-doubling in the atmospheric abundance of the potent greenhouse gas methane (CH 4 ) since preindustrial times, and indications of renewed growth following a leveling off in recent years, questions arise as to future trends and resulting climate and public health impacts from continued growth without mitigation. Changes in atmospheric methane lifetime are determined by factors which regulate the abundance of OH, the primary methane removal mechanism, including changes in CH 4 itself. We investigate the role of emissions of short-lived species and climate in determining the evolution of methane lifetime against loss by tropospheric OH, (τ CH4_OH ), in a suite of historical (1860–2005) and future Representative Concentration Pathway (RCP) simulations (2006–2100), conducted with the Geophysical Fluid Dynamics Laboratory (GFDL) fully coupled chemistry-climate model (CM3). From preindustrial to present, CM3 simulates an overall 5% increase in τ CH4_OH due to a doubling of the methane burden which offsets coincident increases in nitrogen oxide (NO x emissions. Over the last two decades, however, the τ CH4_OH declines steadily, coinciding with the most rapid climate warming and observed slow-down in CH 4 growth rates, reflecting a possible negative feedback through the CH 4 sink. Sensitivity simulations with CM3 suggest that the aerosol indirect effect (aerosol-cloud interactions) plays a significant role in cooling the CM3 climate. The projected decline in aerosols under all RCPs contributes to climate warming over the 21st century, which influences the future evolution of OH concentration and τ CH4_OH . Projected changes in τ CH4_OH from 2006 to 2100 range from −13% to +4%. The only projected increase occurs in the most extreme warming case (RCP8.5) due to the near-doubling of the CH 4 abundance, reflecting a positive feedback on the climate system. The largest decrease occurs in the RCP4.5 scenario due to changes in short-lived climate forcing agents which reinforce climate warming and enhance OH. This decrease is more-than-halved in a sensitivity simulation in which only well-mixed greenhouse gas radiative forcing changes along the RCP4.5 scenario (5% vs. 13%).

Description: Limitations of wind extraction from 4-D-Var assimilation of trace gases Atmospheric Chemistry and Physics Discussions, 12, 32985-33023, 2012 Author(s): D. R. Allen, K. W. Hoppel, G. E. Nedoluha, D. D. Kuhl, N. L. Baker, L. Xu, and T. E. Rosmond Time-dependent variational data assimilation allows the possibility of extracting wind information from observations of long-lived trace gases. Since trace gas observations are not available at sufficient resolution for deriving feature-track winds, they must be combined with model background information to produce an analysis. If done with time-dependent variational assimilation, wind information may be extracted via the adjoint of the linearized tracer continuity equation. This paper presents idealized experiments that illustrate the mechanics of tracer-wind extraction and demonstrate some of the limitations of this procedure. We first examine tracer-wind extraction using a simple one-dimensional advection equation. The analytic solution for a single trace gas observation is discussed along with numerical solutions for multiple observations. The limitations of tracer-wind extraction are then explored using highly idealized ozone experiments performed with a development version of the Navy Global Environmental Model (NAVGEM) in which stratospheric globally-distributed hourly stratospheric ozone profiles are assimilated in a single 6-h update cycle in January 2009. Starting with perfect background ozone conditions, but imperfect dynamical conditions, ozone errors develop over the 6-h background window. Wind increments are introduced in the analysis in order to reduce the differences between background ozone and ozone observations. For "perfect" observations (unbiased and no random error), this results in root mean square (RMS) vector wind error reductions of up to ∼ 3 m s −1 in the winter hemisphere and tropics. Wind extraction is more difficult in the summer hemisphere due to weak ozone gradients and smaller background wind errors. The limitations of wind extraction are also explored for observations with imposed random errors and for limited sampling patterns. As expected, the amount of wind information extracted degrades as observation errors or data voids increase. In the case of poorly specified observation error covariances, assimilation of ozone data with imposed errors may result in erroneous wind increments, since the assimilation is constrained too tightly to the noisy observations.

Description: Mass deposition fluxes of Saharan mineral dust to the tropical northeast Atlantic Ocean: an intercomparison of methods Atmospheric Chemistry and Physics Discussions, 12, 33025-33081, 2012 Author(s): N. Niedermeier, A. Held, T. Müller, B. Heinold, K. Schepanski, I. Tegen, K. Kandler, M. Ebert, S. Weinbruch, K. Read, J. Lee, K. W. Fomba, K. Müller, H. Herrmann, and A. Wiedensohler The aim of this study is to determine the mass deposition flux of mineral dust to the tropical northeast Atlantic Ocean at the Cape Verde Atmospheric Observatory (CVAO) on the island Sao Vicente for January 2009. Five different methods were applied to estimate the deposition flux, using different meteorological and microphysical measurements, remote sensing, and regional dust transport simulations. The set of observations comprises micrometeorological measurements with an ultra-sonic anemometer and profile measurements using 2-D anemometers at two different heights, and microphysical measurements of the size-resolved mass concentrations of mineral dust. In addition, the total mass concentration of mineral dust was derived from absorption photometer observations and passive sampling. The regional dust model COSMO-MUSCAT was used for simulations of dust emission and transport, including dry and wet deposition processes. The four observation-based methods yield a monthly average deposition flux of mineral dust of 12–29 ng m −2 s −1 . The simulation results come close to the upper range of the measurements with an average value of 47 ng m −2 s −1 . It is shown, that the mass deposition flux of mineral dust obtained by the combination of micrometeorological (ultra-sonic anemometer) and microphysical measurements (particle mass size distribution of mineral dust) is within 5% to modeled mass deposition fluxes when the mineral dust is relatively homogenously distributed over the investigated area.

Description: Cloud condensation nuclei production associated with atmospheric nucleation: a synthesis based on existing literature and new results Atmospheric Chemistry and Physics, 12, 12037-12059, 2012 Author(s): V.-M. Kerminen, M. Paramonov, T. Anttila, I. Riipinen, C. Fountoukis, H. Korhonen, E. Asmi, L. Laakso, H. Lihavainen, E. Swietlicki, B. Svenningsson, A. Asmi, S. N. Pandis, M. Kulmala, and T. Petäjä This paper synthesizes the available scientific information connecting atmospheric nucleation with subsequent cloud condensation nuclei (CCN) formation. We review both observations and model studies related to this topic, and discuss the potential climatic implications. We conclude that CCN production associated with atmospheric nucleation is both frequent and widespread phenomenon in many types of continental boundary layers, and probably also over a large fraction of the free troposphere. The contribution of nucleation to the global CCN budget spans a relatively large uncertainty range, which, together with our poor understanding of aerosol-cloud interactions, results in major uncertainties in the radiative forcing by atmospheric aerosols. In order to better quantify the role of atmospheric nucleation in CCN formation and Earth System behavior, more information is needed on (i) the factors controlling atmospheric CCN production and (ii) the properties of both primary and secondary CCN and their interconnections. In future investigations, more emphasis should be put on combining field measurements with regional and large-scale model studies.

Description: Intercomparison of shortwave radiative transfer schemes in global aerosol modeling: results from the AeroCom Radiative Transfer Experiment Atmospheric Chemistry and Physics Discussions, 12, 32631-32706, 2012 Author(s): C. A. Randles, S. Kinne, G. Myhre, M. Schulz, P. Stier, J. Fischer, L. Doppler, E. Highwood, C. Ryder, B. Harris, J. Huttunen, Y. Ma, R. T. Pinker, B. Mayer, D. Neubauer, R. Hitzenberger, L. Oreopoulos, D. Lee, G. Pitari, G. Di Genova, J. Quaas, Fred G. Rose, S. Kato, S. T. Rumbold, I. Vardavas, N. Hatzianastassiou, C. Matsoukas, H. Yu, F. Zhang, H. Zhang, and P. Lu In this study we examine the performance of 31 global model radiative transfer schemes in cloud-free conditions with prescribed gaseous absorbers and no aerosols (Rayleigh atmosphere), with prescribed scattering-only aerosols, and with more absorbing aerosols. Results are compared to benchmark results from high-resolution, multi-angular line-by-line radiation models. For purely scattering aerosols, model bias relative to the line-by-line models in the top-of-the atmosphere aerosol radiative forcing ranges from roughly −10 to 20%, with over- and underestimates of radiative cooling at higher and lower sun elevation, respectively. Inter-model diversity (relative standard deviation) increases from ~10 to 15% as sun elevation increases. Inter-model diversity in atmospheric and surface forcing decreases with increased aerosol absorption, indicating that the treatment of multiple-scattering is more variable than aerosol absorption in the models considered. Aerosol radiative forcing results from multi-stream models are generally in better agreement with the line-by-line results than the simpler two-stream schemes. Considering radiative fluxes, model performance is generally the same or slightly better than results from previous radiation scheme intercomparisons. However, the inter-model diversity in aerosol radiative forcing remains large, primarily as a result of the treatment of multiple-scattering. Results indicate that global models that estimate aerosol radiative forcing with two-stream radiation schemes may be subject to persistent biases introduced by these schemes, particularly for regional aerosol forcing.

Description: Immersion freezing of birch pollen washing water Atmospheric Chemistry and Physics Discussions, 12, 32911-32943, 2012 Author(s): S. Augustin, S. Hartmann, B. Pummer, H. Grothe, D. Niedermeier, T. Clauss, J. Voigtländer, L. Tomsche, H. Wex, and F. Stratmann In the present study, the immersion freezing behavior of birch pollen, i.e. its ice nucleating active (INA) macromolecules, was investigated at the Leipzig Aerosol Cloud Interaction Simulator (LACIS). For that, washing water of two different birch pollen samples with different regional origin (Northern birch and Southern birch) were used. The immersion freezing of droplets generated from the pollen washing water was already observed at temperatures higher than −20 °C, for both samples. Main differences between the Northern birch pollen and the Southern birch pollen were obvious in a temperature range, between −18 °C and −24 °C, where the ice fraction increased with decreasing temperature. There, the Northern birch pollen washing water featured two different slopes, with one being steeper and one being similar to the slope of the Southern birch pollen washing water. As we assume single INA macromolecules being the reason for the ice nucleation, we concluded that the Northern birch pollen are able to produce at least two different types of INA macromolecules. We were able to determine the heterogeneous nucleation rates for both INA macromolecule types and so could explain the ice nucleation behavior of both, the Southern and the Northern birch pollen washing water.

Description: Effect of sampling variation on error of rainfall variables measured by optical disdrometer Atmospheric Measurement Techniques Discussions, 5, 8895-8924, 2012 Author(s): X. C. Liu, T. C. Gao, and L. Liu During the sampling process of precipitation particles by optical disdrometers, the randomness of particles and sampling variability has great impact on the accuracy of precipitation variables. Based on a marked point model of raindrop size distribution, the effect of sampling variation on drop size distribution and velocity distribution measurement using optical disdrometers are analyzed by Monte Carlo simulation. The results show that the samples number, rain rate, drop size distribution, and sampling size have different influences on the accuracy of rainfall variables. The relative errors of rainfall variables caused by sampling variation in a descending order as: water concentration, mean diameter, mass weighed mean diameter, mean volume diameter, radar reflectivity factor, and number density, which are independent with samples number basically; the relative error of rain variables are positively correlated with the margin probability, which is also positively correlated with the rain rate and the mean diameter of raindrops; the sampling size is one of the main factors that influence the margin probability, with the decreasing of sampling area, especially the decreasing of short side of sample size, the probability of margin raindrops is getting greater, hence the error of rain variables are getting greater, and the variables of median size raindrops have the maximum error. To ensure the relative error of rainfall variables measured by optical disdrometer less than 1%, the width of light beam should be at least 40 mm.

Description: Evidence and quantitation of aromatic organosulfates in ambient aerosols in Lahore, Pakistan Atmospheric Chemistry and Physics Discussions, 12, 32795-32823, 2012 Author(s): S. Kundu, T. A. Quraishi, G. Yu, C. Suarez, F. N. Keutsch, and E. A. Stone Organosulfates are important components of atmospheric organic aerosols, yet their structures, abundances, sources and formation processes are not adequately understood. This study presents the identification and quantitation of benzyl sulfate in atmospheric aerosols, which is the first reported atmospheric organosulfate with aromatic carbon backbone. Benzyl sulfate was identified and quantified in fine particulate matter (PM 2.5 ) collected in Lahore, Pakistan during 2007–2008. An authentic standard of benzyl sulfate was synthesized, standardized, and identified in atmospheric aerosols using ultra-performance liquid chromatography (UPLC) coupled with quadrupole time-of-flight (Q-ToF) mass spectrometry (MS). Benzyl sulfate was quantified in aerosol samples using UPLC coupled to negative electrospray ionization triple quadrupole (TQ) MS. The highest benzyl sulfate concentrations were recorded in November and January 2007 (0.50 ± 0.11 ng m −3 ) whereas the lowest concentration was observed in July (0.05 ± 0.02 ng m −3 ). To evaluate matrix effects, benzyl sulfate concentrations were determined using external calibration and the method of standard addition; comparable concentrations were detected by the two methods, which ruled out significant matrix effects in benzyl sulfate quantitation. Three additional organosulfates with m / z 187, 201 and 215 were qualitatively identified as aromatic organosulfates with additional methyl substituents by high-resolution mass measurements and tandem MS. The observed aromatic organosulfates form a homologous series analogous to toluene, xylene, and trimethylbenzene, which are abundant anthropogenic volatile organic compounds (VOC), suggesting that aromatic organosulfates may be formed by secondary reactions. Further studies are needed to elucidate the sources and formation pathways of aromatic organosulfates in the atmosphere.