Simultaneous multislice acquisition with multi-contrast segmented EPI for separation of signal contributions in dynamic contrast-enhanced imaging

We present a method to efficiently separate signal in magnetic resonance imaging (MRI) into a base signal S0, representing the mainly T1-weighted component without T2*-relaxation, and its T2*-weighted counterpart by the rapid acquisition of multiple contrasts for advanced pharmacokinetic modelling. This is achieved by incorporating simultaneous multislice (SMS) imaging into a multi-contrast, segmented echo planar imaging (EPI) sequence to allow extended spatial coverage, which covers larger body regions without time penalty. Simultaneous acquisition of four slices was combined with segmented EPI for fast imaging with three gradient echo times in a preclinical perfusion study. Six female domestic pigs, German-landrace or hybrid-form, were scanned for 11 minutes respectively during administration of gadolinium-based contrast agent. Influences of reconstruction methods and training data were investigated. The separation into T1- and T2*-dependent signal contributions was achieved by fitting a standard analytical model to the acquired multi-echo data. The application of SMS yielded sufficient temporal resolution for the detection of the arterial input function in major vessels, while anatomical coverage allowed perfusion analysis of muscle tissue. The separation of the MR signal into T1- and T2*-dependent components allowed the correction of susceptibility related changes. We demonstrate a novel sequence for dynamic contrast-enhanced MRI that meets the requirements of temporal resolution (Δt < 1.5 s) and image quality. The incorporation of SMS into multi-contrast, segmented EPI can overcome existing limitations of dynamic contrast enhancement and dynamic susceptibility contrast methods, when applied separately. The new approach allows both techniques to be combined in a single acquisition with a large spatial coverage

Process-evaluation of tropospheric humidity simulated by general circulation models using water vapor isotopologues: 1. Comparison between models and observations

The goal of this study is to determine how H_2O and HDO measurements in water vapor can be used to detect and diagnose biases in the representation of processes controlling tropospheric humidity in atmospheric general circulation models (GCMs). We analyze a large number of isotopic data sets (four satellite, sixteen ground-based remote-sensing, five surface in situ and three aircraft data sets) that are sensitive to different altitudes throughout the free troposphere. Despite significant differences between data sets, we identify some observed HDO/H_2O characteristics that are robust across data sets and that can be used to evaluate models. We evaluate the isotopic GCM LMDZ, accounting for the effects of spatiotemporal sampling and instrument sensitivity. We find that LMDZ reproduces the spatial patterns in the lower and mid troposphere remarkably well. However, it underestimates the amplitude of seasonal variations in isotopic composition at all levels in the subtropics and in midlatitudes, and this bias is consistent across all data sets. LMDZ also underestimates the observed meridional isotopic gradient and the contrast between dry and convective tropical regions compared to satellite data sets. Comparison with six other isotope-enabled GCMs from the SWING2 project shows that biases exhibited by LMDZ are common to all models. The SWING2 GCMs show a very large spread in isotopic behavior that is not obviously related to that of humidity, suggesting water vapor isotopic measurements could be used to expose model shortcomings. In a companion paper, the isotopic differences between models are interpreted in terms of biases in the representation of processes controlling humidity

A conjugate gradient algorithm for the astrometric core solution of Gaia

The ESA space astrometry mission Gaia, planned to be launched in 2013, has been designed to make angular measurements on a global scale with micro-arcsecond accuracy. A key component of the data processing for Gaia is the astrometric core solution, which must implement an efficient and accurate numerical algorithm to solve the resulting, extremely large least-squares problem. The Astrometric Global Iterative Solution (AGIS) is a framework that allows to implement a range of different iterative solution schemes suitable for a scanning astrometric satellite. In order to find a computationally efficient and numerically accurate iteration scheme for the astrometric solution, compatible with the AGIS framework, we study an adaptation of the classical conjugate gradient (CG) algorithm, and compare it to the so-called simple iteration (SI) scheme that was previously known to converge for this problem, although very slowly. The different schemes are implemented within a software test bed for AGIS known as AGISLab, which allows to define, simulate and study scaled astrometric core solutions. After successful testing in AGISLab, the CG scheme has been implemented also in AGIS. The two algorithms CG and SI eventually converge to identical solutions, to within the numerical noise (of the order of 0.00001 micro-arcsec). These solutions are independent of the starting values (initial star catalogue), and we conclude that they are equivalent to a rigorous least-squares estimation of the astrometric parameters. The CG scheme converges up to a factor four faster than SI in the tested cases, and in particular spatially correlated truncation errors are much more efficiently damped out with the CG scheme.Comment: 24 pages, 16 figures. Accepted for publication in Astronomy & Astrophysic