ALBERT

Description: On the origin of the occasional spring nitrate peak in Greenland snow Atmospheric Chemistry and Physics, 14, 13361-13376, 2014 Author(s): L. Geng, J. Cole-Dai, B. Alexander, J. Erbland, J. Savarino, A. J. Schauer, E. J. Steig, P. Lin, Q. Fu, and M. C. Zatko Ice core nitrate concentrations peak in the summer in both Greenland and Antarctica. Two nitrate concentration peaks in one annual layer have been observed some years in ice cores in Greenland from samples dating post-1900, with the additional nitrate peak occurring in the spring. The origin of the spring nitrate peak was hypothesized to be pollution transport from the mid-latitudes in the industrial era. We performed a case study on the origin of a spring nitrate peak in 2005 measured from a snowpit at Summit, Greenland, covering 3 years of snow accumulation. The effect of long-range transport of nitrate on this spring peak was excluded by using sulfate as a pollution tracer. The isotopic composition of nitrate (δ 15 N, δ 18 O and Δ 17 O) combined with photochemical calculations suggest that the occurrence of this spring peak is linked to a significantly weakened stratospheric ozone (O 3 ) layer. The weakened O 3 layer resulted in elevated UVB (ultraviolet-B) radiation on the snow surface, where the production of OH and NO x from the photolysis of their precursors was enhanced. Elevated NO x and OH concentrations resulted in enhanced nitrate production mainly through the NO 2 + OH formation pathway, as indicated by decreases in δ 18 O and Δ 17 O of nitrate associated with the spring peak. We further examined the nitrate concentration record from a shallow ice core covering the period from 1772 to 2006 and found 19 years with double nitrate peaks after the 1950s. Out of these 19 years, 14 of the secondary nitrate peaks were accompanied by sulfate peaks, suggesting long-range transport of nitrate as their source. In the other 5 years, low springtime O 3 column density was observed, suggesting enhanced local production of nitrate as their source. The results suggest that, in addition to direct transport of nitrate from polluted regions, enhanced local photochemistry can also lead to a spring nitrate peak. The enhanced local photochemistry is probably associated with the interannual variability of O 3 column density in the Arctic, which leads to elevated surface UV radiation in some years. In this scenario, enhanced photochemistry caused increased local nitrate production under the condition of elevated local NO x abundance in the industrial era.

Description: Spatial and temporal variation of CO over Alberta using measurements from satellite, aircrafts, and ground stations Atmospheric Chemistry and Physics Discussions, 14, 31767-31812, 2014 Author(s): H. S. Marey, Z. Hashisho, L. Fu, and J. Gille Alberta is Canada's largest oil producer and its oil sand deposits comprise 30% of the world's oil reserves. The process of bitumen extraction and upgrading releases trace gases and aerosols to the atmosphere. In this study we present satellite-based analysis to explore, for the first time, various contributing factors that affect tropospheric carbon monoxide (CO) levels over Alberta. The multispectral product that uses both near-infrared (NIR) and the thermal-infrared (TIR) radiances for CO retrieval from the Measurements of Pollution in the Troposphere (MOPITT) are examined for the 12 year period from 2002–2013. Moderate Resolution Imaging Spectroradiometer (MODIS) thermal anomaly product from 2001 to 2013 is employed to investigate the seasonal and temporal variations of forest fires. Additionally, in situ CO measurements at industrial and urban sites are compared to satellite data. Furthermore, the available MOZAIC/IAGOS (Measurement of Ozone, Water Vapor, Carbon Monoxide, Nitrogen Oxide by Airbus In-Service Aircraft/In service Aircraft for Global Observing System) aircraft CO profiles (April 2009–December 2011) are used to validate MOPITT CO data. The climatological time curtain plot and spatial maps for CO over northern Alberta indicate the signatures of transported CO for two distinct biomass burning seasons, summer and spring. Distinct seasonal patterns of CO at the urban site s (Edmonton and Calgary cities) point to the strong influence of traffic. Meteorological parameters play an important role on the CO spatial distribution at various pressure levels. Northern Alberta shows stronger upward lifting motion which leads to larger CO total column values while the poor dispersion in central and south Alberta exacerbate s the surface CO pollution. Inter-annual variations of satellite data depict a slightly decreasing trend for both regions while the decline trend is more evident from ground observations, especially at the urban sites. MOPITT CO vertical averages and MOZAIC/IAGOS aircraft profiles were in good agreement within the standard deviation at all pressure levels. There is consistency between the time evolution of high CO episodes monitored by satellite and ground measurements and the fire frequency peak time which implies that biomass burning has affected the tropospheric CO distribution in northern Alberta. These findings have further demonstrated the potential use of MOPITT V5 multispectral (NIR+TIR) product for assessing a complicated surface process.

Description: Particulate emissions from residential wood combustion in Europe – revised estimates and an evaluation Atmospheric Chemistry and Physics Discussions, 14, 31719-31765, 2014 Author(s): H. A. C. Denier van der Gon, R. Bergström, C. Fountoukis, C. Johansson, S. N. Pandis, D. Simpson, and A. Visschedijk Currently residential wood combustion (RWC) is increasing in Europe because of rising fossil fuel prices but also due to climate change mitigation policies. However, especially in small-scale applications, RWC may cause high emissions of particulate matter (PM). Recently we have developed a new high-resolution (7 km × 7 km) anthropogenic carbonaceous aerosol emission inventory for Europe. The inventory indicated that about half of the total PM 2.5 emission in Europe is carbonaceous aerosol and identified RWC as the largest organic aerosol (OA) source in Europe. The inventory was partly based on national reported PM emissions. Use of this OA inventory as input for two Chemical Transport Models (CTMs), PMCAMx and EMEP MSC-W, revealed major underestimations of OA in winter time, especially for regions dominated by RWC. Interestingly, this was not universal but appeared to differ by country. In the present study we constructed a new bottom-up emission inventory for RWC accounting for the semi-volatile components of the emissions. The new RWC emissions are higher than those in the previous inventory by a factor of 2–3 but with substantial inter-country variation. The new emission inventory served as input for the CTMs and a substantially improved agreement between measured and predicted organic aerosol was found. The new RWC inventory improves the model calculated OA significantly. Comparisons to Scandinavian source apportionment studies also indicate substantial improvements in the modeled wood-burning component of OA. This suggests that primary organic aerosol emission inventories need to be revised to include the semi-volatile OA that is formed almost instantaneously due to cooling of the flue gas or exhaust. Since RWC is a key source of fine PM in Europe, a major revision of the emission estimates as proposed here is likely to influence source-receptor matrices and modelled source apportionment. Since usage of biofuels, such as wood, in small combustion units is a globally significant source, this insight may also dramatically change global estimates of organic aerosol emissions.

Description: Impact of planetary boundary layer turbulence on model climate and tracer transport Atmospheric Chemistry and Physics Discussions, 14, 31627-31674, 2014 Author(s): E. L. McGrath-Spangler, A. Molod, L. E. Ott, and S. Pawson Planetary boundary layer (PBL) processes are important for weather, climate, and tracer transport and concentration. One measure of the strength of these processes is the PBL depth. However, no single PBL depth definition exists and several studies have found that the estimated depth can vary substantially based on the definition used. In the Goddard Earth Observing System (GEOS-5) atmospheric general circulation model, the PBL depth is particularly important because it is used to calculate the turbulent length scale that is used in the estimation of turbulent mixing. This study analyzes the impact of using three different PBL depth definitions in this calculation. Two definitions are based on the scalar eddy diffusion coefficient and the third is based on the bulk Richardson number. Over land, the bulk Richardson number definition estimates shallower nocturnal PBLs than the other estimates while over water this definition generally produces deeper PBLs. The near surface wind velocity, temperature, and specific humidity responses to the change in turbulence are spatially and temporally heterogeneous, resulting in changes to tracer transport and concentrations. Near surface wind speed increases in the bulk Richardson number experiment cause Saharan dust increases on the order of 1 × 10 −4 kg m −2 downwind over the Atlantic Ocean. Carbon monoxide (CO) surface concentrations are modified over Africa during boreal summer, producing differences on the order of 20 ppb, due to the model's treatment of emissions from biomass burning. While differences in carbon dioxide (CO 2 ) are small in the time mean, instantaneous differences are on the order of 10 ppm and these are especially prevalent at high latitude during boreal winter. Understanding the sensitivity of trace gas and aerosol concentration estimates to PBL depth is important for studies seeking to calculate surface fluxes based on near-surface concentrations and to studies projecting future concentrations.

Description: Aerosol hygroscopicity and cloud condensation nuclei activity during the AC 3 Exp campaign: implications for cloud condensation nuclei parameterization Atmospheric Chemistry and Physics, 14, 13423-13437, 2014 Author(s): F. Zhang, Y. Li, Z. Li, L. Sun, R. Li, C. Zhao, P. Wang, Y. Sun, X. Liu, J. Li, P. Li, G. Ren, and T. Fan Aerosol hygroscopicity and cloud condensation nuclei (CCN) activity under background conditions and during pollution events are investigated during the Aerosol-CCN-Cloud Closure Experiment (AC 3 Exp) campaign conducted at Xianghe, China in summer 2013. A gradual increase in size-resolved activation ratio (AR) with particle diameter ( D p ) suggests that aerosol particles have different hygroscopicities. During pollution events, the activation diameter ( D a ) measured at low supersaturation (SS) was significantly increased compared to background conditions. An increase was not observed when SS was 〉 0.4%. The hygroscopicity parameter (κ) was ~ 0.31–0.38 for particles in accumulation mode under background conditions. This range in magnitude of κ was ~ 20%, higher than κ derived under polluted conditions. For particles in nucleation or Aitken mode, κ ranged from 0.20–0.34 for background and polluted cases. Larger particles were on average more hygroscopic than smaller particles. The situation was more complex for heavy pollution particles because of the diversity in particle composition and mixing state. A non-parallel observation CCN closure test showed that uncertainties in CCN number concentration estimates ranged from 30–40%, which are associated with changes in particle composition as well as measurement uncertainties associated with bulk and size-resolved CCN methods. A case study showed that bulk CCN activation ratios increased as total condensation nuclei (CN) number concentrations ( N CN ) increased on background days. The background case also showed that bulk AR correlated well with the hygroscopicity parameter calculated from chemical volume fractions. On the contrary, bulk AR decreased with increasing total N CN during pollution events, but was closely related to the fraction of the total organic mass signal at m/z 44 ( f 44 ), which is usually associated with the particle's organic oxidation level. Our study highlights the importance of chemical composition in determining particle activation properties and underlines the significance of long-term observations of CCN under different atmospheric environments, especially regions with heavy pollution.

Description: Factors affecting flood insurance purchase in residential properties in Johor, Malaysia Natural Hazards and Earth System Science, 14, 3297-3310, 2014 Author(s): U. G. Aliagha, T. E. Jin, W. W. Choong, M. Nadzri Jaafar, and H. M. Ali High-impact floods have become a virtually annual experience in Malaysia, yet flood insurance has remained a grossly neglected part of comprehensive integrated flood risk management. Using discriminant analysis, this study seeks to identify the demand-side variables that best predict flood insurance purchase and risk aversion between two groups of residential homeowners in three districts of Johor State, Malaysia: those who purchased flood insurance and those who did not. Our results revealed an overall 34% purchase rate, with Kota Tinggi district having the highest (44%) and thus the highest degree of flood risk aversion. The Wilks' lambda F test for equality of group means, standardised discriminant function coefficients, structure correlation, and canonical correlation has clearly shown that there are strong significant attribute differences between the two groups of homeowners, based on the measures of objective flood risk exposure, subjective risk perception, and socio-economic cum demographic variables. However, the measures of subjective risk perception were found to be more predictive of flood insurance purchase and flood risk aversion.

Description: A century of sea level data and the UK's 2013/14 storm surges: an assessment of extremes and clustering using the Newlyn tide gauge record Ocean Science, 10, 1031-1045, 2014 Author(s): M. P. Wadey, I. D. Haigh, and J. M. Brown For the UK's longest and most complete sea level record (Newlyn), we assess extreme high waters and their temporal clustering; prompted by the 2013/2014 winter of storms and flooding. These are set into context against this almost 100-year record. We define annual periods for which storm activity and high sea levels can be compared on a year-by-year basis. Amongst the storms and high tides which affected Newlyn, the recent winter produced the largest recorded high water level (3 February 2014) and five other high water events above a 1 in 1-year return period. The large magnitude of tide and mean sea level, and the close inter-event spacings (of large return period high waters), suggests that the 2013/2014 extreme high water level "season" can be considered the most extreme on record. However, storm and sea level events may be classified in different ways. For example, in the context of sea level rise (which we calculate linearly as 1.81 ± 0.1 mm yr −1 from records between 1915 to 2014), a lower probability combination of surge and tide occurred on 29 January 1948, whilst the 1995/1996 storm surge season saw the most high waters of ≥ the 1 in 1-year return period. We provide a basic categorisation of the four types of extreme high water level cluster, ranging from consecutive tidal cycles to multiple years. The assessment is extended to other UK sites (with shorter sea level records and different tide-surge characteristics), which suggests 2013/2014 was particularly unusual. Further work will assess clustering mechanisms and flood system "memory".

Description: Derivation and evaluation of landslide-triggering thresholds by a Monte Carlo approach Hydrology and Earth System Sciences, 18, 4913-4931, 2014 Author(s): D. J. Peres and A. Cancelliere Assessment of landslide-triggering rainfall thresholds is useful for early warning in prone areas. In this paper, it is shown how stochastic rainfall models and hydrological and slope stability physically based models can be advantageously combined in a Monte Carlo simulation framework to generate virtually unlimited-length synthetic rainfall and related slope stability factor of safety data, exploiting the information contained in observed rainfall records and field-measurements of soil hydraulic and geotechnical parameters. The synthetic data set, dichotomized in triggering and non-triggering rainfall events, is analyzed by receiver operating characteristics (ROC) analysis to derive stochastic-input physically based thresholds that optimize the trade-off between correct and wrong predictions. Moreover, the specific modeling framework implemented in this work, based on hourly analysis, enables one to analyze the uncertainty related to variability of rainfall intensity within events and to past rainfall (antecedent rainfall). A specific focus is dedicated to the widely used power-law rainfall intensity–duration ( I – D ) thresholds. Results indicate that variability of intensity during rainfall events influences significantly rainfall intensity and duration associated with landslide triggering. Remarkably, when a time-variable rainfall-rate event is considered, the simulated triggering points may be separated with a very good approximation from the non-triggering ones by a I – D power-law equation, while a representation of rainfall as constant–intensity hyetographs globally leads to non-conservative results. This indicates that the I – D power-law equation is adequate to represent the triggering part due to transient infiltration produced by rainfall events of variable intensity and thus gives a physically based justification for this widely used threshold form, which provides results that are valid when landslide occurrence is mostly due to that part. These conditions are more likely to occur in hillslopes of low specific upslope contributing area, relatively high hydraulic conductivity and high critical wetness ratio. Otherwise, rainfall time history occurring before single rainfall events influences landslide triggering, determining whether a threshold based only on rainfall intensity and duration may be sufficient or it needs to be improved by the introduction of antecedent rainfall variables. Further analyses show that predictability of landslides decreases with soil depth, critical wetness ratio and the increase of vertical basal drainage (leakage) that occurs in the presence of a fractured bedrock.

Description: Climate change and stream temperature projections in the Columbia River basin: habitat implications of spatial variation in hydrologic drivers Hydrology and Earth System Sciences, 18, 4897-4912, 2014 Author(s): D. L. Ficklin, B. L. Barnhart, J. H. Knouft, I. T. Stewart, E. P. Maurer, S. L. Letsinger, and G. W. Whittaker Water temperature is a primary physical factor regulating the persistence and distribution of aquatic taxa. Considering projected increases in air temperature and changes in precipitation in the coming century, accurate assessment of suitable thermal habitats in freshwater systems is critical for predicting aquatic species' responses to changes in climate and for guiding adaptation strategies. We use a hydrologic model coupled with a stream temperature model and downscaled general circulation model outputs to explore the spatially and temporally varying changes in stream temperature for the late 21st century at the subbasin and ecological province scale for the Columbia River basin (CRB). On average, stream temperatures are projected to increase 3.5 °C for the spring, 5.2 °C for the summer, 2.7 °C for the fall, and 1.6 °C for the winter. While results indicate changes in stream temperature are correlated with changes in air temperature, our results also capture the important, and often ignored, influence of hydrological processes on changes in stream temperature. Decreases in future snowcover will result in increased thermal sensitivity within regions that were previously buffered by the cooling effect of flow originating as snowmelt. Other hydrological components, such as precipitation, surface runoff, lateral soil water flow, and groundwater inflow, are negatively correlated to increases in stream temperature depending on the ecological province and season. At the ecological province scale, the largest increase in annual stream temperature was within the Mountain Snake ecological province, which is characterized by migratory coldwater fish species. Stream temperature changes varied seasonally with the largest projected stream temperature increases occurring during the spring and summer for all ecological provinces. Our results indicate that stream temperatures are driven by local processes and ultimately require a physically explicit modeling approach to accurately characterize the habitat regulating the distribution and diversity of aquatic taxa.

Description: Monitoring and modelling of soil–plant interactions: the joint use of ERT, sap flow and Eddy Covariance data to characterize the volume of an orange tree root zone Hydrology and Earth System Sciences Discussions, 11, 13353-13384, 2014 Author(s): G. Cassiani, J. Boaga, D. Vanella, M. T. Perri, and S. Consoli Mass and energy exchanges between soil, plants and atmosphere control a number of key environmental processes involving hydrology, biota and climate. The understanding of these exchanges also play a critical role for practical purposes e.g. in precision agriculture. In this paper we present a methodology based on coupling innovative data collection and models in order to obtain quantitative estimates of the key parameters of such complex flow system. In particular we propose the use of hydro-geophysical monitoring via 4-D Electrical Resistivity Tomography (ERT) in conjunction with measurements of plant transpiration via sap flow and evapotranspiration from Eddy Covariance (EC). This abundance of data is fed to a spatially distributed soil model in order to characterize the distribution of active roots. We conducted experiments in an orange orchard in Eastern Sicily (Italy), characterized by the typical Mediterranean semi-arid climate. The subsoil dynamics, particularly influenced by irrigation and root uptake, were characterized mainly by the ERT setup, consisting of 48 buried electrodes on 4 instrumented micro boreholes (about 1.2 m deep) placed at the corners of a square (about 1.3 m in side) surrounding the orange tree, plus 24 mini-electrodes on the surface spaced 0.1 m on a square grid. During the monitoring, we collected repeated ERT and TDR soil moisture measurements, soil water samples, sap flow measurements from the orange tree and EC data. We conducted a laboratory calibration of the soil electrical properties as a function of moisture content and pore water electrical conductivity. Irrigation, precipitation, sap flow and ET data are available allowing knowledge of the system's long term forcing conditions on the system. This information was used to calibrate a 1-D Richards' equation model representing the dynamics of the volume monitored via 3-D ERT. Information on the soil hydraulic properties was collected from laboratory and field experiments. The successful results of the calibrated modeling exercise allow the quantification of the soil volume interested by root water uptake. This volume is much smaller (with a surface area less than 2 m 2 , and about 40 cm thickness) than expected and assumed in the design of classical drip irrigation schemes that prove to be losing at least half of the irrigated water that is not uptaken by the plants.

Description: Energy exchange and wave action conservation for magnetohydrodynamic (MHD) waves in a general, slowly varying medium Annales Geophysicae, 32, 1495-1510, 2014 Author(s): A. D. M. Walker Magnetohydrodynamic (MHD) waves in the solar wind and magnetosphere are propagated in a medium whose velocity is comparable to or greater than the wave velocity and which varies in both space and time. In the approximation where the scales of the time and space variation are long compared with the period and wavelength, the ray-tracing equations can be generalized and then include an additional first-order differential equation that determines the variation of frequency. In such circumstances the wave can exchange energy with the background: wave energy is not conserved. In such processes the wave action theorem shows that the wave action, defined as the ratio of the wave energy to the frequency in the local rest frame, is conserved. In this paper we discuss ray-tracing techniques and the energy exchange relation for MHD waves. We then provide a unified account of how to deal with energy transport by MHD waves in non-uniform media. The wave action theorem is derived directly from the basic MHD equations for sound waves, transverse Alfvén waves, and the fast and slow magnetosonic waves. The techniques described are applied to a number of illustrative cases. These include a sound wave in a medium undergoing a uniform compression, an isotropic Alfvén wave in a steady-state shear layer, and a transverse Alfvén wave in a simple model of the magnetotail undergoing compression. In each case the nature and magnitude of the energy exchange between wave and background is found.

Description: Examining the major contributors and controlling factors of ozone production in a rural area of the Yangtze River Delta region during harvest season Atmospheric Chemistry and Physics Discussions, 14, 30913-30945, 2014 Author(s): X. Pan, Y. Kanaya, H. Tanimoto, S. Inomata, Z. Wang, S. Kudo, and I. Uno Open biomass burning (OBB) has been reported to emit substantial amounts of non-methane hydrocarbons (NMHCs), and the mixing of OBB with urban plumes could exacerbate regional ozone (O 3 ) pollution. In the present study, an observational field campaign was performed in a rural area at the edge of Yangtze River Delta region (YRDR) during harvest season when intensive open burning of wheat residues was observed. The O 3 production rate at the site was calculated using a photochemical box model (Regional Atmospheric Chemical Mechanism, Version 2) constrained by real-time ambient measurements (e.g., O 3 , volatile organic compounds (VOCs), the sum of NO 2 + NO (NO x ), J values). During the period impacted by OBB, the O 3 concentration frequently exceeded 100 ppbv. Analysis showed that the net O 3 production was pronounced, in particular when the site was characterized by a prevailing southerly wind that also brought substantial amounts of NO x emitted from urban areas. At these times, the maximum rate of O 3 production was 20 ppbv h −1 with potential production rate of 102 ppbv on a daily basis. The O 3 production at the site was typically VOC-sensitive in the morning because NO x dominated the plumes. However, in the afternoon, conditions became NO x -sensitive due to the rapid photochemical consumption of NO x in the production of O 3 . A positive matrix factorization analysis indicated that solvent usage and OBB were the primary contributors to the mass fraction of ambient NMHCs observed at the study site, and were responsible for 35 and 23% of the total O 3 production, respectively. The preferential presence of NO x in the morning inhibited net O 3 production; meanwhile O 3 built up in the afternoon due to a decrease in NO x concentrations. These results indicated that a joint effort in the regulation of solvent (aromatics) usage and OBB, as well as NO x from on-road vehicle exhaust may be effective in eliminating high-O 3 pollution risk in the rural areas of the YRDR.

Description: On the use of the post-closure method uncertainty band to evaluate the performance of land surface models against eddy covariance flux data Biogeosciences Discussions, 11, 16911-16951, 2014 Author(s): J. Ingwersen, K. Imukova, P. Högy, and T. Streck The energy balance of eddy covariance (EC) flux data is normally not closed. Therefore, at least if used for modeling, EC flux data are usually post-closed, i.e. the measured turbulent fluxes are adjusted so as to close the energy balance. At the current state of knowledge, however, it is not clear how to partition the missing energy in the right way. Eddy flux data therefore contain some uncertainty due to the unknown nature of the energy balance gap, which should be considered in model evaluation and the interpretation of simulation results. We propose to construct the post-closure method uncertainty band (PUB), which essentially designates the differences between non-adjusted flux data and flux data adjusted with the three post-closure methods (Bowen ratio, latent heat flux (LE) and sensible heat flux ( H ) method). To demonstrate this approach, simulations with the NOAH-MP land surface model were evaluated based on EC measurements conducted at a winter wheat stand in Southwest Germany in 2011, and the performance of the Jarvis and Ball–Berry stomatal resistance scheme was compared. The width of the PUB of the LE was up to 110 W m –2 (21% of net radiation). Our study shows that it is crucial to account for the uncertainty of EC flux data originating from lacking energy balance closure. Working with only a single post-closing method might result in severe misinterpretations in model-data comparisons.

Description: Carbon, oxygen and biological productivity in the Southern Ocean in and out the Kerguelen plume: CARIOCA drifter results Biogeosciences Discussions, 11, 16877-16909, 2014 Author(s): L. Merlivat, J. Boutin, and F. d'Ovidio The Kerguelen Plateau region in the Indian sector of the Southern Ocean supports annually a large-scale phytoplankton bloom which is naturally fertilized with iron. As part of the second Kerguelen Ocean and Plateau compared Study expedition (KEOPS2) in austral spring (October–November 2011), one Carioca buoy was deployed east of the Kerguelen plateau. It drifted eastward downstream in the Kerguelen plume. Hourly surface measurements of p CO 2 , O 2 and ancillary observations were collected between 1 November 2011 to 12 February 2012 with the aim of characterizing the spatial and temporal variability of the biological Net Community Production (NCP) downstream the Kerguelen plateau, assess the impact of iron-induced productivity on the biological carbon consumption and consequently on the CO 2 flux exchanged at the air–sea interface. The trajectory of the buoy until mid-December was within the longitude range, 72–83° E, close to the polar front and then in the polar frontal zone, PFZ, until 97° E. From 17 November to 16 December, the buoy drifted within the Kerguelen plume following a filament carrying dissolved iron, DFe, for a total distance of 700 km. In the first part of the trajectory, the ocean surface waters are a sink for CO 2 and a source for CO 2 , with fluxes of respective mean values equal to −8 and +38 mmol CO 2 m −2 d −1 . Eastward, as the buoy escapes the iron enriched filament, the fluxes are in opposite direction, with respective mean values of +5 and −48 mmol O 2 m −2 d −1 . These numbers clearly indicate the strong impact of biological processes on the biogeochemistry in the surface waters within the Kerguelen plume in November-mid-December, while it is undetectable eastward in the PFZ from mid-December to mid-February. While the buoy follows the Fe enriched filament, simultaneous observations of dissolved inorganic carbon, DIC, and dissolved oxygen, O 2 , highlight biological events lasting from 2 to 4 days. Stoichiometric ratios, O 2 /C, between 1.1 and 1.4 are observed indicating new and regenerated production regimes. NCP estimates range from 60 to 140 mmol C m −2 d −1 . Based on the relationship between the time a water parcel has left the plateau and its iron content, we have highlighted that the main control on the value of NCP is the availability of iron in the upper water column, with the largest NCP occurring in waters that have recently left the plateau and presented the highest iron concentrations.

Description: Comparison of UV/Vis and FDOM sensors for in situ monitoring of stream DOC concentrations Biogeosciences Discussions, 11, 16855-16876, 2014 Author(s): G.-Y. Yoo, Y. Jeong, E.-J. Lee, J.-H. Park, and N.-H. Oh Optical measurements using ultra-violet/visible (UV/Vis) spectrophotometric sensors and fluorescent dissolved organic matter (FDOM) sensors have recently been used as proxies of dissolved organic carbon (DOC) concentrations of streams and rivers at high temporal resolution. Despite of the merits of the sensors, temperature changes and particulate matter in water can interfere the sensor readings, over- or under-estimating DOC concentrations. However, little efforts have been made to compare responses of the two types of the sensors in natural conditions. We conducted both laboratory experiments and in situ monitoring with a UV/Vis sensor and a FDOM sensor during the three storm events in the fall of 2012 and the spring of 2013 in a forest stream in Korea in order to compare their performance. Laboratory experiments using the Suwannee River natural organic matter, humic acid, and fulvic acid demonstrated strong linear relationships between both the sensor signals and measured DOC concentrations with R 2 ≥ 0.98. Although temperature compensation might not be needed for the UV/Vis sensor, it was sensitive to relativley small changes in turbidity. In contrast, the FDOM sensor was insenstive to relatively low turbidity while the FDOM sensor outputs decreased significantly as temperature increased, requiring temperature compensated FDOM (e.g. FDOM 20 for 20 °C) for in situ monitoring of DOC. The results suggest that both sensors can be employed as a~proxy for stream DOC concentrations after temperature and turbidity compensation in a forest stream where terrestrially derived humic-like materials are dominant components.

Description: Modeling the impact of riverine DON removal by marine bacterioplankton on primary production in the Arctic Ocean Biogeosciences Discussions, 11, 16953-16992, 2014 Author(s): V. Le Fouest, M. Manizza, B. Tremblay, and M. Babin The planktonic and biogeochemical dynamics of the Arctic shelves exhibit a strong variability in response to Arctic warming. In this study, in order to elucidate on the processes regulating the production of phytoplankton (PP) and bacterioplankton (BP) and their interactions, we employ a biogeochemical model coupled to a pan-Arctic ocean-sea ice model (MITgcm) to explicitly simulate and quantify the contribution of usable dissolved organic nitrogen (DON) drained by the major circum-Arctic rivers on PP and BP in a scenario of melting sea ice (1998–2011). Model simulations suggest that on average between 1998 and 2011, the removal of usable RDON by bacterioplankton is responsible of a ~26% increase of the annual BP for the whole Arctic Ocean. With respect to total PP, the model simulates an increase of ~8% on an annual basis and of ~18% in summer. Recycled ammonium is responsible for the PP increase. The recycling of RDON by bacterioplankton promotes higher BP and PP but there is no significant temporal trend in the BP : PP ratio within the ice-free shelves over the 1998–2011 period. This suggests no significant evolution in the balance between autotrophy and heterotrophy in the last decade with a constant annual flux of RDON into the coastal ocean although changes in RDON supply and further reduction in sea ice cover could potentially alter this delicate balance.

Description: The double layers in the plasma sheet boundary layer during magnetic reconnection Nonlinear Processes in Geophysics Discussions, 1, 1657-1671, 2014 Author(s): J. Guo and B. Yu We studied the evolutions of double layers which appear after the magnetic reconnection through two-dimensional electromagnetic particle-in-cell simulation. The simulation results show that the double layers are formed in the plasma sheet boundary layer after magnetic reconnection. At first, the double layers which have unipolar structures are formed. And then the double layers turn into bipolar structures, which will couple with another new weak bipolar structure. Thus a new double layer or tripolar structure comes into being. The double layers found in our work are about several ten Debye lengths, which accords with the observation results. It is suggested that the electron beam formed during the magnetic reconnection is responsible for the production of the double layers.

Description: Atmospheric submicron aerosol composition and particulate organic nitrate formation in a boreal forestland–urban mixed region Atmospheric Chemistry and Physics, 14, 13483-13495, 2014 Author(s): L. Q. Hao, A. Kortelainen, S. Romakkaniemi, H. Portin, A. Jaatinen, A. Leskinen, M. Komppula, P. Miettinen, D. Sueper, A. Pajunoja, J. N. Smith, K. E. J. Lehtinen, D. R. Worsnop, A. Laaksonen, and A. Virtanen The Puijo aerosol–cloud observation station is a unique measurement site for its location in the mixed region between the boreal forestland and the municipality of Kuopio, Finland. A measurement campaign was carried out at the station during fall 2010. An Aerodyne high-resolution time-of-flight aerosol mass spectrometer (HR-Tof-AMS) was deployed to characterize the atmospheric submicron aerosols. Positive matrix factorization (PMF) was applied to the unified high-resolution mass spectra organic species with NO + and NO 2 + ions to discover the intrinsic relationships between the organic and inorganic species and their daily cycles. On average, the submicron aerosols in this study were dominated by organic and sulfate species, composing 48.2 and 28.7% of total observed aerosol mass, respectively, with smaller contributions from ammonium (9.3%), nitrate (4.9%), chloride (0.8%) and BC (8.1%). The sources of these species included the primary emissions originating from the city area, secondary formation from both natural and anthropogenic emissions and regional transport. The PMF analysis succeeded in separating the mixed organic and inorganic spectra into three distinct organic and one inorganic factors. For organic factors, the semi-volatile oxygenated organic aerosol (SVOOA) and low-volatility oxygenated OA (LVOOA) accounted for 54.8 and 36.3% of total organic masses, respectively, while the hydrocarbon-like organic aerosol (HOA) accounted for 8.9% of total organics, with its main source from urban emissions. The inorganic factor is identified as NH 4 NO 3 , comprising 6.9% of the fitted aerosol mass by PMF. Based on the PMF results, the nitrate species were separated into organic and inorganic components, with the organic nitrates contributing one-third of the total nitrate mass. The results highlight both anthropogenic and biogenic emissions as important atmospheric aerosol sources in a forest–urban mixed region.

Description: Formaldehyde (HCHO) in air, snow and interstitial air at Concordia (East Antarctic plateau) in summer Atmospheric Chemistry and Physics Discussions, 14, 32027-32070, 2014 Author(s): S. Preunkert, M. Legrand, M. Frey, A. Kukui, J. Savarino, H. Gallée, M. King, B. Jourdain, W. Vicars, and D. Helmig During the 2011/12 and 2012/13 austral summers HCHO was investigated for the first time in ambient air, snow, and interstitial air at the Concordia site located near Dome C on the East Antarctic plateau by deploying an Aerolaser AL-4021 analyser. Snow emission fluxes were estimated from vertical gradients of mixing ratios observed between 1 cm and 1 m above the snow surface as well as between interstitial air a few cm below the surface and in air just above the snow-pack. Typical flux values range between 1 to 2 × 10 12 molecules m −2 s −1 at night and 3 to 5 × 10 12 molecules m −2 s −1 at noon. Shading experiments suggest that the photochemical HCHO production in the snowpack at Concordia remains negligible compared to temperature-driven air–snow exchanges. At 1 m above the snow surface, the observed mean mixing ratio of 130 pptv and its diurnal cycle characterized by a slight decrease around noon are quite well reproduced by 1-D simulations that include snow emissions and gas phase methane oxidation chemistry.

Description: Sources and fluxes of organic nitrogen in precipitation over the southern East/Japan Sea: potential impacts on marine productivity Atmospheric Chemistry and Physics Discussions, 14, 31987-32025, 2014 Author(s): G. Yan and G. Kim We measured total dissolved reactive nitrogen in precipitation samples collected at Uljin, a Korean coastal site upwind of the southern East/Japan Sea (EJS), selected as a representative study site of atmospheric deposition over the northwestern Pacific margin. NO 3 − was found to be the most abundant nitrogen species, followed by NH 4 + and dissolved organic nitrogen (DON). Air mass back trajectory analysis revealed that a significant fraction of the inorganic nitrogen (NO 3 − and NH 4 + ) originated from mainland Asia, whereas the DON was primarily derived from Korea. Using varimax-rotated factor analysis in combination with major ions as tracers, agricultural activities in Korea were identified as the primary sources of DON in these samples. In addition, a positive correlation was found at Uljin between the size of organic fraction in total reactive nitrogen and nitrogen to carbon atomic ratio in organic matter. This correlation has also been observed at other locations worldwide, implying the utilization potential of atmospheric organic nitrogen might increase with its proportion in total nitrogen. Combining wet deposition measurements in this study with literature values for dry deposition observed at a remote island in the EJS, the total atmospheric depositional flux of reactive nitrogen was estimated to be 115 mmol N m −2 yr −1 over the southern EJS. Our study sheds new light on the potentially significant contribution to primary productivity of the northwestern Pacific Ocean by atmospheric deposition of nitrogen, especially the organic fraction.

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

Description: Intercomparison of stratospheric gravity wave observations with AIRS and IASI Atmospheric Measurement Techniques, 7, 4517-4537, 2014 Author(s): L. Hoffmann, M. J. Alexander, C. Clerbaux, A. W. Grimsdell, C. I. Meyer, T. Rößler, and B. Tournier Gravity waves are an important driver for the atmospheric circulation and have substantial impact on weather and climate. Satellite instruments offer excellent opportunities to study gravity waves on a global scale. This study focuses on observations from the Atmospheric Infrared Sounder (AIRS) onboard the National Aeronautics and Space Administration Aqua satellite and the Infrared Atmospheric Sounding Interferometer (IASI) onboard the European MetOp satellites. The main aim of this study is an intercomparison of stratospheric gravity wave observations of both instruments. In particular, we analyzed AIRS and IASI 4.3 μm brightness temperature measurements, which directly relate to stratospheric temperature. Three case studies showed that AIRS and IASI provide a clear and consistent picture of the temporal development of individual gravity wave events. Statistical comparisons based on a 5-year period of measurements (2008–2012) showed similar spatial and temporal patterns of gravity wave activity. However, the statistical comparisons also revealed systematic differences of variances between AIRS and IASI that we attribute to the different spatial measurement characteristics of both instruments. We also found differences between day- and nighttime data that are partly due to the local time variations of the gravity wave sources. While AIRS has been used successfully in many previous gravity wave studies, IASI data are applied here for the first time for that purpose. Our study shows that gravity wave observations from different hyperspectral infrared sounders such as AIRS and IASI can be directly related to each other, if instrument-specific characteristics such as different noise levels and spatial resolution and sampling are carefully considered. The ability to combine observations from different satellites provides an opportunity to create a long-term record, which is an exciting prospect for future climatological studies of stratospheric gravity wave activity.

Description: Collection efficiency of the soot-particle aerosol mass spectrometer (SP-AMS) for internally mixed particulate black carbon Atmospheric Measurement Techniques, 7, 4507-4516, 2014 Author(s): M. D. Willis, A. K. Y. Lee, T. B. Onasch, E. C. Fortner, L. R. Williams, A. T. Lambe, D. R. Worsnop, and J. P. D. Abbatt The soot-particle aerosol mass spectrometer (SP-AMS) uses an intra-cavity infrared laser to vaporize refractory black carbon (rBC) containing particles, making the particle beam–laser beam overlap critical in determining the collection efficiency (CE) for rBC and associated non-refractory particulate matter (NR-PM). This work evaluates the ability of the SP-AMS to quantify rBC and NR-PM mass in internally mixed particles with different thicknesses of organic coating. Using apparent relative ionization efficiencies for uncoated and thickly coated rBC particles, we report measurements of SP-AMS sensitivity to NR-PM and rBC, for Regal Black, the recommended particulate calibration material. Beam width probe (BWP) measurements are used to illustrate an increase in sensitivity for highly coated particles due to narrowing of the particle beam, which enhances the CE of the SP-AMS by increasing the laser beam–particle beam overlap. Assuming complete overlap for thick coatings, we estimate CE for bare Regal Black particles of 0.6 ± 0.1, which suggests that previously measured SP-AMS sensitivities to Regal Black were underestimated by up to a factor of 2. The efficacy of the BWP measurements is highlighted by studies at a busy road in downtown Toronto and at a non-roadside location, which show particle beam widths similar to, but greater than that of bare Regal Black and coated Regal Black, respectively. Further BWP measurements at field locations will help to constrain the range of CE for fresh and aged rBC-containing particles. The ability of the SP-AMS to quantitatively assess the composition of internally mixed particles is validated through measurements of laboratory-generated organic coated particles, which demonstrate that the SP-AMS can quantify rBC and NR-PM over a wide range of particle compositions and rBC core sizes.

Description: Dynamical analysis of sea-breeze hodograph rotation in Sardinia Atmospheric Chemistry and Physics, 14, 13471-13481, 2014 Author(s): N. Moisseeva and D. G. Steyn This study investigates the diurnal evolution of sea-breeze (SB) rotation over an island at the middle latitudes. Earlier research on sea breezes in Sardinia shows that the onshore winds around various coasts of the island exhibit both the theoretically predicted clockwise rotation as well as seemingly anomalous anticlockwise rotation. A non-hydrostatic fully compressible numerical model (WRF) is used to simulate wind fields on and around the island on previously studied sea-breeze days, and is shown to capture the circulation on all coasts accurately. Diurnal rotation of wind is examined, and patterns of clockwise and anticlockwise rotation are identified. A dynamical analysis is performed by extracting individual forcing terms from the horizontal momentum equations. Analysis of several regions around the island shows that the direction of rotation is a result of a complex interaction between near-surface and synoptic pressure gradient, Coriolis and advection forcings. An idealized simulation is performed over an artificial island with dramatically simplified topography yet similar dimensions and latitude to Sardinia. Dynamical analysis of the idealized case reveals a rather different pattern of hodograph rotation to the real Sardinia, yet similar underlying dynamics. The research provides new insights into the dynamics underlying sea-breeze hodograph rotation, especially in coastal zones with a complex topography and/or coastline.

Description: Validation of first chemistry mode retrieval results from new limb- imaging FTS GLORIA with correlative MIPAS-STR observations Atmospheric Measurement Techniques Discussions, 7, 12691-12717, 2014 Author(s): W. Woiwode, O. Suminska-Ebersoldt, H. Oelhaf, M. Höpfner, G. V. Belyaev, A. Ebersoldt, F. Friedl-Vallon, J.-U. Grooß, T. Gulde, M. Kaufmann, A. Kleinert, M. Krämer, E. Kretschmer, T. Kulessa, G. Maucher, T. Neubert, C. Piesch, P. Preusse, M. Riese, H. Rongen, C. Sartorius, G. Schardt, A. Schönfeld, D. Schuettemeyer, M. K. Sha, F. Stroh, J. Ungermann, C. M. Volk, and J. Orphal We report first chemistry mode retrieval results from the new airborne limb- imaging infrared FTS (Fourier transform spectrometer) GLORIA and comparisons with observations by the conventional airborne limb- scanning infrared FTS MIPAS-STR. For GLORIA, the flights aboard the high-altitude research aircraft M55 Geophysica during the ESSenCe campaign (ESa Sounder Campaign 2011) were the very first in field deployment after several years of development. The simultaneous observations of GLORIA and MIPAS-STR during the flight on 16 December 2011 inside the polar vortex and under the conditions of optically partially transparent polar stratospheric clouds (PSCs) provided us the unique opportunity to compare the observations by two different infrared FTS generations directly. The retrieval results of temperature, HNO 3 , O 3 , H 2 O, CFC-11 and CFC-12 show reasonable agreement of GLORIA with MIPAS-STR and collocated in-situ observations. For the horizontally binned hyperspectral limb-images, the GLORIA sampling outnumbered the horizontal cross-track sampling of MIPAS-STR by up to one order of magnitude. Depending on the target parameter, typical vertical resolutions of 0.5 to 2.0 km were obtained for GLORIA and are typically by factors of 2 to 4 better compared to MIPAS-STR. While the improvement of the performance, characterisation and data processing of GLORIA are subject of ongoing work, the presented first results already demonstrate the considerable gain in sampling and vertical resolution achieved with GLORIA.

Description: Reference quality upper-air measurements: GRUAN data processing for the Vaisala RS92 radiosonde Atmospheric Measurement Techniques, 7, 4463-4490, 2014 Author(s): R. J. Dirksen, M. Sommer, F. J. Immler, D. F. Hurst, R. Kivi, and H. Vömel The GCOS (Global Climate Observing System) Reference Upper-Air Network (GRUAN) data processing for the Vaisala RS92 radiosonde was developed to meet the criteria for reference measurements. These criteria stipulate the collection of metadata, the use of well-documented correction algorithms, and estimates of the measurement uncertainty. An important and novel aspect of the GRUAN processing is that the uncertainty estimates are vertically resolved. This paper describes the algorithms that are applied in version 2 of the GRUAN processing to correct for systematic errors in radiosonde measurements of pressure, temperature, humidity, and wind, as well as how the uncertainties related to these error sources are derived. Currently, the RS92 is launched on a regular basis at 13 out of 15 GRUAN sites. An additional GRUAN requirement for performing reference measurements with the RS92 is that the manufacturer-prescribed procedure for the radiosonde's preparation, i.e. heated reconditioning of the sensors and recalibration during ground check, is followed. In the GRUAN processing however, the recalibration of the humidity sensors that is applied during ground check is removed. For the dominant error source, solar radiation, laboratory experiments were performed to investigate and model its effect on the RS92's temperature and humidity measurements. GRUAN uncertainty estimates are 0.15 K for night-time temperature measurements and approximately 0.6 K at 25 km during daytime. The other uncertainty estimates are up to 6% relative humidity for humidity, 10–50 m for geopotential height, 0.6 hPa for pressure, 0.4–1 m s −1 for wind speed, and 1° for wind direction. Daytime temperature profiles for GRUAN and Vaisala processing are comparable and consistent within the estimated uncertainty. GRUAN daytime humidity profiles are up to 15% moister than Vaisala processed profiles, of which two-thirds is due to the radiation dry bias correction and one-third is due to an additional calibration correction. Redundant measurements with frost point hygrometers (CFH and NOAA FPH) show that GRUAN-processed RS92 humidity profiles and frost point data agree within 15% in the troposphere. No systematic biases occur, apart from a 5% dry bias for GRUAN data around −40 °C at night.

Description: ObsPack: a framework for the preparation, delivery, and attribution of atmospheric greenhouse gas measurements Earth System Science Data, 6, 375-384, 2014 Author(s): K. A. Masarie, W. Peters, A. R. Jacobson, and P. P. Tans Observation Package (ObsPack) is a framework designed to bring together atmospheric greenhouse gas observations from a variety of sampling platforms, prepare them with specific applications in mind, and package and distribute them in a self-consistent and well-documented product. Data products created using the ObsPack framework (called "ObsPack products") are intended to support carbon cycle modeling studies and represent a next generation of value-added greenhouse gas observation products modeled after the cooperative GLOBALVIEW products introduced in 1996. Depending on intended use, ObsPack products may include data in their original form reformatted using the ObsPack framework or may contain derived data consisting of averages, subsets, or smoothed representations of original data. All products include extensive ancillary information (metadata) intended to help ensure the data are used appropriately, their calibration and quality assurance history are clearly described, and that individuals responsible for the measurements (data providers or principal investigators (PIs)) are properly acknowledged for their work. ObsPack products are made freely available using a distribution strategy designed to improve communication between data providers and product users. The strategy includes a data usage policy that requires users to directly communicate with data providers and an automated e-mail notification system triggered when a product is accessed. ObsPack products will be assigned a unique digital object identifier (DOI) to ensure each product can be unambiguously identified in scientific literature. Here we describe the ObsPack framework and its potential role in supporting the evolving needs of both data providers and product users.

Description: Elastic envelope inversion using multicomponent seismic data without low frequency Nonlinear Processes in Geophysics Discussions, 1, 1757-1802, 2014 Author(s): C. Huang, L. Dong, Y. Liu, and B. Chi Low frequency is a key issue to reduce the nonlinearity of elastic full waveform inversion. Hence, the lack of low frequency in recorded seismic data is one of the most challenging problems in elastic full waveform inversion. Theoretical derivations and numerical analysis are presented in this paper to show that envelope operator can retrieve strong low frequency modulation signal demodulated in multicomponent data, no matter what the frequency bands of the data is. With the benefit of such low frequency information, we use elastic envelope of multicomponent data to construct the objective function and present an elastic envelope inversion method to recover the long-wavelength components of the subsurface model, especially for the S-wave velocity model. Numerical tests using synthetic data for the Marmousi-II model prove the effectiveness of the proposed elastic envelope inversion method, especially when low frequency is missing in multicomponent data and when initial model is far from the true model. The elastic envelope can reduce the nonlinearity of inversion and can provide an excellent starting model.

Description: Large eddy simulation of sediment transport over rippled beds Nonlinear Processes in Geophysics, 21, 1169-1184, 2014 Author(s): J. C. Harris and S. T. Grilli Wave-induced boundary layer (BL) flows over sandy rippled bottoms are studied using a numerical model that applies a one-way coupling of a "far-field" inviscid flow model to a "near-field" large eddy simulation (LES) Navier–Stokes (NS) model. The incident inviscid velocity and pressure fields force the LES, in which near-field, wave-induced, turbulent bottom BL flows are simulated. A sediment suspension and transport model is embedded within the coupled flow model. The numerical implementation of the various models has been reported elsewhere, where we showed that the LES was able to accurately simulate both mean flow and turbulent statistics for oscillatory BL flows over a flat, rough bed. Here we show that the model accurately predicts the mean velocity fields and suspended sediment concentration for oscillatory flows over full-scale vortex ripples. Tests show that surface roughness has a significant effect on the results. Beyond increasing our insight into wave-induced oscillatory bottom BL physics, sophisticated coupled models of sediment transport such as that presented have the potential to make quantitative predictions of sediment transport and erosion/accretion around partly buried objects in the bottom, which is important for a vast array of bottom deployed instrumentation and other practical ocean engineering problems.

Description: The transient variation of the complexes of the low latitude ionosphere within the equatorial ionization anomaly region of Nigeria Nonlinear Processes in Geophysics Discussions, 1, 1855-1903, 2014 Author(s): A. B. Rabiu, B. O. Ogunsua, I. A. Fuwape, and J. A. Laoye The quest to find an index for proper characterization and description of the dynamical response of the ionosphere to external influences and its various internal irregularities has led to the study of the day to day variations of the chaoticity and dynamical complexity of the ionosphere. This study was conducted using Global Positioning System (GPS) Total Electron Content (TEC) time series, measured in the year 2011, from 5 GPS receiver stations in Nigeria which lies within the Equatorial Ionization Anomaly region. The nonlinear aspect of the TEC time series were obtained by detrending the data. The detrended TEC time series were subjected to various analyses for phase space reconstruction and to obtain the values of chaotic quantifiers which are Lyapunov exponents LE, correlation dimension, and Tsallis entropy for the study of dynamical complexity. The results show positive Lyapunov exponents for all days which indicate chaoticity of the ionosphere with no definite pattern for both quiet and disturbed days. However values of LE were lower for the storm period compared to its nearest relative quiet periods for all the stations. Considering all the days of the year the daily/transient variations show no definite pattern for each month but day to day values of Lyapunov exponent for the entire year show a wavelike semiannual variation pattern with lower values around March, April, September and October, a change in pattern which demonstrates the self-organized critical phenomenon of the system. This can be seen from the correlation dimension with values between 2.7 and 3.2 with lower values occurring mostly during storm periods demonstrating a phase transition from higher dimension during the quiet periods to lower dimension during storms for most of the stations. The values of Tsallis entropy show similar variation pattern with that of Lyapunov exponent with a lot of agreement in their comparison, with all computed values of Lyapunov exponent correlating with values of Tsallis entropy within the range of 0.79 to 0.82. These results show that Lyapunov quantifiers can be used together as indices in the study of the variations of the dynamical complexity of the ionosphere. The presence of chaos and high variations in the dynamical complexity, even at quiet periods in the ionosphere may be due to the internal dynamics and inherent irregularities of the ionosphere which exhibit non-linear properties. However, this inherent dynamics may be complicated by external factors like geomagnetic storms. This may be the main reason for the drop in the values of Lyapunov exponent and Tsallis entropy during storms. The results also show a strong interplay between determinism and stochasticity, as the ionosphere shows its response to changes in solar activities and in its internal dynamics. The dynamical behavior of the ionosphere throughout the year as described by these quantifiers, were discussed in this work.

Description: Modelling of the anthropogenic tritium transient and its decay product helium-3 in the Mediterranean Sea using a high-resolution regional model Ocean Science Discussions, 11, 2691-2732, 2014 Author(s): M. Ayache, J.-C. Dutay, P. Jean-Baptiste, K. Beranger, T. Arsouze, J. Beuvier, J. Palmieri, B. Le-vu, and W. Roether This numerical study provides the first simulation of the anthropogenic tritium invasion and its decay product helium-3 ( 3 He) in the Mediterranean Sea. The simulation covers the entire tritium ( 3 H) transient generated by the atmospheric nuclear-weapon tests performed in the 1950s and early 1960s and run till 2011. Tritium, helium-3 and their derived age estimates are particularly suitable for studying intermediate and deep-water ventilation and spreading of water masses at intermediate/deep levels. The simulation is made using a high resolution regional model NEMO-MED12 forced at the surface with prescribed tritium evolution derived from observations. The simulation is compared to measurements of tritium and helium-3 performed along large-scale transects in the Mediterranean Sea during the last few decades on cruises of Meteor M5/6, M31/1, M44/4, M51/2, M84/3, and Poseidon 234. The results show that the input function used for the tritium, generates a realistic distribution of the main hydrographic features of the Mediterranean Sea circulation. In the eastern basin, the results highlight the weak formation of Adriatic Deep Water in the model, which explains its weak contribution to the Eastern Mediterranean Deep Water in the Ionian sub-basin. It produces a realistic representation of the Eastern Mediterranean Transient signal, simulating a deep-water formation in the Aegean sub-basin at the beginning of the 1993, with a realistic timing of deep-water renewal in the eastern basin. In the western basin, the unusual intense deep convection event of winter 2005 in the Gulf of Lions during the Western Mediterranean Transition is simulated. However the spreading of the recently ventilated deep water toward the South is too weak. The ventilation and spreading of the Levantine Intermediate Water from the eastern basin toward the western basin is simulated with realistic tracer-age distribution compared to observation-based estimates.

Description: Atmospheric measurements of ratios between CO 2 and co-emitted species from traffic: a tunnel study in the Paris megacity Atmospheric Chemistry and Physics, 14, 12871-12882, 2014 Author(s): L. Ammoura, I. Xueref-Remy, V. Gros, A. Baudic, B. Bonsang, J.-E. Petit, O. Perrussel, N. Bonnaire, J. Sciare, and F. Chevallier Measurements of CO 2 , CO, NO x and selected Volatile Organic Compounds (VOCs) mole fractions were performed continuously during a 10-day period in the Guy Môquet tunnel in Thiais, a peri-urban area about 15 km south of the centre of Paris, between 28 September and 8 October 2012. This data set is used here to identify the characteristics of traffic-emitted CO 2 by evaluating its ratios to co-emitted species for the first time in the Paris region. High coefficients of determination ( r 2 〉 0.7) are observed between CO 2 and certain compounds that are characteristic of the traffic source (CO, NO x , benzene, xylenes and acetylene). Weak correlations ( r 2 〈 0.2) are found with species such as propane, n-butane and i -butane that are associated with fuel evaporation, an insignificant source for CO 2 . To better characterise the traffic signal we focus only on species that are well-correlated with CO 2 and on rush-hour periods characterised by the highest traffic-related mole fractions. From those mole fractions we remove the nighttime-average weekday mole fraction obtained for each species that we infer to be the most appropriate background signal for our study. Then we calculate observed Δspecies / ΔCO 2 ratios, which we compare with the ones provided by the 2010 bottom–up high-resolved regional emission inventory from Airparif (the association in charge of monitoring the air quality in Île-de-France), focusing on local emission data for the specific road of the tunnel. We find an excellent agreement (2%) between the local inventory emission CO / CO 2 ratio and our observed ΔCO / ΔCO 2 ratio. Former tunnel experiments carried out elsewhere in the world provided observed ΔCO / ΔCO 2 ratios that differ from 49 to 592% to ours. This variability can be related to technological improvement of vehicles, differences in driving conditions, and fleet composition. We also find a satisfactory agreement with the Airparif inventory for n-propylbenzene, n-pentane and xylenes to CO 2 ratios. For most of the other species, the ratios obtained from the local emission inventory overestimate the observed ratios to CO 2 by 34 to more than 300%. However, the emission ratios of NO x , o -xylene and i -pentane are underestimated by 30 to 79%. One main cause of such high differences between the inventory and our observations is likely the obsolete feature of the VOCs speciation matrix of the inventory that has not been updated since 1998, although law regulations on some VOCs have occurred since that time. Our study bears important consequences, discussed in the conclusion, for the characterisation of the urban CO 2 plume and for atmospheric inverse modelling of urban CO 2 emissions.

Description: The use of SMILES data to study ozone loss in the Arctic winter 2009/2010 and comparison with Odin/SMR data using assimilation techniques Atmospheric Chemistry and Physics, 14, 12855-12869, 2014 Author(s): K. Sagi, D. Murtagh, J. Urban, H. Sagawa, and Y. Kasai The Superconducting Submillimeter-Wave Limb-Emission Sounder (SMILES) on board the International Space Station observed ozone in the stratosphere with high precision from October 2009 to April 2010. Although SMILES measurements only cover latitudes from 38° S to 65° N, the combination of data assimilation methods and an isentropic advection model allows us to quantify the ozone depletion in the 2009/2010 Arctic polar winter by making use of the instability of the polar vortex in the northern hemisphere. Ozone data from both SMILES and Odin/SMR (Sub-Millimetre Radiometer) for the winter were assimilated into the Dynamical Isentropic Assimilation Model for OdiN Data (DIAMOND). DIAMOND is an off-line wind-driven transport model on isentropic surfaces. Wind data from the operational analyses of the European Centre for Medium- Range Weather Forecasts (ECMWF) were used to drive the model. In this study, particular attention is paid to the cross isentropic transport of the tracer in order to accurately assess the ozone loss. The assimilated SMILES ozone fields agree well with the limitation of noise induced variability within the SMR fields despite the limited latitude coverage of the SMILES observations. Ozone depletion has been derived by comparing the ozone field acquired by sequential assimilation with a passively transported ozone field initialized on 1 December 2009. Significant ozone loss was found in different periods and altitudes from using both SMILES and SMR data: The initial depletion occurred at the end of January below 550 K with an accumulated loss of 0.6–1.0 ppmv (approximately 20%) by 1 April. The ensuing loss started from the end of February between 575 K and 650 K. Our estimation shows that 0.8–1.3 ppmv (20–25 %) of O 3 has been removed at the 600 K isentropic level by 1 April in volume mixing ratio (VMR).

Description: Henry's law constants of polyols Atmospheric Chemistry and Physics, 14, 12815-12837, 2014 Author(s): S. Compernolle and J.-F. Müller Henry's law constants (HLC) are derived for several polyols bearing between 2 and 6 hydroxyl groups, based on literature data for water activity, vapour pressure and/or solubility. While deriving HLC and depending on the case, also infinite dilution activity coefficients (IDACs), solid state vapour pressures or activity coefficient ratios are obtained as intermediate results. An error analysis on the intermediate quantities and the obtained HLC is included. For most compounds, these are the first values reported, while others compare favourably with literature data in most cases. Using these values and those from a previous work (Compernolle and Müller, 2014), an assessment is made on the partitioning of polyols, diacids and hydroxy acids to droplet and aqueous aerosol.

Description: Nitrous oxide emissions from a commercial cornfield ( Zea mays ) measured using the eddy covariance technique Atmospheric Chemistry and Physics, 14, 12839-12854, 2014 Author(s): H. Huang, J. Wang, D. Hui, D. R. Miller, S. Bhattarai, S. Dennis, D. Smart, T. Sammis, and K. C. Reddy Increases in observed atmospheric concentrations of the long-lived greenhouse gas nitrous oxide (N 2 O) have been well documented. However, information on event-related instantaneous emissions during fertilizer applications is lacking. With the development of fast-response N 2 O analyzers, the eddy covariance (EC) technique can be used to gather instantaneous measurements of N 2 O concentrations to quantify the exchange of nitrogen between the soil and atmosphere. The objectives of this study were to evaluate the performance of a new EC system, to measure the N 2 O flux with the system, and finally to examine relationships of the N 2 O flux with soil temperature, soil moisture, precipitation, and fertilization events. An EC system was assembled with a sonic anemometer and a fast-response N 2 O analyzer (quantum cascade laser spectrometer) and applied in a cornfield in Nolensville, Tennessee during the 2012 corn growing season (4 April–8 August). Fertilizer amounts totaling 217 kg N ha −1 were applied to the experimental site. Results showed that this N 2 O EC system provided reliable N 2 O flux measurements. The cumulative emitted N 2 O amount for the entire growing season was 6.87 kg N 2 O-N ha −1 . Seasonal fluxes were highly dependent on soil moisture rather than soil temperature. This study was one of the few experiments that continuously measured instantaneous, high-frequency N 2 O emissions in crop fields over a growing season of more than 100 days.

Description: On the hiatus in the acceleration of tropical upwelling since the beginning of the 21st century Atmospheric Chemistry and Physics, 14, 12803-12814, 2014 Author(s): J. Aschmann, J. P. Burrows, C. Gebhardt, A. Rozanov, R. Hommel, M. Weber, and A. M. Thompson Chemistry–climate models predict an acceleration of the upwelling branch of the Brewer–Dobson circulation as a consequence of increasing global surface temperatures, resulting from elevated levels of atmospheric greenhouse gases. The observed decrease of ozone in the tropical lower stratosphere during the last decades of the 20th century is consistent with the anticipated acceleration of upwelling. However, more recent satellite observations of ozone reveal that this decrease has unexpectedly stopped in the first decade of the 21st century, challenging the implicit assumption of a continuous acceleration of tropical upwelling. In this study we use three decades of chemistry-transport-model simulations (1980–2013) to investigate this phenomenon and resolve this apparent contradiction. Aside from a high-bias between 1985–1990, our model is able to reproduce the observed tropical lower stratosphere ozone record. A regression analysis identifies a significant decrease in the early period followed by a statistically robust trend-change after 2002, in qualitative agreement with the observations. We demonstrate that this trend-change is correlated with structural changes in the vertical transport, represented in the model by diabatic heating rates taken from the reanalysis product Era-Interim. These changes lead to a hiatus in the acceleration of tropical upwelling between 70–30 hPa and a southward shift of the tropical pipe at 30 and 100 hPa during the past decade, which appear to be the primary causes for the observed trend-change in ozone.

Description: Emissions of biogenic volatile organic compounds and subsequent photochemical production of secondary organic aerosol in mesocosm studies of temperate and tropical plant species Atmospheric Chemistry and Physics, 14, 12781-12801, 2014 Author(s): K. P. Wyche, A. C. Ryan, C. N. Hewitt, M. R. Alfarra, G. McFiggans, T. Carr, P. S. Monks, K. L. Smallbone, G. Capes, J. F. Hamilton, T. A. M. Pugh, and A. R. MacKenzie Silver birch ( Betula pendula ) and three Southeast Asian tropical plant species ( Ficus cyathistipula , Ficus benjamina and Caryota millis ) from the pantropical fig and palm genera were grown in a purpose-built and environment-controlled whole-tree chamber. The volatile organic compounds emitted from these trees were characterised and fed into a linked photochemical reaction chamber where they underwent photo-oxidation under a range of controlled conditions (relative humidity or RH ~65–89%, volatile organic compound-to-NO x or VOC / NO x ~3–9 and NO x ~2 ppbV). Both the gas phase and the aerosol phase of the reaction chamber were monitored in detail using a comprehensive suite of on-line and off-line chemical and physical measurement techniques. Silver birch was found to be a high monoterpene and sesquiterpene but low isoprene emitter, and its emissions were observed to produce measurable amounts of secondary organic aerosol (SOA) via both nucleation and condensation onto pre-existing seed aerosol ( Y SOA 26–39%). In contrast, all three tropical species were found to be high isoprene emitters with trace emissions of monoterpenes and sesquiterpenes. In tropical plant experiments without seed aerosol there was no measurable SOA nucleation, but aerosol mass was shown to increase when seed aerosol was present. Although principally isoprene emitting, the aerosol mass produced from tropical fig was mostly consistent (i.e. in 78 out of 120 aerosol mass calculations using plausible parameter sets of various precursor specific yields) with condensation of photo-oxidation products of the minor volatile organic compounds (VOCs) co-emitted; no significant aerosol yield from condensation of isoprene oxidation products was required in the interpretations of the experimental results. This finding is in line with previous reports of organic aerosol loadings consistent with production from minor biogenic VOCs co-emitted with isoprene in principally isoprene-emitting landscapes in Southeast Asia. Moreover, in general the amount of aerosol mass produced from the emissions of the principally isoprene-emitting plants was less than would be expected from published single-VOC experiments, if co-emitted species were solely responsible for the final SOA mass. Interpretation of the results obtained from the fig data sets leaves room for a potential role for isoprene in inhibiting SOA formation under certain ambient atmospheric conditions, although instrumental and experimental constraints impose a level of caution in the interpretation of the results. Concomitant gas- and aerosol-phase composition measurements also provide a detailed overview of numerous key oxidation mechanisms at work within the systems studied, and their combined analysis provides insight into the nature of the SOA formed.

Description: Sources of black carbon aerosols in South Asia and surrounding regions during the Integrated Campaign for Aerosols, Gases and Radiation Budget (ICARB) Atmospheric Chemistry and Physics Discussions, 14, 30727-30759, 2014 Author(s): R. Kumar, M. C. Barth, V. S. Nair, G. G. Pfister, S. S. Babu, S. K. Satheesh, K. K. Moorthy, and G. R. Carmichael The dominant sources of black carbon (BC) in South Asia and surrounding regions are inferred during March–May 2006 (Integrated Campaign for Aerosols, Gases and Radiation Budget, ICARB) period by introducing BC tracers in the Weather Research and Forecasting Model coupled with Chemistry. The model reproduced the magnitude, temporal and spatial variability of BC distribution observed during the ICARB ship-cruise. Average and SD (representing the spatial and temporal variability) in observed and modeled BC mass concentrations along the ship-track are estimated as 755 ± 734 and 732 ± 913 ng m −3 respectively. Average modeled values at most of the inland stations were also found to fall within the range of observed values. Model results show that ICARB measurements were fairly well representative of the Arabian Sea and the Bay of Bengal during the pre-monsoon season. Results show that anthropogenic and biomass burning emissions, respectively, accounted for 70 and 28 % of the average ± SD BC mass concentration (1480 ± 5920 ng m −3 ) in South Asia. BC emissions from residential (49 %) and industrial (37 %) sectors appear to be the major anthropogenic sources, except in the Himalayas where vehicular emissions dominated. We find that, while all parts of continental India contributed to anthropogenic BC over the Bay of Bengal, contribution over the Arabian Sea came mostly from southern Peninsula. We also show that regional-scale transport of anthropogenic emissions contribute up to 30 % of BC mass concentrations in western and eastern India, suggesting that it is important to consider non-local sources along with the local emissions while designing strategies for mitigating BC emissions.

Description: Rapid transport of East Asian pollution to the deep tropics Atmospheric Chemistry and Physics Discussions, 14, 30705-30726, 2014 Author(s): M. J. Ashfold, J. A. Pyle, A. D. Robinson, M. S. M. Nadzir, S. M. Phang, A. A. Samah, S. Ong, H. E. Ung, L. K. Peng, S. E. Yong, and N. R. P. Harris Anthropogenic emissions from East Asia have increased over recent decades, and under the prevailing westerly winds, these increases have led to changes in atmospheric composition as far afield as North America. Here we show that, during Northern Hemisphere (NH) winter, pollution originating in East Asia also directly affects atmospheric composition in the deep tropics. We present observations of marked intra-seasonal variability in the anthropogenic tracer perchloroethene (C 2 Cl 4 ) collected at two locations in Borneo during the NH winter of 2008/09. We use the NAME trajectory model to show that the observed enhancements in C 2 Cl 4 mixing ratio are caused by rapid meridional transport, in the form of "cold surges", from the relatively polluted East Asian land mass. In these events air masses can move across 〉 30° of latitude in 4 days. We then present data from the Monitoring Atmospheric Composition and Climate reanalysis which suggests that air masses high in C 2 Cl 4 may also contain levels of the pollutants carbon monoxide and ozone that are approximately double the typical "background" levels in Borneo. Convection in Southeast Asia can be enhanced by cold surges, and further trajectory calculations indicate that the polluted air masses can subsequently be lifted to the tropical upper troposphere. This suggests a potentially important connection between mid-latitude pollution sources and the very low stratosphere.

Description: Deployment of a sequential two-photon laser-induced fluorescence sensor for the detection of gaseous elemental mercury at ambient levels: fast, specific, ultrasensitive detection with parts-per-quadrillion sensitivity Atmospheric Measurement Techniques, 7, 4251-4265, 2014 Author(s): D. Bauer, S. Everhart, J. Remeika, C. Tatum Ernest, and A. J. Hynes The operation of a laser-based sensor for gas-phase elemental mercury, Hg(0), is described. It utilizes sequential two-photon laser excitation with detection of blue-shifted laser-induced fluorescence (LIF) to provide a highly specific detection scheme that precludes detection of anything other than atomic mercury. It has high sensitivity, fast temporal resolution, and can be deployed for in situ measurements in the open atmosphere with essentially no perturbation of the environment. An ambient sample can also be pulled through a fluorescence cell, allowing for standard addition calibrations of the concentration. No type of preconcentration is required and there appears to be no significant interferences from other atmospheric constituents, including gas-phase oxidized mercury species. As a consequence, it is not necessary to remove oxidized mercury, commonly referred to as reactive gaseous mercury (RGM), from the air sample. The instrument has been deployed as part of an instrument intercomparison and compares well with conventional instrumentation that utilizes preconcentration on gold followed by analysis using cold-vapor atomic fluorescence spectroscopy (CVAFS). Currently, the achievable detection sensitivity is ~ 15 pg m −3 (~ 5 × 10 4 atoms cm −3 , ~ 2 ppq) at a sampling rate of 0.1 Hz, i.e., averaging 100 shots with a 10 Hz laser system. Preliminary results are described for a 50 Hz instrument that utilizes a modified excitation sequence and has monitored ambient elemental mercury with an effective sampling rate of 10 Hz. Additional work is required to produce the precision necessary to perform eddy correlation measurements. Addition of a pyrolysis channel should allow for the measurement of total gaseous mercury (TGM) and hence RGM (by difference) with good sensitivity and time resolution.

Description: Surface-to-mountaintop transport characterised by radon observations at the Jungfraujoch Atmospheric Chemistry and Physics, 14, 12763-12779, 2014 Author(s): A. D. Griffiths, F. Conen, E. Weingartner, L. Zimmermann, S. D. Chambers, A. G. Williams, and M. Steinbacher Atmospheric composition measurements at Jungfraujoch are affected intermittently by boundary-layer air which is brought to the station by processes including thermally driven (anabatic) mountain winds. Using observations of radon-222, and a new objective analysis method, we quantify the land-surface influence at Jungfraujoch hour by hour and detect the presence of anabatic winds on a daily basis. During 2010–2011, anabatic winds occurred on 40% of days, but only from April to September. Anabatic wind days were associated with warmer air temperatures over a large fraction of Europe and with a shift in air-mass properties, even when comparing days with a similar mean radon concentration. Excluding days with anabatic winds, however, did not lead to a better definition of the unperturbed aerosol background than a definition based on radon alone. This implies that a radon threshold reliably excludes local influences from both anabatic and non-anabatic vertical-transport processes.

Description: Prerequisites for application of hyperbolic relaxed eddy accumulation on managed grasslands and alternative net ecosystem exchange flux partitioning Atmospheric Measurement Techniques, 7, 4237-4250, 2014 Author(s): M. Riederer, J. Hübner, J. Ruppert, W. A. Brand, and T. Foken Relaxed eddy accumulation is still applied in ecosystem sciences for measuring trace gas fluxes. On managed grasslands, the length of time between management events and the application of relaxed eddy accumulation has an essential influence on the determination of the proportionality factor b and thus on the resulting flux. In this study this effect is discussed for the first time. Also, scalar similarity between proxy scalars and scalars of interest is affected until the ecosystem has completely recovered. Against this background, CO 2 fluxes were continuously measured and 13 CO 2 isofluxes were determined with a high measurement precision on two representative days in summer 2010. Moreover, a common method for the partitioning of the net ecosystem exchange into assimilation and respiration based on temperature and light response was compared with an isotopic approach directly based on the isotope discrimination of the biosphere. This approach worked well on the grassland site and could enhance flux partitioning results by better reproducing the environmental conditions.

Description: Extending the satellite data record of tropospheric ozone profiles from Aura-TES to MetOp-IASI: characterisation of optimal estimation retrievals Atmospheric Measurement Techniques, 7, 4223-4236, 2014 Author(s): H. Oetjen, V. H. Payne, S. S. Kulawik, A. Eldering, J. Worden, D. P. Edwards, G. L. Francis, H. M. Worden, C. Clerbaux, J. Hadji-Lazaro, and D. Hurtmans We apply the Tropospheric Emission Spectrometer (TES) ozone retrieval algorithm to Infrared Atmospheric Sounding Instrument (IASI) radiances and characterise the uncertainties and information content of the retrieved ozone profiles. This study focuses on mid-latitudes for the year 2008. We validate our results by comparing the IASI ozone profiles to ozone sondes. In the sonde comparisons, we find a negative bias (1–10%) in the IASI profiles in the lower to mid-troposphere and a positive bias (up to 14%) in the upper troposphere/lower stratosphere (UTLS) region. For the described cases, the degrees of freedom for signal are on average 3.2, 0.3, 0.8, and 0.9 for the columns 0 km – top of atmosphere, (0–6), (0–11), and (8–16) km, respectively. We find that our biases with respect to sondes and our degrees of freedom for signal for ozone are comparable to previously published results from other IASI ozone algorithms. In addition to evaluating biases, we validate the retrieval errors by comparing predicted errors to the sample covariance matrix of the IASI observations themselves. For the predicted versus empirical error comparison, we find that these errors are consistent and that the measurement noise and the interference of temperature and water vapour on the retrieval together mostly explain the empirically derived random errors. In general, the precision of the IASI ozone profiles is better than 20%.

Description: Solar irradiances measured using SPN1 radiometers: uncertainties and clues for development Atmospheric Measurement Techniques, 7, 4267-4283, 2014 Author(s): J. Badosa, J. Wood, P. Blanc, C. N. Long, L. Vuilleumier, D. Demengel, and M. Haeffelin The fast development of solar radiation and energy applications, such as photovoltaic and solar thermodynamic systems, has increased the need for solar radiation measurement and monitoring, for not only the global but also the diffuse and direct components. End users look for the best compromise between getting close to state-of-the-art measurements and keeping low capital, maintenance and operating costs. Among the existing commercial options, SPN1 is a relatively low cost solar radiometer that estimates global and diffuse solar irradiances from seven thermopile sensors under a shading mask and without moving parts. This work presents a comprehensive study of SPN1 accuracy and sources of uncertainty, drawing on laboratory experiments, numerical modelling and comparison studies between measurements from this sensor and state-of-the art instruments for six diverse sites. Several clues are provided for improving the SPN1 accuracy and agreement with state-of-the art measurements.

Description: Differences in aerosol absorption Ångström exponents between correction algorithms for a particle soot absorption photometer measured on the South African Highveld Atmospheric Measurement Techniques, 7, 4285-4298, 2014 Author(s): J. Backman, A. Virkkula, V. Vakkari, J. P. Beukes, P. G. Van Zyl, M. Josipovic, S. Piketh, P. Tiitta, K. Chiloane, T. Petäjä, M. Kulmala, and L. Laakso Absorption Ångström exponents (AAEs) calculated from filter-based absorption measurements are often used to give information on the origin of the ambient aerosol, for example, to distinguish between urban pollution and biomass burning aerosol. Filter-based absorption measurements are widely used and are common at aerosol monitoring stations globally. Several correction algorithms are used to account for artefacts associated with filter-based absorption techniques. These algorithms are of profound importance when determining the absolute amount of absorption by the aerosol. However, this study shows that there are substantial differences between the AAEs calculated from these corrections. Depending on the used correction, AAEs can change by as much as 46%. The study also highlights that the difference between AAEs calculated using different corrections can lead to conflicting conclusions on the type of aerosol when using the same data set. The AAE ranged between 1.17 for non-corrected data to 1.96 for the correction that gave the greatest values. Furthermore, the study implies that the AAEs reported for a site depend on at which filter transmittance the filter is changed. In this work, the AAEs were calculated from data measured with a three-wavelength particle soot absorption photometer (PSAP) at Elandsfontein on the South African Highveld for 23 months. The sample air of the PSAP was diluted to prolong filter change intervals, by a factor of 15. The correlation coefficient between the dilution-corrected PSAP and a non-diluted Multi-Angle Absorption Photometer (MAAP) was 0.9. Thus, the study also shows that the applicability of the PSAP can be extended to remote sites that are not often visited or suffer from high levels of pollution.

Description: Tropical troposphere to stratosphere transport of carbon monoxide and long-lived trace species in the Chemical Lagrangian Model of the Stratosphere (CLaMS) Geoscientific Model Development, 7, 2895-2916, 2014 Author(s): R. Pommrich, R. Müller, J.-U. Grooß, P. Konopka, F. Ploeger, B. Vogel, M. Tao, C. M. Hoppe, G. Günther, N. Spelten, L. Hoffmann, H.-C. Pumphrey, S. Viciani, F. D'Amato, C. M. Volk, P. Hoor, H. Schlager, and M. Riese Variations in the mixing ratio of trace gases of tropospheric origin entering the stratosphere in the tropics are of interest for assessing both troposphere to stratosphere transport fluxes in the tropics and the impact of these transport fluxes on the composition of the tropical lower stratosphere. Anomaly patterns of carbon monoxide (CO) and long-lived tracers in the lower tropical stratosphere allow conclusions about the rate and the variability of tropical upwelling to be drawn. Here, we present a simplified chemistry scheme for the Chemical Lagrangian Model of the Stratosphere (CLaMS) for the simulation, at comparatively low numerical cost, of CO, ozone, and long-lived trace substances (CH 4 , N 2 O, CCl 3 F (CFC-11), CCl 2 F 2 (CFC-12), and CO 2 ) in the lower tropical stratosphere. For the long-lived trace substances, the boundary conditions at the surface are prescribed based on ground-based measurements in the lowest model level. The boundary condition for CO in the lower troposphere (below about 4 km) is deduced from MOPITT measurements. Due to the lack of a specific representation of mixing and convective uplift in the troposphere in this model version, enhanced CO values, in particular those resulting from convective outflow are underestimated. However, in the tropical tropopause layer and the lower tropical stratosphere, there is relatively good agreement of simulated CO with in situ measurements (with the exception of the TROCCINOX campaign, where CO in the simulation is biased low ≈10–15 ppbv). Further, the model results (and therefore also the ERA-Interim winds, on which the transport in the model is based) are of sufficient quality to describe large scale anomaly patterns of CO in the lower stratosphere. In particular, the zonally averaged tropical CO anomaly patterns (the so called "tape recorder" patterns) simulated by this model version of CLaMS are in good agreement with observations, although the simulations show a too rapid upwelling compared to observations as a consequence of the overestimated vertical velocities in the ERA-Interim reanalysis data set. Moreover, the simulated tropical anomaly patterns of N 2 O are in good agreement with observations. In the simulations, anomaly patterns of CH 4 and CFC-11 were found to be very similar to those of N 2 O; for all long-lived tracers, positive anomalies are simulated because of the enhanced tropical upwelling in the easterly shear phase of the quasi-biennial oscillation.

Description: Moraine-dammed lake failures in Patagonia and assessment of outburst susceptibility in the Baker Basin Natural Hazards and Earth System Science, 14, 3243-3259, 2014 Author(s): P. Iribarren Anacona, K.P. Norton, and A. Mackintosh Glacier retreat since the Little Ice Age has resulted in the development or expansion of hundreds of glacial lakes in Patagonia. Some of these lakes have produced large (≥ 10 6 m 3 ) Glacial Lake Outburst Floods (GLOFs) damaging inhabited areas. GLOF hazard studies in Patagonia have been mainly based on the analysis of short-term series (≤ 50 years) of flood data and until now no attempt has been made to identify the relative susceptibility of lakes to failure. Power schemes and associated infrastructure are planned for Patagonian basins that have historically been affected by GLOFs, and we now require a thorough understanding of the characteristics of dangerous lakes in order to assist with hazard assessment and planning. In this paper, the conditioning factors of 16 outbursts from moraine-dammed lakes in Patagonia were analysed. These data were used to develop a classification scheme designed to assess outburst susceptibility, based on image classification techniques, flow routine algorithms and the Analytical Hierarchy Process. This scheme was applied to the Baker Basin, Chile, where at least seven moraine-dammed lakes have failed in historic time. We identified 386 moraine-dammed lakes in the Baker Basin of which 28 were classified with high or very high outburst susceptibility. Commonly, lakes with high outburst susceptibility are in contact with glaciers and have moderate (〉 8°) to steep (〉 15°) dam outlet slopes, akin to failed lakes in Patagonia. The proposed classification scheme is suitable for first-order GLOF hazard assessments in this region. However, rapidly changing glaciers in Patagonia make detailed analysis and monitoring of hazardous lakes and glaciated areas upstream from inhabited areas or critical infrastructure necessary, in order to better prepare for hazards emerging from an evolving cryosphere.

Description: Seismology of the Oso-Steelhead landslide Natural Hazards and Earth System Sciences Discussions, 2, 7309-7327, 2014 Author(s): C. Hibert, C. P. Stark, and G. Ekström We carry out a combined analysis of the short- and long-period seismic signals generated by the devastating Oso-Steelhead landslide that occurred on 22 March 2014. The seismic records show that the Oso-Steelhead landslide was not a single slope failure, but a succession of multiple failures distinguished by two major collapses that occurred approximately three minutes apart. The first generated long-period surface waves that were recorded at several proximal stations. We invert these long-period signals for the forces acting at the source, and obtain estimates of the first failure runout and kinematics, as well as its mass after calibration against the mass-center displacement estimated from remote-sensing imagery. Short-period analysis of both events suggests that the source dynamics of the second are more complex than the first. No distinct long-period surface waves were recorded for the second failure, which prevents inversion for its source parameters. However, by comparing the seismic energy of the short-period waves generated by both events we are able to estimate the volume of the second. Our analysis suggests that the volume of the second failure is about 15–30% of the total landslide volume, which is in agreement with ground observations.

Description: Dependence of sandpile avalanche frequency–size distribution on coverage extent and compactness of embedded toppling threshold heterogeneity: implications for the variation of Gutenberg–Richter b value Nonlinear Processes in Geophysics, 21, 1185-1193, 2014 Author(s): L.-Y. Chiao and Q. Liu The effects of the spatiotemporal evolution of failure threshold heterogeneity on the dynamics of fault criticality, and thus on regional seismogenesis, have attracted strong interest in the field of regional seismotectonics. The heterogeneity might be a manifestation of the macroscopic distribution and multiscale strength variation of asperities, the distinct regional stress level, and (microscopically) heterogeneous fault surface roughness or friction regimes. In this study, rather than attempting to mimic the complex microscale slipping physics on a fault surface, sandpile cellular automata were implemented with a straightforward toppling rule. The objective is to examine the influence of distinct configurations of the embedded heterogeneous toppling threshold field on the global system avalanche event statistics. The examination results revealed that increasing the coverage extent and decreasing the compactness of the heterogeneous failure threshold, rather than the magnitude, range of contrast, diversity, or the geometric configuration of the threshold heterogeneity, leads to a systematic increase in the scaling exponent of the avalanche event power law statistics, implying the importance of mutual interaction among toppling sites with distinct thresholds. For tectonic provinces with differing stress regimes evolving spatio temporally, it is postulated that the distinct extent and compactness of the heterogeneous failure threshold are critical factors that manifest in the reported dynamic variations of seismicity scaling.

Description: Evaluation of empirical mode decomposition for quantifying multi-decadal variations in sea level records Nonlinear Processes in Geophysics Discussions, 1, 1833-1854, 2014 Author(s): D. P. Chambers The ability of empirical mode decomposition (EMD) to extract multidecadal variability from sea level records is tested using three simulations: one based on a series of purely sinusoidal modes, one based on scaled climate indices of El Niño and the Pacific Decadal Oscillation (PDO), and the final one including a single month with an extreme sea level event. All simulations include random noise of similar variance to high-frequency variability in the San Francisco tide gauge record. The intrinsic mode functions (IMFs) computed using EMD were compared to the prescribed oscillations. In all cases, the longest-period modes are significantly distorted, with incorrect amplitudes and phases. This affects the estimated acceleration computed from the longest periodic IMF. In these simulations, the acceleration was underestimated in the case with purely sinusoidal modes, and overestimated by nearly 100% in the case with prescribed climate modes. Additionally, in all cases, extra low-frequency modes uncorrelated with the prescribed variability are found. These experiments suggest that using EMD to identify multidecadal variability and accelerations in sea level records should be used with caution.

Description: Benjamin–Feir instability of waves in the presence of current Nonlinear Processes in Geophysics Discussions, 1, 1803-1832, 2014 Author(s): I. V. Shugan, H. H. Hwung, and R. Y. Yang The development of Benjamin–Feir instability of Stokes waves in the presence of variable current is presented. We employ a model of a resonance system having three coexisting nonlinear waves and nonuniform current. The model is free from the narrow-band approximation for surface waves and relatively weak adverse current. The modulation instability of Stokes waves in nonuniform moving media has special properties. Interaction with countercurrent accelerates the growth of sideband modes on a short spatial scale. An increase in initial wave steepness intensifies the wave energy exchange accompanied by wave breaking dissipation, results in asymmetry of sideband modes and a frequency downshift with an energy transfer jump to the lower sideband mode, and depresses the higher sideband and carrier wave. Nonlinear waves may even overpass the blocking barrier produced by strong adverse current. The frequency downshift of the energy peak is permanent and the system does not revert to its initial state. We find reasonable correspondence between the results of model simulations and available experimental results for wave interaction with blocking opposing current. Large transient or freak waves with amplitude and steepness several times those of normal waves may form during temporal nonlinear focusing of the resonant waves accompanied by energy income from sufficiently strong opposing current. We employ the resonance model for the estimation of the maximum amplification of wave amplitudes as a function of gradually increasing opposing current and compare the result obtained with recently published experimental results and modeling results obtained with the nonlinear Schrödinger equation.

Description: A multi-layer land surface energy budget model for implicit coupling with global atmospheric simulations Geoscientific Model Development Discussions, 7, 8649-8701, 2014 Author(s): J. Ryder, J. Polcher, P. Peylin, C. Ottlé, Y. Chen, E. van Gorsel, V. Haverd, M. J. McGrath, K. Naudts, J. Otto, A. Valade, and S. Luyssaert In Earth system modelling, a description of the energy budget of the vegetated surface layer is fundamental as it determines the meteorological conditions in the planetary boundary layer and as such contributes to the atmospheric conditions and its circulation. The energy budget in most Earth system models has long been based on a "big-leaf approach", with averaging schemes that represent in-canopy processes. Such models have difficulties in reproducing consistently the energy balance in field observations. We here outline a newly developed numerical model for energy budget simulation, as a component of the land surface model ORCHIDEE-CAN (Organising Carbon and Hydrology In Dynamic Ecosystems – CANopy). This new model implements techniques from single-site canopy models in a practical way. It includes representation of in-canopy transport, a multilayer longwave radiation budget, height-specific calculation of aerodynamic and stomatal conductance, and interaction with the bare soil flux within the canopy space. Significantly, it avoids iterations over the height of tha canopy and so maintains implicit coupling to the atmospheric model LMDz. As a first test, the model is evaluated against data from both an intensive measurement campaign and longer term eddy covariance measurements for the intensively studied Eucalyptus stand at Tumbarumba, Australia. The model performs well in replicating both diurnal and annual cycles of fluxes, as well as the gradients of sensible heat fluxes. However, the model overestimates sensible heat flux against an underestimate of the radiation budget. Improved performance is expected through the implementation of a more detailed calculation of stand albedo and a more up-to-date stomatal conductance calculation.

Description: The numerical simulation on ionospheric perturbations in electric field before large earthquakes Annales Geophysicae, 32, 1487-1493, 2014 Author(s): S. F. Zhao, X. M. Zhang, Z. Y. Zhao, and X. H. Shen Many observational results have shown electromagnetic abnormality in the ionosphere before large earthquakes. The theoretical simulation can help us to understand the internal mechanism of these anomalous electromagnetic signals resulted from seismic regions. In this paper, the horizontal and vertical components of electric and magnetic field at the topside ionosphere are simulated by using the full wave method that is based on an improved transfer matrix method in the lossy anisotropic horizontally stratified ionosphere. Taken account into two earthquakes with electric field perturbations recorded by the DEMETER satellite, the numerical results reveal that the propagation and penetration of ULF (ultra-low-frequency) electromagnetic waves into the ionosphere is related to the spatial distribution of electron and ion densities at different time and locations, in which the ion density has less effect than electron density on the field intensity. Compared with different frequency signals, the minimum values of electric and magnetic field excited by earthquakes can be detected by satellite in current detection capability have also been calculated, and the lower frequency wave can be detected easier.

Description: Rare temperature histories and cirrus ice number density in a parcel and a one-dimensional model Atmospheric Chemistry and Physics, 14, 13013-13022, 2014 Author(s): D. M. Murphy A parcel and a one-dimensional model are used to investigate the temperature dependence of ice crystal number density. The number of ice crystals initially formed in a cold cirrus cloud is very sensitive to the nucleation mechanism and the detailed history of cooling rates during nucleation. A possible small spread in the homogeneous freezing threshold due to varying particle composition is identified as a sensitive nucleation parameter. In a parcel model, the slow growth rate of ice crystals at low temperatures inherently leads to a strong increase in ice number density at low temperatures. This temperature dependence is not observed. The model temperature dependence occurs for a wide range of assumptions and for either homogeneous or, less strongly, heterogeneous freezing. However, the parcel model also shows that random temperature fluctuations result in an extremely wide range of ice number densities. A one-dimensional model is used to show that the rare temperature trajectories resulting in the lowest number densities are disproportionately important. Low number density ice crystals sediment and influence a large volume of air. When such fall streaks are included, the ice number becomes less sensitive to the details of nucleation than it is in a parcel model. The one-dimensional simulations have a more realistic temperature dependence than the parcel mode. The one-dimensional model also produces layers with vertical dimensions of meters even if the temperature forcing has a much broader vertical wavelength. Unlike warm clouds, cirrus clouds are frequently surrounded by supersaturated air. Sedimentation through supersaturated air increases the importance of any process that produces small numbers of ice crystals. This paper emphasizes the relatively rare temperature trajectories that produce the fewest crystals. Other processes are heterogeneous nucleation, sedimentation from the very bottom of clouds, annealing of disordered to hexagonal ice, and entrainment.

Description: An improved dust emission model – Part 1: Model description and comparison against measurements Atmospheric Chemistry and Physics, 14, 13023-13041, 2014 Author(s): J. F. Kok, N. M. Mahowald, G. Fratini, J. A. Gillies, M. Ishizuka, J. F. Leys, M. Mikami, M.-S. Park, S.-U. Park, R. S. Van Pelt, and T. M. Zobeck Simulations of the dust cycle and its interactions with the changing Earth system are hindered by the empirical nature of dust emission parameterizations in weather and climate models. Here we take a step towards improving dust cycle simulations by using a combination of theory and numerical simulations to derive a physically based dust emission parameterization. Our parameterization is straightforward to implement into large-scale models, as it depends only on the wind friction velocity and the soil's threshold friction velocity. Moreover, it accounts for two processes missing from most existing parameterizations: a soil's increased ability to produce dust under saltation bombardment as it becomes more erodible, and the increased scaling of the dust flux with wind speed as a soil becomes less erodible. Our treatment of both these processes is supported by a compilation of quality-controlled vertical dust flux measurements. Furthermore, our scheme reproduces this measurement compilation with substantially less error than the existing dust flux parameterizations we were able to compare against. A critical insight from both our theory and the measurement compilation is that dust fluxes are substantially more sensitive to the soil's threshold friction velocity than most current schemes account for.

Description: The impact of volcanic aerosol on the Northern Hemisphere stratospheric polar vortex: mechanisms and sensitivity to forcing structure Atmospheric Chemistry and Physics, 14, 13063-13079, 2014 Author(s): M. Toohey, K. Krüger, M. Bittner, C. Timmreck, and H. Schmidt Observations and simple theoretical arguments suggest that the Northern Hemisphere (NH) stratospheric polar vortex is stronger in winters following major volcanic eruptions. However, recent studies show that climate models forced by prescribed volcanic aerosol fields fail to reproduce this effect. We investigate the impact of volcanic aerosol forcing on stratospheric dynamics, including the strength of the NH polar vortex, in ensemble simulations with the Max Planck Institute Earth System Model. The model is forced by four different prescribed forcing sets representing the radiative properties of stratospheric aerosol following the 1991 eruption of Mt. Pinatubo: two forcing sets are based on observations, and are commonly used in climate model simulations, and two forcing sets are constructed based on coupled aerosol–climate model simulations. For all forcings, we find that simulated temperature and zonal wind anomalies in the NH high latitudes are not directly impacted by anomalous volcanic aerosol heating. Instead, high-latitude effects result from enhancements in stratospheric residual circulation, which in turn result, at least in part, from enhanced stratospheric wave activity. High-latitude effects are therefore much less robust than would be expected if they were the direct result of aerosol heating. Both observation-based forcing sets result in insignificant changes in vortex strength. For the model-based forcing sets, the vortex response is found to be sensitive to the structure of the forcing, with one forcing set leading to significant strengthening of the polar vortex in rough agreement with observation-based expectations. Differences in the dynamical response to the forcing sets imply that reproducing the polar vortex responses to past eruptions, or predicting the response to future eruptions, depends on accurate representation of the space–time structure of the volcanic aerosol forcing.

Description: Evaluation of a regional chemistry transport model using a newly developed regional OMI NO 2 retrieval Atmospheric Chemistry and Physics Discussions, 14, 31039-31090, 2014 Author(s): G. Kuhlmann, Y. F. Lam, H. M. Cheung, A. Hartl, J. C. H. Fung, P. W. Chan, and M. O. Wenig In this paper, we evaluate a high-resolution chemistry transport model (CTM) (3 km x 3 km spatial resolution) with the new Hong Kong (HK) NO 2 retrieval developed for the Ozone Monitoring Instrument (OMI) on-board the Aura satellite. The three-dimensional atmospheric chemistry was modelled in the Pearl River Delta (PRD) region in southern China by the Models-3 Community Multiscale Air Quality (CMAQ) modelling system from October 2006 to January 2007. In the HK NO 2 retrieval, tropospheric air mass factors (AMF) were recalculated using high-resolution ancillary parameters of surface reflectance, NO 2 profile shapes and aerosol profiles of which the latter two were taken from the CMAQ simulation. We also tested four different aerosol parametrizations. Ground level measurements by the PRD Regional Air Quality Monitoring (RAQM) network were used as additional independent measurements. The HK NO 2 retrieval increases the NO 2 vertical column densities (VCD) by (+31 ± 38) %, when compared to NASA's standard product (SP2), and reduces the mean bias (MB) between satellite and ground measurements by 26 percentage points from −41 to −15 %. The correlation coefficient r is low for both satellite datasets ( r = 0.35) due to the high spatial variability of NO 2 concentrations. The correlation between CMAQ and the RAQM network is low ( r ≈ 0.3) and the model underestimates the NO 2 concentrations in the north-western model domain (Foshan and Guangzhou). We compared the CMAQ NO 2 time series of the two main plumes with our regional OMI NO 2 product. The model overestimates the NO 2 VCDs by about 15 % in Hong Kong and Shenzhen, while the correlation coefficient is satisfactory ( r = 0.56). In Foshan and Guangzhou, the correlation is low ( r = 0.37) and the model underestimates the VCDs strongly (MB = −40 %). In addition, we estimated that the OMI VCDs are also underestimated by about 10 to 20 % in Foshan and Guangzhou because of the influence of the model parameters on the AMF. In this study, we demonstrate that the HK OMI NO 2 retrieval reduces the bias of the satellite measurements and thus the dataset can be used to study the magnitude of NO 2 concentrations in a regional model. The low bias can be achieved if AMFs are recalculated with more accurate surface reflectance, aerosol profiles and NO 2 profiles; only NO 2 profiles have been replaced in earlier studies. Since unbiased concentrations are important, for example, in air pollution studies, the results of this paper can be very helpful in future model evaluation studies.

Description: Aerosols over Continental Portugal (1978–1993): their sources and an impact on the regional climate Atmospheric Chemistry and Physics Discussions, 14, 31009-31038, 2014 Author(s): A. L. Morozova and I. A. Mironova Understanding of aerosol sources which affect climate is an important problem open in front of scientists as well as policymakers. The role of aerosols in local climate variability depends on a~balance between aerosol absorbing and scattering particles as well as on variability of environmental conditions. In this paper we investigate variability of aerosol content (both absorbing and scattering UV radiation) over Continental Portugal in dependence on aerosol sources (volcanic eruptions, dust events, wildfires and anthropogenic pollution). The effect of the aerosol on the climate is studied analyzing their contribution to variations of temperature, sunshine duration and precipitation over Portuguese regions. The present analysis is based on a developed modern multiple regression technique allowing us to build the statistical correlation models to determine both the main local aerosol sources and aerosol's influence on the climate of the Continental Portugal during 1978–1993 time period. The analysis allows us to conclude that the main sources driving the variations of the aerosol content over studied locations are wildfires, mineral dust intrusions and anthropogenic pollution. The relations between the aerosol content variations and the atmospheric parameters depend on the level of urbanization of the studied region, the type of aerosol and the season. The most significant finding is the decrease of the daily temperature (and diurnal temperature range) related to the decrease of sunshine duration observed during the summer periods of increased content of the absorbing aerosols in the atmosphere.

Description: The role of organic condensation on ultrafine particle growth during nucleation events Atmospheric Chemistry and Physics Discussions, 14, 30761-30798, 2014 Author(s): D. Patoulias, C. Fountoukis, I. Riipinen, and S. N. Pandis A new aerosol dynamics model (DMANx) has been developed that simulates the aerosol size/composition distribution and includes the condensation of organic vapors on nanoparticles through the implementation of the recently developed Volatility Basis Set framework. Simulations were performed for Hyytiala (Finland) and Finokalia (Greece), two locations with different organic sources where detailed measurements were available to constrain the new model. We investigate the effect of condensation of organics and chemical aging reactions of secondary organic aerosol (OA) on ultrafine particle growth and particle number concentration. This work highlights the importance of the pathways of oxidation of biogenic volatile organic compounds and the production of extremely low-volatility organics. At Hyytiala, organic condensation dominates the growth process of new particles. The low-volatility secondary OA contributes to particle growth during the early growth stage, but after a few hours most of the growth is due to semi-volatile secondary OA. At Finokalia, simulations show that organics have a complementary role to new particle growth contributing 45% to the total mass of new particles. Condensation of organics increases the number concentration of particles that can act as CCN ( N 100 ) by 13% at Finokalia and 25% at Hyytiala. The sensitivity of our results to the surface tension used is discussed.

Description: Comparison of aeolian snow transport events and snow mass fluxes between observations and simulations made by the regional climate model MAR in Adélie Land, East Antarctica The Cryosphere Discussions, 8, 6007-6032, 2014 Author(s): A. Trouvilliez, H. Gallée, F. Naaim-Bouvet, C. Genthon, C. Amory, V. Favier, C. Agosta, L. Piard, and H. Bellot The regional climate model MAR including a coupled snow pack/aeolian snow transport parameterisation is compared with aeolian snow mass fluxes at a fine spatial resolution (5 km horizontally and 2 m vertically) and at a fine temporal resolution (30 min) over 1 month in Antarctica. Numerous feedbacks are taken into account in the MAR including the drag partitioning caused by the roughness elements. Wind speed is correctly simulated with a positive value of the Nash test (0.60 and 0.37) but the wind speeds above 10 m s −1 are underestimated. The aeolian snow transport events are correctly reproduced with a good temporal resolution except for the aeolian snow transport events with a particles' maximum height below 1 m. The simulated threshold friction velocity, calculated without snowfall, is overestimated. The simulated aeolian snow mass fluxes between 0 to 2 m have the same variations but are underestimated compared to the second-generation FlowCapt values and so is the simulated relative humidity at 2 m. This underestimation is not entirely due to the underestimation of the simulated wind speed. The MAR underestimates the aeolian snow quantity that pass through the first two meters by a factor ten compared to the second-generation FlowCapt value (13 990 kg m −1 and 151 509 kg m −1 respectively). It will conduct the MAR, with this parametrisation, to underestimate the effect of the aeolian snow transport on the Antarctic surface mass balance.

Description: Influence of weak layer heterogeneity and slab properties on slab tensile failure propensity and avalanche release area The Cryosphere Discussions, 8, 6033-6057, 2014 Author(s): J. Gaume, G. Chambon, N. Eckert, M. Naaim, and J. Schweizer Dry-snow slab avalanches are generally caused by a sequence of fracture processes including failure initiation in a weak snow layer underlying a cohesive slab followed by crack propagation within the weak layer (WL) and tensile fracture through the slab. During past decades, theoretical and experimental work has gradually improved our knowledge of the fracture process in snow. However, our limited understanding of crack propagation and fracture arrest propensity prevents the evaluation of avalanche release sizes and thus impedes hazard assessment. To address this issue, slab tensile failure propensity is examined using a mechanically-based statistical model of the slab–WL system based on the finite element method. This model accounts for WL heterogeneity, stress redistribution by elasticity of the slab and the slab possible tensile failure. Two types of avalanche release are distinguished in the simulations: (1) full-slope release if the heterogeneity is not sufficient to stop crack propagation and to trigger a tensile failure within the slab, (2) partial-slope release if fracture arrest and slab tensile failure occurs due to the WL heterogeneity. The probability of these two release types is presented as a function of the characteristics of WL heterogeneity and of the slab. One of the main outcomes is that, for realistic values of the parameters, the tensile failure propensity is mainly influenced by slab properties. Hard and thick snow slabs are more prone to wide-scale crack propagation and thus lead to larger avalanches (full-slope release). In this case, the avalanche size is mainly influenced by topographical and morphological features such as rocks, trees, slope curvature and the spatial variability of the snow depth as it is often claimed in the literature.

Description: Biogeochemical indicators of peatland degradation – a case study of a temperate bog in northern Germany Biogeosciences Discussions, 11, 16825-16854, 2014 Author(s): J. P. Krüger, J. Leifeld, S. Glatzel, S. Szidat, and C. Alewell Peatlands store a great proportion of the global soil carbon pool and can loose carbon via the atmosphere due to degradation. In Germany, most of the greenhouse gas emissions from organic soils are attributed to sites managed as grassland. Here we investigated a land-use gradient from near-natural wetland (NW) to an extensively managed (GE) to an intensively managed grassland site (GI), all formed in the same bog complex in northern Germany. Vertical depth profiles of δ 13 C, δ 15 N, ash content, C/N ratio, bulk density, as well as radiocarbon ages were studied to identify peat degradation and to calculate carbon loss. At all sites, including the near-natural site, δ 13 C depth profiles indicate aerobic decomposition in the upper horizons. Depth profiles of δ 15 Ndiffered significantly between sites with increasing δ 15 N values in the top layers with increasing intensity of use, indicating that the peat is more decomposed. At both grassland sites, the ash content peaked within the first centimeter. In the near-natural site, ash contents were highest in 10–60 cm depth. This indicates that not only the managed grasslands, but also the near-natural site, is influenced by anthropogenic activities, most likely due to the drainage of the surrounding area. However, we found very young peat material in the first centimeter of the NW, indicating recent peat growth. The NW site accumulates carbon today even though it is and probably was influenced by anthropogenic activities in the past indicated by δ 13 C and ash content depth profiles. Based on the enrichment of ash content and changes in bulk density, we calculated carbon loss from these sites in retrograde. As expected land use intensification leads to a higher carbon loss which is supported by the higher peat ages at the intensive managed grassland site. All investigated biogeochemical parameters together indicate degradation of peat due to (i) conversion to grassland, (ii) historical drainage as well as recent development and (iii) land use intensification.

Description: A model study on changes of European and Swiss particulate matter, ozone and nitrogen deposition between 1990 and 2020 due to the revised Gothenburg protocol Atmospheric Chemistry and Physics, 14, 13081-13095, 2014 Author(s): S. Aksoyoglu, J. Keller, G. Ciarelli, A. S. H. Prévôt, and U. Baltensperger We report a study of changes in air quality due to emission reductions using the chemical transport model CAMx. The model domain includes all of Europe with a nested domain over Switzerland. The model simulations were performed with emissions for 1990 (the reference year for the Gothenburg Protocol), 2005 (the reference year for the revised Gothenburg Protocol), 2006 (for model validation) and 2020 (the target year for the revised Gothenburg Protocol) using three emission scenarios prepared by IIASA/GAINS. Changes in ozone, particulate matter and nitrogen deposition are the central theme of the study. The modelled relative changes in the annual average PM 2.5 concentrations between 1990 and 2005 look reasonable based on various PM 10 and PM 2.5 observations in the past. The results obtained in this study suggest that annual mean concentrations of PM 2.5 decreased by about 20–50% in Europe. Simulations using the baseline scenario (BL 2020) suggest that PM 2.5 concentrations in 2020 will be about 30% lower than those in 2005. The largest predicted decrease in PM 2.5 , based on the MTFR (maximum technically feasible reduction) scenario, was about 60% and was located mainly in the eastern part of Europe. In the case of ozone, both model results and measurements show an increase in the mean ozone mixing ratios between 1990 and 2005. The observations, however, suggest a larger increase, indicating the importance of background ozone levels. Although emission reductions caused a decrease in peak ozone values, average ozone levels in polluted regions increased due to reduced titration with nitric oxide (NO). This caused a change in the frequency distribution of ozone. Model simulations using emission scenarios for 2020 suggest that annual average ozone mixing ratios will continue to increase. Changes in the levels of the damage indicators AOT40 for forests and SOMO35 are reported as well. The model results suggest that nitrogen deposition has decreased by 10–30% in the eastern part of Europe since 1990, while it has increased by about 20% in the Iberian Peninsula. The decrease is mainly due to the deposition of oxidized nitrogen species, whereas deposition of reduced nitrogen compounds increased. In Switzerland, nitrogen deposition is larger in the northern part of the Alps, where ammonia emissions are the highest. Applying the baseline scenario, we found that the deposition of oxidized nitrogen compounds will have decreased by a further 40% by 2020, whereas deposition of reduced species will continue to increase. This will lead to a 10–20% decrease in the total nitrogen deposition in most of the model domain, with a 10% increase in the eastern part of Europe.

Description: An improved dust emission model – Part 2: Evaluation in the Community Earth System Model, with implications for the use of dust source functions Atmospheric Chemistry and Physics, 14, 13043-13061, 2014 Author(s): J. F. Kok, S. Albani, N. M. Mahowald, and D. S. Ward The complex nature of mineral dust aerosol emission makes it a difficult process to represent accurately in weather and climate models. Indeed, results in the companion paper indicate that many large-scale models underestimate the dust flux's sensitivity to the soil's threshold friction velocity for erosion. We hypothesize that this finding explains why many dust cycle simulations are improved by using an empirical dust source function that shifts emissions towards the world's most erodible regions. Here, we both test this hypothesis and evaluate the performance of the new dust emission parameterization presented in the companion paper. We do so by implementing the new emission scheme into the Community Earth System Model (CESM) and comparing the resulting dust cycle simulations against an array of measurements. We find that the new scheme shifts emissions towards the world's most erodible regions in a manner that is strikingly similar to the effect of implementing a widely used source function based on satellite observations of dust source regions. Furthermore, model comparisons against aerosol optical depth measurements show that the new physically based scheme produces a statistically significant improvement in CESM's representation of dust emission, which exceeds the improvement produced by implementing a source function. These results indicate that the need to use an empirical source function is eliminated, at least in CESM, by the additional physics in the new scheme, and in particular by its increased sensitivity to the soil's threshold friction velocity. Since the threshold friction velocity is affected by climate changes, our results further suggest that many large-scale models underestimate the global dust cycle's climate sensitivity.

Description: Looking through the haze: evaluating the CALIPSO level 2 aerosol optical depth using airborne high spectral resolution lidar data Atmospheric Measurement Techniques, 7, 4317-4340, 2014 Author(s): R. R. Rogers, M. A. Vaughan, C. A. Hostetler, S. P. Burton, R. A. Ferrare, S. A. Young, J. W. Hair, M. D. Obland, D. B. Harper, A. L. Cook, and D. M. Winker The Cloud–Aerosol Lidar with Orthogonal Polarization (CALIOP) instrument onboard the Cloud–Aerosol Lidar and Pathfinder Satellite Observations (CALIPSO) spacecraft has provided over 8 yr of nearly continuous vertical profiling of Earth's atmosphere. In this paper we investigate the V3.01 and V3.02 CALIOP 532 nm aerosol layer optical depth (AOD) product (i.e the AOD of individual layers) and the column AOD product (i.e., the sum AOD of the complete column) using an extensive database of coincident measurements. The CALIOP AOD measurements and AOD uncertainty estimates are compared with collocated AOD measurements collected with the NASA High Spectral Resolution Lidar (HSRL) in the North American and Caribbean regions. In addition, the CALIOP aerosol lidar ratios are investigated using the HSRL measurements. In general, compared with the HSRL values, the CALIOP layer AOD are biased high by less than 50% for AOD 〈 0.3 with higher errors for higher AOD. Less than 60% of the HSRL AOD measurements are encompassed within the CALIOP layer 1 SD uncertainty range (around the CALIOP layer AOD), so an error estimate is created to encompass 68% of the HSRL data. Using this new metric, the CALIOP layer AOD error is estimated using the HSRL layer AOD as ±0.035 ± 0.05 · (HSRL layer AOD) at night and ±0.05 ± 0.05 · (HSRL layer AOD) during the daytime. Furthermore, the CALIOP layer AOD error is found to correlate with aerosol loading as well as aerosol subtype, with the AODs in marine and dust layers agreeing most closely with the HSRL values. The lidar ratios used by CALIOP for polluted dust, polluted continental, and biomass burning layers are larger than the values measured by the HSRL in the CALIOP layers, and therefore the AODs for these types retrieved by CALIOP were generally too large. We estimated the CALIOP column AOD error can be expressed as ±0.05 ± 0.07 · (HSRL column AOD) at night and ±0.08 ± 0.1 · (HSRL column AOD) during the daytime. Multiple sources of error contribute to both positive and negative errors in the CALIOP column AOD, including multiple layers in the column of different aerosol types, lidar ratio errors, cloud misclassification, and undetected aerosol layers. The undetected layers were further investigated and we found that the layer detection algorithm works well at night, although undetected aerosols in the free troposphere introduce a mean underestimate of 0.02 in the column AOD in the data set examined. The decreased signal-to-noise ratio (SNR) during the daytime led to poorer performance of the layer detection. This caused the daytime CALIOP column AOD to be less accurate than during the nighttime, because CALIOP frequently does not detect optically thin aerosol layers with AOD 〈 0.1. Given that the median vertical extent of aerosol detected within any column was 1.6 km during the nighttime and 1.5 km during the daytime, we can estimate the minimum extinction detection threshold to be 0.012 km −1 at night and 0.067 km −1 during the daytime in a layer median sense. This extensive validation of level 2 CALIOP AOD products extends previous validation studies to nighttime lighting conditions and provides independent measurements of the lidar ratio; thus, allowing the assessment of the effect on the CALIOP AOD of using inappropriate lidar ratio values in the extinction retrieval.

Description: The use of NO 2 absorption cross section temperature sensitivity to derive NO 2 profile temperature and stratospheric–tropospheric column partitioning from visible direct-sun DOAS measurements Atmospheric Measurement Techniques, 7, 4299-4316, 2014 Author(s): E. Spinei, A. Cede, W. H. Swartz, J. Herman, and G. H. Mount This paper presents a temperature sensitivity method (TESEM) to accurately calculate total vertical NO 2 column, atmospheric slant NO 2 profile-weighted temperature ( T ), and to separate stratospheric and tropospheric columns from direct-sun (DS), ground-based measurements using the retrieved T . TESEM is based on differential optical absorption spectroscopy (DOAS) fitting of the linear temperature-dependent NO 2 absorption cross section, σ ( T ), regression model (Vandaele et al., 2003). Separation between stratospheric and tropospheric columns is based on the primarily bimodal vertical distribution of NO 2 and an assumption that stratospheric effective temperature can be represented by temperature at 27 km ± 3 K, and tropospheric effective temperature is equal to surface temperature within 3–5 K. These assumptions were derived from the Global Modeling Initiative (GMI) chemistry-transport model (CTM) simulations over two northern midlatitude sites in 2011. TESEM was applied to the Washington State University Multi-Function DOAS instrument (MFDOAS) measurements at four midlatitude locations with low and moderate NO 2 anthropogenic emissions: (1) the Jet Propulsion Laboratory's Table Mountain Facility (JPL-TMF), CA, USA (34.38° N/117.68° W); (2) Pullman, WA, USA (46.73° N/117.17° W); (3) Greenbelt, MD, USA (38.99° N/76.84° W); and (4) Cabauw, the Netherlands (51.97° N/4.93° E) during July 2007, June–July 2009, July–August and October 2011, November 2012–May 2013, respectively. NO 2 T and total, stratospheric, and tropospheric NO 2 vertical columns were determined over each site.

Description: On the link between hygroscopicity, volatility, and oxidation state of ambient and water-soluble aerosol in the Southeastern United States Atmospheric Chemistry and Physics Discussions, 14, 30835-30877, 2014 Author(s): K. M. Cerully, A. Bougiatioti, J. R. Hite Jr., H. Guo, L. Xu, N. L. Ng, R. Weber, and A. Nenes The formation of secondary organic aerosol (SOA) combined with the partitioning of semi-volatile organic components can impact numerous aerosol properties including cloud condensation nuclei (CCN) activity, hygroscopicity and volatility. During the summer 2013 Southern Oxidant and Aerosol Study (SOAS) field campaign in a rural site in the Southeastern United States, a suite of instruments including a CCN counter, a thermodenuder (TD) and a high resolution time-of-flight aerosol mass spectrometer (AMS) were used to measure CCN activity, aerosol volatility, composition and oxidation state. Particles were either sampled directly from ambient or through a Particle Into Liquid Sampler (PILS), allowing the investigation of the water-soluble aerosol component. Ambient aerosol exhibited size-dependent composition with larger particles being more hygroscopic. The hygroscopicity of thermally-denuded aerosol was similar between ambient and PILS-generated aerosol and showed limited dependence on volatilization. Results of AMS 3-factor Positive Matrix Factorization (PMF) analysis for the PILS-generated aerosol showed that the most hygroscopic components are most likely the most and the least volatile features of the aerosol. No clear relationship was found between organic hygroscopicity and oxygen-to-carbon ratio; in fact, Isoprene organic aerosol (Isoprene-OA) was found to be the most hygroscopic factor, while at the same time being the least oxidized and likely most volatile of all PMF factors. Considering the diurnal variation of each PMF factor and its associated hygroscopicity, Isoprene-OA and More Oxidized – Oxidized Oxygenated Organic Aerosol (MO-OOA) are the prime contributors to hygroscopicity and covary with Less Oxidized – Oxidized Oxygenated Organic Aerosol (LO-OOA) in a way that induces the observed diurnal invariance in total organic hygroscopicity. Biomass Burning Organic Aerosol (BBOA) contributed little to aerosol hygroscopicity, which is expected since there was little biomass burning activity during the sampling period examined.

Description: Elemental composition and clustering of α-pinene oxidation products for different oxidation conditions Atmospheric Chemistry and Physics Discussions, 14, 30799-30833, 2014 Author(s): A. P. Praplan, S. Schobesberger, F. Bianchi, M. P. Rissanen, M. Ehn, T. Jokinen, H. Junninen, A. Adamov, A. Amorim, J. Dommen, J. Duplissy, J. Hakala, A. Hansel, M. Heinritzi, J. Kangasluoma, J. Kirkby, M. Krapf, A. Kürten, K. Lehtipalo, F. Riccobono, L. Rondo, N. Sarnela, M. Simon, A. Tomé, J. Tröstl, P. M. Winkler, C. Williamson, P. Ye, J. Curtius, U. Baltensperger, N. M. Donahue, M. Kulmala, and D. R. Worsnop This study presents the difference between oxidised organic compounds formed by α-pinene ozonolysis and hydroxyl radical (OH) oxidation in the CLOUD environmental chamber. The results from three Atmospheric Pressure interface Time-Of-Flight (APi-TOF) mass spectrometers measuring simultaneously the composition of naturally charged, as well as neutral species (via chemical ionisation with nitrate) are discussed. Natural chemical ionisation takes place in the CLOUD chamber and organic oxidised compounds form clusters with nitrate, bisulphate, bisulphate/sulphuric acid clusters, ammonium, and dimethylaminium, or get protonated. This process is selective towards various oxidised organic compounds, so that in order to get a comprehensive picture of the elemental composition of oxidation products, several instruments must be used. A comparison between oxidation products containing 10 and 20 carbon atoms is presented. Oxidation products from ozonolysis showed a higher oxidation state than the ones from OH oxidation. Also, highly oxidised organic compounds are shown to be formed in the early stages of the oxidation, for low α-pinene levels.

Description: Solar cycle in current reanalyses: (non)linear attribution study Atmospheric Chemistry and Physics Discussions, 14, 30879-30912, 2014 Author(s): A. Kuchar, P. Sacha, J. Miksovsky, and P. Pisoft This study focusses on the variability of temperature, ozone and circulation characteristics in the stratosphere and lower mesosphere with regard to the influence of the 11 year solar cycle. It is based on attribution analysis using multiple nonlinear techniques (Support Vector Regression, Neural Networks) besides the traditional linear approach. The analysis was applied to several current reanalysis datasets for the 1979–2013 period, including MERRA, ERA-Interim and JRA-55, with the aim to compare how this type of data resolves especially the double-peaked solar response in temperature and ozone variables and the consequent changes induced by these anomalies. Equatorial temperature signals in the lower and upper stratosphere were found to be sufficiently robust and in qualitative agreement with previous observational studies. The analysis also pointed to the solar signal in the ozone datasets (i.e. MERRA and ERA-Interim) not being consistent with the observed double-peaked ozone anomaly extracted from satellite measurements. Consequently the results obtained by linear regression were confirmed by the nonlinear approach through all datasets, suggesting that linear regression is a relevant tool to sufficiently resolve the solar signal in the middle atmosphere. Furthermore, the seasonal dependence of the solar response was also discussed, mainly as a source of dynamical causalities in the wave propagation characteristics in the zonal wind and the induced meridional circulation in the winter hemispheres. The hypothetical mechanism of a weaker Brewer Dobson circulation was reviewed together with discussion of polar vortex stability.

Description: The oxidation capacity of the boreal forest: first simulated reactivities of O 3 and NO 3 Atmospheric Chemistry and Physics Discussions, 14, 30947-31007, 2014 Author(s): D. Mogensen, R. Gierens, J. N. Crowley, P. Keronen, S. Smolander, A. Sogachev, A. C. Nölscher, L. Zhou, M. Kulmala, M. J. Tang, J. Williams, and M. Boy Using the 1D atmospheric chemistry–transport model SOSAA, we have investigated the atmospheric reactivity of a boreal forest ecosystem during the HUMPPA-COPEC-10 campaign (summer 2010, at SMEAR II in Southern Finland). For the very first time, we present vertically resolved model simulations of the NO 3 - and O 3 -reactivity ( R ) together with the modelled and measured reactivity of OH. We find that OH is the most reactive oxidant ( R ~3 s −1 ) followed by NO 3 ( R ~0.07 s −1 ) and O 3 ( R ~2 × 10 −5 s −1 ). The missing OH-reactivity was found to be large in accordance with measurements (~65%) as would be expected from the chemical subset described in the model. The accounted OH radical sinks were inorganic compounds (~41%, mainly due to reaction with CO), emitted monoterpenes (~14%) and oxidised biogenic volatile organic compounds (~44%). The missing reactivity is expected to be due to unknown biogenic volatile organic compounds and their photoproducts, indicating that the true main sink of OH is not expected to be inorganic compounds. The NO 3 radical was found to react mainly with primary emitted monoterpenes (~60%) and inorganic compounds (~37%, including NO 2 ). NO 2 is, however, only a temporary sink of NO 3 under the conditions of the campaign and does not affect the NO 3 concentration. We discuss the difference between instantaneous and steady state reactivity and present the first boreal forest steady state lifetime of NO 3 (113 s). O 3 almost exclusively reacts with inorganic compounds (~91%, mainly NO, but also NO 2 during night) and less with primary emitted sesquiterpenes (~6%) and monoterpenes (~3%). When considering the concentration of the oxidants investigated, we find that O 3 is the oxidant that is capable of removing pollutants fastest. As part of this study, we developed a simple empirical parameterisation for conversion of measured spectral irradiance into actinic flux. Further, the meteorological conditions were evaluated using radiosonde observations and ground based measurements. The overall vertical structure of the boundary layer is discussed, together with validation of the surface energy balance and turbulent fluxes. The sensible heat and momentum fluxes above the canopy were on average overestimated, while the latent heat flux was underestimated.

Description: Sensitivity of the air–sea CO 2 exchange in the Baltic Sea and Danish inner waters to atmospheric short term variability Biogeosciences Discussions, 11, 16993-17042, 2014 Author(s): A. S. Lansø, J. Bendtsen, J. H. Christensen, L. L. Sørensen, H. Chen, H. A. J. Meijer, and C. Geels Minimising the uncertainties in estimates of air–sea CO 2 exchange is an important step toward increasing the confidence in assessments of the CO 2 cycle. Using an atmospheric transport model makes it possible to investigate the direct impact of atmospheric parameters on the air–sea CO 2 flux along with its sensitivity to e.g. short-term temporal variability in wind speed, atmospheric mixing height and the atmospheric CO 2 concentration. With this study the importance of high spatiotemporal resolution of atmospheric parameters for the air–sea CO 2 flux is assessed for six sub-basins within the Baltic Sea and Danish inner waters. A new climatology of surface water partial pressure of CO 2 ( p CO 2 ) has been developed for this coastal area based on available data from monitoring stations and underway p CO 2 measuring systems. Parameterisations depending on wind speed were applied for the transfer velocity to calculate the air–sea CO 2 flux. Two model simulations were conducted – one including short term variability in atmospheric CO 2 (VAT), and one where it was not included (CAT). A seasonal cycle in the air–sea CO 2 flux was found for both simulations for all sub-basins with uptake of CO 2 in summer and release of CO 2 to the atmosphere in winter. During the simulated period 2005–2010 the average annual net uptake of atmospheric CO 2 for the Baltic Sea, Danish Straits and Kattegat was 287 and 471 Gg C yr -1 for the VAT and CAT simulations, respectively. The obtained difference of 184 Gg C yr -1 was found to be significant, and thus ignoring short term variability in atmospheric CO 2 does have a sizeable effect on the air–sea CO 2 exchange. The combination of the atmospheric model and the new p CO 2 fields has also made it possible to make an estimate of the marine part of the Danish CO 2 budget for the first time. A net annual uptake of 2613 Gg C yr -1 was found for the Danish waters. A large uncertainty is connected to the air–sea CO 2 flux in particular caused by the transfer velocity parameterisation and the applied p CO 2 climatology. However, the present study underlines the importance of including short term variability in the atmospheric CO 2 concentration in future model studies of the air–sea exchange in order to minimise the uncertainty.

Description: Epic landslide erosion from mountain roads in Yunnan, China – challenges for sustainable development Natural Hazards and Earth System Science, 14, 3093-3104, 2014 Author(s): R. C. Sidle, M. Ghestem, and A. Stokes Expanding systems of mountain roads in developing countries have significantly increased the risk of landslides and sedimentation, and have created vulnerabilities for residents and aquatic resources. We measured landslide erosion along seven road segments in steep terrain in the upper Salween River basin, Yunnan, China and estimated sediment delivery to channels. Landslide erosion rates along the roads ranged from 2780 to 48 235 Mg ha −1 yr −1 , the upper end of this range being the highest rate ever reported along mountain roads. The two roads with the highest landslide erosion (FG1 = 12 966 Mg ha −1 yr −1 ; DXD = 48 235 Mg ha −1 yr −1 ) had some of the highest sediment delivery rates to channels (about 80 and 86%, respectively). Overall, 3 times more landslides occurred along cut slopes compared to fill slopes, but fill slope failures had a combined mass 〉 1.3 times that of cut slope failures. Many small landslides occurred along road cuts, but these were often trapped on the road surface. Given the magnitude of the landslide problem and the lack of attention to this issue, a more sustainable approach for mountain road development is outlined based on an analysis of landslide susceptibility and how thresholds for landslide trigger mechanisms would be modified by road location and different construction techniques.

Description: Evaluation of forest fire models on a large observation database Natural Hazards and Earth System Science, 14, 3077-3091, 2014 Author(s): J. B. Filippi, V. Mallet, and B. Nader This paper presents the evaluation of several fire propagation models using a large set of observed fires. The observation base is composed of 80 Mediterranean fire cases of different sizes, which come with the limited information available in an operational context (burned surface and approximative ignition point). Simulations for all cases are carried out with four different front velocity models. The results are compared with several error scoring methods applied to each of the 320 simulations. All tasks are performed in a fully automated manner, with simulations run as first guesses with no tuning for any of the models or cases. This approach leads to a wide range of simulation performance, including some of the bad simulation results to be expected in an operational context. Disregarding the quality of the input data, it is found that the models can be ranked based on their performance and that the most complex models outperform the more empirical ones. Data and source codes used for this paper are freely available to the community.

Description: Monitoring and prediction in Early Warning Systems (EWS) for rapid mass movements Natural Hazards and Earth System Sciences Discussions, 2, 7149-7179, 2014 Author(s): M. Stähli, M. Sättele, C. Huggel, B. W. McArdell, P. Lehmann, A. Van Herwijnen, A. Berne, M. Schleiss, A. Ferrari, A. Kos, D. Or, and S. M. Springman Rapid mass movements (RMM) pose a substantial risk to people and infrastructure. Reliable and cost-efficient measures have to be taken to reduce this risk. One of these measures includes establishing and advancing the State of Practice in the application of Early Warning Systems (EWS). EWS have been developed during the past decades and are rapidly increasing. In this document, we focus on the technical part of EWS, i.e. the prediction and timely recognition of imminent hazards, as well as on monitoring slopes at risk and released mass movements. Recent innovations in assessing spatial precipitation, as well as monitoring and modelling precursors, the triggering and deformation of RMM offer new opportunities for next-generation EWS. However, technical advancement can only be transferred into more reliable, operational EWS with an intense dialog between scientists, engineers and those in charge of warning. To this end, further experience with new comprehensive prototype systems jointly operated by scientists and practitioners will be essential.

Description: A new approach to mapping permafrost and change incorporating uncertainties in ground conditions and climate projections The Cryosphere, 8, 2177-2194, 2014 Author(s): Y. Zhang, I. Olthof, R. Fraser, and S. A. Wolfe Spatially detailed information on permafrost distribution and change with climate is important for land use planning, infrastructure development, and environmental assessments. However, the required soil and surficial geology maps in the North are coarse, and projected climate scenarios vary widely. Considering these uncertainties, we propose a new approach to mapping permafrost distribution and change by integrating remote sensing data, field measurements, and a process-based model. Land cover types from satellite imagery are used to capture the general land conditions and to improve the resolution of existing permafrost maps. For each land cover type, field observations are used to estimate the probabilities of different ground conditions. A process-based model is used to quantify the evolution of permafrost for each ground condition under three representative climate scenarios (low, medium, and high warming). From the model results, the probability of permafrost occurrence and the most likely permafrost conditions are determined. We apply this approach at 20 m resolution to a large area in Northwest Territories, Canada. Mapped permafrost conditions are in agreement with field observations and other studies. The data requirements, model robustness, and computation time are reasonable, and this approach may serve as a practical means to mapping permafrost and changes at high resolution in other regions.

Description: Biogeography of jellyfish in the North Atlantic, by traditional and genomic methods Earth System Science Data Discussions, 7, 629-667, 2014 Author(s): P. Licandro, M. Blackett, A. Fischer, A. Hosia, J. Kennedy, R. R. Kirby, K. Raab, R. Stern, and P. Tranter Scientific debate on whether the recent increase in reports of jellyfish outbreaks is related to a true rise in their abundance, have outlined the lack of reliable records of Cnidaria and Ctenophora. Here we describe different data sets produced within the EU program EUROBASIN, which have been assembled with the aim of presenting an up to date overview of the diversity and standing stocks of jellyfish in the North Atlantic region. Using a net adapted to sample gelatinous zooplankton quantitatively, Cnidaria and Ctenophora were collected in the epipelagic layer during spring-summer 2010–2013, in inshore and offshore waters between 59–68° N Lat and 62° W–5° E Long. Jellyfish were also identified and counted in samples opportunistically collected by other sampling equipment in the same region and at two coastal stations in the Bay of Biscay and in the Gulf of Cadiz. Continuous Plankton Recorder (CPR) samples collected in 2009–2012 were re-analysed with the aim of identifying the time and location of Cnidarian blooms across the North Atlantic basin. Overall the data show high variability in jellyfish abundance and diversity, mainly in relation with different water masses and with the bathymetry. Higher densities were generally recorded on the shelves, where populations tend to be more diversified due to the presence of meropelagic medusae. Comparisons of net records from the G.O. Sars transatlantic cruise show that information on jellyfish diversity differs significantly depending on the sampling gear utilised. Indeed, the big trawls mostly collect relatively large scyphozoan and hydrozoan species, while small hydrozoans and early stages of ctenophora are only caught by smaller nets. Based on CPR data from 2009–2012, blooms of Cnidarians occurred in all seasons across the whole North Atlantic basin. Molecular analysis revealed that, in contrast with what was previously hypothesized, the CPR is able to detect blooms of meroplanktonic and holoplanktonic hydrozoans and scyphozoans. Combining different types of data, key jellyfish taxa for the spring-summer period were identified in the northern North Atlantic regions. Key species for the central and southern North Atlantic could be inferred based on Cnidarian blooms identified by the CPR survey, although this should be confirmed further by comparison with quantitative data. The identification by DNA barcoding of 23 jellyfish specimens collected during the EUROBASIN cruises contributes to increasing the still very limited number of jellyfish sequences available on GenBank. All observations presented here can be downloaded from PANGAEA ( http://doi.pangaea.de/10.1594/PANGAEA.835732 ).

Description: A comparison between VEGA 1, 2 and Giotto flybys of comet 1P/Halley: implications for Rosetta Annales Geophysicae, 32, 1441-1453, 2014 Author(s): M. Volwerk, K.-H. Glassmeier, M. Delva, D. Schmid, C. Koenders, I. Richter, and K. Szegö Three flybys of comet 1P/Halley, by VEGA 1, 2 and Giotto, are investigated with respect to the occurrence of mirror mode waves in the cometosheath and field line draping in the magnetic pile-up region around the nucleus. The time interval covered by these flybys is approximately 8 days, which is also the approximate length of an orbit or flyby of Rosetta around comet 67P/Churyumov–Gerasimenko. Thus any significant changes observed around Halley are changes that might occur for Rosetta during one pass of 67P/CG. It is found that the occurrence of mirror mode waves in the cometosheath is strongly influenced by the dynamical pressure of the solar wind and the outgassing rate of the comet. Field line draping happens in the magnetic pile-up region. Changes in nested draping regions (i.e. regions with different B x directions) can occur within a few days, possibly influenced by changes in the outgassing rate of the comet and thereby the conductivity of the cometary ionosphere.

Description: A Fabry–Perot interferometer-based camera for two-dimensional mapping of SO 2 distributions Atmospheric Measurement Techniques, 7, 3705-3715, 2014 Author(s): J. Kuhn, N. Bobrowski, P. Lübcke, L. Vogel, and U. Platt We examine a new imaging method for the remote sensing of volcanic gases, which relies on the regularly spaced narrow-band absorption structures in the UV–VIS of many molecules. A Fabry–Perot interferometer (FPI) is used to compare the scattered sunlight radiance at wavelengths corresponding to absorption bands with the radiance at wavelengths in between the bands, thereby identifying and quantifying the gas. In this first theoretical study, we present sample calculations for the detection of sulfur dioxide (SO 2 ). Optimum values for the FPI setup parameters are proposed. Furthermore, the performance of the FPI method is compared to SO 2 cameras. We show that camera systems using an FPI are far less influenced by changes in atmospheric radiative transfer (e.g., due to aerosol) and have a great potential as a future technique for examining emissions of SO 2 (or other gases) from volcanic sources and other point sources.

Description: Comparison of nitric oxide measurements in the mesosphere and lower thermosphere from ACE-FTS, MIPAS, SCIAMACHY, and SMR Atmospheric Measurement Techniques Discussions, 7, 12735-12794, 2014 Author(s): S. Bender, M. Sinnhuber, T. von Clarmann, G. Stiller, B. Funke, M. López-Puertas, J. Urban, K. Pérot, K. A. Walker, and J. P. Burrows We compare the nitric oxide measurements in the mesosphere and lower thermosphere (60 to 150 km) from four instruments: ACE-FTS, MIPAS, SCIAMACHY, and SMR. We use the daily zonal mean data in that altitude range for the years 2004–2010 (ACE-FTS), 2005–2012 (MIPAS), 2008–2012 (SCIAMACHY), and 2003–2012 (SMR). We first compare the data qualitatively with respect to the morphology, focussing on the major features, and then compare the time series directly and quantitatively. In three geographical regions, we compare the vertical density profiles on coincident measurement days. Since none of the instruments delivers continuous daily measurements in this altitude region, we carried out a multi-linear regression analysis. This regression analysis considers annual and semi-annual variability in form of harmonic terms and inter-annual variability by responding linearly to the solar Lyman-α radiation index and the geomagnetic Kp index. This analysis helps to find similarities and differences in the individual data sets with respect to the inter-annual variations caused by geomagnetic and solar variability. We find that the data sets are consistent and that they only disagree on minor aspects. SMR and ACE-FTS deliver the longest time series in the mesosphere and they both agree remarkably well. The shorter time series from MIPAS and SCIAMACHY also agree with them where they overlap. The data agree within ten to twenty percent when the number densities are large, but they can differ by 50 to 100% in some cases.

Description: Middle-atmospheric zonal and meridional wind profiles from polar, tropical and midlatitudes with the ground-based microwave Doppler wind radiometer WIRA Atmospheric Measurement Techniques, 7, 4491-4505, 2014 Author(s): R. Rüfenacht, A. Murk, N. Kämpfer, P. Eriksson, and S. A. Buehler WIRA is a ground-based microwave Doppler spectroradiometer specifically designed for the measurement of profiles of horizontal wind in the upper stratosphere and lower mesosphere region where no other continuously running measurement technique exists. A proof of principle has been delivered in a previous publication. A technical upgrade including a new high-frequency amplifier and sideband filter has improved the signal to noise ratio by a factor of 2.4. Since this upgrade the full horizontal wind field comprising zonal and meridional wind profiles is continuously measured. A completely new retrieval based on optimal estimation has been set up. Its characteristics are detailed in the present paper. Since the start of the routine operation of the first prototype in September 2010, WIRA has been measuring at four different locations at polar, mid- and tropical latitudes (67°22' N/26°38' E, 46°57' N/7°26' E, 43°56' N/5°43' E and 21°04' S/55°23' E) for time periods between 5.5 and 11 months. The data presented in this paper are daily average wind profiles with typical uncertainties and resolutions of 10 to 20 m s −1 and 10 to 16 km, respectively. A comparison between the data series from WIRA and European Centre for Medium-Range Weather Forecasts (ECMWF) model data revealed agreement within 10% in the stratospheric zonal wind. The meridional wind profiles agree within their error bars over the entire sensitive altitude range of WIRA. However, significant differences in the mesospheric zonal wind speed of up to 50% have been found.

Description: Do climate factors govern soil microbial community composition and biomass at a regional scale? Biogeosciences Discussions, 11, 17729-17756, 2014 Author(s): L. Ma, C. Guo, X. Lü, S. Yuan, and R. Wang Soil microbial communities play important role in organic matter decomposition, nutrient cycling and vegetation dynamic. However, little is known about factors driving soil microbial community composition at large scales. The objective of this study was to determine whether climate dominates among environmental factors governing microbial community composition and biomass at a regional scale. Here, we compared soil microbial communities using phospholipid fatty acid method across 7 land use types from 23 locations in North-East China Transect (850 km x 50 km). The results showed that soil water availability and land use changes exhibited the dominant effects on soil microbial community composition and biomass at the regional scale, while climate factors (expressed as a function of large-scale spatial variation) did not show strong relationships with distribution of microbial community composition. Likewise, factors such as spatial structure, soil texture, nutrient availability and vegetation types were not important. Wetter soils had higher contributions of gram-positive bacteria, whereas drier soils had higher contributions of gram-negative bacteria and fungi. Heavily disturbed soils had lower contributions of gram-negative bacteria and fungi than historically disturbed and undisturbed soils. The lowest microbial biomass appeared in the wettest and driest soils. In conclusion, dominant climate factors, commonly known to structure distribution of macroorganisms, were not the most important drivers governing regional pattern of microbial communities because of inclusion of irrigated and managed practices. In comparison, soil water regime and land use types appear to be primary determinants of microbial community composition and biomass.

Description: Fossilized bioelectric wire – the trace fossil Trichichnus Biogeosciences Discussions, 11, 17707-17728, 2014 Author(s): M. Kędzierski, A. Uchman, Z. Sawlowicz, and A. Briguglio The trace fossil Trichichnus is proposed as an indicator of fossil bioelectric bacterial activity at the interface oxic – anoxic zone of marine sediments. This fulfils the idea that such processes, commonly found in the modern realm, should be also present in the geological past. Trichichnus is an exceptional trace fossil due to its very thin diameter (mostly less than 1 mm) and common pyritic filling. It is ubiquitous in some fine-grained sediments, where it has been interpreted as a burrow formed deeper than any other trace fossils, below the redox boundary. Trichichnus formerly referred to as deeply burrowed invertebrates, has been found as remnant of a fossilized intrasediment bacterial mat that is pyritized. As visualized in 3-D by means of X-ray computed microtomography scanner, Trichichnus forms dense filamentous fabric, which reflects that produced by modern large, mat-forming, sulphide-oxidizing bacteria, belonging mostly to Trichichnus -related taxa, which are able to house a complex bacterial consortium. Several stages of Trichichnus formation, including filamentous, bacterial mat and its pyritization, are proposed to explain an electron exchange between oxic and suboxic/anoxic layers in the sediment. Therefore, Trichichnus can be considered a fossilized "electric wire".

Description: Nitrate limitation and ocean acidification interact with UV-B to reduce photosynthetic performance in the diatom Phaeodactylum tricornutum Biogeosciences Discussions, 11, 17675-17706, 2014 Author(s): W. Li, K. Gao, and J. Beardall It has been proposed that ocean acidification (OA) will interact with other environmental factors to influence the overall impact of global change on biological systems. Accordingly we investigated the influence of nitrogen limitation and OA on the physiology of diatoms by growing the diatom Phaeodactylum tricornutum Bohlin under elevated (1000 μatm, HC) or ambient (390 μatm, LC) levels of CO 2 with replete (110 μmol L -1 , HN) or reduced (10 μmol L -1 , LN) levels of NO 3 - and subjecting the cells to solar radiation with or without UV irradiance to determine their susceptibility to UV radiation (280–400 nm). Our results indicate that OA and UVB induced significantly higher inhibition of both the photosynthetic rate and quantum yield under LN than under HN conditions. UVA or/and UVB increased the cells' non-photochemical quenching (NPQ) regardless of the CO 2 levels. Under LN and OA conditions, activity of superoxide dismutase and catalase activities were enhanced, along with the highest sensitivity to UVB and the lowest ratio of repair to damage of PSII. HC-grown cells showed a faster recovery rate of yield under HN but not under LN conditions. The finding that nitrate limitation and ocean acidification interact with UV-B to reduce photosynthetic performance of the diatom P. tricornutum implies that ocean primary production and the marine biological C pump will be affected by the OA under multiple stressors.

Description: Vegetation heterogeneity and landscape position exert strong controls on soil CO 2 efflux in a moist, Appalachian watershed Biogeosciences Discussions, 11, 17631-17673, 2014 Author(s): J. W. Atkins, H. E. Epstein, and D. L. Welsch In topographically complex watersheds, landscape position and vegetation heterogeneity can alter the soil water regime through both lateral and vertical redistribution, respectively. These alterations of soil moisture may have significant impacts on the spatial heterogeneity of biogeochemical cycles throughout the watershed. To evaluate how landscape position and vegetation heterogeneity affect soil CO 2 efflux ( F SOIL ) we conducted observations across the Weimer Run watershed (373 ha), located near Davis, West Virginia, for three growing seasons with varying precipitation (2010 – 1042 mm; 2011 – 1739 mm; 2012 – 1244 mm; precipitation data from BDKW2 station, MesoWest, University of Utah). An apparent soil temperature threshold of 11 °C at 12 cm depth on F SOIL was observed in our data – where F SOIL rates greatly increase in variance above this threshold. For analysis, F SOIL values above this threshold were isolated and examined. Differences in F SOIL among years were apparent by elevation ( F 4,633 = 3.17; p = 0.013) and by vegetation cover ( F 4, 633 = 2.96; p = 0.019). For the Weimer Run watershed, vegetation exerts the major control on soil CO 2 efflux ( F SOIL ), with the plots beneath shrubs at all elevations for all years showing the greatest mean rates of F SOIL (6.07 μmol CO 2 m -2 s -1 ) compared to plots beneath closed-forest canopy (4.69 μmol CO 2 m -2 s -1 ) and plots located in open, forest gaps (4.09 μmol CO 2 m -2 s -1 ) plots. During periods of high soil moisture, we find that CO 2 efflux rates are constrained and that maximum efflux rates in this system occur during periods of average to below average soil water availability. These findings offer valuable insight into the processes occurring within these topographically complex, temperate and humid systems, and the interactions of abiotic and biotic factors mediating biogeochemical cycles. With possible changing rainfall patterns as predicted by climate models, it is important to understand the couplings between water and carbon cycling at the watershed and landscape scales, and their potential dynamics under global change scenarios.

Description: Daily burned area and carbon emissions from boreal fires in Alaska Biogeosciences Discussions, 11, 17579-17629, 2014 Author(s): S. Veraverbeke, B. M. Rogers, and J. T. Randerson Boreal fires burn carbon-rich organic soils, thereby releasing large quantities of trace gases and aerosols that influence atmospheric composition and climate. To better understand the factors regulating boreal fire emissions, we developed a statistical model of carbon consumption by fire for Alaska with a spatial resolution of 500 m and a temporal resolution of one day. We used the model to estimate variability in carbon emissions between 2001 and 2012. Daily burned area was mapped using imagery from the Moderate Resolution Imaging Spectroradiometer combined with perimeters from the Alaska Large Fire Database. Carbon consumption was calibrated using available field measurements from black spruce forests in Alaska. We built two nonlinear multiplicative models to separately predict above- and belowground carbon consumption by fire in response to environmental variables including elevation, day of burning within the fire season, pre-fire tree cover and the differenced normalized burn ratio (dNBR). Higher belowground consumption occurred later in the season and for mid-elevation regions. Aboveground and belowground consumption also increased as a function of tree cover and the dNBR, suggesting a causal link between the processes regulating these two components of consumption. Between 2001 and 2012, the median fuel consumption was 2.48 kg C m -2 and the median pixel-based uncertainty (SD of prediction error) was 0.38 kg C m -2 . There were considerable amounts of burning in other cover types than black spruce and consumption in pure black spruce stands was generally higher. Fuel consumption originated primarily from the belowground fraction (median = 2.30 kg C m -2 for all cover types and 2.63 kg C m -2 for pure black spruce stands). Total carbon emissions varied considerably from year to year, with the highest emissions occurring during 2004 (67 Tg C), 2005 (44 Tg C), 2009 (25 Tg C), and 2002 (16 Tg C) and a mean of 14 Tg C per year between 2001 and 2012. Our analysis highlights the importance of accounting for the spatial heterogeneity within fuels and consumption when extrapolating emissions in space and time. This data on daily burned area and emissions may be useful for in understanding controls and limits on fire growth, and predicting potential feedbacks of changing fire regimes.

Description: GPS radio occultation with TerraSAR-X and TanDEM-X: sensitivity of lower troposphere sounding to the Open-Loop Doppler model Atmospheric Measurement Techniques Discussions, 7, 12719-12733, 2014 Author(s): F. Zus, G. Beyerle, S. Heise, T. Schmidt, and J. Wickert The Global Positioning System (GPS) radio occultation (RO) technique provides valuable input for numerical weather prediction and is considered as a data source for climate related research. Numerous studies outline the high precision and accuracy of RO atmospheric soundings in the upper troposphere and lower stratosphere. In this altitude region (8–25 km) RO atmospheric soundings are considered to be free of any systematic error. In the tropical (30° S–30° N) Lower (

Description: The sharp decline of East Asian summer monsoon at mid-Holocene indicated by the lake-wetland transition in the Sanjiang Plain, northeastern China Climate of the Past Discussions, 10, 4595-4622, 2014 Author(s): Z. Q. Zhang, G. P. Wang, X. G. Lv, H. J. Jia, and Q. H. Xu The timing of the waxing and wining of the East Asian summer monsoon during the Holocene is still under debate. In present study, we present the high-resolution grain-size and LOI records from a well-dated mud/peat profile to reveal the lake-wetland transition in the Sanjiang Plain and discuss its significance to Holocene monsoon evolutions. The results show that the shallow-water lakes have developed in low-lying areas of the plain before 4600 yr BP, corresponding to the Holocene monsoon maximum. Thereafter, the wetlands began to initiate with the extinction of the paleolakes, marking a lake-shrinking stage with the relative dry climate. Considering the prevalent monsoon climate in the Sanjiang Plain, we suggest the lake-wetland transition at 4600 yr BP indicate a sharp decline of the summer monsoon rather than the basin infilling process. Such a remarkable monsoon weakening event has been widely documented in northern China, and we associated it with the ocean–atmosphere interacting processes in low-latitude regions.

Description: Large-eddy simulations of contrails in a turbulent atmosphere Atmospheric Chemistry and Physics Discussions, 14, 29499-29546, 2014 Author(s): J. Picot, R. Paoli, O. Thouron, and D. Cariolle In this work, the evolution of contrails in the vortex and dissipation regimes is studied by means of fully three-dimensional large-eddy simulation (LES) coupled to a Lagrangian particle tracking method to treat the ice phase. This is the first paper where fine-scale atmospheric turbulence is generated and sustained by means of a stochastic forcing that mimics the properties of stably stratified turbulent flows as those occurring in the upper troposphere lower stratosphere. The initial flow-field is composed by the turbulent background flow and a wake flow obtained from separate LES of the jet regime. Atmospheric turbulence is the main driver of the wake instability and the structure of the resulting wake is sensitive to the intensity of the perturbations, primarily in the vertical direction. A stronger turbulence accelerates the onset of the instability, which results in shorter contrail decent and more effective mixing in the interior of the plume. However, the self-induced turbulence that is produced in the wake after the vortex break-up dominates over background turbulence at the end of the vortex regime and dominates the mixing with ambient air. This results in global microphysical characteristics such as ice mass and optical depth that are be slightly affected by the intensity of atmospheric turbulence. On the other hand, the background humidity and temperature have a first order effect on the survival of ice crystals and particle size distribution, which is in line with recent and ongoing studies in the literature.

Description: Influence of frontal cyclone evolution on the 2009 (Ekman) and 2010 (Franklin) Loop Current eddy detachment events Ocean Science, 10, 947-965, 2014 Author(s): Y. S. Androulidakis, V. H. Kourafalou, and M. Le Hénaff The anticyclonic Loop Current Eddy (LCE) shedding events are strongly associated with the evolution of Loop Current Frontal Eddies (LCFEs) over the eastern Gulf of Mexico (GoM). A numerical simulation, in tandem with in situ measurements and satellite data, was used to investigate the Loop Current (LC) evolution and the surrounding LCFE formation, structure, growth and migration during the Eddy Ekman and Eddy Franklin shedding events in the summers of 2009 and 2010, respectively. During both events, northern GoM LCFEs appeared vertically coherent to at least 1500 m in temperature observations. They propagated towards the base of the LC, where, together with the migration of Campeche Bank (southwest GoM shelf) eddies from south of the LC, contributed to its "necking-down". Growth of Campeche Bank LCFEs involved in Eddy Franklin was partially attributed to Campeche Bank waters following upwelling events. Slope processes associated with such upwelling included offshore exports of high positive potential vorticity that may trigger cyclone formation and growth. The advection and growth of LCFEs, originating from the northern and southern GoM, and their interaction with the LC over the LCE detachment area favor shedding conditions and may contribute to the final separation of the LCE.

Description: A two thousand year annual record of snow accumulation rates for Law Dome, East Antarctica Climate of the Past Discussions, 10, 4469-4497, 2014 Author(s): J. Roberts, C. Plummer, T. Vance, T. van Ommen, A. Moy, S. Poynter, A. Treverrow, M. Curran, and S. George Accurate high resolution records of snow accumulation rates in Antarctica are crucial for estimating ice sheet mass balance and subsequent sea level change. Snowfall rates at Law Dome, East Antarctica, have been linked with regional atmospheric circulation to mid-latitudes as well as regional Antarctic snowfall. Here, we extend the Law Dome accumulation record from 750 to 2035 years, using recent annual layer dating that extends to AD −22. Accumulation rates were calculated as the ratio of measured to modelled layer thicknesses, multiplied by the long term mean accumulation rate. The modelled layer thicknesses were based on a power law vertical strain rate profile fitted to observed annual layer thickness. The periods AD 380–442, AD 727–783 and AD 1970–2009 have above average snow accumulation rates, while AD 663–704, AD 933–975 and AD 1429–1468 were below average. The calculated snow accumulation rates show good correlation with atmospheric reanalysis estimates, and significant spatial correlation over a wide expanse of East Antarctica, demonstrating that the Law Dome record captures larger scale variability across a large region of East Antarctica well beyond the immediate vicinity of the Law Dome summit. Spectral analysis reveals periodicities in the snow accumulation record which may be related to ENSO and Interdecadal Pacific Oscillation frequencies.

Description: Sensitivity of East African savannah vegetation to historical moisture-balance variation Climate of the Past, 10, 2067-2080, 2014 Author(s): I. Ssemmanda, V. Gelorini, and D. Verschuren Fossil pollen records provide key insight into the sensitivity of terrestrial ecosystems to climate change. However, tracing vegetation response to relatively modest historical climate fluctuations is often complicated by the overriding signature of anthropogenic landscape disturbance. Here we use high-resolution pollen data from a ~200-year lake-sediment record in open wooded savannah of Queen Elizabeth National Park (southwestern Uganda) to assess the sensitivity of the tropical lowland grassland–forest transition to historical, decade-scale moisture-balance fluctuations. Specifically we trace vegetation response to three episodes of higher average rainfall dated to the 1820s–1830s, ca. 1865–1890 and from 1962 to around 2000. Our pollen data indeed reveal a sequence of three wet periods, separated by two drier periods. During the inferred wetter episodes we find increases in the percent pollen abundance of trees and shrubs from moist semi-deciduous forest ( Allophylus , Macaranga , Alchornea , Celtis ), riparian forest ( Phoenix reclinata ) and wooded savannah ( Acalypha , Rhus -type vulgaris , Combretaceae/Melastomataceae) as well as taxa common in the local rift-valley grasslands ( Acacia , Ficus ), together creating strong temporary reductions in Poaceae pollen (to 45–55% of the terrestrial pollen sum). During intervening dry periods, Poaceae pollen attained values of 65–75%, and dryland herbs such as Commelina , Justicia -type odora and Chenopodiaceae expanded at the expense of Asteraceae, Solanum -type, Swertia usambarensis -type, and (modestly so) Urticaceae. Noting that the overall richness of arboreal taxa remained high but their combined abundance low, we conclude that the landscape surrounding Lake Chibwera has been an open wooded savannah throughout the past 200 years, with historical moisture-balance variation exerting modest effects on local tree cover (mostly the abundance of Acacia and Ficus ) and the occurrence of damp soil areas promoting Phoenix reclinata . The strong apparent expansion of true forest trees during wet episodes can be explained partly by enhanced pollen input via a temporarily activated upland stream. Pollen from exotic trees and cultural indicators appears from the 1970s onwards, but their combined influence fails to mask the signature of natural vegetation dynamics in the pollen record.

Description: CREST (Climate REconstruction SofTware): a probability density function (PDF)-based quantitative climate reconstruction method Climate of the Past, 10, 2081-2098, 2014 Author(s): M. Chevalier, R. Cheddadi, and B. M. Chase Several methods currently exist to quantitatively reconstruct palaeoclimatic variables from fossil botanical data. Of these, probability density function (PDF)-based methods have proven valuable as they can be applied to a wide range of plant assemblages. Most commonly applied to fossil pollen data, their performance, however, can be limited by the taxonomic resolution of the pollen data, as many species may belong to a given pollen type. Consequently, the climate information associated with different species cannot always be precisely identified, resulting in less-accurate reconstructions. This can become particularly problematic in regions of high biodiversity. In this paper, we propose a novel PDF-based method that takes into account the different climatic requirements of each species constituting the broader pollen type. PDFs are fitted in two successive steps, with parametric PDFs fitted first for each species and then a combination of those individual species PDFs into a broader single PDF to represent the pollen type as a unit. A climate value for the pollen assemblage is estimated from the likelihood function obtained after the multiplication of the pollen-type PDFs, with each being weighted according to its pollen percentage. To test its performance, we have applied the method to southern Africa as a regional case study and reconstructed a suite of climatic variables (e.g. winter and summer temperature and precipitation, mean annual aridity, rainfall seasonality). The reconstructions are shown to be accurate for both temperature and precipitation. Predictable exceptions were areas that experience conditions at the extremes of the regional climatic spectra. Importantly, the accuracy of the reconstructed values is independent of the vegetation type where the method is applied or the number of species used. The method used in this study is publicly available in a software package entitled CREST (Climate REconstruction SofTware) and will provide the opportunity to reconstruct quantitative estimates of climatic variables even in areas with high geographical and botanical diversity.

Description: Results from a full coupling of the HIRHAM regional climate model and the MIKE SHE hydrological model for a Danish catchment Hydrology and Earth System Sciences, 18, 4733-4749, 2014 Author(s): M. A. D. Larsen, J. C. Refsgaard, M. Drews, M. B. Butts, K. H. Jensen, J. H. Christensen, and O. B. Christensen A major challenge in the emerging research field of coupling of existing regional climate models (RCMs) and hydrology/land-surface models is the computational interaction between the models. Here we present results from a full two-way coupling of the HIRHAM RCM over a 4000 km × 2800 km domain at 11 km resolution and the combined MIKE SHE-SWET hydrology and land-surface models over the 2500 km 2 Skjern River catchment. A total of 26 one-year runs were performed to assess the influence of the data transfer interval (DTI) between the two models and the internal HIRHAM model variability of 10 variables. DTI frequencies between 12 and 120 min were assessed, where the computational overhead was found to increase substantially with increasing exchange frequency. In terms of hourly and daily performance statistics the coupled model simulations performed less accurately than the uncoupled simulations, whereas for longer-term cumulative precipitation the opposite was found, especially for more frequent DTI rates. Four of six output variables from HIRHAM, precipitation, relative humidity, wind speed and air temperature, showed statistically significant improvements in root-mean-square error (RMSE) by reducing the DTI. For these four variables, the HIRHAM RMSE variability corresponded to approximately half of the influence from the DTI frequency and the variability resulted in a large spread in simulated precipitation. Conversely, DTI was found to have only a limited impact on the energy fluxes and discharge simulated by MIKE SHE.

Description: Seasonal cycle and long-term trend of solar energy fluxes through Arctic sea ice The Cryosphere, 8, 2219-2233, 2014 Author(s): S. Arndt and M. Nicolaus Arctic sea ice has not only decreased in volume during the last decades, but has also changed in its physical properties towards a thinner and more seasonal ice cover. These changes strongly impact the energy budget, and might affect the ice-associated ecosystems. In this study, we quantify solar shortwave fluxes through sea ice for the entire Arctic during all seasons. To focus on sea-ice-related processes, we exclude fluxes through open water, scaling linearly with sea ice concentration. We present a new parameterization of light transmittance through sea ice for all seasons as a function of variable sea ice properties. The maximum monthly mean solar heat flux under the ice of 30 × 10 5 Jm −2 occurs in June, enough heat to melt 0.3 m of sea ice. Furthermore, our results suggest that 96% of the annual solar heat input through sea ice occurs during only a 4-month period from May to August. Applying the new parameterization to remote sensing and reanalysis data from 1979 to 2011, we find an increase in transmitted light of 1.5% yr −1 for all regions. This corresponds to an increase in potential sea ice bottom melt of 63% over the 33-year study period. Sensitivity studies reveal that the results depend strongly on the timing of melt onset and the correct classification of ice types. Assuming 2 weeks earlier melt onset, the annual transmitted solar radiation to the upper ocean increases by 20%. Continuing the observed transition from a mixed multi-year/first-year sea ice cover to a seasonal ice cover results in an increase in light transmittance by an additional 18%.

Description: Instability and change detection in exponential families and generalized linear models, with a study of Atlantic tropical storms Nonlinear Processes in Geophysics, 21, 1133-1143, 2014 Author(s): Y. Lu and S. Chatterjee Exponential family statistical distributions, including the well-known normal, binomial, Poisson, and exponential distributions, are overwhelmingly used in data analysis. In the presence of covariates, an exponential family distributional assumption for the response random variables results in a generalized linear model. However, it is rarely ensured that the parameters of the assumed distributions are stable through the entire duration of the data collection process. A failure of stability leads to nonsmoothness and nonlinearity in the physical processes that result in the data. In this paper, we propose testing for stability of parameters of exponential family distributions and generalized linear models. A rejection of the hypothesis of stable parameters leads to change detection. We derive the related likelihood ratio test statistic. We compare the performance of this test statistic to the popular normal distributional assumption dependent cumulative sum (Gaussian CUSUM) statistic in change detection problems. We study Atlantic tropical storms using the techniques developed here, so to understand whether the nature of these tropical storms has remained stable over the last few decades.

Description: The Rossby radius in the Arctic Ocean Ocean Science, 10, 967-975, 2014 Author(s): A. J. G. Nurser and S. Bacon The first (and second) baroclinic deformation (or Rossby) radii are presented north of ~60° N, focusing on deep basins and shelf seas in the high Arctic Ocean, the Nordic seas, Baffin Bay, Hudson Bay and the Canadian Arctic Archipelago, derived from climatological ocean data. In the high Arctic Ocean, the first Rossby radius increases from ~5 km in the Nansen Basin to ~15 km in the central Canadian Basin. In the shelf seas and elsewhere, values are low (1–7 km), reflecting weak density stratification, shallow water, or both. Seasonality strongly impacts the Rossby radius only in shallow seas, where winter homogenization of the water column can reduce it to below 1 km. Greater detail is seen in the output from an ice–ocean general circulation model, of higher resolution than the climatology. To assess the impact of secular variability, 10 years (2003–2012) of hydrographic stations along 150° W in the Beaufort Gyre are also analysed. The first-mode Rossby radius increases over this period by ~20%. Finally, we review the observed scales of Arctic Ocean eddies.

Description: Interpretative modelling of a geological cross section from boreholes: sources of uncertainty and their quantification Solid Earth, 5, 1189-1203, 2014 Author(s): R. M. Lark, S. Thorpe, H. Kessler, and S. J. Mathers We conducted a designed experiment to quantify sources of uncertainty in geologists' interpretations of a geological cross section. A group of 28 geologists participated in the experiment. Each interpreted borehole record included up to three Palaeogene bedrock units, including the target unit for the experiment: the London Clay. The set of boreholes was divided into batches from which validation boreholes had been withheld; as a result, we obtained 129 point comparisons between the interpreted elevation of the base of the London Clay and its observed elevation in a borehole not used for that particular interpretation. Analysis of the results showed good general agreement between the observed and interpreted elevations, with no evidence of systematic bias. Between-site variation of the interpretation error was spatially correlated, and the variance appeared to be stationary. The between-geologist component of variance was smaller overall, and depended on the distance to the nearest borehole. There was also evidence that the between-geologist variance depends on the degree of experience of the individual. We used the statistical model of interpretation error to compute confidence intervals for any one interpretation of the base of the London Clay on the cross section, and to provide uncertainty measures for decision support in a hypothetical route-planning process. The statistical model could also be used to quantify error propagation in a full 3-D geological model produced from interpreted cross sections.

Description: Using cloud ice flux to parametrise large-scale lightning Atmospheric Chemistry and Physics, 14, 12665-12682, 2014 Author(s): D. L. Finney, R. M. Doherty, O. Wild, H. Huntrieser, H. C. Pumphrey, and A. M. Blyth Lightning is an important natural source of nitrogen oxide especially in the middle and upper troposphere. Hence, it is essential to represent lightning in chemistry transport and coupled chemistry–climate models. Using ERA-Interim meteorological reanalysis data we compare the lightning flash density distributions produced using several existing lightning parametrisations, as well as a new parametrisation developed on the basis of upward cloud ice flux at 440 hPa. The use of ice flux forms a link to the non-inductive charging mechanism of thunderstorms. Spatial and temporal distributions of lightning flash density are compared to tropical and subtropical observations for 2007–2011 from the Lightning Imaging Sensor (LIS) on the Tropical Rainfall Measuring Mission (TRMM) satellite. The well-used lightning flash parametrisation based on cloud-top height has large biases but the derived annual total flash density has a better spatial correlation with the LIS observations than other existing parametrisations. A comparison of flash density simulated by the different schemes shows that the cloud-top height parametrisation has many more instances of moderate flash densities and fewer low and high extremes compared to the other parametrisations. Other studies in the literature have shown that this feature of the cloud-top height parametrisation is in contrast to lightning observations over certain regions. Our new ice flux parametrisation shows a clear improvement over all the existing parametrisations with lower root mean square errors (RMSEs) and better spatial correlations with the observations for distributions of annual total, and seasonal and interannual variations. The greatest improvement with the new parametrisation is a more realistic representation of the zonal distribution with a better balance between tropical and subtropical lightning flash estimates. The new parametrisation is appropriate for testing in chemistry transport and chemistry–climate models that use a lightning parametrisation.

Description: Tropospheric carbon monoxide over the Pacific during HIPPO: two-way coupled simulation of GEOS-Chem and its multiple nested models Atmospheric Chemistry and Physics, 14, 12649-12663, 2014 Author(s): Y.-Y. Yan, J.-T. Lin, Y. Kuang, D. Yang, and L. Zhang Global chemical transport models (CTMs) are used extensively to study air pollution and transport at a global scale. These models are limited by coarse horizontal resolutions that do not allow for a detailed representation of small-scale nonlinear processes over the pollutant source regions. Here we couple the global GEOS-Chem CTM and its three high-resolution nested models to simulate the tropospheric carbon monoxide (CO) over the Pacific Ocean during five High-performance Instrumented Airborne Platform for Environmental Research (HIAPER) Pole-to-Pole Observations (HIPPO) campaigns between 2009 and 2011. We develop a two-way coupler, the PeKing University CouPLer (PKUCPL), allowing for the exchange and interaction of chemical constituents between the global model (at 2.5° long. × 2° lat.) and the three nested models (at 0.667° long. × 0.5° lat.) covering Asia, North America, and Europe. The coupler obtains nested model results to modify the global model simulation within the respective nested domains, and simultaneously acquires global model results to provide lateral boundary conditions (LBCs) for the nested models. Compared to the global model alone, the two-way coupled simulation results in enhanced CO concentrations in the nested domains. Sensitivity tests suggest the enhancement to be a result of improved representation of the spatial distributions of CO, nitrogen oxides, and non-methane volatile organic compounds, the meteorological dependence of natural emissions, and other resolution-dependent processes. The relatively long lifetime of CO allows for the enhancement to be accumulated and carried across the globe. We found that the two-way coupled simulation increased the global tropospheric mean CO concentrations in 2009 by 10.4%, with a greater enhancement at 13.3% in the Northern Hemisphere. Coincidently, the global tropospheric mean hydroxyl radical (OH) was reduced by 4.2%, resulting in a 4.2% enhancement in the methyl chloroform lifetime (MCF; via reaction with tropospheric OH). The resulting CO and OH contents and MCF lifetime are closer to observation-based estimates. Both the global and the two-way coupled models capture the general spatiotemporal patterns of HIPPO CO over the Pacific. The two-way coupled simulation is much closer to HIPPO CO, with a mean bias of 1.1 ppb (1.4%) below 9 km compared to the bias at −7.2 ppb (−9.2%) for the global model alone. The improvement is most apparent over the North Pacific. Our test simulations show that the global model alone could resemble the two-way coupled simulation (especially below 4 km) by increasing its global CO emissions by 15% for HIPPO-1 and HIPPO-3, by 25% for HIPPO-2 and HIPPO-4, and by 35% for HIPPO-5. This has important implications for using the global model alone to constrain CO emissions. Thus, the two-way coupled simulation is a significantly improved model tool for studying the global impacts of air pollutants from major anthropogenic source regions.

Description: Testing gridded land precipitation data and precipitation and runoff reanalyses (1982–2010) between 45° S and 45° N with Normalized Difference Vegetation Index data Hydrology and Earth System Sciences Discussions, 11, 13175-13205, 2014 Author(s): S. O. Los The realistic simulation of key components of the land-surface hydrological cycle – precipitation, runoff, evaporation and transpiration – in general circulation models of the atmosphere is crucial to assess adverse weather impacts on environment and society. Here, gridded precipitation data from observations and precipitation and runoff fields from reanalyses were tested with satellite-derived global vegetation index data for 1982–2010 and latitudes between 45° S and 45° N. Data were obtained from the Climate Research Unit (CRU), the Global Precipitation Climatology Project (GPCP) and Tropical Rainfall Monitoring Mission (TRMM; analysed for 1998–2010 only) and (precipitation and runoff) reanalyses were obtained from the National Centers for Environmental Prediction/National Center for Atmospheric Research (NCEP/NCAR), the European Centre for Medium-Range Weather Forecasts (ECMWF) and the NASA Global Modelling and Assimilation Office (GMAO). Annual land-surface precipitation was converted to annual potential vegetation net primary productivity (NPP) and was compared to mean annual Normalized Difference Vegetation Index data measured by the Advanced Very High Resolution Radiometer (1982–1999) and MODIS (2001–2010). The effect of spatial resolution on the agreement between NPP and NDVI was investigated as well. The CRU and TRMM derived NPP agreed most closely with the NDVI data. The GPCP data showed weaker spatial agreement, largely because of their lower spatial resolution, but similar temporal agreement. MERRA Land and ERA Interim precipitation reanalyses showed similar spatial agreement as the GPCP data and good temporal agreement in semi-arid regions of the Americas, Asia, Australia and southern Africa. The NCEP/NCAR reanalysis showed the lowest spatial agreement which could only in part be explained by its lower spatial resolution. No reanalysis showed realistic interannual precipitation variations for northern tropical Africa. Inclusion of runoff in the NPP prediction resulted only in (marginally) better agreement for the MERRA Land reanalysis and worse agreement for the NCEP/NCAR and ERA Interim reanalyses.