ALBERT

Description: Predicting future thaw slump activity requires a sound understanding of the atmospheric drivers and geomorphic controls on mass wasting across a range of time scales. On sub-seasonal time scales, sparse measurements indicate that mass wasting at active slumps is often limited by the energy available for melting ground ice, but other factors such as rainfall or the formation of an insulating veneer may also be relevant. To study the sub-seasonal drivers, we derive topographic changes from single-pass radar interferometric data acquired by the TanDEM-X satellite (12 m resolution). The high vertical precision (around 30 cm), frequent observations (11 days) and large coverage (5000 km2) allow us to track volume losses as drivers such as the available energy change during summer in two study regions. We find that thaw slumps in the Tuktoyaktuk coastlands, Canada, are not energy limited in June, as they undergo limited mass wasting (height loss of around 0 cm/day) despite the ample available energy, indicating the widespread presence of an insulating snow or debris veneer. Later in summer, height losses generally increase (around 3 cm/day), but they do so in distinct ways. For many slumps, mass wasting tracks the available energy, a temporal pattern that is also observed at coastal yedoma cliffs on the Bykovsky Peninsula, Russia. However, the other two common temporal trajectories are asynchronous with the available energy, as they track strong precipitation events or show a sudden speed-up in late August, respectively. The observed temporal patterns are poorly related to slump characteristics like the slump area. The contrasting temporal behaviour of nearby thaw slumps highlights the importance of complex local and temporally varying controls on mass wasting.

Description: A suite of oxygenated volatile organic compounds (OVOCs – acetaldehyde, acetone, propanal, butanal and butanone) were measured concurrently in the surface water and atmosphere of the South China Sea and Sulu Sea in November 2011. A strong correlation was observed between all OVOC concentrations in the surface seawater along the entire cruise track, except for acetaldehyde, suggesting similar sources and sinks in the surface ocean. Additionally, several phytoplankton groups, such as haptophytes or pelagophytes, were also correlated to all OVOCs indicating that phytoplankton may be an important source for marine OVOCs in the South China and Sulu Seas. Humic and protein like fluorescent dissolved organic matter (FDOM) components seemed to be additional precursors for butanone and acetaldehyde. The atmospheric OVOC mixing ratios were relative high compared with literature values, suggesting the coastal region of North Borneo as a local hot spot for atmospheric OVOCs. The flux of atmospheric OVOCs was largely into the ocean for all 5 gases, with a few important exceptions near the coast of Borneo. The calculated amount of OVOCs entrained into the ocean seemed to be an important source of OVOCs to the surface ocean. When the fluxes were out of the ocean, marine OVOCs were found to be enough to control the local measured OVOC distribution in the atmosphere. Based on our model calculations, at least 0.4 ppb of marine derived acetone and butanone can reach the upper troposphere, where they may have an important influence on hydrogen oxide radical formation over the western Pacific Ocean.

Description: The Extrapolar SWIFT model is a fast ozone chemistry scheme for interactive calculation of the extrapolar stratospheric ozone layer in coupled general circulation models (GCMs). In contrast to the widely used prescribed ozone, the SWIFT ozone layer interacts with the model dynamics and can respond to atmospheric variability or climatological trends. The Extrapolar SWIFT model employs a repro-modelling approach, where algebraic functions are used to approximate the numerical output of a full stratospheric chemistry and transport model (ATLAS). The full model solves a coupled chemical differential equations system with 55 initial and boundary conditions (mixing ratio of various chemical species and atmospheric parameters). Hence the rate of change of ozone over 24  h is a function of 55 variables. Using covariances between these variables, we can find linear combinations in order to reduce the parameter space to the following nine basic variables: latitude, pressure altitude, temperature, local ozone column, mixing ratio of ozone and of the ozone depleting families (Cly, Bry, NOy and HOy). We will show that these 9 variables are sufficient to characterize the rate of change of ozone. An automated procedure fits a polynomial function of fourth degree to the rate of change of ozone obtained from several simulations with the ATLAS model. One polynomial function is determined per month which yields the rate of change of ozone over 24 h. A key aspect for the robustness of the Extrapolar SWIFT model is to include a wide range of stratospheric variability in the numerical output of the ATLAS model, also covering atmospheric states that will occur in a future climate (e.g. temperature and meridional circulation changes or reduction of stratospheric chlorine loading). For validation purposes, the Extrapolar SWIFT model has been integrated into the ATLAS model replacing the full stratospheric chemistry scheme. Simulations with SWIFT in ATLAS have proven that the systematic error is small and does not accumulate during the course of a simulation. In the context of a 10 year simulation, the ozone layer, simulated by SWIFT, shows a stable annual cycle, with inter-annual variations comparable to the ATLAS model. The application of Extrapolar SWIFT requires the evaluation of polynomial functions with 30–100 terms. Nowadays, computers can calculate such polynomial functions at thousands of model grid points in seconds. SWIFT provides the desired numerical efficiency and computes the ozone layer 104 times faster than the chemistry scheme in the ATLAS CTM.

Description: Ice retreat in the eastern Eurasian Arctic is a consequence of atmospheric and oceanic processes and regional feedback mechanisms acting on the ice cover, both in winter and summer. A correct representation of these processes in numerical models is important, since it will improve predictions of sea ice anomalies along the Northeast Passage and beyond. In this study, we highlight the importance of winter ice dynamics for local summer sea ice anomalies in thickness, volume and extent. By means of airborne sea ice thickness surveys made over pack ice areas in the south-eastern Laptev Sea, we show that years of offshore-directed sea ice transport have a thinning effect on the late-winter sea ice cover. To confirm the preconditioning effect of enhanced offshore advection in late winter on the summer sea ice cover, we perform a sensitivity study using a numerical model. Results verify that the preconditioning effect plays a bigger role for the regional ice extent. Furthermore, they indicate an increase in volume export from the Laptev Sea as a consequence of enhanced offshore advection, which has far-reaching consequences for the entire Arctic sea ice mass balance. Moreover we show that ice dynamics in winter not only preconditions local summer ice extent, but also accelerate fast-ice decay.

Description: Permafrost temperatures are increasing in Alaska due to climate change and in some cases permafrost is thawing and degrading. In areas where degradation has already occurred the effects can be dramatic, resulting in changing ecosystems, carbon release, and damage to infrastructure. However, in many areas we lack baseline data, such as subsurface temperatures, needed to assess future changes and potential risk areas. Besides climate, the physical properties of the vegetation cover and subsurface material have a major influence on the thermal state of permafrost. These properties are often directly related to the type of ecosystem overlaying permafrost. In this paper we demonstrate that classifying the landscape into general ecotypes is an effective way to scale up permafrost thermal data collected from field monitoring sites. Additionally, we find that within some ecotypes the absence of a moss layer is indicative of the absence of near-surface permafrost. As a proof of concept, we used the ground temperature data collected from the field sites to recode an ecotype land cover map into a map of mean annual ground temperature ranges at 1 m depth based on analysis and clustering of observed thermal regimes. The map should be useful for decision making with respect to land use and understanding how the landscape might change under future climate scenarios.

Description: Climate trends in the Antarctic region remain poorly characterized, owing to the brevity and scarcity of direct climate observations and the large magnitude of interannual to decadal-scale climate variability. Here, within the framework of the PAGES Antarctica2k working group, we build an enlarged database of ice core water stable isotope records from Antarctica, consisting of 112 records. We produce both unweighted and weighted isotopic (δ18O) composites and temperature reconstructions since 0 CE, binned at 5- and 10-year resolution, for seven climatically distinct regions covering the Antarctic continent. Following earlier work of the Antarctica2k working group, we also produce composites and reconstructions for the broader regions of East Antarctica, West Antarctica and the whole continent. We use three methods for our temperature reconstructions: (i) a temperature scaling based on the δ18O–temperature relationship output from an ECHAM5-wiso model simulation nudged to ERA-Interim atmospheric reanalyses from 1979 to 2013, and adjusted for the West Antarctic Ice Sheet region to borehole temperature data, (ii) a temperature scaling of the isotopic normalized anomalies to the variance of the regional reanalysis temperature and (iii) a composite-plus-scaling approach used in a previous continent-scale reconstruction of Antarctic temperature since 1 CE but applied to the new Antarctic ice core database. Our new reconstructions confirm a significant cooling trend from 0 to 1900 CE across all Antarctic regions where records extend back into the 1st millennium, with the exception of the Wilkes Land coast and Weddell Sea coast regions. Within this long-term cooling trend from 0 to 1900 CE, we find that the warmest period occurs between 300 and 1000 CE, and the coldest interval occurs from 1200 to 1900 CE. Since 1900 CE, significant warming trends are identified for the West Antarctic Ice Sheet, the Dronning Maud Land coast and the Antarctic Peninsula regions, and these trends are robust across the distribution of records that contribute to the unweighted isotopic composites and also significant in the weighted temperature reconstructions. Only for the Antarctic Peninsula is this most recent century-scale trend unusual in the context of natural variability over the last 2000 years. However, projected warming of the Antarctic continent during the 21st century may soon see significant and unusual warming develop across other parts of the Antarctic continent. The extended Antarctica2k ice core isotope database developed by this working group opens up many avenues for developing a deeper understanding of the response of Antarctic climate to natural and anthropogenic climate forcings. The first long-term quantification of regional climate in Antarctica presented herein is a basis for data–model comparison and assessments of past, present and future driving factors of Antarctic climate.

Description: Permafrost presence is determined by a complex interaction of climatic, topographic, and ecological conditions operating over long time scales. In particular, vegetation and organic layer characteristics may act to protect permafrost in regions with a mean annual air temperature (MAAT) above 0°C. In this study, we document the presence of residual permafrost plateaus in the western Kenai Peninsula lowlands of south-central Alaska, a region with a MAAT of 1.5+/-1 °C (1981–2010). Continuous ground temperature measurements between 16 September 2012 and 15 September 2015, using calibrated thermistor strings, documented the presence of warm permafrost (-0.04 to -0.08 °C). Field measurements (probing) on several plateau features during the fall of 2015 showed that the depth to the permafrost table averaged 1.48m but at some locations was as shallow as 0.53 m. Late winter surveys (augering, coring, and GPR) in 2016 showed that the average seasonally frozen ground thickness was 0.45 m, overlying a talik above the permafrost table. Measured permafrost thickness ranged from 0.33 to 〉6.90 m. Manual interpretation of historic aerial photography acquired in 1950 indicates that residual permafrost plateaus covered 920 ha as mapped across portions of four wetland complexes encompassing 4810 ha. However, between 1950 and ca. 2010, permafrost plateau extent decreased by 60.0 %, with lateral feature degradation accounting for 85.0% of the reduction in area. Permafrost loss on the Kenai Peninsula is likely associated with a warming climate, wildfires that remove the protective forest and organic layer cover, groundwater flow at depth, and lateral heat transfer from wetland surface waters in the summer. Better understanding the resilience and vulnerability of ecosystem-protected permafrost is critical for mapping and predicting future permafrost extent and degradation across all permafrost regions that are currently warming. Further work should focus on reconstructing permafrost history in south-central Alaska as well as additional contemporary observations of these ecosystem-protected permafrost sites south of the regions with relatively stable permafrost.

Description: Permafrost presence is determined by a complex interaction of climatic, topographic, and ecological conditions operating over long time scales. In particular, vegetation and organic layer characteristics may act to protect permafrost in regions with a mean annual air temperature (MAAT) above 0 °C. In this study, we document the presence of residual permafrost plateaus on the western Kenai Peninsula lowlands of southcentral Alaska, a region with a MAAT of 1.5 ± 1 °C (1981 to 2010). Continuous ground temperature measurements between 16 September 2012 and 15 September 2015, using calibrated thermistor strings, documented the presence of warm permafrost (−0.04 to −0.08 °C). Field measurements (probing) on several plateau features during the fall of 2015 showed that the depth to the permafrost table averaged 1.48 m but was as shallow as 0.53 m. Late winter surveys (drilling, coring, and GPR) in 2016 showed that the average seasonally frozen ground thickness was 0.45 m, overlying a talik above the permafrost table. Measured permafrost thickness ranged from 0.33 to 〉 6.90 m. Manual interpretation of historic aerial photography acquired in 1950 indicates that residual permafrost plateaus covered 920 ha as mapped across portions of four wetland complexes encompassing 4810 ha. However, between 1950 and ca. 2010, permafrost plateau extent decreased by 60 %, with lateral feature degradation accounting for 85 % of the reduction in area. Permafrost loss on the Kenai Peninsula is likely associated with a warming climate, wildfires that remove the protective forest and organic layer cover, groundwater flow at depth, and lateral heat transfer from wetland surface waters in the summer. Better understanding the resilience and vulnerability of ecosystem-protected permafrost is critical for mapping and predicting future permafrost extent and degradation across all permafrost regions that are currently warming. Further work should focus on reconstructing permafrost history in southcentral Alaska as well as additional contemporary observations of these ecosystem-protected permafrost sites lying south of the regions with relatively stable permafrost.

Description: Coastal erosion and flooding transform terrestrial landscapes into marine environments. In the Arctic, these processes inundate terrestrial permafrost with seawater and create submarine permafrost. Permafrost begins to warm under marine conditions, which can destabilize the sea floor and may release greenhouse gases. We report on the transition of terrestrial to submarine permafrost at a site where the timing of inundation can be inferred from the rate of coastline retreat. On Muostakh Island in the central Laptev Sea, East Siberia, changes in annual coastline position have been measured for decades and vary highly spatially. We hypothesize that these rates are inversely related to the inclination of the upper surface of submarine ice-bonded permafrost (IBP) based on the consequent duration of inundation with increasing distance from the shoreline. We compared rapidly eroding and stable coastal sections of Muostakh Island and find permafrost-table inclinations, determined using direct current resistivity, of 1 and 5 %, respectively. Determinations of submarine IBP depth from a drilling transect in the early 1980s were compared to resistivity profiles from 2011. Based on borehole observations, the thickness of unfrozen sediment overlying the IBP increased from 0 to 14m below sea level with increasing distance from the shoreline. The geoelectrical profiles showed thickening of the unfrozen sediment overlying ice-bonded permafrost over the 28 years since drilling took place. We use geoelectrical estimates of IBP depth to estimate permafrost degradation rates since inundation. Degradation rates decreased from over 0.4ma-1 following inundation to around 0.1ma-1 at the latest after 60 to 110 years and remained constant at this level as the duration of inundation increased to 250 years. We suggest that long-term rates are lower than these values, as the depth to the IBP increases and thermal and porewater solute concentration gradients over depth decrease. For the study region, recent increases in coastal erosion rate and changes in benthic temperature and salinity regimes are expected to affect the depth to submarine permafrost, leading to coastal regions with shallower IBP.

Description: Sources and Transformations of Anthropogenic Nitrogen along an Urban River-Estuarine Continuum Michael J. Pennino, Sujay S. Kaushal, Sudhir Murthy, Joel Blomquist, Jeff Cornwell, and Lora Harris Biogeosciences Discuss., doi:10.5194/bg-2016-264,2016 Manuscript under review for BG (discussion: open, 0 comments) The results of this manuscript report the analysis of the fate and transport of wastewater and anthropogenic nitrogen along the Potomac River Estuary, from Washington D.C. to the Chesapeake Bay. In conjunction with a mass balance approach, nitrate isotopes were used to estimate fluxes and trace the sources and transformations N along the estuary. This study shows that estuaries have a large capacity to transform N inputs, but with large seasonal variability due to hydrological extremes.

Description: Boreal fire records in Northern Hemisphere ice cores: A review Michel Legrand, Joseph McConnell, Hubertus Fischer, Eric W. Wolff, Susanne Preunkert, Nathan Chellman, Daiana Leuenberger, Olivia Maselli, Michael Sigl, Simon Schüpbach, and Mike Flannigan Clim. Past Discuss., doi:10.5194/cp-2016-70,2016 Manuscript under review for CP (discussion: open, 0 comments) Here we review previous attempts made to reconstruct past forest fire using chemical signals recorded in Greenland ice. We showed that the Greenland ice records of ammonium, found to be a good fire proxy, consistently indicate changing fire activity in Canada in response to past climatic conditions that occurred since the last 15 000 years including the little ice age and the last large climatic transition.

Description: Tectonothermal evolution in the core of an arcuate fold and thrust belt: the south-eastern sector of the Cantabrian Zone (Variscan belt, north-western Spain) María Luz Valín, Susana García-López, Covadonga Brime, Fernando Bastida, and Jesús Aller Solid Earth, 7, 1003-1022, doi:10.5194/se-7-1003-2016, 2016 The tectonothermal evolution of an area located in the core of the Ibero-Armorican Arc (Variscan belt) has been determined by using the conodont colour alteration index (CAI), Kübler index of illite (KI), the Árkai index of chlorite (AI) and the analysis of clay minerals and rock cleavage. The area is part of the Cantabrian Zone (CZ), which represents the foreland fold and thrust belt of the orogen. It has been thrust by several large units of the CZ, what resulted in the generation of a large number of synorogenic Carboniferous sediments. CAI, KI and AI values show an irregular distribution of metamorphic grade, independent of stratigraphic position. Two tectonothermal events have been distinguished in the area. The first one, poorly defined, is mainly located in the northern part. It gave rise to very-low-grade metamorphism in some areas and it was associated with a deformation event that resulted in the emplacement of the last large thrust unit and development of upright folds and associated cleavage ( S 1 ). The second tectonothermal event gave rise to low-grade metamorphism and cleavage ( S 2 ) crosscutting earlier upright folds in the central, western and southern parts of the study area. The event continued with the intrusion of small igneous rock bodies, which gave rise to contact metamorphism and hydrothermal alteration. This event was linked to an extensional episode due to a gravitational instability at the end of the Variscan deformation. This tectonothermal evolution occurred during the Gzhelian–Sakmarian. Subsequently, several hydrothermal episodes took place and local crenulation cleavage developed during the Alpine deformation.

Description: Microstructures and deformation mechanisms in Opalinus Clay: insights from scaly clay from the Main Fault in the Mont Terri Rock Laboratory (CH) Ben Laurich, Janos L. Urai, and Christophe Nussbaum Solid Earth Discuss., doi:10.5194/se-2016-94,2016 Manuscript under review for SE (discussion: open, 0 comments) Scaly clay is a well-known rock fabric that can develop in tectonic systems and that can alter the physical rock properties of a formation. However, internal microstructure and evolution of this fabric remains poorly understood. We examined the scaly microstructure of progressively faulted Opalinus Clay using optical as well as scanning electron microscopy. We present an evolutionary model for scaly clay that is of interest in the evaluation of Opalinus Clay to host radioactive waste.

Description: Microstructures and deformation mechanisms in Opalinus Clay: insights from scaly clay from the Main Fault in the Mont Terri Rock Laboratory (CH) Ben Laurich, Janos L. Urai, and Christophe Nussbaum Solid Earth Discuss., doi:10.5194/se-2016-94,2016 Manuscript under review for SE (discussion: open, 0 comments) Scaly clay is a well-known rock fabric that can develop in tectonic systems and that can alter the physical rock properties of a formation. However, internal microstructure and evolution of this fabric remains poorly understood. We examined the scaly microstructure of progressively faulted Opalinus Clay using optical as well as scanning electron microscopy. We present an evolutionary model for scaly clay that is of interest in the evaluation of Opalinus Clay to host radioactive waste.

Description: Fixation kinetics of chelated and non-chelated zinc in semi-arid alkaline soils: application to zinc management Theophilus K. Udeigwe, Madeleine Eichmann, and Matthew C. Menkiti Solid Earth, 7, 1023-1031, doi:10.5194/se-7-1023-2016, 2016 Micronutrient fixation leads to the reduction of plant-available portions. This study examined the fixation kinetics of chelated and non-chelated zinc (Zn) in semi-arid soils. About 30 % more Zn was fixed in the non-chelated system in the first 14 days. Slope of change of Zn to Cu (i.e., Cu–Zn) was higher than other micronutrients. Zn fixation kinetics were better described by a power-function model. Tools developed from this study will be used to advance micronutrient management.

Description: Improved MODIS Dark Target Aerosol Optical Depth algorithm over land: Angular effect Correction Yerong Wu, Martin de Graaf, and Massimo Menenti Atmos. Meas. Tech. Discuss., doi:10.5194/amt-2016-185,2016 Manuscript under review for AMT (discussion: open, 0 comments) In this paper, we developed a new algorithm to improve the retrieval of Aerosol Optical Depth (AOD) over land with satellite measurement, by refining the shape of the spectral surface reflectance to improve the estimated background reflectance in the channel used for AOD retrieval. The results show that the angular effects of the retrievals are largely reduced, including fewer occurrences of negative retrievals. This implies that the users can get more accurate data without angular bias.

Description: Investigation of processes controlling summertime gaseous elemental mercury oxidation at midlatitudinal marine, coastal, and inland sites Zhuyun Ye, Huiting Mao, Che-Jen Lin, and Su Youn Kim Atmos. Chem. Phys., 16, 8461-8478, doi:10.5194/acp-16-8461-2016, 2016 In this study, a state-of-the-art chemical mechanism was incorporated into a box model to investigate the atmospheric Hg cycling in different environments. As a result, for each of the three environments, GOM diurnal cycles of over half the selected cases were reasonably represented by the box model. A realistic model can be a powerful tool, providing important information on atmospheric Hg cycling and implications for policy makers.

Description: Effects of polar stratospheric clouds in the Nimbus 7 LIMS Version 6 data set Ellis Remsberg and V. Lynn Harvey Atmos. Meas. Tech., 9, 2927-2946, doi:10.5194/amt-9-2927-2016, 2016 Emissions from polar stratospheric cloud (PSC) particles affect the retrieved ozone and water vapor from the Limb Infrared Monitor of the Stratosphere (LIMS) satellite experiment. Threshold criteria are applied to the retrieved ozone for the detection and screening of those effects. The PSC effects correlate very well with regions of coldest temperatures (

Description: The 1816 ‘year without a summer’ in an atmospheric reanalysis Philip Brohan, Gilbert P. Compo, Stefan Brönnimann, Robert J. Allan, Renate Auchmann, Yuri Brugnara, Prashant D. Sardeshmukh, and Jeffrey S. Whitaker Clim. Past Discuss., doi:10.5194/cp-2016-78,2016 Manuscript under review for CP (discussion: open, 0 comments) We have used modern weather forecasting tools to reconstruct the dreadful European weather of 200 years ago – 1816 was the ‘year without a summer’; harvests failed, and people starved. We can show that 1816’s extreme climate was caused by the eruption of the Tambora volcano the previous year. This means we have some chance of predicting such extreme summers if they occur in future, though this is still a challenge to today’s forecast models.

Description: Palynological evidence for late Miocene stepwise aridification on the northeastern Tibetan Plateau Jia Liu, Ji Jun Li, Chun Hui Song, Hao Yu, Ting Jiang Peng, Zheng Chuang Hui, and Xi Yan Ye Clim. Past, 12, 1473-1484, doi:10.5194/cp-12-1473-2016, 2016 The late Cenozoic basins in the northeastern Tibetan Plateau document both the tectonic uplift process and its associated environmental changes. Here, we investigated a late Miocene sporopollen record from the Tianshui Basin in the northeastern Tibetan Plateau. The results show that a persistent aridification trend parallels the global cooling of the late Miocene, and the stepwise vegetation succession is consistent with the major uplift events of the Tibetan Plateau.

Description: Gone or just out of sight? The apparent disappearance of aromatic litter components in soils Thimo Klotzbücher, Karsten Kalbitz, Chiara Cerli, Peter J. Hernes, and Klaus Kaiser SOIL, 2, 325-335, doi:10.5194/soil-2-325-2016, 2016 Uncertainties concerning stabilization of organic compounds in soil limit our basic understanding on soil organic matter (SOM) formation and our ability to model and manage effects of global change on SOM stocks. One controversially debated aspect is the contribution of aromatic litter components, such as lignin and tannins, to stable SOM forms. Here, we summarize and discuss the inconsistencies and propose research options to clear them.

Description: Structure of Suasselkä Postglacial Fault in northern Finland obtained by analysis of local events and ambient seismic noise Nikita Afonin, Elena Kozlovskaya, Ilmo Kukkonen, and DAFNE/FINLAND Working Group Solid Earth Discuss., doi:10.5194/se-2016-90,2016 Manuscript under review for SE (discussion: open, 0 comments) Understanding inner structure of seismogenic faults and their ability to reactivate is particularly important in investigating continental intraplate seismicity regime. In our study we address this problem using analysis of local seismic events and ambient seismic noise recorded by the temporary DAFNE array in northern Fennoscandian Shield. The main purpose of the DAFNE/FINLAND passive seismic array experiment was to characterize the present-day seismicity of the Suasselkä post-glacial fault (SPGF) that was proposed as one potential target for the DAFNE (Drilling Active Faults in Northern Europe) project. The DAFNE/FINLAND array comprised the area of about 20 to 100 km and consisted of 8 short-period and 4 broad-band 3-component autonomous seismic stations installed in the close vicinity of the fault area. The array recorded continuous seismic data during September, 2011–May, 2013. Recordings of the array have being analyzed in order to identify and locate natural earthquakes from the fault area and to discriminate them from the blasts in the Kittilä Gold Mine. As a result, we found several dozens of natural seismic events originating from the fault area, which proves that the fault is still seismically active. In order to study the inner structure of the SPGF we use cross-correlation of ambient seismic noise recorded by the array. Analysis of azimuthal distribution of noise sources demonstrated that during the time interval under consideration the distribution of noise sources is close to the uniform one. The continuous data were processed in several steps including single station data analysis, instrument response removal and time-domain stacking. The data were used to estimate empirical Green’s functions between pairs of stations in the frequency band of 0.1–1 Hz and to calculate correspondent surface wave dispersion curves. The S-wave velocity models were obtained as a result of dispersion curves inversion. The results suggest that the area of the SPGF corresponds to a narrow region of low S-wave velocities surrounded by rocks with high S-wave velocities. We interpret this low velocity region as a non-healed mechanically weak fault damage zone (FDZ) that remained after the last major earthquake that occurred after the last glaciation.

Description: Structure of Suasselkä Postglacial Fault in northern Finland obtained by analysis of local events and ambient seismic noise Nikita Afonin, Elena Kozlovskaya, Ilmo Kukkonen, and DAFNE/FINLAND Working Group Solid Earth Discuss., doi:10.5194/se-2016-90,2016 Manuscript under review for SE (discussion: open, 0 comments) Understanding inner structure of seismogenic faults and their ability to reactivate is particularly important in investigating continental intraplate seismicity regime. In our study we address this problem using analysis of local seismic events and ambient seismic noise recorded by the temporary DAFNE array in northern Fennoscandian Shield. The main purpose of the DAFNE/FINLAND passive seismic array experiment was to characterize the present-day seismicity of the Suasselkä post-glacial fault (SPGF) that was proposed as one potential target for the DAFNE (Drilling Active Faults in Northern Europe) project. The DAFNE/FINLAND array comprised the area of about 20 to 100 km and consisted of 8 short-period and 4 broad-band 3-component autonomous seismic stations installed in the close vicinity of the fault area. The array recorded continuous seismic data during September, 2011–May, 2013. Recordings of the array have being analyzed in order to identify and locate natural earthquakes from the fault area and to discriminate them from the blasts in the Kittilä Gold Mine. As a result, we found several dozens of natural seismic events originating from the fault area, which proves that the fault is still seismically active. In order to study the inner structure of the SPGF we use cross-correlation of ambient seismic noise recorded by the array. Analysis of azimuthal distribution of noise sources demonstrated that during the time interval under consideration the distribution of noise sources is close to the uniform one. The continuous data were processed in several steps including single station data analysis, instrument response removal and time-domain stacking. The data were used to estimate empirical Green’s functions between pairs of stations in the frequency band of 0.1–1 Hz and to calculate correspondent surface wave dispersion curves. The S-wave velocity models were obtained as a result of dispersion curves inversion. The results suggest that the area of the SPGF corresponds to a narrow region of low S-wave velocities surrounded by rocks with high S-wave velocities. We interpret this low velocity region as a non-healed mechanically weak fault damage zone (FDZ) that remained after the last major earthquake that occurred after the last glaciation.

Description: Comparison of box counting and correlation dimension methods in well logging data analysis associate with the texture of volcanic rocks D. Mou and Z. W. Wang Nonlin. Processes Geophys. Discuss., doi:10.5194/npg-2014-85,2016 Manuscript under review for NPG (discussion: open, 0 comments) We have developed a fractal analysis method to estimate the dimension of well logging curves in Liaohe oil field, China. The box counting and correlation dimension are methods that can be applied to predict the texture of volcanic rocks with calculation the fractal dimension of well logging curves. The well logging curves are composed of gamma ray (GR), compensated neutron logs (CNL), acoustic (AC), density (DEN), Resistivity lateral log deep ( R LLD ), every curve contains a total of 6000 logging data. The dimension of well logging curves are calculated using box counting and correlation algorithms respectively. It is shown that two types of dimension of CNL, DEN and AC have the same average value. The box counting dimension of volcanic lava is lower than the pyroclastic rock obviously. The majority of correlation dimension of volcanic lava is lower than the pyroclastic rock, but a small amount of correlation dimension of volcanic lava is equal to the pyroclastic rock. It is demonstrated that the box counting dimension is more suitable for predicting the texture of volcanic rocks. Applications to logging data, A well show the relationship between the fractal dimension and the texture of volcanic rock in certain depth.

Description: First year of the new Arctic AWIPEV-COSYNA cabled Underwater Observatory in Kongsfjorden, Spitsbergen Philipp Fischer, Max Schwanitz, Reiner Loth, Uwe Posner, Markus Brand, and Friedhelm Schröder Ocean Sci. Discuss., doi:10.5194/os-2016-52,2016 Manuscript under review for OS (discussion: open, 0 comments) We observed oceanographic and community data from October 2013 to November 2014 in the shallow waters of Kongsfjorden on the west coast of Svalbard (Norway) using remote controlled hydrographic and optic sensors. Daily vertical profiles of temperature, salinity and turbidity were sampled together with stereo images of the macrobiotic community, including fish. A distinct seasonal cycle in total species abundances was found with surprisingly high animal counts during the polar winter.

Description: Landfast ice thickness in the Canadian Arctic Archipelago from observations and models Stephen E. L. Howell, Frédéric Laliberté, Ron Kwok, Chris Derksen, and Joshua King The Cryosphere, 10, 1463-1475, doi:10.5194/tc-10-1463-2016, 2016 The Canadian Ice Service record of observed landfast ice and snow thickness represents one of the longest in the Arctic that spans over 5 decades. We analyze this record to report on long-term trends and variability of ice and snow thickness within the Canadian Arctic Archipelago (CAA). Results indicate a thinning of ice at several sites in the CAA. State-of-the-art climate models still have difficultly capturing observed ice thickness values in the CAA and should be used with caution.

Description: Coastal dynamics and submarine permafrost in shallow water of the central Laptev Sea, East Siberia Pier Paul Overduin, Sebastian Wetterich, Frank Günther, Mikhail N. Grigoriev, Guido Grosse, Lutz Schirrmeister, Hans-Wolfgang Hubberten, and Aleksandr Makarov The Cryosphere, 10, 1449-1462, doi:10.5194/tc-10-1449-2016, 2016 How fast does permafrost warm up and thaw after it is covered by the sea? Ice-rich permafrost in the Laptev Sea, Siberia, is rapidly eroded by warm air and waves. We used a floating electrical technique to measure the depth of permafrost thaw below the sea, and compared it to 60 years of coastline retreat and permafrost depths from drilling 30 years ago. Thaw is rapid right after flooding of the land and slows over time. The depth of permafrost is related to how fast the coast retreats.

Description: Age of the Mt. Ortles ice cores, the Tyrolean Iceman and glaciation of the highest summit of South Tyrol since the Northern Hemisphere Climatic Optimum Paolo Gabrielli, Carlo Barbante, Giuliano Bertagna, Michele Bertó, Daniel Binder, Alberto Carton, Luca Carturan, Federico Cazorzi, Giulio Cozzi, Giancarlo Dalla Fontana, Mary Davis, Fabrizio De Blasi, Roberto Dinale, Gianfranco Dragà, Giuliano Dreossi, Daniela Festi, Massimo Frezzotti, Jacopo Gabrieli, Stephan Galos, Patrick Ginot, Petra Heidenwolf, Theo M. Jenk, Natalie Kehrwald, Donald Kenny, Olivier Magand, Volkmar Mair, Vladimir Mikhalenko, Ping Nan Lin, Klaus Oeggl, Gianni Piffer, Mirko Rinaldi, Ulrich Schotterer, Margit Schwikowski, Roberto Seppi, Andrea Spolaor, Barbara Stenni, Davdi Tonidandel, CChiara Uglietti, Victor Zagorodnov, Thomas Zanoner, and Piero Zennaro The Cryosphere Discuss., doi:10.5194/tc-2016-159,2016 Manuscript under review for TC (discussion: open, 0 comments) New ice cores were extracted from Alto dell'Ortles, the highest glacier of South Tyrol in the Italian Alps, to check whether prehistoric ice, that is coeval to the famous 5200 years old Tyrolean Iceman, is still preserved in this region. Dating of the ice cores confirms the hypothesis and indicates that the drilling site was glaciated since the end of the Northern Hemisphere Climatic Optimum (7000 BP). We also infer that an unprecedented acceleration of the glacier flow has just recently begun.

Description: Combined diurnal variations of discharge and hydrochemistry of the Isunnguata Sermia outlet of the Greenland Ice Sheet give in sight on sub glacial conditions Joseph Graly, Joel Harrington, and Neil Humphrey The Cryosphere Discuss., doi:10.5194/tc-2016-137,2016 Manuscript under review for TC (discussion: open, 0 comments) At a major outlet of the Greenland Ice Sheet in West Greenland, we find that the chemical solutes is the emerging subglacial waters are out of phase with water discharge and can spike in concentration during waning flow. This suggests that the subglacial waters are spreading out across a large area of the glacial bed throughout the day, stimulating chemical weathering beyond the major water distribution channels.

Description: Why did the storm ex-Gaston (2010) fail to redevelop during the PREDICT experiment? Thomas M. Freismuth, Blake Rutherford, Mark A. Boothe, and Michael T. Montgomery Atmos. Chem. Phys., 16, 8511-8519, doi:10.5194/acp-16-8511-2016, 2016 Numerical model analyses are used to investigate the role of dry, environmental air in the failed redevelopment of a tropical cyclone (ex-Gaston, 2010). As early as 12:00 UTC 2 September 2010, a dry layer at and above 600 hPa results in a decrease in the vertical mass flux and vertical, relative vorticity. The intrusion of dry air led to a reduction in vorticity and a compromised pouch at these middle levels. This study supports work looking at the role of dry air in moist convection.

Description: A comparative analysis of UV nadir-backscatter and infrared limb-emission ozone data assimilation Rossana Dragani Atmos. Chem. Phys., 16, 8539-8557, doi:10.5194/acp-16-8539-2016, 2016 We present an assessment of ultraviolet nadir-backscatter and infrared limb-emission ozone data assimilation. Understanding the differences in the impact produced by the assimilation of limb and nadir ozone data is relevant to many applications, including climate reanalysis and air quality. The results show the potential and limitations of each dataset and support the need for a more balanced long-term availability of both types of sensors than currently envisaged.

Description: Continuous measurements of isotopic composition of water vapour on the East Antarctic Plateau Mathieu Casado, Amaelle Landais, Valérie Masson-Delmotte, Christophe Genthon, Erik Kerstel, Samir Kassi, Laurent Arnaud, Ghislain Picard, Frederic Prie, Olivier Cattani, Hans-Christian Steen-Larsen, Etienne Vignon, and Peter Cermak Atmos. Chem. Phys., 16, 8521-8538, doi:10.5194/acp-16-8521-2016, 2016 Climatic conditions in Concordia are very cold (−55 °C in average) and very dry, imposing difficult conditions to measure the water vapour isotopic composition. New developments in infrared spectroscopy enable now the measurement of isotopic composition in water vapour traces (down to 20 ppmv). Here we present the results results of a first campaign of measurement of isotopic composition of water vapour in Concordia, the site where the 800 000 years long ice core was drilled.

Description: Surface Renewal as a Significant Mechanism for Dust Emission Jie Zhang, Zhenjiao Teng, Ning Huang, Lei Guo, and Yaping Shao Atmos. Chem. Phys. Discuss., doi:10.5194/acp-2016-421,2016 Manuscript under review for ACP (discussion: open, 0 comments) In spite of the tremendous efforts, many questions remain unanswered regarding dust emission mechanisms. A series of wind-tunnel experiments are carried out on dust emissions from different soil surfaces to better understand relevant mechanisms. Here are some interesting results that demonstrate the importance of surface renewal mechanism which was normally neglected in previous researches and is strongly recommended to be considered in future dust models.

Description: Variations in O 3 , CO, and CH 4 over the Bay of Bengal during the summer monsoon season: Ship-borne measurements and model simulations Imran A. Girach, Narendra Ojha, Prabha R. Nair, Andrea Pozzer, Yogesh K. Tiwari, K. Ravi Kumar, and Jos Lelieveld Atmos. Chem. Phys. Discuss., doi:10.5194/acp-2016-595,2016 Manuscript under review for ACP (discussion: open, 0 comments) This study presents first ship-borne measurements of trace gases over the Bay of Bengal during summer monsoon. The observed variations in trace gases are shown to be due to dynamics/transport and en-route photochemistry. Analysis of meteorological and chemical fields shows that significantly lower ozone during rainfall is associated with the downdrafts. A regional model reproduces the observed variations and revealed the rapid transport of ozone across the Bay of Bengal during an event.

Description: Chemical Characteristics of Marine Fine Aerosols over Sea and at Offshore Islands during Three Cruise Sampling Campaigns in the Taiwan Strait– Sea Salts and Anthropogenic Particles Tsung-Chang Li, Chung-Shin Yuan, Chung-Hsuang Hung, Hsun-Yu Lin, Hu-Ching Huang, and Chon-Lin Lee Atmos. Chem. Phys. Discuss., doi:10.5194/acp-2016-384,2016 Manuscript under review for ACP (discussion: open, 0 comments) This manuscript is significant because it is important to establish data on the background physicochemical characteristics of atmospheric fine particles, as well as their spatial distributions and seasonal variations given various meteorological conditions, and the results of this study will be very much valuable for clarifying the transportation PM 2.5 over sea and at the offshore islands at the Taiwan Strait, which can further fill the the important blank for the global atmospheric research.

Description: Continuous measurements of isotopic composition of water vapour on the East Antarctic Plateau Mathieu Casado, Amaelle Landais, Valérie Masson-Delmotte, Christophe Genthon, Erik Kerstel, Samir Kassi, Laurent Arnaud, Ghislain Picard, Frederic Prie, Olivier Cattani, Hans-Christian Steen-Larsen, Etienne Vignon, and Peter Cermak Atmos. Chem. Phys., 16, 8521-8538, doi:10.5194/acp-16-8521-2016, 2016 Climatic conditions in Concordia are very cold (−55 °C in average) and very dry, imposing difficult conditions to measure the water vapour isotopic composition. New developments in infrared spectroscopy enable now the measurement of isotopic composition in water vapour traces (down to 20 ppmv). Here we present the results results of a first campaign of measurement of isotopic composition of water vapour in Concordia, the site where the 800 000 years long ice core was drilled.

Description: On the role of tropopause folds in summertime tropospheric ozone over the eastern Mediterranean and the Middle East Dimitris Akritidis, Andrea Pozzer, Prodromos Zanis, Evangelos Tyrlis, Bojan Škerlak, Michael Sprenger, and Jos Lelieveld Atmos. Chem. Phys. Discuss., doi:10.5194/acp-2016-547,2016 Manuscript under review for ACP (discussion: open, 0 comments) We investigate the contribution of tropopause folds in the summertime tropospheric ozone pool over the eastern Mediterranean and the Middle East. For this purpose we use the EMAC atmospheric chemistry climate model and a fold identification algorithm. A clear increase of ozone is found in the middle troposphere due to fold activity. The interannual variability of near surface ozone over the eastern Mediterranean is related to that of both tropopause folds and ozone in the free troposphere.

Description: A Polarimetric Scattering Database for Non-spherical Ice Particles at MicrowaveWavelengths Yinghui Lu, Zhiyuan Jiang, Kultegin Aydin, Johannes Verlinde, Eugene E. Clothiaux, and Giovanni Botta Atmos. Meas. Tech. Discuss., doi:10.5194/amt-2016-228,2016 Manuscript under review for AMT (discussion: open, 0 comments) The database contains the complete(polarimatric) scattering information for different types of ice particles at different incident and scattered radiation directions at 4 microwave wavelengths. These results are useful for understanding the dependence of ice-particle scattering properties on ice-particle orientation with respect to the incident and scattered radiation. It is also useful in ice property retrievals, radar forward simulation.

Description: Variations in O 3 , CO, and CH 4 over the Bay of Bengal during the summer monsoon season: Ship-borne measurements and model simulations Imran A. Girach, Narendra Ojha, Prabha R. Nair, Andrea Pozzer, Yogesh K. Tiwari, K. Ravi Kumar, and Jos Lelieveld Atmos. Chem. Phys. Discuss., doi:10.5194/acp-2016-595,2016 Manuscript under review for ACP (discussion: open, 0 comments) This study presents first ship-borne measurements of trace gases over the Bay of Bengal during summer monsoon. The observed variations in trace gases are shown to be due to dynamics/transport and en-route photochemistry. Analysis of meteorological and chemical fields shows that significantly lower ozone during rainfall is associated with the downdrafts. A regional model reproduces the observed variations and revealed the rapid transport of ozone across the Bay of Bengal during an event.

Description: A comparative analysis of UV nadir-backscatter and infrared limb-emission ozone data assimilation Rossana Dragani Atmos. Chem. Phys., 16, 8539-8557, doi:10.5194/acp-16-8539-2016, 2016 We present an assessment of ultraviolet nadir-backscatter and infrared limb-emission ozone data assimilation. Understanding the differences in the impact produced by the assimilation of limb and nadir ozone data is relevant to many applications, including climate reanalysis and air quality. The results show the potential and limitations of each dataset and support the need for a more balanced long-term availability of both types of sensors than currently envisaged.

Description: Why did the storm ex-Gaston (2010) fail to redevelop during the PREDICT experiment? Thomas M. Freismuth, Blake Rutherford, Mark A. Boothe, and Michael T. Montgomery Atmos. Chem. Phys., 16, 8511-8519, doi:10.5194/acp-16-8511-2016, 2016 Numerical model analyses are used to investigate the role of dry, environmental air in the failed redevelopment of a tropical cyclone (ex-Gaston, 2010). As early as 12:00 UTC 2 September 2010, a dry layer at and above 600 hPa results in a decrease in the vertical mass flux and vertical, relative vorticity. The intrusion of dry air led to a reduction in vorticity and a compromised pouch at these middle levels. This study supports work looking at the role of dry air in moist convection.

Description: Thermal infrared laser heterodyne spectro-radiometry for solar occultation atmospheric CO 2 measurements Alex Hoffmann, Neil A. Macleod, Marko Huebner, and Damien Weidmann Atmos. Meas. Tech. Discuss., doi:10.5194/amt-2016-142,2016 Manuscript under review for AMT (discussion: open, 0 comments) This paper focuses on the demonstration and assessment of thermal infrared laser heterodyne spectro-radiometry for the remote sensing of carbon dioxide (CO 2 ). A research instrument has been developed and operated from the ground using direct sunlight to measure CO 2 to a high precision. This technology would enable the development of high-performance miniature ground-based sounders to complement existing measurement networks and contribute to the improvement of global carbon emission assessment.

Description: A Polarimetric Scattering Database for Non-spherical Ice Particles at MicrowaveWavelengths Yinghui Lu, Zhiyuan Jiang, Kultegin Aydin, Johannes Verlinde, Eugene E. Clothiaux, and Giovanni Botta Atmos. Meas. Tech. Discuss., doi:10.5194/amt-2016-228,2016 Manuscript under review for AMT (discussion: open, 0 comments) The database contains the complete(polarimatric) scattering information for different types of ice particles at different incident and scattered radiation directions at 4 microwave wavelengths. These results are useful for understanding the dependence of ice-particle scattering properties on ice-particle orientation with respect to the incident and scattered radiation. It is also useful in ice property retrievals, radar forward simulation.

Description: Chemical Characteristics of Marine Fine Aerosols over Sea and at Offshore Islands during Three Cruise Sampling Campaigns in the Taiwan Strait– Sea Salts and Anthropogenic Particles Tsung-Chang Li, Chung-Shin Yuan, Chung-Hsuang Hung, Hsun-Yu Lin, Hu-Ching Huang, and Chon-Lin Lee Atmos. Chem. Phys. Discuss., doi:10.5194/acp-2016-384,2016 Manuscript under review for ACP (discussion: open, 0 comments) This manuscript is significant because it is important to establish data on the background physicochemical characteristics of atmospheric fine particles, as well as their spatial distributions and seasonal variations given various meteorological conditions, and the results of this study will be very much valuable for clarifying the transportation PM 2.5 over sea and at the offshore islands at the Taiwan Strait, which can further fill the the important blank for the global atmospheric research.

Description: Evaluation of Machine Learning Algorithms for Classification of Primary Biological Aerosol using a new UV-LIF spectrometer Simon Ruske, David O. Topping, Virginia E. Foot, Paul H. Kaye, Warren R. Stanley, Ian Crawford, Andrew P. Morse, and Martin W. Gallagher Atmos. Meas. Tech. Discuss., doi:10.5194/amt-2016-214,2016 Manuscript under review for AMT (discussion: open, 0 comments) Particles such as Bacteria, Pollen and Fungal spores have important implications within the environment and public health sectors. Here we evaluate the performance of various different methods for distinguishing between these different types of particles using a new instrument. We demonstrate that there may be better alternatives to the currently used methods which can be further investigated in future research.

Description: Identification of the cloud base height over the central Himalayan region: Intercomparison of Ceilometer and Doppler Lidar K. K. Shukla, K. Niranjan Kumar, D. V. Phanikumar, R. K. Newsom, V. R. Kotamarthi, T. B. M. J. Ouarda, and M. V. Ratnam Atmos. Meas. Tech. Discuss., doi:10.5194/amt-2016-162,2016 Manuscript under review for AMT (discussion: open, 0 comments) Estimation of Cloud base height was carried out by using various ground based instruments (Doppler Lidar and Ceilometer) and satellite datasets (MODIS) over central Himalayan region for the first time. The present study demonstrates the potential of Doppler Lidar in precise estimation of cloud base height and updraft velocities. More such deployments will be invaluable inputs for regional weather prediction models over complex Himalayan terrains.

Description: On the role of tropopause folds in summertime tropospheric ozone over the eastern Mediterranean and the Middle East Dimitris Akritidis, Andrea Pozzer, Prodromos Zanis, Evangelos Tyrlis, Bojan Škerlak, Michael Sprenger, and Jos Lelieveld Atmos. Chem. Phys. Discuss., doi:10.5194/acp-2016-547,2016 Manuscript under review for ACP (discussion: open, 0 comments) We investigate the contribution of tropopause folds in the summertime tropospheric ozone pool over the eastern Mediterranean and the Middle East. For this purpose we use the EMAC atmospheric chemistry climate model and a fold identification algorithm. A clear increase of ozone is found in the middle troposphere due to fold activity. The interannual variability of near surface ozone over the eastern Mediterranean is related to that of both tropopause folds and ozone in the free troposphere.

Description: Increased nitrate and decreased δ 15 N–NO 3 − in the Greenland Arctic after 1940 attributed to North American oil burning Nathan J. Chellman, Meredith G. Hastings, and Joseph R. McConnell The Cryosphere Discuss., doi:10.5194/tc-2016-163,2016 Manuscript under review for TC (discussion: open, 0 comments) This manuscript analyzes the changing sources of nitrate deposition to Greenland since 1760 CE using a dataset consisting of sub-seasonally resolved nitrogen isotopes of nitrate and source tracers. Correlations amongst ion concentration, source tracers, and the δ 15 N–NO 3 − provide evidence of the impact of biomass burning and fossil fuel combustion emissions of nitrogen oxides and suggest that oil combustion is the likely driver of increased nitrate concentration in Greenland ice since 1940 CE.

Description: Inversion of geothermal heat flux in a thermomechanically coupled nonlinear Stokes ice sheet model Hongyu Zhu, Noemi Petra, Georg Stadler, Tobin Isaac, Thomas J. R. Hughes, and Omar Ghattas The Cryosphere, 10, 1477-1494, doi:10.5194/tc-10-1477-2016, 2016 We study how well the basal geothermal heat flux can be inferred from surface velocity observations using a thermomechanically coupled nonlinear Stokes ice sheet model. The prospects and limitations of this inversion is studied in two and three dimensional model problems. We also argue that a one-way coupled approach for the adjoint equations motivated by staggered solvers for forward multiphysics problems can lead to an incorrect gradient and premature termination of the optimization iteration.

Description: Increased nitrate and decreased δ 15 N–NO 3 − in the Greenland Arctic after 1940 attributed to North American oil burning Nathan J. Chellman, Meredith G. Hastings, and Joseph R. McConnell The Cryosphere Discuss., doi:10.5194/tc-2016-163,2016 Manuscript under review for TC (discussion: open, 0 comments) This manuscript analyzes the changing sources of nitrate deposition to Greenland since 1760 CE using a dataset consisting of sub-seasonally resolved nitrogen isotopes of nitrate and source tracers. Correlations amongst ion concentration, source tracers, and the δ 15 N–NO 3 − provide evidence of the impact of biomass burning and fossil fuel combustion emissions of nitrogen oxides and suggest that oil combustion is the likely driver of increased nitrate concentration in Greenland ice since 1940 CE.

Description: Improving the inter-hemispheric gradient of total column atmospheric CO 2 and CH 4 in simulations with the ECMWF semi-Lagrangian atmospheric global model Anna Agusti-Panareda, Michail Diamantakis, Victor Bayona, Friedrich Klappenbach, and Andre Butz Geosci. Model Dev. Discuss., doi:10.5194/gmd-2016-143,2016 Manuscript under review for GMD (discussion: open, 0 comments) This paper demonstrates how important mass fixers can be in the simulation of long-lived greenhouse gases using transport models based on the highly efficient semi-lagrangian advection scheme. Mass fixers can have a large impact on the representation of the inter-hemispheric gradient of CO 2 and CH 4 , a crucial feature of their distribution. This work is relevant for models simulating atmospheric composition that use semi-lagrangian advection schemes both for climate and air quality applications.

Description: sUAS and their application in observing geomorphological processes Jozef Gallik and Lenka Bolešová Solid Earth, 7, 1033-1042, doi:10.5194/se-7-1033-2016, 2016 Technology is moving ahead very fast, and so researchers have new possibilities for their research. We tried to demonstrate benefits of using remote-sensing technology (Phantom 1 drone) such as its accuracy in the terrain, easy access to hardly accessible areas, and the possibility to collect data even during unfavourable weather conditions. The high mountainous environment provided us great conditions for testing the drone as a device for very easy and accurate mapping of natural phenomena.

Description: Response of seasonal soil freeze depth to climate change across China Xiaoqing Peng, Oliver W. Frauenfeld, Tingjun Zhang, Kang Wang, Bin Cao, Xinyue Zhong, Hang Su, and Cuicui Mu The Cryosphere Discuss., doi:10.5194/tc-2016-129,2016 Manuscript under review for TC (discussion: open, 0 comments) Previous researches paid more attention on the permafrost, e.g. active layer thickness, soil temperature, permafrost area extent, and associated with permafrost degradation leading other changes. However, seasonally frozen ground, vast area extent, did not focus by so much attention. Here, we combined more than 800 observation station data and gridded data to investigate soil freeze depth across China. The results indicate that soil freeze depth increase with climate warming.

Description: Slight glacier reduction over the northwestern Tibetan Plateau despite significant recent warming Yetang Wang, Shugui Hou, Wenling An, Hongxi Pang, and Yaping Liu The Cryosphere Discuss., doi:10.5194/tc-2016-165,2016 Manuscript under review for TC (discussion: open, 0 comments) This study further confirms "Pamir–Karakoram–Western-Kunlun-Mountain (northwestern Tibetan Plateau) Glacier Anomaly". Slight glacier reduction over the northwestern Tibetan Plateau may result from more accumulation from increased precipitation in winter which to great extent protects it from mass reductions under climate warming during 1961–2000. Warming slowdown since 2000 happening at this region may further mitigate glacier mass reduction.

Description: Response of seasonal soil freeze depth to climate change across China Xiaoqing Peng, Oliver W. Frauenfeld, Tingjun Zhang, Kang Wang, Bin Cao, Xinyue Zhong, Hang Su, and Cuicui Mu The Cryosphere Discuss., doi:10.5194/tc-2016-129,2016 Manuscript under review for TC (discussion: open, 0 comments) Previous researches paid more attention on the permafrost, e.g. active layer thickness, soil temperature, permafrost area extent, and associated with permafrost degradation leading other changes. However, seasonally frozen ground, vast area extent, did not focus by so much attention. Here, we combined more than 800 observation station data and gridded data to investigate soil freeze depth across China. The results indicate that soil freeze depth increase with climate warming.

Description: Effect of tropical cyclones on the stratosphere–troposphere exchange observed using satellite observations over the north Indian Ocean M. Venkat Ratnam, S. Ravindra Babu, S. S. Das, G. Basha, B. V. Krishnamurthy, and B. Venkateswararao Atmos. Chem. Phys., 16, 8581-8591, doi:10.5194/acp-16-8581-2016, 2016 The impact of cyclones that occurred over the north Indian Ocean during 2007–2013 on the STE process is quantified using satellite observations. It is shown that cyclones have a significant impact on the tropopause structure, ozone and water vapour budget, and consequentially STE in the UTLS region. The cross-tropopause mass flux from the stratosphere to the troposphere for cyclonic storms is found to be 0.05 ± 0.29 × 10 −3  kg m −2 , and for very severe cyclonic storms it is 0.5 ± 1.07 × 10 −3  kg m −2 .

Description: Variation of CCN activity during new particle formation events in the North China Plain Nan Ma, Chunsheng Zhao, Jiangchuan Tao, Zhijun Wu, Simonas Kecorius, Zhibin Wang, Johannes Größ, Hongjian Liu, Yuxuan Bian, Ye Kuang, Monique Teich, Gerald Spindler, Konrad Müller, Dominik van Pinxteren, Hartmut Herrmann, Min Hu, and Alfred Wiedensohler Atmos. Chem. Phys., 16, 8593-8607, doi:10.5194/acp-16-8593-2016, 2016 New particle formation (NPF) is one of main sources of cloud condensation nuclei (CCN) in the atmosphere. Based on in situ measurements, we found that CCN activity of newly formed particles largely differs in different NPF events. It is therefore difficult to find a simple parameterization of CCN activity for NPF events. Using a fixed size-resolved activation ratio curve or critical diameter is very likely to result in large biases up to 50 % in the calculated N CCN during NPF events.

Description: Real time retrieval of volcanic cloud particles and SO 2 by satellite using an improved simplified approach Sergio Pugnaghi, Lorenzo Guerrieri, Stefano Corradini, and Luca Merucci Atmos. Meas. Tech., 9, 3053-3062, doi:10.5194/amt-9-3053-2016, 2016 Volcanic plume removal (VPR) is a procedure developed to retrieve the ash optical depth, effective radius and mass, and sulfur dioxide mass contained in a volcanic cloud from the thermal radiance at 8.7, 11, and 12 µm. It is based on an estimation of a virtual image representing what the sensor would have seen in a multispectral thermal image if the volcanic cloud were not present. Ash and sulfur dioxide were retrieved by the first version of the VPR using a very simple atmospheric model that ignored the layer above the volcanic cloud. This new version takes into account the layer of atmosphere above the cloud as well as thermal radiance scattering along the line of sight of the sensor. In addition to improved results, the new version also offers an easier and faster preliminary preparation and includes other types of volcanic particles (andesite, obsidian, pumice, ice crystals, and water droplets). As in the previous version, a set of parameters regarding the volcanic area, particle types, and sensor is required to run the procedure. However, in the new version, only the mean plume temperature is required as input data. In this work, a set of parameters to compute the volcanic cloud transmittance in the three quoted bands, for all the aforementioned particles, for both Mt. Etna (Italy) and Eyjafjallajökull (Iceland) volcanoes, and for the Terra and Aqua MODIS instruments is presented. Three types of tests are carried out to verify the results of the improved VPR. The first uses all the radiative transfer simulations performed to estimate the above mentioned parameters. The second one makes use of two synthetic images, one for Mt. Etna and one for Eyjafjallajökull volcanoes. The third one compares VPR and Look-Up Table (LUT) retrievals analyzing the true image of Eyjafjallajökull volcano acquired by MODIS aboard the Aqua satellite on 11 May 2010 at 14:05 GMT.

Description: Temperature-dependence of aerosol optical depth over the southeastern US Tero Mielonen, Anca Hienola, Thomas Kühn, Joonas Merikanto, Antti Lipponen, Tommi Bergman, Hannele Korhonen, Pekka Kolmonen, Larisa Sogacheva, Darren Ghent, Antti Arola, Gerrit de Leeuw, and Harri Kokkola Atmos. Chem. Phys. Discuss., doi:10.5194/acp-2016-625,2016 Manuscript under review for ACP (discussion: open, 0 comments) We studied the temperature dependence of AOD and its radiative effects over the southeastern US. We used spaceborne observations of AOD, LST and tropospheric NO 2 with simulations of ECHAM-HAMMOZ. The level of AOD in this region is governed by anthropogenic emissions but the temperature dependency is most likely caused by BVOC emissions. According to the observations and simulations, the regional clear-sky DRE for biogenic aerosols is −0.43 ± 0.88 W/m 2 /K and −0.86 ± 0.06 W/m 2 /K, respectively.

Description: Stratospheric Air Sub-sampler (SAS) and its application to analysis of Δ 17 O(CO 2 ) from small air samples collected with an AirCore Dorota Janina Mrozek, Carina van der Veen, Magdalena E. G. Hofmann, Huilin Chen, Rigel Kivi, Pauli Heikkinen, and Thomas Röckmann Atmos. Meas. Tech. Discuss., doi:10.5194/amt-2016-124,2016 Manuscript under review for AMT (discussion: open, 0 comments) Stratospheric Air Sub-sampler (SAS) is a device to collect and to store the stratospheric profile of air collected with an AirCore (Karion et al, 2010) in numerous sub-samples. The sub-samples (each of 25 mL at ambient temperature and pressure) can be later introduced to the continuous flow systems to measure for example the isotopic composition of CO 2 . The performance of the coupled system is demonstrated for a set of air samples from an AirCore flight in November 2014 near Sodankylä, Finland.

Description: Real time retrieval of volcanic cloud particles and SO 2 by satellite using an improved simplified approach Sergio Pugnaghi, Lorenzo Guerrieri, Stefano Corradini, and Luca Merucci Atmos. Meas. Tech., 9, 3053-3062, doi:10.5194/amt-9-3053-2016, 2016 Volcanic plume removal (VPR) is a procedure developed to retrieve the ash optical depth, effective radius and mass, and sulfur dioxide mass contained in a volcanic cloud from the thermal radiance at 8.7, 11, and 12 µm. It is based on an estimation of a virtual image representing what the sensor would have seen in a multispectral thermal image if the volcanic cloud were not present. Ash and sulfur dioxide were retrieved by the first version of the VPR using a very simple atmospheric model that ignored the layer above the volcanic cloud. This new version takes into account the layer of atmosphere above the cloud as well as thermal radiance scattering along the line of sight of the sensor. In addition to improved results, the new version also offers an easier and faster preliminary preparation and includes other types of volcanic particles (andesite, obsidian, pumice, ice crystals, and water droplets). As in the previous version, a set of parameters regarding the volcanic area, particle types, and sensor is required to run the procedure. However, in the new version, only the mean plume temperature is required as input data. In this work, a set of parameters to compute the volcanic cloud transmittance in the three quoted bands, for all the aforementioned particles, for both Mt. Etna (Italy) and Eyjafjallajökull (Iceland) volcanoes, and for the Terra and Aqua MODIS instruments is presented. Three types of tests are carried out to verify the results of the improved VPR. The first uses all the radiative transfer simulations performed to estimate the above mentioned parameters. The second one makes use of two synthetic images, one for Mt. Etna and one for Eyjafjallajökull volcanoes. The third one compares VPR and Look-Up Table (LUT) retrievals analyzing the true image of Eyjafjallajökull volcano acquired by MODIS aboard the Aqua satellite on 11 May 2010 at 14:05 GMT.

Description: Evaluation of IWV from the numerical weather prediction WRF model with PPP GNSS processing for Bulgaria Tzvetan Simeonov, Dmitry Sidorov, Felix Norman Teferle, Georgi Milev, and Guergana Guerova Atmos. Meas. Tech. Discuss., doi:10.5194/amt-2016-152,2016 Manuscript under review for AMT (discussion: open, 0 comments) Global Navigation Satellite Systems (GNSS) meteorology is an established operational service providing hourly updated GNSS tropospheric products to the National Meteorologic Services (NMS) in Europe. In the last decade through the ground-based GNSS network densification and new processing strategies like Precise Point Positioning (PPP) it has become possible to obtain sub-hourly tropospheric products for monitoring severe weather events. In this work one year (January–December 2013) of sub-hourly GNSS tropospheric products (Zenith Total Delay) are computed using the PPP strategy for seven stations in Bulgaria. In order to take advantage of the sub-hourly GNSS data to derive Integrated Water Vapour (IWV) surface pressure and temperature with similar temporal resolution is required. As the surface observations are on 3 hourly basis the first step is to compare the surface pressure and temperature from numerical weather prediction model Weather Forecasting and Research (WRF) with observations at three synoptic stations in Bulgaria. The mean difference between the two data-sets for 1) surface pressure is less than 0.5 hPa and the correlation is over 0.989, 2) temperature the largest mean difference is 1.1 °C and the correlation coefficient is over 0.957 and 3) IWV mean difference is in range 0.1–1.1 mm. The evaluation of WRF on annual bases shows IWV underestimation between 0.5 and 1.5 mm at five stations and overestimation at Varna and Rozhen. Varna and Rozhen have also much smaller correlation 0.9 and 0.76. The study of the monthly IWV variation shows that at those locations the GNSS IWV has unexpected drop in April and March respectively. The reason for this drop is likely problems with station raw data. At the remaining 5 stations a very good agreement between GNSS and WRF is observed with high correlation during the cold part of 2013 i.e. March, October and December (0.95) and low correlation during the warm part of 2013 i.e. April to August (below 0.9). The diurnal cycle of the WRF model shows a dry bias in the range of 0.5-1.5 mm. Between 00 and 01 UTC the GNSS IWV tends to be underestimate IWV which is likely due to the processing window used. The precipitation efficiency from GNSS and WRF show very good agreement on monthly bases with a maximum in May-June and minimum in August–September. The annual precipitation efficiency in 2013 at Lovech and Burgas is about 6 %.

Description: Contrasting radiation and soil heat fluxes in Arctic shrub and wet sedge tundra Inge Juszak, Werner Eugster, Monique M. P. D. Heijmans, and Gabriela Schaepman-Strub Biogeosciences, 13, 4049-4064, doi:10.1594/PANGAEA.860561, 2016 Changes in Arctic vegetation composition and structure feed back to climate and permafrost. Using field observations at a Siberian tundra site, we find that dwarf shrubs absorb more solar radiation than wet sedges and thus amplify surface warming, especially during snow melt. On the other hand, permafrost thaw was enhanced below sedges as a consequence of high soil moisture. Standing dead sedge leaves affected the radiation budget strongly and deserve more scientific attention.

Description: Holocene hydrological changes in the Rhône River (NW Mediterranean) as recorded in the marine mud belt Maria-Angela Bassetti, Serge Berné, Marie-Alexandrine Sicre, Bernard Dennielou, Yoann Alonso, Roselyne Buscail, Bassem Jalali, Bertil Hebert, and Christophe Menniti Clim. Past, 12, 1539-1553, doi:10.5194/cp-12-1539-2016, 2016 This work represents the first attempt to decipher the linkages between rapid climate changes and continental Holocene paleohydrology in the NW Mediterranean shallow marine setting. Between 11 and 4 ka cal BP, terrigenous input increased and reached a maximum at 7 ka cal BP, probably as a result of a humid phase. From ca. 4 ka cal BP to the present, enhanced variability in the land-derived material is possibly due to large-scale atmospheric circulation and rainfall patterns in western Europe.

Description: Effects of Multiple Doppler Radar data assimilation on the numerical simulation of a Flash Flood Event during the HyMeX campaign Ida Maiello, Sabrina Gentile, Rossella Ferretti, Luca Baldini, Nicoletta Roberto, Errico Picciotti, Pier Paolo Alberoni, and Frank S. Marzano Hydrol. Earth Syst. Sci. Discuss., doi:10.5194/hess-2016-320,2016 Manuscript under review for HESS (discussion: open, 0 comments) An analysis to evaluate the impact of assimilating multiple radar data with a three dimensional variational (3D-Var) system on a heavy precipitation event is presented. The main goal is to establish a general methodology to quantitatively assess the performance of flash-flood numerical weather prediction at mesoscale. In this respect, during the first Special Observation Period (SOP1) of HyMeX (Hydrological cycle in the Mediterranean Experiment) campaign several Intensive Observing Periods (IOPs) were launched and nine occurred in Italy. Among them IOP4 is chosen for this study because of its low predictability. This event hit central Italy on 14 September 2012 producing heavy precipitation and causing several damages. Data taken from three C-band radars running operationally during the event are assimilated to improve high resolution initial conditions. In order to evaluate the impact of the assimilation procedure at different horizontal resolution and to assess the impact of assimilating multiple radars data, several experiments using Weather Research and Forecasting (WRF) model are performed. Finally, the statistical indexes as accuracy, equitable threat score, false alarm ratio and frequency bias are used to objectively compare the experiments, using rain gauges data as benchmark.

Description: Model-based study of the role of rainfall and land use–land cover in the changes in the occurrence and intensity of Niger red floods in Niamey between 1953 and 2012 Claire Casse, Marielle Gosset, Théo Vischel, Guillaume Quantin, and Bachir Alkali Tanimoun Hydrol. Earth Syst. Sci., 20, 2841-2859, doi:10.5194/hess-20-2841-2016, 2016 Since 1950, the Niger River basin has overcome drastic changes. In Niamey city, the highest river levels and the longest flooded period ever recorded occurred in 2003, 2010, 2012 and 2013, with heavy casualties and property damage. The reasons for these changes, and the relative role of climate versus Land Use Land Cover changes are still debated and are investigated in this paper using observations and modelling. We conclude on the successive role of cover and then rainfall variability.

Description: Upscaling instantaneous to daily evapotranspiration using modelled daily shortwave radiation for remote sensing applications: an Artificial Neural Network approach Loise Wandera, Kaninska Mallick, Gerard Kiely, Olivier Roupsard, Matthias Peichl, and Vicenzo Magliulo Hydrol. Earth Syst. Sci. Discuss., doi:10.5194/hess-2016-344,2016 Manuscript under review for HESS (discussion: open, 0 comments) Upscaling instantaneous to daily evapotranspiration (ET i to ET d ) is one of the central challenge in regional vegetation water use mapping using polar orbiting satellites. Here we developed a robust ET i upscaling for global studies using the ratio between daily and instantaneous global radiation (R Sd /R Si ). Using data from 126 FLUXNET tower sites this study demonstrated R Sd /R si ratio to be the most robust factor explaining ET d /ET i variability across variable sky conditions and multiple biomes.

Description: 4-D imaging of sub-second dynamics in pore-scale processes using real-time synchrotron X-ray tomography Katherine J. Dobson, Sophia B. Coban, Samuel A. McDonald, Joanna N. Walsh, Robert C. Atwood, and Philip J. Withers Solid Earth, 7, 1059-1073, doi:10.5194/se-7-1059-2016, 2016 State-of-the-art synchrotron x-ray imaging was used to observe micron scale transport processes in real time. The 20 Hz 3-D image acquisition rates give experimental data free from motion artefacts, and suitable for detailed quantitative analysis of the dynamic fluid distribution, flow pathways and processes. The method marks a major breakthrough in our ability to capture both sub-second and lower frequency non-equilibrium process in many geological or engineering systems.

Description: Alpine tectonic wedging and crustal delamination in the Cantabrian Mountains (NW Spain) Jorge Gallastegui, Javier A. Pulgar, and Josep Gallart Solid Earth, 7, 1043-1057, doi:10.5194/se-7-1043-2016, 2016 We present the interpretation of deep seismic reflection profile ESCIN-2 across the boundary between the Cantabrian Mountains and the Cenozoic Duero Basin (NW Spain). 2-D seismic modelling allowed us to construct a N-S geological cross section which is compatible with the rest of the geophysical data. This profile shows the importance of the N-S Alpine deformation. A wedge of Cantabrian margin crust indented southwards into the delaminated Iberian crust, forcing its northward subduction.

Description: Recent advances in smart devices have sustained them as a better alternative for the design of human–machine interaction (HMI), because they are equipped with accelerometer sensor, gyroscope sensor, and an advanced operating system. This paper presents a continuous hand gestures recognition technique that is capable of continuous recognition of hand gestures using three-axis accelerometer and gyroscope sensors in a smart device. To reduce the influence of unstableness of a hand making the gesture and compress the data, a gesture coding algorithm is developed. An automatic gesture spotting algorithm is developed to detect the start and end points of meaningful gesture segments. Finally, a gesture is recognized by comparing the gesture code with gesture database using dynamic time warping algorithm. In addition, a prototype system is developed to recognize the continuous hand gestures-based HMI. With the smartphone, the user is able to perform the predefined gestures and control smart appliances using the Samsung AllShare protocol.

Description: As nanometer technology advances, conventional optical proximity correction (OPC) that minimizes the edge placement error (EPE) at the nominal process condition alone often leads to poor process windows. To improve the mask printability across various process corners, process-window OPC optimizes EPE for multiple process corners, but often suffers long runtime, due to repeated lithographic simulations. This paper presents an efficient process variation (PV)-aware mask optimization framework, namely PVOPC, to simultaneously minimize EPE and PV band with fast convergence. The PVOPC framework includes EPE-sensitivity-driven dynamic fragmentation, PV-aware EPE modeling, and correction with three new EPE-converging techniques and a systematic subresolution-assisted feature insertion algorithm. Experimental results show that our approach efficiently achieves high-quality EPE and PV band results.

Description: Advertisement, IEEE.

Description: Prospective authors are requested to submit new, unpublished manuscripts for inclusion in the upcoming event described in this call for papers.

Description: Binning after volume production is a widely accepted technique to classify fabricated integrated circuits (ICs) into different clusters depending on different degrees of specification compliance. This allows the manufacturer to sell nonoptimal devices at lower rates, so adapting to customer’s quality-price requirements. The binning procedure can be carried out by measuring every single circuit performances, but this approach is costly and time-consuming. On the contrary, if alternate measurements are used to characterize the bins, the procedure is considerably enhanced. In such a case, the specification bin boundaries become arbitrary shape regions due to the highly nonlinear mappings between the specifications space and the alternate measurements space. The binning strategy proposed in this paper functions with the same efficiency regardless of these shapes. The digital encoding of the bins in the alternate measurements space using octrees is the key idea of the proposal. The strategy has two phases: 1) the training phase and 2) the binning phase. In the training phase, the specification bins are encoded using octrees. This first phase requires sufficient samples of each class to generate the octree under realistic variations, but it only needs to be performed once. The binning phase corresponds to the actual production binning of the fabricated ICs. This is achieved by evaluating the alternate measurements in the previously generated octree. The binning phase is fast due to the inherent sparsity of the octree data structure. In order to illustrate the proposal, the method has been applied to a band-pass Butterworth filter considering three specification bins as a proof of concept. Successful simulation results are reported showing considerable advantages as compared to a support vector machine (SVM)-based classifier. Similar bin misclassifications are obtained with both methods, 1.68% using octrees and 1.83% using SVM, while binning time is $5times $ times faster using octrees than using the SVM-based classifier.

Description: Presents information on forthcoming Engineering Management society conferences.

Description: Provides a listing of the editors, board members, and current staff for this issue of the publication.

Description: The small-signal vibration model of Tb x Dy 1– x Fe y plate was developed by mechanical-electrical analogy method. Based on the model, we demonstrated the reasonableness of measuring the small-signal magnetomechanical behaviors by a laser Doppler vibrometer. The strain coefficient of the Tb x Dy 1– x Fe y plates was measured as a function of frequency and bias field ( $H_{mathrm{ dc}})$ , and Young’s modulus, mechanical quality factor ( $Q_{m})$ , relative permeability, and magnetomechanical coupling coefficient were investigated as a function of $H_{mathrm{ dc}}$ . Many novel characteristics were observed under the drive of a small-signal field (7.96 A/m). The change tendency of the strain coefficient at resonance differs from that at low frequency, and the resonant strain coefficients are significantly high (>85 nm/A) in a wide range of bias field from 4.78 to 55.72 kA/m. Both the negative- $Delta E$ and positive- $Delta E$ effects are observed, and the negative- $Delta E$ effect in the low field range is also remarkable. In particular, $Q_{m}$ sharply decreases from the initial value of 104 to a minimum value of 11.4 and, then, increases slowly, and the ratio of the maximum variation of $Q_{m}$ over $H_{mathrm{ dc}}$ to the minimum value of $Q_{m}$ exceeds $sim 812.3$ %. This is an important systematic investigation on the small-signal dynamic magnetomechanical behavior of Tb x Dy 1– x Fe y , and the results are highly beneficial to the designing of magnetostrictive devices.

Description: The presence of particles, which can intrude into the air bearing, is one of the most common factors in the failure of hard disk drives (HDDs). Previous studies have investigated particles trajectories with the assumption of ideal trapping or reflecting boundary conditions in air-filled drives. However, only the colliding particle with insufficient energy to escape the potential well will be trapped by the surface. In this paper, considering the particle-surface energy during the collision, the trapping criterion of the incident normal critical velocity ( $V_{text {ni}}^{mathrm {ast }})$ for Al 2 O 3 particles is developed as the boundary conditions for different colliding surfaces inside a 2.5 in drive. Then, trapping status for Al 2 O 3 particles and particles trajectories inside the drive are simulated by using the commercial computational fluid dynamics solver FLUENT with user-defined functions. The results reveal that the particles will travel longer distances until trapped by HDD components when considering the trapping criterion. In addition, smaller particles will more likely degrade the head–disk interface reliability, since they easily stick on the disk surface.

Description: Benefiting from its simple switching scheme (only a bidirectional current source), high-speed and low-power spin-transfer torque (STT) has been regarded as one of the most promising switching mechanisms for a magnetic tunnel junction (MTJ)-based non-volatile memory and logic circuits. However, it suffers from a number of reliability issues like write error induced by its intrinsic stochasticity, process variation, and so on. In order to reduce the write error rate, the mainstream solution is to enlarge the write pulse duration at the expense of write energy dissipation. Some self-terminated write circuits have been proposed to avoid the wasted write energy. But the hardware cost of these write circuits is especially large. In this paper, a novel cost-efficient self-terminated write circuit is proposed using two simple built-in sensing circuits. The proposed write circuit is simulated with a physics-based STT-MTJ compact model and a commercial CMOS 40 nm design kit. The simulation result shows about 35% reduction of circuit area and 10% lower energy consumption in comparison with that in prior work. In addition, the Error-Free write operation under process variation of both the CMOS transistor and the STT-MTJ is achieved due to its large sense margin ( $sim 320$ mV).

Description: Provides a listing of board members, committee members, editors, and society officers.

Description: Domain wall nanomagnet (DWNM)-based devices have been extensively studied as a promising alternative to the conventional CMOS technology in both the memory and logic implementations due to their non-volatility, near-zero standby power, and high integration density characteristics. In this paper, we leverage a physics-based model of a DWNM device to design a highly scalable current-mode majority gate to achieve a novel one bit full-adder (FA) circuit. The modeled DWNM specifications are calibrated with the experimentally measured data. The functionality of the proposed DWNM-based FA (DWNM-FA) is verified using a SPICE circuit simulator. The detailed analysis and the calculations have been performed to realize the proposed DWNM-FA delay and power consumption corresponding to the various induced input currents at different operating temperatures. The power-delay product of DWNM-FA is examined to tune the operation within the optimum induced input current region to obtain desired power-delay requirements over a range of 200 $mu text{A}$ to 1 mA at temperatures from 298 to 378 K. Finally, the comparison results exhibit 52% and 49% area improvement as well as 41% and 31% improvement in device count complexity over CMOS-based and magnetic tunnel junction-based FA designs, respectively.

Description: Bulk CoCr 2 O 4 undergoes a transition from paramagnetic to long-range ferrimagnetic phase at $T_{c}$ (94 K) to a long-range and/or short-range spiral order at $T_{s}$ ( $sim 24$ K), and finally shows a lock-in-transition below 15 K. The spiral component induces an electric polarization and also a spontaneous magnetization for which it is said to be multiferroic. Reducing the size of a CoCr 2 O 4 multiferroic material to $sim 50$ nm by a coprecipitation method, we obtain a pure cubic phase with space group, Fd3m and lattice parameter (8.334 ± 0.003 °A). A rich sequence of magnetic transitions are examined by measuring temperature and field-dependent magnetization and diffused neutron scattering (DNS) using polarized neutron at different temperatures. While paramagnetic to ferrimagnetic transition is enhanced from 97 K in bulk to 99 K at 0.5 kOe field, followed by a decrease in lock-in-transition ( $T_{L}$ ) from 15 K in bulk to 8 K, spiral ordering temperature does not show a significant change. A strong disagreement between paramagnetic moment obtained from the fitting of $chi ^{-1}=({T}/{C})+({1}/{chi _{o}})-({b}/{T-theta })$ and ferrimagnetic moment obtained from the $M$ versus $H$ loop taken at 2 K, nonsaturated magnetization at 50–100 kOe field, two order of magnitude higher coercivity ( $H_{c}$ ), and splitting of ac susceptibly confirm the core–shell structure of the particles. Furthermore, a magnetic scattering analysis clearly shows that while the paramagnetic to ferrimagnetic transition is continuous, the spiral ordering is sharp, short range, and commensurate in contrast to incommensurate spiral order observed single crystal of CoCr 2 O 4 .

Description: The finite-element analysis for the simulation of magnetic fields in electrical machines leads to an index-1 differential-algebraic equation (DAE) (as opposed to a conventional ordinary differential equation), because the electrical conductivity can be zero in certain regions. First, we construct a DAE-compatible time integration scheme which is energy-balanced, meaning that in a linear system, the input stored and lost powers sum exactly to zero. Second, we use a method based on the energy balance to compute torque. We show that the energy balance method approaches the virtual work principle applied at remeshing layer, as the time step is refined. A similar result also holds if the rotation of the rotor is implemented by Nitsche’s method, which is an instance of the so-called mortar methods.

Description: Investigations about induction sensors, electromagnetic launchers, shields, transformers, and power line-induced currents address increased number of low-frequency research and industrial applications. In general, a magneto-quasi-static (MQS) approximation is considered for the solutions of low-frequency problems in electromagnetics. This approximation leads to a diffusion process when displacement currents are neglected. However, keeping the displacement currents, Maxwell’s equations are valid at low frequencies. In this manner, the finite difference time domain (FDTD) method must be modified by the slowing down propagation velocity at low-frequency regime. In this paper, important and crucial points of the MQS approximation and its application in the FDTD method are clarified in the sense of analytical and numerical aspects. A material scaling technique of dielectric permittivity for the QS FDTD application is analyzed within comprehensive investigations. Furthermore, a criterion for choosing a proper value of scaling parameter will be revealed. Finally, effects of proper and improper values of the scaling parameter are presented with validated analytical and numerical results.

Description: The ability to manipulate the relative magnetization alignment between ferromagnetic source and drain electrodes attached to a molecule or small quantum dot is a prerequisite for a number of spintronic device applications. The influence of electrode shape and field orientation on pairwise magnetization reversal mechanisms in nanogap and point-contact structures is investigated here using micromagnetic simulations. A favorable device geometry and setup are identified for enabling planar, monodomain source and drain electrodes with a magnetization alignment that may be controllably switched between a parallel and antiparallel configuration.

Description: In typical inductive power transfer (IPT) applications, soft magnetic materials are employed in both transmitter and receiver coils as low reluctance paths in order to improve the inductive coupling and support the guidance of magnetic flux. Considering their critical role in the overall performance, this paper investigates the impact of magnetic material properties on the basic features of inductive powering systems. To this end, a group of Ni–Zn and Ni–Cu–Zn spinel ferrite samples with a wide value range of magnetic properties is prepared. The materials are characterized in terms of their structural and magnetic properties, and they are used to build various substrate assemblies for the IPT transmitter (Tx) coil. Conducting Al and Cu sheets and biasing permanent magnets are combined with ferrite disks to simulate common Tx coil designs. This allows the identification of the correlation between the material properties, especially power losses, and the measured power transfer efficiency and shielding effectiveness.

Description: This paper presents a method to estimate steady wind and a maneuvering strategy for loitering under a strong, steady wind. The wind estimation uses Global Positioning System velocity only through a novel filter design without airspeed measurement. The wind estimate is then used to guide the aircraft to crab in a direction perpendicular to the wind, thereby avoiding large changes in flight variables if an orbit-type maneuver is attempted. The proposed method is demonstrated through flight tests.

Description: We propose biased estimators to find the direction of arrival of emitters present in the mainlobe of a spinning antenna-based electronic intelligence system. The proposed estimators were constructed by using Bayesian techniques and by performing a linear transformation and an affine transformation on the maximum likelihood estimator. From a Monte Carlo simulation and experimental results, we demonstrate that the proposed estimators outperform the limit set by the popular performance benchmark, the Cramer-Rao lower bound.

Description: Objective: We present the framework for wearable joint rehabilitation assessment following musculoskeletal injury. We propose a multimodal sensing (i.e., contact based and airborne measurement of joint acoustic emission) system for at-home monitoring. Methods: We used three types of microphones—electret, MEMS, and piezoelectric film microphones—to obtain joint sounds in healthy collegiate athletes during unloaded flexion/extension, and we evaluated the robustness of each microphone's measurements via: 1) signal quality and 2) within-day consistency. Results: First, air microphones acquired higher quality signals than contact microphones (signal-to-noise-and-interference ratio of 11.7 and 12.4 dB for electret and MEMS, respectively, versus 8.4 dB for piezoelectric). Furthermore, air microphones measured similar acoustic signatures on the skin and 5 cm off the skin (∼4.5× smaller amplitude). Second, the main acoustic event during repetitive motions occurred at consistent joint angles (intra-class correlation coefficient ICC(1, 1) = 0.94 and ICC(1, k) = 0.99). Additionally, we found that this angular location was similar between right and left legs, with asymmetry observed in only a few individuals. Conclusion: We recommend using air microphones for wearable joint sound sensing; for practical implementation of contact microphones in a wearable device, interface noise must be reduced. Importantly, we show that airborne signals can be measured consistently and that healthy left and right knees often produce a similar pattern in acoustic emissions. Significance: These proposed methods have the potential for enabling knee joint acoustics measurement outside the clinic/lab and permitting long-term monitoring of knee health for patients rehabilitating an acute knee joint injury.

Description: A new method was proposed for chirp signal detection and estimation built on the frame-based fast Fourier transform (FFT). The proposed method uses the peak frequency difference between FFT frames to detect a chirp signal and estimate chirp rate. This approach differs from conventional methods and is easy to implement. It generates more accurate chirp rate estimation especially under a low signal-to-noise ratio. Simulation and experimental data are used to verify the proposed methods.

Description: As telehealth applications emerge, the need for accurate and reliable biosignal quality indices has increased. One typical modality used in remote patient monitoring is the electrocardiogram (ECG), which is inherently susceptible to several different noise sources, including environmental (e.g., powerline interference), experimental (e.g., movement artifacts), and physiological (e.g., muscle and breathing artifacts). Accurate measurement of ECG quality can allow for automated decision support systems to make intelligent decisions about patient conditions. This is particularly true for in-home monitoring applications, where the patient is mobile and the ECG signal can be severely corrupted by movement artifacts. In this paper, we propose an innovative ECG quality index based on the so-called modulation spectral signal representation. The representation quantifies the rate of change of ECG spectral components, which are shown to be different from the rate of change of typical ECG noise sources. The proposed modulation spectral-based quality index, MS-QI, was tested on 1) synthetic ECG signals corrupted by varying levels of noise, 2) single-lead recorded data using the Hexoskin garment during three activity levels (sitting, walking, running), 3) 12-lead recorded data using conventional ECG machines (Computing in Cardiology 2011 dataset), and 4) two-lead ambulatory ECG recorded from arrhythmia patients (MIT-BIH Arrhythmia Database). Experimental results showed the proposed index outperforming two conventional benchmark quality measures, particularly in the scenarios involving recorded data in real-world environments.

Description: Platelet-rich plasma (PRP) is a volume of autologous plasma that has a higher platelet concentration above baseline. It has already been approved as a new therapeutic modality and investigated in clinics, such as bone repair and regeneration, and oral surgery, with low cost-effectiveness ratio. At present, PRP is mostly prepared using a centrifuge. However, this method has several shortcomings, such as long preparation time (30 min), complexity in operation, and contamination of red blood cells (RBCs). In this paper, a new PRP preparation approach was proposed and tested. Ultrasound waves (4.5 MHz) generated from piezoelectric ceramics can establish standing waves inside a syringe filled with the whole blood. Subsequently, RBCs would accumulate at the locations of pressure nodes in response to acoustic radiation force, and the formed clusters would have a high speed of sedimentation. It is found that the PRP prepared by the proposed device can achieve higher platelet concentration and less RBCs contamination than a commercial centrifugal device, but similar growth factor (i.e., PDGF-ββ). In addition, the sedimentation process under centrifugation and sonication was simulated using the Mason–Weaver equation and compared with each other to illustrate the differences between these two technologies and to optimize the design in the future. Altogether, ultrasound method is an effective method of PRP preparation with comparable outcomes as the commercially available centrifugal products.

Description: The quasimaximum likelihood (QML) estimator of polynomial-phase signals (PPSs) is based on the maximization of the short-time Fourier transform and suffers from aliasing when signals are sampled below the Nyquist sampling rate. In this paper, a phase unwrapping procedure has been proposed as an additional step in the QML to estimate parameters of such signals. Statistical study has shown excellent performance of the proposed approach.

Description: A new wireless sensor was designed, fabricated, and applied for in situ monitoring of tensile force at a wound site. The sensor was comprised of a thin strip of magnetoelastic material with its two ends connected to suture threads for securing the sensor across a wound repair site. Since the sensor was remotely interrogated by applying an ac magnetic field and capturing the resulting magnetic field, it did not require direct wire connections to an external device or internal battery for long-term use. Due to its magnetoelastic property, the application of a tensile force changed the magnetic permeability of the sensor, altering the amplitude of the measured magnetic field. This study presents two sensor designs: one for high and one for low-force ranges. A sensor was fabricated by directly adhering the magnetoelastic strip to the suture. This sensor showed good sensitivity at low force, but its response saturated at about 1.5 N. To monitor high tensile force, the magnetoelastic strip was attached to a metal strip for load sharing. The suture thread was attached to the both ends of the metal strip so only a fraction of the applied force was directed to the sensor, allowing it to exhibit good sensitivity even at 44.5 N. The sensor was applied to two ex vivo models: a sutured section of porcine skin and a whitetail deer Achilles tendon. The results demonstrate the potential for in vivo force monitoring at a wound repair site.

Description: This study presents a precise way to detect the third ( $S_{3}$ ) heart sound, which is recognized as an important indication of heart failure, based on nonlinear single decomposition and time–frequency localization. The detection of the $S_{3}$ is obscured due to its significantly low energy and frequency. Even more, the detected $S_{3}$ may be misunderstood as an abnormal second heart sound with a fixed split, which was not addressed in the literature. To detect such $S_{3}$ , the Hilbert vibration decomposition method is applied to decompose the heart sound into a certain number of subcomponents while intactly preserving the phase information. Thus, the time information of all of the decomposed components are unchanged, which further expedites the identification and localization of any module/section of a signal properly. Next, the proposed localization step is applied to the decomposed subcomponents by using smoothed pseudo Wigner–Ville distribution followed by the reassignment method. Finally, based on the positional information, the $S_{3}$ is distinguished and confirmed by measuring time delays between the $S_{2}$ and $S_{3}$ . In total, 82 sets of cardiac cycles collected from different databases including Texas Heart Institute database are examined for evaluation of the proposed method. The result analysis shows that the proposed method can detect the $S_{3}$ correctly, even when the - ormalized temporal energy of $S_{3}$ is larger than 0.16, and the frequency of those is larger than 34 Hz. In a performance analysis, the proposed method demonstrates that the accuracy rate of $S_{3}$ detection is as high as 93.9%, which is significantly higher compared with the other methods. Such findings prove the robustness of the proposed idea for detecting substantially low-energized $S_{3}$ .

Description: Electromagnetic (EM) tracking systems are highly susceptible to field distortion. The interference can cause measurement errors up to a few centimeters in clinical environments, which limits the reliability of these systems. Unless corrected for, this measurement error imperils the success of clinical procedures. It is therefore fundamental to dynamically calibrate EM tracking systems and compensate for measurement error caused by field distorting objects commonly present in clinical environments. We propose to combine a motion model with observations of redundant EM sensors and compensate for field distortions in real time. We employ a simultaneous localization and mapping technique to accurately estimate the pose of the tracked instrument while creating the field distortion map. We conducted experiments with six degrees-of-freedom motions in the presence of field distorting objects in research and clinical environments. We applied our approach to improve the EM tracking accuracy and compared our results to a conventional sensor fusion technique. Using our approach, the maximum tracking error was reduced by 67% for position measurements and by 64% for orientation measurements. Currently, clinical applications of EM trackers are hampered by the adverse distortion effects. Our approach introduces a novel method for dynamic field distortion compensation, independent from preoperative calibrations or external tracking devices, and enables reliable EM navigation for potential applications.