Determination of respiratory gas flow by electrical impedance tomography in an animal model of mechanical ventilation

Background A recent method determines regional gas flow of the lung by electrical impedance tomography (EIT). The aim of this study is to show the applicability of this method in a porcine model of mechanical ventilation in healthy and diseased lungs. Our primary hypothesis is that global gas flow measured by EIT can be correlated with spirometry. Our secondary hypothesis is that regional analysis of respiratory gas flow delivers physiologically meaningful results. Methods In two sets of experiments n = 7 healthy pigs and n = 6 pigs before and after induction of lavage lung injury were investigated. EIT of the lung and spirometry were registered synchronously during ongoing mechanical ventilation. In-vivo aeration of the lung was analysed in four regions-of-interest (ROI) by EIT: 1) global, 2) ventral (non-dependent), 3) middle and 4) dorsal (dependent) ROI. Respiratory gas flow was calculated by the first derivative of the regional aeration curve. Four phases of the respiratory cycle were discriminated. They delivered peak and late inspiratory and expiratory gas flow (PIF, LIF, PEF, LEF) characterizing early or late inspiration or expiration. Results Linear regression analysis of EIT and spirometry in healthy pigs revealed a very good correlation measuring peak flow and a good correlation detecting late flow. PIFEIT = 0.702 · PIFspiro + 117.4, r2 = 0.809; PEFEIT = 0.690 · PEFspiro-124.2, r2 = 0.760; LIFEIT = 0.909 · LIFspiro + 27.32, r2 = 0.572 and LEFEIT = 0.858 · LEFspiro-10.94, r2 = 0.647. EIT derived absolute gas flow was generally smaller than data from spirometry. Regional gas flow was distributed heterogeneously during different phases of the respiratory cycle. But, the regional distribution of gas flow stayed stable during different ventilator settings. Moderate lung injury changed the regional pattern of gas flow. Conclusions We conclude that the presented method is able to determine global respiratory gas flow of the lung in different phases of the respiratory cycle. Additionally, it delivers meaningful insight into regional pulmonary characteristics, i.e. the regional ability of the lung to take up and to release air.

The spatial distribution of water in the inner coma of Comet 9P/Tempel 1: Comparison between models and observations

The near nucleus coma of Comet 9P/Tempel 1 has been simulated with the 3D Direct Simulation Monte Carlo (DSMC) code PDSC++ (Su, C.-C. [2013]. Parallel Direct Simulation Monte Carlo (DSMC) Methods for Modeling Rarefied Gas Dynamics. PhD Thesis, National Chiao Tung University, Taiwan) and the derived column densities have been compared to observations of the water vapour distribution found by using infrared imaging spectrometer on the Deep Impact spacecraft (Feaga, L.M., A’Hearn, M.F., Sunshine, J.M., Groussin, O., Farnham, T.L. [2007]. Icarus 191(2), 134–145. http://dx.doi.org/10.1016/j.icarus.2007.04.038). Modelled total production rates are also compared to various observations made at the time of the Deep Impact encounter. Three different models were tested. For all models, the shape model constructed from the Deep Impact observations by Thomas et al. (Thomas, P.C., Veverka, J., Belton, M.J.S., Hidy, A., A’Hearn, M.F., Farnham, T.L., et al. [2007]. Icarus, 187(1), 4–15. http://dx.doi.org/10.1016/j.icarus.2006.12.013) was used. Outgassing depending only on the cosine of the solar insolation angle on each shape model facet is shown to provide an unsatisfactory model. Models constructed on the basis of active areas suggested by Kossacki and Szutowicz (Kossacki, K., Szutowicz, S. [2008]. Icarus, 195(2), 705–724. http://dx.doi.org/10.1016/j.icarus.2007.12.014) are shown to be superior. The Kossacki and Szutowicz model, however, also shows deficits which we have sought to improve upon. For the best model we investigate the properties of the outflow.

Measurement of flow harmonics with multi-particle cumulants in Pb+Pb collisions at √sNN=2.76 TeV with the ATLAS detector

ATLAS measurements of the azimuthal anisotropy in lead–lead collisions at √sNN = 2.76 TeV are shown using a dataset of approximately 7μb−1 collected at the LHC in 2010. The measurements are performed for charged particles with transversemomenta 0.5 < pT < 20 GeV and in the pseudorapidity range |η| < 2.5. The anisotropy is characterized by the Fourier coefficients, vn, of the charged-particle azimuthal angle distribution for n = 2–4. The Fourier coefficients are evaluated using multi-particle cumulants calculated with the generating function method. Results on the transverse momentum, pseudorapidity and centrality dependence of the vn coefficients are presented. The elliptic flow, v2, is obtained from the two-, four-, six- and eight-particle cumulants while higher-order coefficients, v3 and v4, are determined with two- and four-particle cumulants. Flow harmonics vn measured with four-particle cumulants are significantly reduced compared to the measurement involving two-particle cumulants. A comparison to vn measurements obtained using different analysis methods and previously reported by the LHC experiments is also shown. Results of measurements of flow fluctuations evaluated with multiparticle cumulants are shown as a function of transverse momentum and the collision centrality. Models of the initial spatial geometry and its fluctuations fail to describe the flow fluctuations measurements.

Identification, quantification, spatiotemporal distribution and genetic variation of major latex secondary metabolites in the common dandelion (Taraxacum officinale agg.)

The secondary metabolites in the roots, leaves and flowers of the common dandelion (Taraxacum officinale agg.) have been studied in detail. However, little is known about the specific constituents of the plant’s highly specialized laticifer cells. Using a combination of liquid and gas chromatography, mass spectrometry and nuclear magnetic resonance spectrometry, we identified and quantified the major secondary metabolites in the latex of different organs across different growth stages in three genotypes, and tested the activity of the metabolites against the generalist root herbivore Diabrotica balteata. We found that common dandelion latex is dominated by three classes of secondary metabolites: phenolic inositol esters (PIEs), triterpene acetates (TritAc) and the sesquiterpene lactone taraxinic acid β-d-glucopyranosyl ester (TA-G). Purification and absolute quantification revealed concentrations in the upper mg g−1 range for all compound classes with up to 6% PIEs, 5% TritAc and 7% TA-G per gram latex fresh weight. Contrary to typical secondary metabolite patterns, concentrations of all three classes increased with plant age. The highest concentrations were measured in the main root. PIE profiles differed both quantitatively and qualitatively between plant genotypes, whereas TritAc and TA-G differed only quantitatively. Metabolite concentrations were positively correlated within and between the different compound classes, indicating tight biosynthetic co-regulation. Latex metabolite extracts strongly repelled D. balteata larvae, suggesting that the latex constituents are biologically active.

Estimating Particle Shape and Orientation Using Volume Tensors

We develop statistical procedures for estimating shape and orientation of arbitrary three-dimensional particles. We focus on the case where particles cannot be observed directly, but only via sections. Volume tensors are used for describing particle shape and orientation, and we derive stereological estimators of the tensors. These estimators are combined to provide consistent estimators of the moments of the so-called particle cover density. The covariance structure associated with the particle cover density depends on the orientation and shape of the particles. For instance, if the distribution of the typical particle is invariant under rotations, then the covariance matrix is proportional to the identity matrix. We develop a non-parametric test for such isotropy. A flexible Lévy-based particle model is proposed, which may be analysed using a generalized method of moments in which the volume tensors enter. The developed methods are used to study the cell organization in the human brain cortex.

Gas transport in firn: multiple-tracer characterisation and model intercomparison for NEEM, Northern Greenland

Air was sampled from the porous firn layer at the NEEM site in Northern Greenland. We use an ensemble of ten reference tracers of known atmospheric history to characterise the transport properties of the site. By analysing uncertainties in both data and the reference gas atmospheric histories, we can objectively assign weights to each of the gases used for the depth-diffusivity reconstruction. We define an objective root mean square criterion that is minimised in the model tuning procedure. Each tracer constrains the firn profile differently through its unique atmospheric history and free air diffusivity, making our multiple-tracer characterisation method a clear improvement over the commonly used single-tracer tuning. Six firn air transport models are tuned to the NEEM site; all models successfully reproduce the data within a 1σ Gaussian distribution. A comparison between two replicate boreholes drilled 64 m apart shows differences in measured mixing ratio profiles that exceed the experimental error. We find evidence that diffusivity does not vanish completely in the lock-in zone, as is commonly assumed. The ice age- gas age difference (1 age) at the firn-ice transition is calculated to be 182+3−9 yr. We further present the first intercomparison study of firn air models, where we introduce diagnostic scenarios designed to probe specific aspects of the model physics. Our results show that there are major differences in the way the models handle advective transport. Furthermore, diffusive fractionation of isotopes in the firn is poorly constrained by the models, which has consequences for attempts to reconstruct the isotopic composition of trace gases back in time using firn air and ice core records.

Distribution of cutin and suberin biomarkers under forest trees with different root systems

Background and aims Differences in chemical composition of root compounds and root systems among tree species may affect organic matter (OM) distribution, source and composition in forest soils. The objective of this study was to elucidate the contribution of species specific cutin and suberin biomarkers as proxies for shoot- and root-derived organic carbon (OC) to soil OM at different depths with increasing distance to the stems of four different tree species. Methods The contribution of cutin- and suberin-derived lipids to OM in a Cutanic Alisol was analyzed with increasing soil depth and distance to the stems of Fagus sylvatica L., Picea abies (L.) Karst., Quercus robur L. and Pseudotsuga menziesii (Mirb.) Franco. Cutin and suberin monomers of plants and soils were analyzed by alkaline hydrolysis and subsequent gas chromatography–mass spectrometry. Results The amount and distribution of suberin-derived lipids in soil clearly reflected the specific root system of the different tree species. The amount of cutin-derived lipids decreased strongly with soil depth, indicating that the input of leaf/needle material is restricted to the topsoil. In contrast to the suberin-derived lipids, the spatial pattern of cutin monomer contribution to soil OM did not depend on tree species. Conclusions Our results document the importance of tree species as a main factor controlling the composition and distribution of OM in forest soils. They reveal the impact of tree species on root-derived OM distribution and the necessity to distinguish among different zones when studying soil OM storage in forests.

Size parameters, size-class distribution and area-number relationship of microscopic charcoal: relevance for fire reconstruction

Charcoal analysis was conducted on sediment cores from three lakes to assess the relationship between the area and number of charcoal particles. Three charcoal-size parameters (maximum breadth, maximum length and area) were measured on sediment samples representing various vegetation types, including shrub tundra, boreal forest and temperate forest. These parameters and charcoal size-class distributions do not differ statistically between two sites where the same preparation technique (glycerine pollen slides) was used, but they differ for the same core when different techniques were applied. Results suggest that differences in charcoal size and size-class distribution are mainly caused by different preparation techniques and are not related to vegetation-type variation. At all three sites, the area and number concentrations of charcoal particles are highly correlated in standard pollen slides; 82–83% of the variability of the charcoal-area concentration can be explained by the particle-number concentration. Comparisons between predicted and measured area concentrations show that regression equations linking charcoal number and area concentrations can be used across sites as long as the same pollen-preparation technique is used. Thus it is concluded that it is unnecessary to measure charcoal areas in standard pollen slides – a time-consuming and tedious process.

Kinetic simulation of neutral∕ionized gas and electrically charged dust in the coma of comet 67P∕Churyumov-Gerasimenko

The cometary coma is a unique phenomenon in the solar system being a planetary atmosphere influenced by little or no gravity. As a comet approaches the sun, the water vapor with some fraction of other gases sublimate, generating a cloud of gas, ice and other refractory materials (rocky and organic dust) ejected from the surface of the nucleus. Sublimating gas molecules undergo frequent collisions and photochemical processes in the near‐nucleus region. Owing to its negligible gravity, comets produce a large and highly variable extensive dusty coma with a size much larger than the characteristic size of the cometary nucleus. The Rosetta spacecraft is en route to comet 67P/Churyumov‐Gerasimenko for a rendezvous, landing, and extensive orbital phase beginning in 2014. Both, interpretation of measurements and safety consideration of the spacecraft require modeling of the comet’s dusty gas environment. In this work we present results of a numerical study of multispecies gaseous and electrically charged dust environment of comet Chyuryumov‐Gerasimenko. Both, gas and dust phases of the coma are simulated kinetically. Photolytic reactions are taken into account. Parameters of the ambient plasma as well as the distribution of electric/magnetic fields are obtained from an MHD simulation [1] of the coma connected to the solar wind. Trajectories of ions and electrically charged dust grains are simulated by accounting for the Lorentz force and the nucleus gravity.

Narrow dust jets in a diffuse gas coma: a natural product of small active regions on comets

Comets often display narrow dust jets but more diffuse gas comae when their eccentric orbits bring them into the inner solar system and sunlight sublimates the ice on the nucleus. Comets are also understood to have one or more active areas covering only a fraction of the total surface active with sublimating volatile ices. Calculations of the gas and dust distribution from a small active area on a comet’s nucleus show that as the gas moves out radially into the vacuum of space it expands tangentially, filling much of the hemisphere centered on the active region. The dust dragged by the gas remains more concentrated over the active area. This explains some puzzling appearances of comets having collimated dust jets but more diffuse gaseous atmospheres. Our test case is 67P/Churyumov–Gerasimenko, the Rosetta mission target comet, whose activity is dominated by a single area covering only 4% of its surface.

The calculation of Afρ and mass loss rate for comets

Ab initio calculations of Afρ are presented using Mie scattering theory and a Direct Simulation Monte Carlo (DSMC) dust outflow model in support of the Rosetta mission and its target 67P/Churyumov-Gerasimenko (CG). These calculations are performed for particle sizes ranging from 0.010 μm to 1.0 cm. The present status of our knowledge of various differential particle size distributions is reviewed and a variety of particle size distributions is used to explore their effect on Afρ , and the dust mass production View the MathML sourcem˙. A new simple two parameter particle size distribution that curtails the effect of particles below 1 μm is developed. The contributions of all particle sizes are summed to get a resulting overall Afρ. The resultant Afρ could not easily be predicted a priori and turned out to be considerably more constraining regarding the mass loss rate than expected. It is found that a proper calculation of Afρ combined with a good Afρ measurement can constrain the dust/gas ratio in the coma of comets as well as other methods presently available. Phase curves of Afρ versus scattering angle are calculated and produce good agreement with observational data. The major conclusions of our calculations are: – The original definition of A in Afρ is problematical and Afρ should be: qsca(n,λ)×p(g)×f×ρqsca(n,λ)×p(g)×f×ρ. Nevertheless, we keep the present nomenclature of Afρ as a measured quantity for an ensemble of coma particles.– The ratio between Afρ and the dust mass loss rate View the MathML sourcem˙ is dominated by the particle size distribution. – For most particle size distributions presently in use, small particles in the range from 0.10 to 1.0 μm contribute a large fraction to Afρ. – Simplifying the calculation of Afρ by considering only large particles and approximating qsca does not represent a realistic model. Mie scattering theory or if necessary, more complex scattering calculations must be used. – For the commonly used particle size distribution, dn/da ∼ a−3.5 to a−4, there is a natural cut off in Afρ contribution for both small and large particles. – The scattering phase function must be taken into account for each particle size; otherwise the contribution of large particles can be over-estimated by a factor of 10. – Using an imaginary index of refraction of i = 0.10 does not produce sufficient backscattering to match observational data. – A mixture of dark particles with i ⩾ 0.10 and brighter silicate particles with i ⩽ 0.04 matches the observed phase curves quite well. – Using current observational constraints, we find the dust/gas mass-production ratio of CG at 1.3 AU is confined to a range of 0.03–0.5 with a reasonably likely value around 0.1.

Effect of the Tiger Stripes on the water vapor distribution in Enceladus' exosphere

The jet activity emanating from Enceladus' exosphere south pole region observed by Cassini is a subject of intensive study. The in situ and remote sensing observations performed since 2005 triggered an active modeling campaign. Such modeling is essential for better understanding of the measurements performed by individual instruments as well as to link them for a more complete picture of the volatile and ice grain distribution in Enceladus' exosphere. This paper is focused on the investigation of the effect that diffuse gas sources along the Tiger Stripes have on distribution of the water vapor in Enceladus' exosphere using the updated version of our multiplume model. We have found that accounting for the gas production by Tiger Stripes is critical for interpretation of the Cassini data. According to our calculations, sources along the Tiger Stripes (apart from those originally identified by Spitale and Porco (2007)) must contribute about 23–32% to the total plume source rate, which varies in the range of (6.4–29) ×1027 s−1. The effect of the previously unidentified source suggested in the paper is found to be critical for explaining the Ultraviolet Imaging Spectrograph 2007 and 2010 observations in the whole range of the elapsed times.

Size-frequency distribution of boulders ≥7 m on comet 67P/Churyumov-Gerasimenko

Aims. We derive for the first time the size-frequency distribution of boulders on a comet, 67P/Churyumov-Gerasimenko (67P), computed from the images taken by the Rosetta/OSIRIS imaging system. We highlight the possible physical processes that lead to these boulder size distributions. Methods. We used images acquired by the OSIRIS Narrow Angle Camera, NAC, on 5 and 6 August 2014. The scale of these images (2.44−2.03 m/px) is such that boulders ≥7 m can be identified and manually extracted from the datasets with the software ArcGIS. We derived both global and localized size-frequency distributions. The three-pixel sampling detection, coupled with the favorable shadowing of the surface (observation phase angle ranging from 48° to 53°), enables unequivocally detecting boulders scattered all over the illuminated side of 67P. Results. We identify 3546 boulders larger than 7 m on the imaged surface (36.4 km2), with a global number density of nearly 100/km2 and a cumulative size-frequency distribution represented by a power-law with index of −3.6 +0.2/−0.3. The two lobes of 67P appear to have slightly different distributions, with an index of −3.5 +0.2/−0.3 for the main lobe (body) and −4.0 +0.3/−0.2 for the small lobe (head). The steeper distribution of the small lobe might be due to a more pervasive fracturing. The difference of the distribution for the connecting region (neck) is much more significant, with an index value of −2.2 +0.2/−0.2. We propose that the boulder field located in the neck area is the result of blocks falling from the contiguous Hathor cliff. The lower slope of the size-frequency distribution we see today in the neck area might be due to the concurrent processes acting on the smallest boulders, such as i) disintegration or fragmentation and vanishing through sublimation; ii) uplifting by gas drag and consequent redistribution; and iii) burial beneath a debris blanket. We also derived the cumulative size-frequency distribution per km2 of localized areas on 67P. By comparing the cumulative size-frequency distributions of similar geomorphological settings, we derived similar power-law index values. This suggests that despite the selected locations on different and often opposite sides of the comet, similar sublimation or activity processes, pit formation or collapses, as well as thermal stresses or fracturing events occurred on multiple areas of the comet, shaping its surface into the appearance we see today.

3D Direct Simulation Monte Carlo Modelling of the Inner Gas Coma of Comet 67P/Churyumov-Gerasimenko: A Parameter Study

Direct Simulation Monte Carlo (DSMC) is a powerful numerical method to study rarefied gas flows such as cometary comae and has been used by several authors over the past decade to study cometary outflow. However, the investigation of the parameter space in simulations can be time consuming since 3D DSMC is computationally highly intensive. For the target of ESA's Rosetta mission, comet 67P/Churyumov-Gerasimenko, we have identified to what extent modification of several parameters influence the 3D flow and gas temperature fields and have attempted to establish the reliability of inferences about the initial conditions from in situ and remote sensing measurements. A large number of DSMC runs have been completed with varying input parameters. In this work, we present the simulation results and conclude on the sensitivity of solutions to certain inputs. It is found that among cases of water outgassing, the surface production rate distribution is the most influential variable to the flow field.

Alteration of natural ³⁷Ar activity concentration in the subsurface by gas transport and water infiltration

High ³⁷Ar activity concentration in soil gas is proposed as a key evidence for the detection of underground nuclear explosion by the Comprehensive Nuclear Test-Ban Treaty. However, such a detection is challenged by the natural background of ³⁷Ar in the subsurface, mainly due to Ca activation by cosmic rays. A better understanding and improved capability to predict ³⁷Ar activity concentration in the subsurface and its spatial and temporal variability is thus required. A numerical model integrating ³⁷Ar production and transport in the subsurface is developed, including variable soil water content and water infiltration at the surface. A parameterized equation for ³⁷Ar production in the first 15 m below the surface is studied, taking into account the major production reactions and the moderation effect of soil water content. Using sensitivity analysis and uncertainty quantification, a realistic and comprehensive probability distribution of natural ³⁷Ar activity concentrations in soil gas is proposed, including the effects of water infiltration. Site location and soil composition are identified as the parameters allowing for a most effective reduction of the possible range of ³⁷Ar activity concentrations. The influence of soil water content on ³⁷Ar production is shown to be negligible to first order, while ³⁷Ar activity concentration in soil gas and its temporal variability appear to be strongly influenced by transient water infiltration events. These results will be used as a basis for practical CTBTO concepts of operation during an OSI.