Computed Ultrasound Tomography in Echo mode (CUTE) of speed of sound for diagnosis and for aberration correction in pulse-echo sonography

Sound speed as a diagnostic marker for various diseases of human tissue has been of interest for a while. Up to now, mostly transmission ultrasound computed tomography (UCT) was able to detect spatially resolved sound speed, and its promise as a diagnostic tool has been demonstrated. However, UCT is limited to acoustically transparent samples such as the breast. We present a novel technique where spatially resolved detection of sound speed can be achieved using conventional pulse-echo equipment in reflection mode. For this purpose, pulse-echo images are acquired under various transmit beam directions and a two-dimensional map of the sound speed is reconstructed from the changing phase of local echoes using a direct reconstruction method. Phantom results demonstrate that a high spatial resolution (1 mm) and contrast (0.5 % of average sound speed) can be achieved suitable for diagnostic purposes. In comparison to previous reflection-mode based methods, CUTE works also in a situation with only diffuse echoes, and its direct reconstruction algorithm enables real-time application. This makes it suitable as an addition to conventional clinical ultrasound where it has the potential to benefit diagnosis in a multimodal approach. In addition, knowledge of the spatial distribution of sound speed allows full aberration correction and thus improved spatial resolution and contrast of conventional B-mode ultrasound. © (2014) COPYRIGHT Society of Photo-Optical Instrumentation Engineers (SPIE). Downloading of the abstract is permitted for personal use only.

The geomorphology and morphometry of the banded terrain in Hellas basin, Mars

Hellas basin is a large impact basin situated in the southern highlands of Mars. The north-western part of the basin has the lowest elevation (-7.5 km) on the planet and contains a possibly unique terrain type, which we informally call “banded terrain”. The banded terrain is made up of smooth-looking banded deposits that display signs of viscous behavior and a paucity of superimposed impact craters. In this study, we use newly acquired high spatial resolution images from the High Resolution Imaging Science Experiment (HiRISE) in addition to existing datasets to characterize the geomorphology, the morphometry and the architecture of the banded terrain. The banded terrain is generally confined to the NW edge of the Alpheus Colles plateau. The individual bands are ~3–15 km-long, ~0.3 km-wide and are separated by narrow inter-band depressions, which are ~65 m-wide and ~10 m-deep. The bands display several morphologies that vary from linear to concentric forms. Morphometric analysis reveals that the slopes along a given linear or lobate band ranges from 0.5° to 15° (average~6°), whereas the concentric bands are located on flatter terrain (average slope~2–3°). Crater-size frequency analysis yields an Amazonian-Hesperian boundary crater retention age for the terrain (~3 Gyr), which together, with the presence of very few degraded craters, either implies a recent emplacement, resurfacing, or intense erosion. The apparent sensitivity to local topography and preference for concentrating in localized depressions is compatible with deformation as a viscous fluid. In addition, the bands display clear signs of degradation and slumping at their margins along with a suite of other features that include fractured mounds, polygonal cracks at variable size-scales, and knobby/hummocky textures. Together, these features suggest an ice-rich composition for at least the upper layers of the terrain, which is currently being heavily modified through loss of ice and intense weathering, possibly by wind.

The influence of physical state on shikimic acid ozonolysis: a case for in situ microspectroscopy

Atmospheric soluble organic aerosol material can become solid or semi-solid. Due to increasing viscosity and decreasing diffusivity, this can impact important processes such as gas uptake and reactivity within aerosols containing such substances. This work explores the dependence of shikimic acid ozonolysis on humidity and thereby viscosity. Shikimic acid, a proxy for oxygenated reactive organic material, reacts with O3 in a Criegee-type reaction. We used an environmental microreactor embedded in a scanning transmission X-ray microscope (STXM) to probe this oxidation process. This technique facilitates in situ measurements with single micron-sized particles and allows to obtain near-edge X-ray absorption fine structure (NEXAFS) spectra with high spatial resolution. Thus, the chemical evolution of the interior of the particles can be followed under reaction conditions. The experiments show that the overall degradation rate of shikimic acid is depending on the relative humidity in a way that is controlled by the decreasing diffusivity of ozone with decreasing humidity. This decreasing diffusivity is most likely linked to the increasing viscosity of the shikimic acid–water mixture. The degradation rate was also depending on particle size, most congruent with a reacto-diffusion limited kinetic case where the reaction progresses only in a shallow layer within the bulk. No gradient in the shikimic acid concentration was observed within the bulk material at any humidity indicating that the diffusivity of shikimic acid is still high enough to allow its equilibration throughout the particles on the timescale of hours at higher humidity and that the thickness of the oxidized layer under dry conditions, where the particles are solid, is beyond the resolution of STXM.

Measuring the gravitational free-fall of antihydrogen

Antihydrogen holds the promise to test, for the first time, the universality of freefall with a system composed entirely of antiparticles. The AEgIS experiment at CERN’s antiproton decelerator aims to measure the gravitational interaction between matter and antimatter by measuring the deflection of a beam of antihydrogen in the Earths gravitational field (g). The principle of the experiment is as follows: cold antihydrogen atoms are synthesized in a Penning-Malberg trap and are Stark accelerated towards a moir´e deflectometer, the classical counterpart of an atom interferometer, and annihilate on a position sensitive detector. Crucial to the success of the experiment is the spatial precision of the position sensitive detector.We propose a novel free-fall detector based on a hybrid of two technologies: emulsion detectors, which have an intrinsic spatial resolution of 50 nm but no temporal information, and a silicon strip / scintillating fiber tracker to provide timing and positional information. In 2012 we tested emulsion films in vacuum with antiprotons from CERN’s antiproton decelerator. The annihilation vertices could be observed directly on the emulsion surface using the microscope facility available at the University of Bern. The annihilation vertices were successfully reconstructed with a resolution of 1–2 μmon the impact parameter. If such a precision can be realized in the final detector, Monte Carlo simulations suggest of order 500 antihydrogen annihilations will be sufficient to determine gwith a 1 % accuracy. This paper presents current research towards the development of this technology for use in the AEgIS apparatus and prospects for the realization of the final detector.

Tracking individual membrane proteins and their biochemistry: The power of direct observation

The advent of single molecule fluorescence microscopy has allowed experimental molecular biophysics and biochemistry to transcend traditional ensemble measurements, where the behavior of individual proteins could not be precisely sampled. The recent explosion in popularity of new super-resolution and super-localization techniques coupled with technical advances in optical designs and fast highly sensitive cameras with single photon sensitivity and millisecond time resolution have made it possible to track key motions, reactions, and interactions of individual proteins with high temporal resolution and spatial resolution well beyond the diffraction limit. Within the purview of membrane proteins and ligand gated ion channels (LGICs), these outstanding advances in single molecule microscopy allow for the direct observation of discrete biochemical states and their fluctuation dynamics. Such observations are fundamentally important for understanding molecular-level mechanisms governing these systems. Examples reviewed here include the effects of allostery on the stoichiometry of ligand binding in the presence of fluorescent ligands; the observation of subdomain partitioning of membrane proteins due to microenvironment effects; and the use of single particle tracking experiments to elucidate characteristics of membrane protein diffusion and the direct measurement of thermodynamic properties, which govern the free energy landscape of protein dimerization. The review of such characteristic topics represents a snapshot of efforts to push the boundaries of fluorescence microscopy of membrane proteins to the absolute limit.

Time variable Earth’s gravity field from SLR satellites

The time variable Earth’s gravity field contains information about the mass transport within the system Earth, i.e., the relationship between mass variations in the atmosphere, oceans, land hydrology, and ice sheets. For many years, satellite laser ranging (SLR) observations to geodetic satellites have provided valuable information of the low-degree coefficients of the Earth’s gravity field. Today, the Gravity Recovery and Climate Experiment (GRACE) mission is the major source of information for the time variable field of a high spatial resolution. We recover the low-degree coefficients of the time variable Earth’s gravity field using SLR observations up to nine geodetic satellites: LAGEOS-1, LAGEOS-2, Starlette, Stella, AJISAI, LARES, Larets, BLITS, and Beacon-C. We estimate monthly gravity field coefficients up to degree and order 10/10 for the time span 2003–2013 and we compare the results with the GRACE-derived gravity field coefficients. We show that not only degree-2 gravity field coefficients can be well determined from SLR, but also other coefficients up to degree 10 using the combination of short 1-day arcs for low orbiting satellites and 10-day arcs for LAGEOS-1/2. In this way, LAGEOS-1/2 allow recovering zonal terms, which are associated with long-term satellite orbit perturbations, whereas the tesseral and sectorial terms benefit most from low orbiting satellites, whose orbit modeling deficiencies are minimized due to short 1-day arcs. The amplitudes of the annual signal in the low-degree gravity field coefficients derived from SLR agree with GRACE K-band results at a level of 77 %. This implies that SLR has a great potential to fill the gap between the current GRACE and the future GRACE Follow-On mission for recovering of the seasonal variations and secular trends of the longest wavelengths in gravity field, which are associated with the large-scale mass transport in the system Earth.

A 49 year hindcast of surface winds over the Iberian Peninsula

The use of hindcast climatic data is quite extended for multiple applications. However, this approach needs the support of a validation process to allow its drawbacks and, therefore, confidence levels to be assessed. In this work, the strategy relies on an hourly wind database resulting from a dynamical downscaling experiment, with a spatial resolution of 10 km, covering the Iberian Peninsula (IP), driven by the ERA40 reanalysis (1959–2001) extended by European Centre for Medium-Range Weather Forecast (ECMWF) analysis (2002–2007) and comprising two main steps. Initially, the skill of the simulation is evaluated comparing the quality-tested observational database (Lorente-Plazas et al., 2014) at local and regional scales. The results show that the model is able to portray the main features of the wind over the IP: annual cycles, wind roses, spatial and temporal variability, as well as the response to different circulation types. In addition, there is a significant added value of the simulation with respect to driving conditions, especially in regions with a complex orography. However, some problems are evident, the major drawback being the systematic overestimation of the wind speed, which is mainly attributed to a missrepresentation of frictional forces. The model skill is also lower along the Mediterranean coast and for the Pyrenees. In a second phase, the high spatio-temporal resolution of the pseudo-real wind database is used to explore the limitations of the observational database. It is shown that missing values do not affect the characterisation of the wind climate over the IP, while the length of the observational period (6 years) is sufficient for most regions, with only a few exceptions. The spatial distribution of the observational sampling schemes should be enhanced to improve the correct assessment of all IP wind regimes, particularly in some mountainous areas.

Using MDEFT MRI Sequences to Target the GPi in DBS Surgery.

OBJECTIVE Recent advances in different MRI sequences have enabled direct visualization and targeting of the Globus pallidus internus (GPi) for DBS surgery. Modified Driven Equilibrium Fourier Transform (MDEFT) MRI sequences provide high spatial resolution and an excellent contrast of the basal ganglia with low distortion. In this study, we investigate if MDEFT sequences yield accurate and reliable targeting of the GPi and compare direct targeting based on MDEFT sequences with atlas-based targeting. METHODS 13 consecutive patients considered for bilateral GPi-DBS for dystonia or PD were included in this study. Preoperative targeting of the GPi was performed visually based on MDEFT sequences as well as by using standard atlas coordinates. Postoperative CT imaging was performed to calculate the location of the implanted leads as well as the active electrode(s). The coordinates of both visual and atlas based targets were compared. The stereotactic coordinates of the lead and active electrode(s) were calculated and projected on the segmented GPi. RESULTS On MDEFT sequences the GPi was well demarcated in most patients. Compared to atlas-based planning the mean target coordinates were located significantly more posterior. Subgroup analysis showed a significant difference in the lateral coordinate between dystonia (LAT = 19.33 ± 0.90) and PD patients (LAT = 20.67 ± 1.69). Projected on the segmented preoperative GPi the active contacts of the DBS electrode in both dystonia and PD patients were located in the inferior and posterior part of the structure corresponding to the motor part of the GPi. CONCLUSIONS MDEFT MRI sequences provide high spatial resolution and an excellent contrast enabling precise identification and direct visual targeting of the GPi. Compared to atlas-based targeting, it resulted in a significantly different mean location of our target. Furthermore, we observed a significant variability of the target among the PD and dystonia subpopulation suggesting accurate targeting for each individual patient.

Sensitivity of the WRF model to PBL parametrisations and nesting techniques: evaluation of wind storms over complex terrain

Simulating surface wind over complex terrain is a challenge in regional climate modelling. Therefore, this study aims at identifying a set-up of the Weather Research and Forecasting Model (WRF) model that minimises system- atic errors of surface winds in hindcast simulations. Major factors of the model configuration are tested to find a suitable set-up: the horizontal resolution, the planetary boundary layer (PBL) parameterisation scheme and the way the WRF is nested to the driving data set. Hence, a number of sensitivity simulations at a spatial resolution of 2 km are carried out and compared to observations. Given the importance of wind storms, the analysis is based on case studies of 24 historical wind storms that caused great economic damage in Switzerland. Each of these events is downscaled using eight different model set-ups, but sharing the same driving data set. The results show that the lack of representation of the unresolved topography leads to a general overestimation of wind speed in WRF. However, this bias can be substantially reduced by using a PBL scheme that explicitly considers the effects of non-resolved topography, which also improves the spatial structure of wind speed over Switzerland. The wind direction, although generally well reproduced, is not very sensitive to the PBL scheme. Further sensitivity tests include four types of nesting methods: nesting only at the boundaries of the outermost domain, analysis nudging, spectral nudging, and the so-called re-forecast method, where the simulation is frequently restarted. These simulations show that restricting the freedom of the model to develop large-scale disturbances slightly increases the temporal agreement with the observations, at the same time that it further reduces the overestimation of wind speed, especially for maximum wind peaks. The model performance is also evaluated in the outermost domains, where the resolution is coarser. The results demonstrate the important role of horizontal resolution, where the step from 6 to 2 km significantly improves model performance. In summary, the combination of a grid size of 2 km, the non-local PBL scheme modified to explicitly account for non-resolved orography, as well as analysis or spectral nudging, is a superior combination when dynamical downscaling is aimed at reproducing real wind fields.

Understanding measured water rotational temperatures and column densities in the very innermost coma of Comet 73P/Schwassmann–Wachmann 3 B

Direct sublimation of a comet nucleus surface is usually considered to be the main source of gas in the coma of a comet. However, evidence from a number of comets including the recent spectacular images of Comet 103P/Hartley 2 by the EPOXI mission indicates that the nucleus alone may not be responsible for all, or possibly at times even most, of the total amount of gas seen in the coma. Indeed, the sublimation of icy grains, which have been injected into the coma, appears to constitute an important source. We use the fully-kinetic Direct Simulation Monte Carlo model of Tenishev et al. (Tenishev, V.M., Combi, M.R., Davidsson, B. [2008]. Astrophys. J., 685, 659−677; Tenishev, V.M., Combi, M.R., Rubin, M. [2011]. Astrophys. J., 732) to reproduce the measurements of column density and rotational temperature of water in Comet 73P-B/Schwassmann–Wachmann 3 obtained with a very high spatial resolution of ∼30 km using IRCS/Subaru in May 2006 (Bonev, B.P., Mumma, M.J., Kawakita, H., Kobayashi, H., Villanueva, G.L. [2008]. Icarus, 196, 241−248). For gas released solely from the cometary nucleus at a heliocentric distance of 1 AU, modeled rotational temperatures start at 110 K close to the surface and decrease to only several tens of degrees by 10–20 nucleus radii. However, the measured decay of both rotational temperature and column density with distance from the nucleus is much slower than predicted by this simple model. The addition of a substantial (distributed) source of gas from icy grains in the model slows the decay in rotational temperature and provides a more gradual drop in column density profiles. Together with a contribution of rotational heating of water molecules by electrons, the combined effects allow a much better match to the IRCS/Subaru observations. From the spatial distributions of water abundance and temperature measured in 73P/SW3-B, we have identified and quantified multiple mechanisms of release. The application of this tool to other comets may permit such studies over a range of heliocentric and geocentric distances.

Computed Ultrasound Tomography in Echo Mode for Imaging Speed of Sound Using Pulse-Echo Sonography: Proof of Principle

The limitations of diagnostic echo ultrasound have motivated research into novel modalities that complement ultrasound in a multimodal device. One promising candidate is speed of sound imaging, which has been found to reveal structural changes in diseased tissue. Transmission ultrasound tomography shows speed of sound spatially resolved, but is limited to the acoustically transparent breast. We present a novel method by which speed-of-sound imaging is possible using classic pulse-echo equipment, facilitating new clinical applications and the combination with state-of-the art diagnostic ultrasound. Pulse-echo images are reconstructed while scanning the tissue under various angles using transmit beam steering. Differences in average sound speed along different transmit directions are reflected in the local echo phase, which allows a 2-D reconstruction of the sound speed. In the present proof-of-principle study, we describe a contrast resolution of 0.6% of average sound speed and a spatial resolution of 1 mm (laterally) × 3 mm (axially), suitable for diagnostic applications.

Towards clinical computed ultrasound tomography in echo-mode: Dynamic range artefact reduction

Computed ultrasound tomography in echo-mode (CUTE) allows imaging the speed of sound inside tissue using hand-held pulse-echo ultrasound. This technique is based on measuring the changing local phase of beamformed echoes when changing the transmit beam steering angle. Phantom results have shown a spatial resolution and contrast that could qualify CUTE as a promising novel diagnostic modality in combination with B-mode ultrasound. Unfortunately, the large intensity range of several tens of dB that is encountered in clinical images poses difficulties to echo phase tracking and results in severe artefacts. In this paper we propose a modification to the original technique by which more robust echo tracking can be achieved, and we demonstrate in phantom experiments that dynamic range artefacts are largely eliminated. Dynamic range artefact reduction also allowed for the first time a clinical implementation of CUTE with sufficient contrast to reproducibly distinguish the different speed of sound in different tissue layers of the abdominal wall and the neck.

An automated method for mapping physical soil and water conservation structures on cultivated land using GIS and remote sensing techniques

An efficient and reliable automated model that can map physical Soil and Water Conservation (SWC) structures on cultivated land was developed using very high spatial resolution imagery obtained from Google Earth and ArcGIS, ERDAS IMAGINE, and SDC Morphology Toolbox for MATLAB and statistical techniques. The model was developed using the following procedures: (1) a high-pass spatial filter algorithm was applied to detect linear features, (2) morphological processing was used to remove unwanted linear features, (3) the raster format was vectorized, (4) the vectorized linear features were split per hectare (ha) and each line was then classified according to its compass direction, and (5) the sum of all vector lengths per class of direction per ha was calculated. Finally, the direction class with the greatest length was selected from each ha to predict the physical SWC structures. The model was calibrated and validated on the Ethiopian Highlands. The model correctly mapped 80% of the existing structures. The developed model was then tested at different sites with different topography. The results show that the developed model is feasible for automated mapping of physical SWC structures. Therefore, the model is useful for predicting and mapping physical SWC structures areas across diverse areas.

Comparison of Routine Brain Imaging at 3 T and 7 T.

OBJECTIVE The aim of this study was to compare quantitative and semiquantitative parameters (signal-to-noise ratio [SNR], contrast-to-noise ratio [CNR], image quality, diagnostic confidence) from a standard brain magnetic resonance imaging examination encompassing common neurological disorders such as demyelinating disease, gliomas, cerebrovascular disease, and epilepsy, with comparable sequence protocols and acquisition times at 3 T and at 7 T. MATERIALS AND METHODS Ten healthy volunteers and 4 subgroups of 40 patients in total underwent comparable magnetic resonance protocols with standard diffusion-weighted imaging, 2D and 3D turbo spin echo, 2D and 3D gradient echo and susceptibility-weighted imaging of the brain (10 sequences) at 3 T and 7 T. The subgroups comprised patients with either lesional (n = 5) or nonlesional (n = 4) epilepsy, intracerebral tumors (n = 11), demyelinating disease (n = 11) (relapsing-remitting multiple sclerosis [MS, n = 9], secondary progressive MS [n = 1], demyelinating disease not further specified [n = 1]), or chronic cerebrovascular disorders [n = 9]). For quantitative analysis, SNR and CNR were determined. For a semiquantitative assessment of the diagnostic confidence, a 10-point scale diagnostic confidence score (DCS) was applied. Two experienced radiologists with additional qualification in neuroradiology independently assessed, blinded to the field strength, 3 pathology-specific imaging criteria in each of the 4 disease groups and rated their diagnostic confidence. The overall image quality was semiquantitatively assessed using a 4-point scale taking into account whether diagnostic decision making was hampered by artifacts or not. RESULTS Without correction for spatial resolution, SNR was higher at 3 T except in the T2 SPACE 3D, DWI single shot, and DIR SPACE 3D sequences. The SNR corrected by the ratio of 3 T/7 T voxel sizes was higher at 7 T than at 3 T in 10 of 11 sequences (all except for T1 MP2RAGE 3D).In CNR, there was a wide variation between sequences and patient cohorts, but average CNR values were broadly similar at 3 T and 7 T.DCS values for all 4 pathologic entities were higher at 7 T than at 3 T. The DCS was significantly higher at 7 T for diagnosis and exclusion of cortical lesions in vascular disease. A tendency to higher DCS at 7 T for cortical lesions in MS was observed, and for the depiction of a central vein and iron deposits within MS lesions. Despite motion artifacts, DCS values were higher at 7 T for the diagnosis and exclusion of hippocampal sclerosis in mesial temporal lobe epilepsy (improved detection of the hippocampal subunits). Interrater agreement was 69.7% at 3 T and 93.3% at 7 T. There was no significant difference in the overall image quality score between 3 T and 7 T taking into account whether diagnostic decision making was hampered by artifacts or not. CONCLUSIONS Ultra-high-field magnetic resonance imaging at 7 T compared with 3 T yielded an improved diagnostic confidence in the most frequently encountered neurologic disorders. Higher spatial resolution and contrast were identified as the main contributory factors.

Geomorphology of the Imhotep region on comet 67P/Churyumov-Gerasimenko from OSIRIS observations

Context. Since August 2014, the OSIRIS Narrow Angle Camera (NAC) onboard the Rosetta spacecraft has acquired high spatial resolution images of the nucleus of comet 67P/Churyumov-Gerasimenko, down to the decimeter scale. This paper focuses on the Imhotep region, located on the largest lobe of the nucleus, near the equator. Aims. We map, inventory, and describe the geomorphology of the Imhotep region. We propose and discuss some processes to explain the formation and ongoing evolution of this region. Methods. We used OSIRIS NAC images, gravitational heights and slopes, and digital terrain models to map and measure the morphologies of Imhotep. Results. The Imhotep region presents a wide variety of terrains and morphologies: smooth and rocky terrains, bright areas, linear features, roundish features, and boulders. Gravity processes such as mass wasting and collapse play a significant role in the geomorphological evolution of this region. Cometary processes initiate erosion and are responsible for the formation of degassing conduits that are revealed by elevated roundish features on the surface. We also propose a scenario for the formation and evolution of the Imhotep region; this implies the presence of large primordial voids inside the nucleus, resulting from its formation process.