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.

MOREMATA: Stata module (Mata) to provide various functions

This package includes various Mata functions. kern(): various kernel functions; kint(): kernel integral functions; kdel0(): canonical bandwidth of kernel; quantile(): quantile function; median(): median; iqrange(): inter-quartile range; ecdf(): cumulative distribution function; relrank(): grade transformation; ranks(): ranks/cumulative frequencies; freq(): compute frequency counts; histogram(): produce histogram data; mgof(): multinomial goodness-of-fit tests; collapse(): summary statistics by subgroups; _collapse(): summary statistics by subgroups; gini(): Gini coefficient; sample(): draw random sample; srswr(): SRS with replacement; srswor(): SRS without replacement; upswr(): UPS with replacement; upswor(): UPS without replacement; bs(): bootstrap estimation; bs2(): bootstrap estimation; bs_report(): report bootstrap results; jk(): jackknife estimation; jk_report(): report jackknife results; subset(): obtain subsets, one at a time; composition(): obtain compositions, one by one; ncompositions(): determine number of compositions; partition(): obtain partitions, one at a time; npartitionss(): determine number of partitions; rsubset(): draw random subset; rcomposition(): draw random composition; colvar(): variance, by column; meancolvar(): mean and variance, by column; variance0(): population variance; meanvariance0(): mean and population variance; mse(): mean squared error; colmse(): mean squared error, by column; sse(): sum of squared errors; colsse(): sum of squared errors, by column; benford(): Benford distribution; cauchy(): cumulative Cauchy-Lorentz dist.; cauchyden(): Cauchy-Lorentz density; cauchytail(): reverse cumulative Cauchy-Lorentz; invcauchy(): inverse cumulative Cauchy-Lorentz; rbinomial(): generate binomial random numbers; cebinomial(): cond. expect. of binomial r.v.; root(): Brent's univariate zero finder; nrroot(): Newton-Raphson zero finder; finvert(): univariate function inverter; integrate_sr(): univariate function integration (Simpson's rule); integrate_38(): univariate function integration (Simpson's 3/8 rule); ipolate(): linear interpolation; polint(): polynomial inter-/extrapolation; plot(): Draw twoway plot; _plot(): Draw twoway plot; panels(): identify nested panel structure; _panels(): identify panel sizes; npanels(): identify number of panels; nunique(): count number of distinct values; nuniqrows(): count number of unique rows; isconstant(): whether matrix is constant; nobs(): number of observations; colrunsum(): running sum of each column; linbin(): linear binning; fastlinbin(): fast linear binning; exactbin(): exact binning; makegrid(): equally spaced grid points; cut(): categorize data vector; posof(): find element in vector; which(): positions of nonzero elements; locate(): search an ordered vector; hunt(): consecutive search; cond(): matrix conditional operator; expand(): duplicate single rows/columns; _expand(): duplicate rows/columns in place; repeat(): duplicate contents as a whole; _repeat(): duplicate contents in place; unorder2(): stable version of unorder(); jumble2(): stable version of jumble(); _jumble2(): stable version of _jumble(); pieces(): break string into pieces; npieces(): count number of pieces; _npieces(): count number of pieces; invtokens(): reverse of tokens(); realofstr(): convert string into real; strexpand(): expand string argument; matlist(): display a (real) matrix; insheet(): read spreadsheet file; infile(): read free-format file; outsheet(): write spreadsheet file; callf(): pass optional args to function; callf_setup(): setup for mm_callf().

Single fluoroscopic image based 2D/3D reconstruction of a scaled, patient-specific vertebral model

Purpose Accurate three-dimensional (3D) models of lumbar vertebrae can enable image-based 3D kinematic analysis. The common approach to derive 3D models is by direct segmentation of CT or MRI datasets. However, these have the disadvantages that they are expensive, timeconsuming and/or induce high-radiation doses to the patient. In this study, we present a technique to automatically reconstruct a scaled 3D lumbar vertebral model from a single two-dimensional (2D) lateral fluoroscopic image. Methods Our technique is based on a hybrid 2D/3D deformable registration strategy combining a landmark-to-ray registration with a statistical shape model-based 2D/3D reconstruction scheme. Fig. 1 shows different stages of the reconstruction process. Four cadaveric lumbar spine segments (total twelve lumbar vertebrae) were used to validate the technique. To evaluate the reconstruction accuracy, the surface models reconstructed from the lateral fluoroscopic images were compared to the associated ground truth data derived from a 3D CT-scan reconstruction technique. For each case, a surface-based matching was first used to recover the scale and the rigid transformation between the reconstructed surface model Results Our technique could successfully reconstruct 3D surface models of all twelve vertebrae. After recovering the scale and the rigid transformation between the reconstructed surface models and the ground truth models, the average error of the 2D/3D surface model reconstruction over the twelve lumbar vertebrae was found to be 1.0 mm. The errors of reconstructing surface models of all twelve vertebrae are shown in Fig. 2. It was found that the mean errors of the reconstructed surface models in comparison to their associated ground truths after iterative scaled rigid registrations ranged from 0.7 mm to 1.3 mm and the rootmean squared (RMS) errors ranged from 1.0 mm to 1.7 mm. The average mean reconstruction error was found to be 1.0 mm. Conclusion An accurate, scaled 3D reconstruction of the lumbar vertebra can be obtained from a single lateral fluoroscopic image using a statistical shape model based 2D/3D reconstruction technique. Future work will focus on applying the reconstructed model for 3D kinematic analysis of lumbar vertebrae, an extension of our previously-reported imagebased kinematic analysis. The developed method also has potential applications in surgical planning and navigation.

ZeroG: overground gait and balance training system

A new overground body-weight support system called ZeroG has been developed that allows patients with severe gait impairments to practice gait and balance activities in a safe, controlled manner. The unloading system is capable of providing up to 300 lb of static support and 150 lb of dynamic (or constant force) support using a custom-series elastic actuator. The unloading system is mounted to a driven trolley, which rides along an overhead rail. We evaluated the performance of ZeroG's unloading system, as well as the trolley tracking system, using benchtop and human-subject testing. Average root-mean-square and peak errors in unloading were 2.2 and 7.2 percent, respectively, over the range of forces tested while trolley tracking errors were less than 3 degrees, indicating the system was able to maintain its position above the subject. We believe training with ZeroG will allow patients to practice activities that are critical to achieving functional independence at home and in the community.

Kernel density correlation based non-rigid point set matching for statistical model-based 2D/3D reconstruction

This paper presents a kernel density correlation based nonrigid point set matching method and shows its application in statistical model based 2D/3D reconstruction of a scaled, patient-specific model from an un-calibrated x-ray radiograph. In this method, both the reference point set and the floating point set are first represented using kernel density estimates. A correlation measure between these two kernel density estimates is then optimized to find a displacement field such that the floating point set is moved to the reference point set. Regularizations based on the overall deformation energy and the motion smoothness energy are used to constraint the displacement field for a robust point set matching. Incorporating this non-rigid point set matching method into a statistical model based 2D/3D reconstruction framework, we can reconstruct a scaled, patient-specific model from noisy edge points that are extracted directly from the x-ray radiograph by an edge detector. Our experiment conducted on datasets of two patients and six cadavers demonstrates a mean reconstruction error of 1.9 mm

A Portable Image Overlay Projection Device for Computer-Aided Open Liver Surgery

Image overlay projection is a form of augmented reality that allows surgeons to view underlying anatomical structures directly on the patient surface. It improves intuitiveness of computer-aided surgery by removing the need for sight diversion between the patient and a display screen and has been reported to assist in 3-D understanding of anatomical structures and the identification of target and critical structures. Challenges in the development of image overlay technologies for surgery remain in the projection setup. Calibration, patient registration, view direction, and projection obstruction remain unsolved limitations to image overlay techniques. In this paper, we propose a novel, portable, and handheld-navigated image overlay device based on miniature laser projection technology that allows images of 3-D patient-specific models to be projected directly onto the organ surface intraoperatively without the need for intrusive hardware around the surgical site. The device can be integrated into a navigation system, thereby exploiting existing patient registration and model generation solutions. The position of the device is tracked by the navigation system’s position sensor and used to project geometrically correct images from any position within the workspace of the navigation system. The projector was calibrated using modified camera calibration techniques and images for projection are rendered using a virtual camera defined by the projectors extrinsic parameters. Verification of the device’s projection accuracy concluded a mean projection error of 1.3 mm. Visibility testing of the projection performed on pig liver tissue found the device suitable for the display of anatomical structures on the organ surface. The feasibility of use within the surgical workflow was assessed during open liver surgery. We show that the device could be quickly and unobtrusively deployed within the sterile environment.

Conserving indigenous animal genetic resources as a coping strategy to adapt to climate change: the Azikheli Buffalo in Northern Mountains of Pakistan

There are clear signs that the agro-pastoralists in the Himalayan and Hindu-Kush mountain ranges will have less cropping opportunities due to reduced possibilities for irrigated agriculture as a result of climate change. The importance of extensive livestock production based on well adapted livestock species may once again increase. This calls for a better documentation and understanding of the adaptation capabilities of indigenous breeds considering a changing environment. The current study investigates the adaptive traits of the Azikheli buffalo to mountain environments through calculating mean, standard error and percentages for different variables. Results from this study suggest that the brown coat color, the small body size and the high fertility are adaptive traits of the Azikheli buffalo that may well suit harsh mountainous environment conditions with greater climate variability. Local farmers find it hard to sustain the Azikheli buffalo’s key adaptive traits because of a low bull to buffalo ratio, possibility of insemination with semen from imported breeds and a lack of institutional support to conserve the Azikheli breed. The breed is crucial for sustaining custodian communities in these mountains and thus needs to be conserved.

Physical activity intensity and surrogate markers for cardiovascular health in adolescents

We examined the impact of physical activity (PA) on surrogate markers of cardiovascular health in adolescents. 52 healthy students (28 females, mean age 14.5 ± 0.7 years) were investigated. Microvascular endothelial function was assessed by peripheral arterial tonometry to determine reactive hyperemic index (RHI). Vagal activity was measured using 24 h analysis of heart rate variability [root mean square of successive normal-to-normal intervals (rMSSD)]. Exercise testing was performed to determine peak oxygen uptake ([Formula: see text]) and maximum power output. PA was assessed by accelerometry. Linear regression models were performed and adjusted for age, sex, skinfolds, and pubertal status. The cohort was dichotomized into two equally sized activity groups (low vs. high) based on the daily time spent in moderate-to-vigorous PA (MVPA, 3,000-5,200 counts(.)min(-1), model 1) and vigorous PA (VPA, >5,200 counts(.)min(-1), model 2). MVPA was an independent predictor for rMSSD (β = 0.448, P = 0.010), and VPA was associated with maximum power output (β = 0.248, P = 0.016). In model 1, the high MVPA group exhibited a higher vagal tone (rMSSD 49.2 ± 13.6 vs. 38.1 ± 11.7 ms, P = 0.006) and a lower systolic blood pressure (107.3 ± 9.9 vs. 112.9 ± 8.1 mmHg, P = 0.046). In model 2, the high VPA group had higher maximum power output values (3.9 ± 0.5 vs. 3.4 ± 0.5 W kg(-1), P = 0.012). In both models, no significant differences were observed for RHI and [Formula: see text]. In conclusion, in healthy adolescents, PA was associated with beneficial intensity-dependent effects on vagal tone, systolic blood pressure, and exercise capacity, but not on microvascular endothelial function.

A realistic validation study of a new nitrogen multiple-breath washout system

Background For reliable assessment of ventilation inhomogeneity, multiple-breath washout (MBW) systems should be realistically validated. We describe a new lung model for in vitro validation under physiological conditions and the assessment of a new nitrogen (N2)MBW system. Methods The N2MBW setup indirectly measures the N2 fraction (FN2) from main-stream carbon dioxide (CO2) and side-stream oxygen (O2) signals: FN2 = 1−FO2−FCO2−FArgon. For in vitro N2MBW, a double chamber plastic lung model was filled with water, heated to 37°C, and ventilated at various lung volumes, respiratory rates, and FCO2. In vivo N2MBW was undertaken in triplets on two occasions in 30 healthy adults. Primary N2MBW outcome was functional residual capacity (FRC). We assessed in vitro error (√[difference]2) between measured and model FRC (100–4174 mL), and error between tests of in vivo FRC, lung clearance index (LCI), and normalized phase III slope indices (Sacin and Scond). Results The model generated 145 FRCs under BTPS conditions and various breathing patterns. Mean (SD) error was 2.3 (1.7)%. In 500 to 4174 mL FRCs, 121 (98%) of FRCs were within 5%. In 100 to 400 mL FRCs, the error was better than 7%. In vivo FRC error between tests was 10.1 (8.2)%. LCI was the most reproducible ventilation inhomogeneity index. Conclusion The lung model generates lung volumes under the conditions encountered during clinical MBW testing and enables realistic validation of MBW systems. The new N2MBW system reliably measures lung volumes and delivers reproducible LCI values.

Reconstruction of patient-specific 3D bone surface from 2D calibrated fluoroscopic images and point distribution model

Reconstruction of patient-specific 3D bone surface from 2D calibrated fluoroscopic images and a point distribution model is discussed. We present a 2D/3D reconstruction scheme combining statistical extrapolation and regularized shape deformation with an iterative image-to-model correspondence establishing algorithm, and show its application to reconstruct the surface of proximal femur. The image-to-model correspondence is established using a non-rigid 2D point matching process, which iteratively uses a symmetric injective nearest-neighbor mapping operator and 2D thin-plate splines based deformation to find a fraction of best matched 2D point pairs between features detected from the fluoroscopic images and those extracted from the 3D model. The obtained 2D point pairs are then used to set up a set of 3D point pairs such that we turn a 2D/3D reconstruction problem to a 3D/3D one. We designed and conducted experiments on 11 cadaveric femurs to validate the present reconstruction scheme. An average mean reconstruction error of 1.2 mm was found when two fluoroscopic images were used for each bone. It decreased to 1.0 mm when three fluoroscopic images were used.

Oral misoprostol for third stage of labor: a randomized placebo-controlled trial

OBJECTIVE: To investigate whether orally administered misoprostol during the third stage of labor is efficient in reducing postpartum blood loss. METHODS: In a double-masked trial, during vaginal delivery women were randomly assigned to receive a single oral dose of misoprostol (600 microg) or placebo in third stage of labor, immediately after cord clamping. The third stage of labor was managed routinely by early cord clamping and controlled cord traction; oxytocin was administered only if blood loss seemed more than usual. Blood loss was estimated by the delivering physician and differences in hematocrit were measured before and after delivery. RESULTS: Mean (+/- standard error of the mean) estimated blood loss (345 +/- 19.5 mL versus 417 +/- 25.9 mL, P = .031) and hematocrit difference (4.5 +/- 0.9% versus 7.9 +/- 1.2%, P = .014) were significantly lower in women who received misoprostol than those who received placebo. Fewer women in the misoprostol group had postpartum hemorrhage (blood loss of at least 500 mL), but that difference was not statistically significant (7% versus 15%, P = .43). Additional oxytocin before or after placental separation was used less often in the misoprostol group (16% versus 38%, P = .047). There were no differences in the length of third stage of labor (8 +/- 0.9 minutes versus 9 +/- 1 minutes, P = .947). There were no differences in pain during third stage of labor, postpartum fever, or diarrhea, but shivering was more frequent in the misoprostol group. CONCLUSION: Oral misoprostol administered in the third stage of labor reduced postpartum blood loss and might be effective in reducing incidence of postpartum hemorrhage.

Comparison of Doppler-derived aortic velocities obtained from various transducer sites in healthy dogs and cats

In normal dogs and dogs with subaortic stenosis, it is known that the subcostal transducer site provides higher left ventricular ejection velocities than does the left apical site. We hypothesized that aortic flow velocities could also be obtained from the right parasternal long-axis view, optimized for the placement of the Doppler cursor as parallel as possible into the aortic root. In 15 healthy dogs and 13 healthy cats, high-pulsed repetition frequency Doppler flow velocity measurements in the proximal aorta were performed using two-dimensional echocardiographic guidance. The mean [ +/- standard error of the mean (SEM)] peak aortic flow velocities in healthy dogs were as follows: subcostal site 1.46 +/- 0.05 m/s; apical site 1.12 +/- 0.06 m/s; right parasternal long-axis site 1.09 +/- 0.05 m/s. In healthy cats, the following peak aortic flow velocities were observed: apical site 0.87 +/- 0.03m/s; right parasternal long-axis site 0.87 +/- 0.03 m/s. Aortic flow velocities obtained from the subcostal site were significantly higher in healthy dogs than those obtained from the left apical and right parasternal long-axis site (P< 0.001). There was no statistical difference between the peak aortic flow velocities obtained from right parasternal long-axis and left apical transducer position in all groups. We conclude therefore that right parasternal long-axis and left apical-derived aortic flow velocities are similar and may be used interchangeably in healthy dogs and cats.

Methodical study of nitrous oxide eddy covariance measurements using quantum cascade laser spectrometery over a Swiss forest

Nitrous oxide fluxes were measured at the Lägeren CarboEurope IP flux site over the multi-species mixed forest dominated by European beech and Norway spruce. Measurements were carried out during a four-week period in October–November 2005 during leaf senescence. Fluxes were measured with a standard ultrasonic anemometer in combination with a quantum cascade laser absorption spectrometer that measured N2O, CO2, and H2O mixing ratios simultaneously at 5 Hz time resolution. To distinguish insignificant fluxes from significant ones it is proposed to use a new approach based on the significance of the correlation coefficient between vertical wind speed and mixing ratio fluctuations. This procedure eliminated roughly 56% of our half-hourly fluxes. Based on the remaining, quality checked N2O fluxes we quantified the mean efflux at 0.8±0.4 μmol m−2 h−1 (mean ± standard error). Most of the contribution to the N2O flux occurred during a 6.5-h period starting 4.5 h before each precipitation event. No relation with precipitation amount could be found. Visibility data representing fog density and duration at the site indicate that wetting of the canopy may have as strong an effect on N2O effluxes as does below-ground microbial activity. It is speculated that above-ground N2O production from the senescing leaves at high moisture (fog, drizzle, onset of precipitation event) may be responsible for part of the measured flux.

High-resolution morphological and biochemical imaging of articular cartilage of the ankle joint at 3.0 T using a new dedicated phased array coil: in vivo reproducibility study

OBJECTIVE: The objective of this study was to evaluate the feasibility and reproducibility of high-resolution magnetic resonance imaging (MRI) and quantitative T2 mapping of the talocrural cartilage within a clinically applicable scan time using a new dedicated ankle coil and high-field MRI. MATERIALS AND METHODS: Ten healthy volunteers (mean age 32.4 years) underwent MRI of the ankle. As morphological sequences, proton density fat-suppressed turbo spin echo (PD-FS-TSE), as a reference, was compared with 3D true fast imaging with steady-state precession (TrueFISP). Furthermore, biochemical quantitative T2 imaging was prepared using a multi-echo spin-echo T2 approach. Data analysis was performed three times each by three different observers on sagittal slices, planned on the isotropic 3D-TrueFISP; as a morphological parameter, cartilage thickness was assessed and for T2 relaxation times, region-of-interest (ROI) evaluation was done. Reproducibility was determined as a coefficient of variation (CV) for each volunteer; averaged as root mean square (RMSA) given as a percentage; statistical evaluation was done using analysis of variance. RESULTS: Cartilage thickness of the talocrural joint showed significantly higher values for the 3D-TrueFISP (ranging from 1.07 to 1.14 mm) compared with the PD-FS-TSE (ranging from 0.74 to 0.99 mm); however, both morphological sequences showed comparable good results with RMSA of 7.1 to 8.5%. Regarding quantitative T2 mapping, measurements showed T2 relaxation times of about 54 ms with an excellent reproducibility (RMSA) ranging from 3.2 to 4.7%. CONCLUSION: In our study the assessment of cartilage thickness and T2 relaxation times could be performed with high reproducibility in a clinically realizable scan time, demonstrating new possibilities for further investigations into patient groups.

Coagulation activity before and after acute psychosocial stress increases with age

OBJECTIVE: To assess whether stress further increases hypercoagulation in older individuals. We investigated whether acute stress-induced changes in coagulation parameters differ with age. It is known that hypercoagulation occurs in response to acute stress and that a shift in hemostasis toward a hypercoagulability state occurs with age. However, it is not yet known whether acute stress further increases hypercoagulation in older individuals, and thus may increase their risk for cardiovascular disease (CVD). METHODS: A total of 63 medication-free nonsmoking men, aged between 20 and 65 years (mean +/- standard error of the mean = 36.7 +/- 1.7 years), underwent an acute standardized psychosocial stress task combining public speaking and mental arithmetic in front of an audience. We measured plasma clotting factor VII activity (FVII:C), fibrinogen, and D-dimer at rest, immediately, and 20 minutes after stress. RESULTS: Increased age predicted greater increases in fibrinogen (beta = 0.26, p = 0.041; DeltaR(2) = 0.05), FVII:C (beta = 0.40, p = .006; DeltaR(2) = 0.11), and D-dimer (beta = 0.51, p < .001; DeltaR(2) = 0.18) from rest to 20 minutes after stress independent of body mass index and mean arterial blood pressure. General linear models revealed significant effects of age and stress on fibrinogen, FVII:C, and D-dimer (main effects: p < .04), and greater D-dimer stress reactivity with older age (interaction age-by-stress: F(1.5/90.4) = 4.36, p = .024; f = 0.33). CONCLUSIONS: Our results suggest that acute stress might increase vulnerability in the elderly for hypercoagulability and subsequent hemostasis-associated diseases like CVD.