Inter-site and inter-scanner diffusion MRI data harmonization.

We propose a novel method to harmonize diffusion MRI data acquired from multiple sites and scanners, which is imperative for joint analysis of the data to significantly increase sample size and statistical power of neuroimaging studies. Our method incorporates the following main novelties: i) we take into account the scanner-dependent spatial variability of the diffusion signal in different parts of the brain; ii) our method is independent of compartmental modeling of diffusion (e.g., tensor, and intra/extra cellular compartments) and the acquired signal itself is corrected for scanner related differences; and iii) inter-subject variability as measured by the coefficient of variation is maintained at each site. We represent the signal in a basis of spherical harmonics and compute several rotation invariant spherical harmonic features to estimate a region and tissue specific linear mapping between the signal from different sites (and scanners). We validate our method on diffusion data acquired from seven different sites (including two GE, three Philips, and two Siemens scanners) on a group of age-matched healthy subjects. Since the extracted rotation invariant spherical harmonic features depend on the accuracy of the brain parcellation provided by Freesurfer, we propose a feature based refinement of the original parcellation such that it better characterizes the anatomy and provides robust linear mappings to harmonize the dMRI data. We demonstrate the efficacy of our method by statistically comparing diffusion measures such as fractional anisotropy, mean diffusivity and generalized fractional anisotropy across multiple sites before and after data harmonization. We also show results using tract-based spatial statistics before and after harmonization for independent validation of the proposed methodology. Our experimental results demonstrate that, for nearly identical acquisition protocol across sites, scanner-specific differences can be accurately removed using the proposed method.

Size-partitioned phytoplankton carbon concentrations retrieved from ocean color data, links to data in netCDF format

Owing to their important roles in biogeochemical cycles, phytoplankton functional types (PFTs) have been the aim of an increasing number of ocean color algorithms. Yet, none of the existing methods are based on phytoplankton carbon (C) biomass, which is a fundamental biogeochemical and ecological variable and the "unit of accounting" in Earth system models. We present a novel bio-optical algorithm to retrieve size-partitioned phytoplankton carbon from ocean color satellite data. The algorithm is based on existing methods to estimate particle volume from a power-law particle size distribution (PSD). Volume is converted to carbon concentrations using a compilation of allometric relationships. We quantify absolute and fractional biomass in three PFTs based on size - picophytoplankton (0.5-2 µm in diameter), nanophytoplankton (2-20 µm) and microphytoplankton (20-50 µm). The mean spatial distributions of total phytoplankton C biomass and individual PFTs, derived from global SeaWiFS monthly ocean color data, are consistent with current understanding of oceanic ecosystems, i.e., oligotrophic regions are characterized by low biomass and dominance of picoplankton, whereas eutrophic regions have high biomass to which nanoplankton and microplankton contribute relatively larger fractions. Global climatological, spatially integrated phytoplankton carbon biomass standing stock estimates using our PSD-based approach yield - 0.25 Gt of C, consistent with analogous estimates from two other ocean color algorithms and several state-of-the-art Earth system models. Satisfactory in situ closure observed between PSD and POC measurements lends support to the theoretical basis of the PSD-based algorithm. Uncertainty budget analyses indicate that absolute carbon concentration uncertainties are driven by the PSD parameter No which determines particle number concentration to first order, while uncertainties in PFTs' fractional contributions to total C biomass are mostly due to the allometric coefficients. The C algorithm presented here, which is not empirically constrained a priori, partitions biomass in size classes and introduces improvement over the assumptions of the other approaches. However, the range of phytoplankton C biomass spatial variability globally is larger than estimated by any other models considered here, which suggests an empirical correction to the No parameter is needed, based on PSD validation statistics. These corrected absolute carbon biomass concentrations validate well against in situ POC observations.

Fractional abundances of isoprenoid and branched glycerol dialkyl glycerol tetraethers (GDGT) of cores SO201-2-12KL and SO201-2-114KL

It has been proposed that North Pacific sea surface temperature (SST) evolution was intimately linked to North Atlantic climate oscillations during the last glacial-interglacial transition. However, during the early deglaciation and the Last Glacial Maximum, the SST development in the subarctic northwest Pacific and the Bering Sea is poorly constrained as most existing deglacial SST records are based on alkenone paleothermometry, which is limited prior to 15 ka B.P. in the subarctic North Pacific realm. By applying the TEXL86 temperature proxy we obtain glacial-Holocene-SST records for the marginal northwest Pacific and the Western Bering Sea. Our TEXL86-based records and existing alkenone data suggest that during the past 15.5 ka, SSTs in the northwest Pacific and the Western Bering Sea closely followed millennial-scale climate fluctuations known from Greenland ice cores, indicating rapid atmospheric teleconnections with abrupt climate changes in the North Atlantic. Our SST reconstructions indicate that in the Western Bering Sea SSTs drop significantly during Heinrich Stadial 1 (HS1), similar to the known North Atlantic climate history. In contrast, progressively rising SST in the northwest Pacific is different to the North Atlantic climate development during HS1. Similarities between the northwest Pacific SST and climate records from the Gulf of Alaska point to a stronger influence of Alaskan Stream waters connecting the eastern and western basin of the North Pacific during this time. During the Holocene, dissimilar climate trends point to reduced influence of the Alaskan Stream in the northwest Pacific.

On independent component analysis applied to unmixing hyperspectral data

One of the most challenging task underlying many hyperspectral imagery applications is the spectral unmixing, which decomposes a mixed pixel into a collection of reectance spectra, called endmember signatures, and their corresponding fractional abundances. Independent Component Analysis (ICA) have recently been proposed as a tool to unmix hyperspectral data. The basic goal of ICA is to nd a linear transformation to recover independent sources (abundance fractions) given only sensor observations that are unknown linear mixtures of the unobserved independent sources. In hyperspectral imagery the sum of abundance fractions associated to each pixel is constant due to physical constraints in the data acquisition process. Thus, sources cannot be independent. This paper address hyperspectral data source dependence and its impact on ICA performance. The study consider simulated and real data. In simulated scenarios hyperspectral observations are described by a generative model that takes into account the degradation mechanisms normally found in hyperspectral applications. We conclude that ICA does not unmix correctly all sources. This conclusion is based on the a study of the mutual information. Nevertheless, some sources might be well separated mainly if the number of sources is large and the signal-to-noise ratio (SNR) is high.

Modeling some real phenomena by fractional differential equations

This paper deals with fractional differential equations, with dependence on a Caputo fractional derivative of real order. The goal is to show, based on concrete examples and experimental data from several experiments, that fractional differential equations may model more efficiently certain problems than ordinary differential equations. A numerical optimization approach based on least squares approximation is used to determine the order of the fractional operator that better describes real data, as well as other related parameters.

Part-time and fractional staff in the learning and skills sector

This Universities and College Union Launch Event presentation reported on the findings of Learning and Skills Research Network (LSRN) London and South East (LSE) Regional Research Project. The presentation reflected on research carried out during 2002-06 on the development and deployment of part-time staff in the Learning and Skills Sector. Although the lifelong learning sector is the largest UK education sector, little attention has as yet been paid to the role of LSC sector part-time staff. Worrying trends of an increasing casualisation of staffing have been reported. The role of part-timers as highly committed (philanthropic) but generally underpaid and exploited staff (ragged-trousered) emerged from the data collected by this investigation, which examined the role of part-timers in several colleges and adult education institutions in London and the South East. The metaphor of the 'ragged-trousered philanthropist' was consciously selected to investigate the interactivity between philantrophy, employment practices for PT staff, and education as social action, in addressing the need for good practice to achieve quality outcomes in learning and teaching. The results are to some extent transferable to other education and training sectors employing part-time staff, e.g. higher education institutions and work-based training organisations.

Detailed analysis of a conservative difference approximation for the time fractional diffusion equation

Diffusion equations that use time fractional derivatives are attractive because they describe a wealth of problems involving non-Markovian Random walks. The time fractional diffusion equation (TFDE) is obtained from the standard diffusion equation by replacing the first-order time derivative with a fractional derivative of order α ∈ (0, 1). Developing numerical methods for solving fractional partial differential equations is a new research field and the theoretical analysis of the numerical methods associated with them is not fully developed. In this paper an explicit conservative difference approximation (ECDA) for TFDE is proposed. We give a detailed analysis for this ECDA and generate discrete models of random walk suitable for simulating random variables whose spatial probability density evolves in time according to this fractional diffusion equation. The stability and convergence of the ECDA for TFDE in a bounded domain are discussed. Finally, some numerical examples are presented to show the application of the present technique.