Informatics for cross-sample analysis with comprehensive two-dimensional gas chromatography and high-resolution mass spectrometry (GCxGC-HRMS)

This paper describes informatics for cross-sample analysis with comprehensive two-dimensional gas chromatography (GCxGC) and high-resolution mass spectrometry (HRMS). GCxGC-HRMS analysis produces large data sets that are rich with information, but highly complex. The size of the data and volume of information requires automated processing for comprehensive cross-sample analysis, but the complexity poses a challenge for developing robust methods. The approach developed here analyzes GCxGC-HRMS data from multiple samples to extract a feature template that comprehensively captures the pattern of peaks detected in the retention-times plane. Then, for each sample chromatogram, the template is geometrically transformed to align with the detected peak pattern and generate a set of feature measurements for cross-sample analyses such as sample classification and biomarker discovery. The approach avoids the intractable problem of comprehensive peak matching by using a few reliable peaks for alignment and peak-based retention-plane windows to define comprehensive features that can be reliably matched for cross-sample analysis. The informatics are demonstrated with a set of 18 samples from breast-cancer tumors, each from different individuals, six each for Grades 1-3. The features allow classification that matches grading by a cancer pathologist with 78% success in leave-one-out cross-validation experiments. The HRMS signatures of the features of interest can be examined for determining elemental compositions and identifying compounds.

Two-phase change in CO2, Antarctic temperature and global climate during Termination II

The end of the Last Glacial Maximum (Termination I), roughly 20 thousand years ago (ka), was marked by cooling in the Northern Hemisphere, a weakening of the Asian monsoon, a rise in atmospheric CO2 concentrations and warming over Antarctica. The sequence of events associated with the previous glacial–interglacial transition (Termination II), roughly 136 ka, is less well constrained. Here we present high-resolution records of atmospheric CO2 concentrations and isotopic composition of N2—an atmospheric temperature proxy—from air bubbles in the EPICA Dome C ice core that span Termination II. We find that atmospheric CO2 concentrations and Antarctic temperature started increasing in phase around 136 ka, but in a second phase of Termination II, from 130.5 to 129 ka, the rise in atmospheric CO2 concentrations lagged that of Antarctic temperature unequivocally. We suggest that during this second phase, the intensification of the low-latitude hydrological cycle resulted in the development of a CO2 sink, which counteracted the CO2 outgassing from the Southern Hemisphere oceans over this period.

Grain coarsening in contact metamorphic carbonates: effects of second-phase particles, fluid flow and thermal perturbations

Under contact metamorphic conditions, carbonate rocks in the direct vicinity of the Adamello pluton reflect a temperature-induced grain coarsening. Despite this large-scale trend, a considerable grain size scatter occurs on the outcrop-scale indicating local influence of second-order effects such as thermal perturbations, fluid flow and second-phase particles. Second-phase particles, whose sizes range from nano- to the micron-scale, induce the most pronounced data scatter resulting in grain sizes too small by up to a factor of 10, compared with theoretical grain growth in a pure system. Such values are restricted to relatively impure samples consisting of up to 10 vol.% micron-scale second-phase particles, or to samples containing a large number of nano-scale particles. The obtained data set suggests that the second phases induce a temperature-controlled reduction on calcite grain growth. The mean calcite grain size can therefore be expressed in the form D 1⁄4 C2 eQ*/RT(dp/fp)m*, where C2 is a constant, Q* is an activation energy, T the temperature and m* the exponent of the ratio dp/fp, i.e. of the average size of the second phases divided by their volume fraction. However, more data are needed to obtain reliable values for C2 and Q*. Besides variations in the average grain size, the presence of second-phase particles generates crystal size distribution (CSD) shapes characterized by lognormal distributions, which differ from the Gaussian-type distributions of the pure samples. In contrast, fluid-enhanced grain growth does not change the shape of the CSDs, but due to enhanced transport properties, the average grain sizes increase by a factor of 2 and the variance of the distribution increases. Stable d18O and d13C isotope ratios in fluid-affected zones only deviate slightly from the host rock values, suggesting low fluid/rock ratios. Grain growth modelling indicates that the fluid-induced grain size variations can develop within several ka. As inferred from a combination of thermal and grain growth modelling, dykes with widths of up to 1 m have only a restricted influence on grain size deviations smaller than a factor of 1.1.To summarize, considerable grain size variations of up to one order of magnitude can locally result from second-order effects. Such effects require special attention when comparing experimentally derived grain growth kinetics with field studies.

Assessment of splanchnic blood flow using magnetic resonance imaging

BACKGROUND AND AIMS: The splanchnic circulation has an important function in the body under both physiological and pathophysiological conditions. Despite its importance, no reliable noninvasive procedures for estimating splanchnic circulation have been established. The aim of this study was to evaluate MRI as a tool for assessing intra-abdominal blood flows of the aorta, portal vein (VPO) and the major intestinal and hepatic vessels. METHODS: In nine healthy volunteers, the proximal aorta (AOP) and distal abdominal aorta (AOD), superior mesenteric artery (SAM), celiac trunk (CTR), hepatic arteries (common and proper hepatic arteries, AHC and AHP, respectively), and VPO were localized on contrast-enhanced magnetic resonance angiography images. Volumetric flow was measured using a two-dimensional cine echocardiogram-gated phase contrast technique. Measurements were taken before and 30 min after continuous intravenous infusion of somatostatin (250 microg/h) and were independently evaluated by two investigators. RESULTS: Blood flow measured by MRI in the VPO, SAM, AOP, AHP, and CTR significantly decreased after drug infusion. Flows in the AOD and AHC showed a tendency to decrease (P>0.05). Interrater agreement on flows in MRI was very good for large vessels (VPO, AOP, and AOD), with a concordance correlation coefficient of 0.94, as well as for smaller vessels such as the CTR, AHC, AHP, and SAM (concordance correlation coefficient =0.78). CONCLUSION: Somatostatin-induced blood flow changes in the splanchnic region were reliably detected by MRI. MRI may be useful for the noninvasive assessment of blood flow changes in the splanchnic region.

Time-lapse pressure tomography for characterizing CO2 plume evolution in a deep saline aquifer

A time-lapse pressure tomography inversion approach is applied to characterize the CO2 plume development in a virtual deep saline aquifer. Deep CO2 injection leads to flow properties of the mixed-phase, which vary depending on the CO2 saturation. Analogous to the crossed ray paths of a seismic tomographic experiment, pressure tomography creates streamline patterns by injecting brine prior to CO2 injection or by injecting small amounts of CO2 into the two-phase (brine and CO2) system at different depths. In a first step, the introduced pressure responses at observation locations are utilized for a computationally rapid and efficient eikonal equation based inversion to reconstruct the heterogeneity of the subsurface with diffusivity (D) tomograms. Information about the plume shape can be derived by comparing D-tomograms of the aquifer at different times. In a second step, the aquifer is subdivided into two zones of constant values of hydraulic conductivity (K) and specific storage (Ss) through a clustering approach. For the CO2 plume, mixed-phase K and Ss values are estimated by minimizing the difference between calculated and “true” pressure responses using a single-phase flow simulator to reduce the computing complexity. Finally, the estimated flow property is converted to gas saturation by a single-phase proxy, which represents an integrated value of the plume. This novel approach is tested first with a doublet well configuration, and it reveals a great potential of pressure tomography based concepts for characterizing and monitoring deep aquifers, as well as the evolution of a CO2 plume. Still, field-testing will be required for better assessing the applicability of this approach.

Planet formation models: the interplay with the planetesimal disc

Context. According to the sequential accretion model (or core-nucleated accretion model), giant planet formation is based first on the formation of a solid core which, when massive enough, can gravitationally bind gas from the nebula to form the envelope. The most critical part of the model is the formation time of the core: to trigger the accretion of gas, the core has to grow up to several Earth masses before the gas component of the protoplanetary disc dissipates. Aims: We calculate planetary formation models including a detailed description of the dynamics of the planetesimal disc, taking into account both gas drag and excitation of forming planets. Methods: We computed the formation of planets, considering the oligarchic regime for the growth of the solid core. Embryos growing in the disc stir their neighbour planetesimals, exciting their relative velocities, which makes accretion more difficult. Here we introduce a more realistic treatment for the evolution of planetesimals' relative velocities, which directly impact on the formation timescale. For this, we computed the excitation state of planetesimals, as a result of stirring by forming planets, and gas-solid interactions. Results: We find that the formation of giant planets is favoured by the accretion of small planetesimals, as their random velocities are more easily damped by the gas drag of the nebula. Moreover, the capture radius of a protoplanet with a (tiny) envelope is also larger for small planetesimals. However, planets migrate as a result of disc-planet angular momentum exchange, with important consequences for their survival: due to the slow growth of a protoplanet in the oligarchic regime, rapid inward type I migration has important implications on intermediate-mass planets that have not yet started their runaway accretion phase of gas. Most of these planets are lost in the central star. Surviving planets have masses either below 10 M⊕ or above several Jupiter masses. Conclusions: To form giant planets before the dissipation of the disc, small planetesimals (~0.1 km) have to be the major contributors of the solid accretion process. However, the combination of oligarchic growth and fast inward migration leads to the absence of intermediate-mass planets. Other processes must therefore be at work to explain the population of extrasolar planets that are presently known.

Prototype of the gas chromatograph-mass spectrometer to investigate volatile species in the lunar soil for the Luna-Resurs mission

In preparation for the Russian Luna-Resurs mission we combined our compact time-of-flight mass spectrometer (TOF-MS) with a chemical pre-separation of the species by gas chromatography (GC). Coupled measurements with both instruments were successfully performed with the prototype of the mass spectrometer and a flight-like gas chromatograph. The system was tested with two test gas mixtures, a mixture of hydrocarbons and a mixture of noble gases. Due to its capability to record mass spectra over the full mass range at once with high sensitivity and a dynamic range of up to 10(6) within 1 s, the TOF-MS system is a valuable extension of the GC analytical system. Based on the measurements with calibration gases performed with the combined GC-MS prototype and under assumption of mean characteristics for the Moon's regolith, the detection limit for volatile species in a soil sample is estimated to 2.10(-10) by mass for hydrocarbons and 2.10(-9) by mass for noble gases. (C) 2015 Elsevier Ltd. All rights reserved.

Determining gypsum growth temperatures using monophase fluid inclusions-Application to the giant gypsum crystals of Naica, Mexico

Determining the formation temperature of minerals using fluid inclusions is a crucial step in understanding rock-forming scenarios. Unfortunately, fluid inclusions in minerals formed at low temperature, such as gypsum, are commonly in a metastable monophase liquid state. To overcome this problem, ultra-short laser pulses can be used to induce vapor bubble nucleation, thus creating a stable two-phase fluid inclusion appropriate for subsequent measurements of the liquid-vapor homogenization temperature, T-h. In this study we evaluate the applicability of T-h data to accurately determine gypsum formation temperatures. We used fluid inclusions in synthetic gypsum crystals grown in the laboratory at different temperatures between 40 degrees C and 80 degrees C under atmospheric pressure conditions. We found an asymmetric distribution of the T-h values, which are systematically lower than the actual crystal growth temperatures, T-g; this is due to (1) the effect of surface tension on liquid-vapor homogenization, and (2) plastic deformation of the inclusion walls due to internal tensile stress occurring in the metastable state of the inclusions. Based on this understanding, we have determined growth temperatures of natural giant gypsum crystals from Naica (Mexico), yielding 47 +/- 1.5 degrees C for crystals grown in the Cave of Swords (120 m below surface) and 54.5 +/- 2 degrees C for giant crystals grown in the Cave of Crystals (290 m below surface). These results support the earlier hypothesis that the population and the size of the Naica crystals were controlled by temperature. In addition, this experimental method opens a door to determining the growth temperature of minerals forming in low-temperature environments.

Stereoselective biotransformation of ketamine in equine liver and lung microsomes

Stereoselectivity has to be considered for pharmacodynamic and pharmacokinetic features of ketamine. Stereoselective biotransformation of ketamine was investigated in equine microsomes in vitro. Concentration curves were constructed over time, and enzyme activity was determined for different substrate concentrations using equine liver and lung microsomes. The concentrations of R/S-ketamine and R/S-norketamine were determined by enantioselective capillary electrophoresis. A two-phase model based on Hill kinetics was used to analyze the biotransformation of R/S-ketamine into R/S-norketamine and, in a second step, into R/S-downstream metabolites. In liver and lung microsomes, levels of R-ketamine exceeded those of S-ketamine at all time points and S-norketamine exceeded R-norketamine at time points below the maximum concentration. In liver and lung microsomes, significant differences in the enzyme velocity (V(max)) were observed between S- and R-norketamine formation and between V(max) of S-norketamine formation when S-ketamine was compared to S-ketamine of the racemate. Our investigations in microsomal reactions in vitro suggest that stereoselective ketamine biotransformation in horses occurs in the liver and the lung with a slower elimination of S-ketamine in the presence of R-ketamine. Scaling of the in vitro parameters to liver and lung organ clearances provided an excellent fit with previously published in vivo data and confirmed a lung first-pass effect.

Concentration-dependent isoflurane effects on withdrawal reflexes in pigs and the role of the stimulation paradigm

In this prospective two-phase experimental trial, 10 pigs were anaesthetized twice with isoflurane only. In the first phase, the individual minimum alveolar concentration (MAC) was determined and in the second phase the effects on withdrawal reflexes of increasing end-tidal isoflurane concentrations (from 1.6% to 2.8%) were assessed. Single, 10 and 60 repeated electrical stimulations were used to evoke withdrawal reflexes which were recorded and quantified by electromyography. Recruitment curves for reflex amplitude for increasing stimulation intensities and isoflurane concentrations were constructed. Isoflurane MAC was 1.9+/-0.3%. Reflexes evoked by repeated stimulation were suppressed at isoflurane concentrations significantly higher than those which suppressed complex movements during MAC determination (P=0.014 and P=0.006 for 10 and 60 repeated stimuli respectively). Isoflurane up to 2.8% was still not able to abolish reflex activity evoked by repeated stimulations in all pigs. Single stimulation reflexes were suppressed at significantly lower concentrations than repeated stimulation reflexes (P=0.008 and P=0.004 for 10 and 60 repeated stimuli, respectively). Reflex amplitude was significantly correlated with isoflurane concentration (P<0.001, r=-0.85) independent of the individual MAC. The findings indicate that the level at which isoflurane suppresses withdrawal reflexes is dependent on the stimulation paradigm (single vs. repeated electrical stimulation), and there is limited value in expressing reflex withdrawal suppression in terms of MAC as purposeful and reflex movements are independently affected by isoflurane in individual animals.

Solute transport in crystalline rocks at Äspö — II: Blind predictions, inverse modelling and lessons learnt from test STT1

Based on the results from detailed structural and petrological characterisation and on up-scaled laboratory values for sorption and diffusion, blind predictions were made for the STT1 dipole tracer test performed in the Swedish A¨ spo¨ Hard Rock Laboratory. The tracers used were nonsorbing, such as uranine and tritiated water, weakly sorbing 22Na+, 85Sr2 +, 47Ca2 +and more strongly sorbing 86Rb+, 133Ba2 +, 137Cs+. Our model consists of two parts: (1) a flow part based on a 2D-streamtube formalism accounting for the natural background flow field and with an underlying homogeneous and isotropic transmissivity field and (2) a transport part in terms of the dual porosity medium approach which is linked to the flow part by the flow porosity. The calibration of the model was done using the data from one single uranine breakthrough (PDT3). The study clearly showed that matrix diffusion into a highly porous material, fault gouge, had to be included in our model evidenced by the characteristic shape of the breakthrough curve and in line with geological observations. After the disclosure of the measurements, it turned out that, in spite of the simplicity of our model, the prediction for the nonsorbing and weakly sorbing tracers was fairly good. The blind prediction for the more strongly sorbing tracers was in general less accurate. The reason for the good predictions is deemed to be the result of the choice of a model structure strongly based on geological observation. The breakthrough curves were inversely modelled to determine in situ values for the transport parameters and to draw consequences on the model structure applied. For good fits, only one additional fracture family in contact with cataclasite had to be taken into account, but no new transport mechanisms had to be invoked. The in situ values for the effective diffusion coefficient for fault gouge are a factor of 2–15 larger than the laboratory data. For cataclasite, both data sets have values comparable to laboratory data. The extracted Kd values for the weakly sorbing tracers are larger than Swedish laboratory data by a factor of 25–60, but agree within a factor of 3–5 for the more strongly sorbing nuclides. The reason for the inconsistency concerning Kds is the use of fresh granite in the laboratory studies, whereas tracers in the field experiments interact only with fracture fault gouge and to a lesser extent with cataclasite both being mineralogically very different (e.g. clay-bearing) from the intact wall rock.

Safety and tolerability of slow-release oral morphine versus methadone in the treatment of opioid dependence

Opioid substitution treatment (OST) for opioid dependence may be limited by adverse events (AEs). Increasing the range of therapeutic options optimizes outcomes and facilitates patient management. An international, multi-center, two-phase study investigated the efficacy and safety of slow-release oral morphine (SROM) versus methadone in patients receiving methadone therapy for opioid dependence. In phase 1 (two way cross-over, 11 weeks each period) patients were randomized to SROM or methadone oral solution. In phase 2 (25 weeks), patients continued treatment with SROM (group A) or switched from methadone to SROM (group B). In total, 211 out of 276 completed phase 1 and 198 entered phase 2 (n = 95 group A, n = 103 group B). Treatment with both SROM and methadone was well tolerated. However, the mean QTc-interval associated with methadone was significantly longer than that under SROM. Higher treatment satisfaction, fewer cravings for heroin, and lower mental stress were reported with SROM. This study adds a significant further weight of evidence that SROM is an effective and well tolerated long-term maintenance treatment for opioid dependence with a beneficial risk profile compared to methadone regarding cardiac effects and supports its clinical utility.

Diagenesis of a light, tight-oil chert reservoir at the Ediacaran/Cambrian boundary, Sultanate of Oman

The Al Shomou Silicilyte Member (Athel Formation) in the South Oman Salt Basin shares many of the characteristics of a light, tight-oil (LTO) reservoir: it is a prolifi c source rock mature for light oil, it produces light oil from a very tight matrix and reservoir, and hydraulic fracking technology is required to produce the oil. What is intriguing about the Al Shomou Silicilyte, and different from other LTO reservoirs, is its position related to the Precambrian/Cambrian Boundary (PCB) and the fact that it is a ‘laminated chert‘ rather than a shale. In an integrated diagenetic study we applied microstructural analyses (SEM, BSE) combined with state-of-the-art stable isotope and trace element analysis of the silicilyte matrix and fractures. Fluid inclusion microthermometry was applied to record the salinity and minimum trapping temperatures. The microstructural investigations reveal a fi ne lamination of the silicilyte matrix with a mean lamina thickness of ca. 20 μm consisting of predominantly organic matter-rich and fi nely crystalline quartz-rich layers, respectively. Authigenic, micron-sized idiomorphic quartz crystals are the main matrix components of the silicilyte. Other diagenetic phases are pyrite, apatite, dolomite, magnesite and barite cements. Porosity values based on neutron density logs and core plug data indicate porosity in the silicilyte ranges from less than 2% to almost to 40%. The majority of the pore space in the silicilyte is related to (primary) inter-crystalline pores, with locally important oversized secondary pores. Pore casts of the silica matrix show that pores are extremely irregular in three dimensions, and are generally interconnected by a complex web or meshwork of fi ne elongate pore throats. Mercury injection capillary data are in line with the microstructural observations suggesting two populations of pore throats, with an effective average modal diameter of 0.4 μm. The acquired geochemical data support the interpretation that the primary source of the silica is the ambient seawater rather than hydrothermal or biogenic. A maximum temperature of ca. 45°C for the formation of microcrystalline quartz in the silicilyte is good evidence that the lithifi cation and crystallization of quartz occurred in the fi rst 5 Ma after deposition. Several phases of brittle fracturing and mineralization occurred in response to salt tectonics during burial. The sequences of fracture-fi lling mineral phases (dolomite - layered chalcedony – quartz – apatite - magnesite I+II - barite – halite) indicates a complex fl uid evolution after silicilyte lithifi cation. Primary, all-liquid fl uid inclusions in the fracturefi lling quartz are good evidence of growth beginning at low temperatures, i.e. ≤ 50ºC. Continuous precipitation during increasing temperature and burial is documented by primary two-phase fl uid inclusions in quartz cements that show brines at 50°C and fi rst hydrocarbons at ca. 70°C. The absolute timing of each mineral phase can be constrained based on U-Pb geochronometry, and basin modelling. Secondary fl uid inclusions in quartz, magnesite and barite indicate reactivation of the fracture system after peak burial temperature during the major cooling event, i.e. uplift, between 450 and 310 Ma. A number of fi rst-order trends in porosity and reservoir-quality distribution are observed which are strongly related to the diagenetic and fl uid history of the reservoir: the early in-situ generation of hydrocarbons and overpressure development arrests diagenesis and preserves matrix porosity. Chemical compaction by pressure dissolution in the fl ank areas could be a valid hypothesis to explain the porosity variations in the silicilitye slabs resulting in lower porosity and poorer connectivity on the fl anks of the reservoir. Most of the hydrocarbon storage and production comes from intervals characterized by Amthor et al. 114488 preserved micropores, not hydrocarbon storage in a fracture system. The absence of oil expulsion results in present-day high oil saturations. The main diagenetic modifi cations of the silicilyte occurred and were completed relatively early in its history, i.e. before 300 Ma. An instrumental factor for preserving matrix porosity is the diffi culty for a given slab to evacuate all the fl uids (water and hydrocarbons), or in other words, the very good sealing capacity of the salt embedding the slab.