Ultrasound measurements of the distance between acupuncture needle tip at P6 and the median nerve

BACKGROUND: Pericard 6 (P6) is one of the most frequently used acupuncture points, especially in preventing nausea and vomiting. At this point, the median nerve is located very superficially. OBJECTIVES: To investigate the distance between the needle tip and the median nerve during acupuncture at P6, we conducted a prospective observational ultrasound (US) imaging study. We tested the hypothesis that de qi (a sensation that is typical of acupuncture needling) is evoked when the needle comes into contact with the epineural tissue and thereby prevents nerve penetration. SETTINGS/LOCATION: The outpatient pain clinic of the Medical University of Vienna, Austria. SUBJECTS: Fifty (50) patients receiving acupuncture treatment including P6 bilaterally. INTERVENTIONS: Patients were examined at both forearms using US (a 10-MHz linear transducer) after insertion of the needle at P6. OUTCOME MEASURES: The distance between the needle tip and the median nerve, the number of nerve contacts and nerve penetrations, as well as the number of successfully elicited de qi sensations were recorded. RESULTS: Complete data could be obtained from 97 cases. The mean distance from the needle tip to the nerve was 1.8 mm (standard deviation 2.2; range 0-11.3). Nerve contacts were recorded in 52 cases, in 14 of which the nerve was penetrated by the needle. De qi was elicited in 85 cases. We found no association between the number of nerve contacts and de qi. The 1-week follow-up showed no complications or neurologic problems. CONCLUSIONS: This is the first investigation demonstrating the relationship between acupuncture needle placement and adjacent neural structures using US technology. The rate of median nerve penetrations by the acupuncture needle at P6 was surprisingly high, but these seemed to carry no risk of neurologic sequelae. De qi at P6 does not depend on median nerve contact, nor does it prevent median nerve penetration.

Comparisons between OMI SO2 data and ground-based SO2 measurements at Turrialba volcano

Turrialba is one of the largest and most active stratovolcanoes in the Central Cordillera of Costa Rica and an excellent target for validation of satellite data using ground based measurements due to its high elevation, relative ease of access, and persistent elevated SO2 degassing. The Ozone Monitoring Instrument (OMI) aboard the Aura satellite makes daily global observations of atmospheric trace gases and it is used in this investigation to obtain volcanic SO2 retrievals in the Turrialba volcanic plume. We present and evaluate the relative accuracy of two OMI SO2 data analysis procedures, the automatic Band Residual Index (BRI) technique and the manual Normalized Cloud-mass (NCM) method. We find a linear correlation and good quantitative agreement between SO2 burdens derived from the BRI and NCM techniques, with an improved correlation when wet season data are excluded. We also present the first comparisons between volcanic SO2 emission rates obtained from ground-based mini-DOAS measurements at Turrialba and three new OMI SO2 data analysis techniques: the MODIS smoke estimation, OMI SO2 lifetime, and OMI SO2 transect techniques. A robust validation of OMI SO2 retrievals was made, with both qualitative and quantitative agreements under specific atmospheric conditions, proving the utility of satellite measurements for estimating accurate SO2 emission rates and monitoring passively degassing volcanoes.

Algorithms for spacecraft formation flying navigation based on wireless positioning system measurements

Spacecraft formation flying navigation continues to receive a great deal of interest. The research presented in this dissertation focuses on developing methods for estimating spacecraft absolute and relative positions, assuming measurements of only relative positions using wireless sensors. The implementation of the extended Kalman filter to the spacecraft formation navigation problem results in high estimation errors and instabilities in state estimation at times. This is due tp the high nonlinearities in the system dynamic model. Several approaches are attempted in this dissertation aiming at increasing the estimation stability and improving the estimation accuracy. A differential geometric filter is implemented for spacecraft positions estimation. The differential geometric filter avoids the linearization step (which is always carried out in the extended Kalman filter) through a mathematical transformation that converts the nonlinear system into a linear system. A linear estimator is designed in the linear domain, and then transformed back to the physical domain. This approach demonstrated better estimation stability for spacecraft formation positions estimation, as detailed in this dissertation. The constrained Kalman filter is also implemented for spacecraft formation flying absolute positions estimation. The orbital motion of a spacecraft is characterized by two range extrema (perigee and apogee). At the extremum, the rate of change of a spacecraft’s range vanishes. This motion constraint can be used to improve the position estimation accuracy. The application of the constrained Kalman filter at only two points in the orbit causes filter instability. Two variables are introduced into the constrained Kalman filter to maintain the stability and improve the estimation accuracy. An extended Kalman filter is implemented as a benchmark for comparison with the constrained Kalman filter. Simulation results show that the constrained Kalman filter provides better estimation accuracy as compared with the extended Kalman filter. A Weighted Measurement Fusion Kalman Filter (WMFKF) is proposed in this dissertation. In wireless localizing sensors, a measurement error is proportional to the distance of the signal travels and sensor noise. In this proposed Weighted Measurement Fusion Kalman Filter, the signal traveling time delay is not modeled; however, each measurement is weighted based on the measured signal travel distance. The obtained estimation performance is compared to the standard Kalman filter in two scenarios. The first scenario assumes using a wireless local positioning system in a GPS denied environment. The second scenario assumes the availability of both the wireless local positioning system and GPS measurements. The simulation results show that the WMFKF has similar accuracy performance as the standard Kalman Filter (KF) in the GPS denied environment. However, the WMFKF maintains the position estimation error within its expected error boundary when the WLPS detection range limit is above 30km. In addition, the WMFKF has a better accuracy and stability performance when GPS is available. Also, the computational cost analysis shows that the WMFKF has less computational cost than the standard KF, and the WMFKF has higher ellipsoid error probable percentage than the standard Measurement Fusion method. A method to determine the relative attitudes between three spacecraft is developed. The method requires four direction measurements between the three spacecraft. The simulation results and covariance analysis show that the method’s error falls within a three sigma boundary without exhibiting any singularity issues. A study of the accuracy of the proposed method with respect to the shape of the spacecraft formation is also presented.

Comparison of arterial blood pressure measurements and hypertension scores obtained by use of three indirect measurement devices in hospitalized dogs

OBJECTIVE: To evaluate the agreement of blood pressure measurements and hypertension scores obtained by use of 3 indirect arterial blood pressure measurement devices in hospitalized dogs. Design-Diagnostic test evaluation. ANIMALS: 29 client-owned dogs. PROCEDURES: 5 to 7 consecutive blood pressure readings were obtained from each dog on each of 3 occasions with a Doppler ultrasonic flow detector, a standard oscillometric device (STO), and a high-definition oscillometric device (HDO). RESULTS: When the individual sets of 5 to 7 readings were evaluated, the coefficient of variation for systolic arterial blood pressure (SAP) exceeded 20% for 0% (Doppler), 11 % (STO), and 28% (HDO) of the sets of readings. After readings that exceeded a 20% coefficient of variation were discarded, repeatability was within 25 (Doppler), 37 (STO), and 39 (HDO) mm Hg for SAP. Correlation of mean values among the devices was between 0.47 and 0.63. Compared with Doppler readings, STO underestimated and HDO overestimated SAP. Limits of agreement between mean readings of any 2 devices were wide. With the hypertension scale used to score SAP, the intraclass correlation of scores was 0.48. Linear-weighted inter-rater reliability between scores was 0.40 (Doppler vs STO), 0.38 (Doppler vs HDO), and 0.29 (STO vs HDO). CONCLUSIONS AND CLINICAL RELEVANCE: Results of this study suggested that no meaningful clinical comparison can be made between blood pressure readings obtained from the same dog with different indirect blood pressure measurement devices.

Bone mineral density (BMD) and computer tomographic measurements of the equine proximal phalanx in correlation with breaking strength

Despite the fact that bone mineral density (BMD) is an important fracture risk predictor in human medicine, studies in equine orthopedic research are still lacking. We hypothesized that BMD correlates with bone failure and fatigue fractures of this bone. Thus, the objectives of this study were to measure the structural and mechanical properties of the proximal phalanx with dual energy X-ray absorptiometry (DXA), to correlate the data obtained from DXA and computer tomography (CT) measurements to those obtained by loading pressure examination and to establish representative region of interest (ROI) for in vitro BMD measurements of the equine proximal phalanx for predicting bone failure force. DXA was used to measure the whole bone BMD and additional three ROI sites in 14 equine proximal phalanges. Following evaluation of the bone density, whole bone, cortical width and area in the mid-diaphyseal plane were measured on CT images. Bones were broken using a manually controlled universal bone crusher to measure bone failure force and reevaluated for the site of fractures on follow-up CT images. Compressive load was applied at a constant displacement rate of 2 mm/min until failure, defined as the first clear drop in the load measurement. The lowest BMD was measured at the trabecular region (mean +/- SD: 1.52 +/- 0.12 g/cm2; median: 1.48 g/cm2; range: 1.38-1.83 g/cm2). There was a significant positive linear correlation between trabelcular BMD and the breaking strength (P = 0.023, r = 0.62). The trabecular region of the proximal phalanx appears to be the only significant indicator of failure of strength in vitro. This finding should be reassessed to further reveal the prognostic value of trabecular BMD in an in vivo fracture risk model.

Quantification of pleural effusion from single area measurements on CT

The objective of this study was to determine if area measurements of pleural fluid on computed tomography (CT) reflect the actual pleural fluid volume (PEvol) as measured at autopsy, to establish a formula to estimate the volume of pleural effusion (PEest), and to test the accuracy and observer reliability of PEest.132 human cadavers, with pleural effusion were divided into phase 1 (n = 32) and phase 2 (n = 100). In phase 1, PEvol was compared to area measurements on axial (axA), sagittal (sagA), and coronal (corA) CT images. Linear regression analysis was used to create a formula to calculate PEest. In phase 2, intra-class correlation (ICC) was used to assess inter-reader reliability and determine the agreement between PEest and PEvol. PEvol correlated to a higher degree to axA (r s mean = 0.738; p < 0.001) than to sagA (r s mean = 0.679, p < 0.001) and corA (r s mean = 0.709; p < 0.001). PEest can be established with the following formula: axA × 0.1 = PEest. Mean difference between PEest and PEvol was less than 40 mL (ICC = 0.837-0.874; p < 0.001). Inter-reader reliability was higher between two experienced readers (ICC = 0.984-0.987; p < 0.001) than between an inexperienced reader and both experienced readers (ICC = 0.660-0.698; p < 0.001). Pleural effusions may be quantified in a rapid, reliable, and reasonably accurate fashion using single area measurements on CT.

Relating Crevassing to Non-Linear Strain in the Floating Part of Jakobshavn Isbrae, West Greenland

Jakobshavn Isbrae is a major ice stream that drains the west-central Greenland ice sheet and becomes afloat in Jakobshavn Isfiord (69degreesN, 49degreesW), where it has maintained the world's fastest-known sustained velocity and calving rate (7 km a(-1)) for at least four decades. The floating portion is approximately 12 km long and 6 km wide. Surface elevations and motion vectors were determined photogrammetrically for about 500 crevasses on the floating ice, and adjacent grounded ice, using aerial photographs obtained 2 weeks apart in July 1985. Surface strain rates were computed from a mesh of 399 quadrilateral elements having velocity measurements at each corner. It is shown that heavy crevassing of floating ice invalidates the assumptions of linear strain theory that (i) surface strain in the floating ice is homogeneous in both space and time, (ii) the squares and products of strain components are nil, and (iii) first- and second-order rotation components are small compared to strain components. Therefore, strain rates and rotation rates were also computed using non-linear strain theory. The percentage difference between computed linear and non-linear second invariants of strain rate per element were greatest (mostly in the range 40-70%) where crevassing is greatest. Isopleths of strain rate parallel and transverse to flow and elevation isopleths relate crevassing to known and inferred pinning points.

Relating in situ gas measurements to the surface outgassing properties of cometary nuclei

The sensitivity of the gas flow field to changes in different initial conditions has been studied for the case of a highly simplified cometary nucleus model. The nucleus model simulated a homogeneously outgassing sphere with a more active ring around an axis of symmetry. The varied initial conditions were the number density of the homogeneous region, the surface temperature, and the composition of the flow (varying amounts of H2O and CO2) from the active ring. The sensitivity analysis was performed using the Polynomial Chaos Expansion (PCE) method. Direct Simulation Monte Carlo (DSMC) was used for the flow, thereby allowing strong deviations from local thermal equilibrium. The PCE approach can be used to produce a sensitivity analysis with only four runs per modified input parameter and allows one to study and quantify non-linear responses of measurable parameters to linear changes in the input over a wide range. Hence the PCE allows one to obtain a functional relationship between the flow field properties at every point in the inner coma and the input conditions. It is for example shown that the velocity and the temperature of the background gas are not simply linear functions of the initial number density at the source. As probably expected, the main influence on the resulting flow field parameter is the corresponding initial parameter (i.e. the initial number density determines the background number density, the temperature of the surface determines the flow field temperature, etc.). However, the velocity of the flow field is also influenced by the surface temperature while the number density is not sensitive to the surface temperature at all in our model set-up. Another example is the change in the composition of the flow over the active area. Such changes can be seen in the velocity but again not in the number density. Although this study uses only a simple test case, we suggest that the approach, when applied to a real case in 3D, should assist in identifying the sensitivity of gas parameters measured in situ by, for example, the Rosetta spacecraft to the surface boundary conditions and vice versa.

Hygroscopic growth of atmospheric aerosol particles based on active remote sensing and radiosounding measurements: selected cases in southeastern Spain

A new methodology based on combining active and passive remote sensing and simultaneous and collocated radiosounding data to study the aerosol hygroscopic growth effects on the particle optical and microphysical properties is presented. The identification of hygroscopic growth situations combines the analysis of multispectral aerosol particle backscatter coefficient and particle linear depolarization ratio with thermodynamic profiling of the atmospheric column. We analyzed the hygroscopic growth effects on aerosol properties, namely the aerosol particle backscatter coefficient and the volume concentration profiles, using data gathered at Granada EARLINET station. Two study cases, corresponding to different aerosol loads and different aerosol types, are used for illustrating the potential of this methodology. Values of the aerosol particle backscatter coefficient enhancement factors range from 2.1 ± 0.8 to 3.9 ± 1.5, in the ranges of relative humidity 60–90 and 40–83%, being similar to those previously reported in the literature. Differences in the enhancement factor are directly linked to the composition of the atmospheric aerosol. The largest value of the aerosol particle backscatter coefficient enhancement factor corresponds to the presence of sulphate and marine particles that are more affected by hygroscopic growth. On the contrary, the lowest value of the enhancement factor corresponds to an aerosol mixture containing sulphates and slight traces of mineral dust. The Hänel parameterization is applied to these case studies, obtaining results within the range of values reported in previous studies, with values of the γ exponent of 0.56 ± 0.01 (for anthropogenic particles slightly influenced by mineral dust) and 1.07 ± 0.01 (for the situation dominated by anthropogenic particles), showing the convenience of this remote sensing approach for the study of hygroscopic effects of the atmospheric aerosol under ambient unperturbed conditions. For the first time, the retrieval of the volume concentration profiles for these cases using the Lidar Radiometer Inversion Code (LIRIC) allows us to analyze the aerosol hygroscopic growth effects on aerosol volume concentration, observing a stronger increase of the fine mode volume concentration with increasing relative humidity.

Evaluation of Polymer Gel Dosimeters for Measurements of Dose and LET in Proton Beams

This project assessed the effectiveness of polymer gel dosimeters as tools for measuring the dose deposited by and LET of a proton beam. A total of three BANG® dosimeter formulations were evaluated: BANG®-3-Pro-2 BANGkits™ for dose measurement and two BANG®-3 variants, the LET-Baseline and LET-Meter dosimeters, for LET measurement. All dosimeters were read out using an OCT scanner. The basic characteristics of the BANGkits™ were assessed in a series of photon and electron irradiations. The dose-response relationship was found to be sigmoidal with a threshold for response of approximately 15 cGy. The active region of the dosimeter, the volume in which dosimeter response is not inhibited by oxygen, was found to make up roughly one fourth of the total dosimeter volume. Delivering a dose across multiple fractions was found to yield a greater response than delivering the same dose in a single irradiation. The dosimeter was found to accurately measure a dose distribution produced by overlapping photon fields, yielding gamma pass rates of 95.4% and 93.1% from two planar gamma analyses. Proton irradiations were performed for measurements of proton dose and LET. Initial irradiations performed through the side of a dosimeter led to OCT artifacts. Gamma pass rates of 85.7% and 89.9% were observed in two planar gamma analyses. In irradiations performed through the base of a dosimeter, gel response was found to increase with height in the dosimeter, even in areas of constant dose. After a correction was applied, gamma pass rates of 94.6% and 99.3% were observed in two planar gamma analyses. Absolute dose measurements were substantially higher (33%-100%) than the delivered doses for proton irradiations. Issues encountered while calibrating the LET-Meter gel restricted analysis of the LET measurement data to the SOBP of a proton beam. LET-Meter overresponse was found to increase linearly with track-average LET across the LET range that could be investigated (1.5 keV/micron – 3.5 keV/micron).

Application of the general linear model to assess the effect of missing data on bone marrow mononuclear cell therapy with infusion timing and follow-up MRI

With most clinical trials, missing data presents a statistical problem in evaluating a treatment's efficacy. There are many methods commonly used to assess missing data; however, these methods leave room for bias to enter the study. This thesis was a secondary analysis on data taken from TIME, a phase 2 randomized clinical trial conducted to evaluate the safety and effect of the administration timing of bone marrow mononuclear cells (BMMNC) for subjects with acute myocardial infarction (AMI).^ We evaluated the effect of missing data by comparing the variance inflation factor (VIF) of the effect of therapy between all subjects and only subjects with complete data. Through the general linear model, an unbiased solution was made for the VIF of the treatment's efficacy using the weighted least squares method to incorporate missing data. Two groups were identified from the TIME data: 1) all subjects and 2) subjects with complete data (baseline and follow-up measurements). After the general solution was found for the VIF, it was migrated Excel 2010 to evaluate data from TIME. The resulting numerical value from the two groups was compared to assess the effect of missing data.^ The VIF values from the TIME study were considerably less in the group with missing data. By design, we varied the correlation factor in order to evaluate the VIFs of both groups. As the correlation factor increased, the VIF values increased at a faster rate in the group with only complete data. Furthermore, while varying the correlation factor, the number of subjects with missing data was also varied to see how missing data affects the VIF. When subjects with only baseline data was increased, we saw a significant rate increase in VIF values in the group with only complete data while the group with missing data saw a steady and consistent increase in the VIF. The same was seen when we varied the group with follow-up only data. This essentially showed that the VIFs steadily increased when missing data is not ignored. When missing data is ignored as with our comparison group, the VIF values sharply increase as correlation increases.^

Thorium and radiocarbon measurements on sediment cores from the Northseast Atlantic

Radiocarbon and 230Thexcess data from six NE Atlantic box cores are considered. The cores form a transect from the Porcupine Abyssal Plain over the East Thulean Rise to the southern end of Feni Drift. The chronology for the cores is established from bulk sediment carbonate radiocarbon data and reveals that sections exhibiting constant accumulation rates can be identified in all the cores, with rates of 3.0-3.5 cm/kyr on the plain through the Holocene and late Holocene rates of 4.3-6.6 cm/kyr elsewhere. Five out of the six cores show accumulations of more 230Thexcess than is produced in the overlying water column, with the greatest inventories (up to 225% of production) in the cores from the rise and drift. A size fraction comparison between two cores from the plain and rise reveals that the higher overall accumulation rates and 230Thexcess inventories in the off-plain cores are due to an increased fine (<5 µm) component fraction, whereas the flux of coarser material is similar to that received on the plain. This suggests that the higher fluxes of materials observed are physically (rather than biogeochemically) driven and also that drift formation has been continuously active in the late Holocene. Sections of all the cores where regular accumulation is defined by the radiocarbon data are modeled first by a linear radiocarbon age/depth model and second by a constant rain (230Thexcess)0 model prorated for the observed core inventories. These modeling approaches yield historical mass accumulation rate estimates which are generally in reasonable agreement (±30%), but the differences observed appear to be well organized in time rather than random.

Inorganic carbon and geochemistry measurements from landfast sea ice in Resolute Passage, Nunavut, Canada, as part of the Arctic-ICE project, May - June 2010

We conducted a six-week investigation of the sea ice inorganic carbon system during the winter-spring transition in the Canadian Arctic Archipelago. Samples for the determination of sea ice geochemistry were collected in conjunction with physical and biological parameters as part of the 2010 Arctic-ICE (Arctic - Ice-Covered Ecosystem in a Rapidly Changing Environment) program, a sea ice-based process study in Resolute Passage, Nunavut. The goal of Arctic-ICE was to determine the physical-biological processes controlling the timing of primary production in Arctic landfast sea ice and to better understand the influence of these processes on the drawdown and release of climatically active gases. The field study was conducted from 1 May to 21 June, 2010.

Snow-pit measurements on Atka Bay landfast sea ice during ANT-Land 2012/2013 season

Up to now, snow cover on Antarctic sea ice and its impact on radar backscatter, particularly after the onset of freeze/thaw processes, are not well understood. Here we present a combined analysis of in situ observations of snow properties from the landfast sea ice in Atka Bay, Antarctica, and high-resolution TerraSAR-X backscatter data, for the transition from austral spring (November 2012) to summer (January 2013). The physical changes in the seasonal snow cover during that time are reflected in the evolution of TerraSAR-X backscatter. We are able to explain 76-93% of the spatio-temporal variability of the TerraSAR-X backscatter signal with up to four snowpack parameters with a root-mean-squared error of 0.87-1.62 dB, using a simple multiple linear model. Over the complete study, and especially after the onset of early-melt processes and freeze/thaw cycles, the majority of variability in the backscatter is influenced by changes in snow/ice interface temperature, snow depth and top-layer grain size. This suggests it may be possible to retrieve snow physical properties over Antarctic sea ice from X-band SAR backscatter.

Measurements of gastropods shell sizes from 23 fossil and extant long-lived lakes

Aim: To investigate shell size variation among gastropod faunas of fossil and recent long-lived European lakes and discuss potential underlying processes. Location: 23 long-lived lakes of the Miocene to Recent of Europe. Methods: Based on a dataset of 1412 species of both fossil and extant lacustrine gastropods, we assessed differences in shell size in terms of characteristics of the faunas (species richness, degree of endemism, differences in family composition) and the lakes (surface area, latitude and longitude of lake centroid, distance to closest neighbouring lake) using multiple and linear regression models. Because of a strong species-area relationship, we used resampling to determine whether any observed correlation is driven by that relationship. Results: The regression models indicated size range expansion rather than unidirectional increase or decrease as the dominant pattern of size evolution. The multiple regression models for size range and maximum and minimum size were statistically significant, while the model with mean size was not. Individual contributions and linear regressions indicated species richness and lake surface area as best predictors for size changes. Resampling analysis revealed no significant effects of species richness on the observed patterns. The correlations are comparable across families of different size classes, suggesting a general pattern. Main conclusions: Among the chosen variables, species richness and lake surface area are the most robust predictors of shell size in long-lived lake gastropods. Although the most outstanding and attractive examples for size evolution in lacustrine gastropods derive from lakes with extensive durations, shell size appears to be independent of the duration of the lake as well as longevity of a species. The analogue of long-lived lakes as 'evolutionary islands' does not hold for developments of shell size because different sets of parameters predict size changes.