Cardiac motion estimation by joint alignment of tagged MRI sequences

Image registration has been proposed as an automatic method for recovering cardiac displacement fields from Tagged Magnetic Resonance Imaging (tMRI) sequences. Initially performed as a set of pairwise registrations, these techniques have evolved to the use of 3D+t deformation models, requiring metrics of joint image alignment (JA). However, only linear combinations of cost functions defined with respect to the first frame have been used. In this paper, we have applied k-Nearest Neighbors Graphs (kNNG) estimators of the -entropy (H ) to measure the joint similarity between frames, and to combine the information provided by different cardiac views in an unified metric. Experiments performed on six subjects showed a significantly higher accuracy (p < 0.05) with respect to a standard pairwise alignment (PA) approach in terms of mean positional error and variance with respect to manually placed landmarks. The developed method was used to study strains in patients with myocardial infarction, showing a consistency between strain, infarction location, and coronary occlusion. This paper also presentsan interesting clinical application of graph-based metric estimators, showing their value for solving practical problems found in medical imaging.

Landslides cause tsunami waves: insights from Aysén Fjord, Chile

On 21 April 2007, an Mw 6.2 earthquake produced an unforeseen chain of events in the Aysén fjord (Chilean Patagonia, 45.5°S). The earthquake triggered hundreds of subaerial landslides along the fjord flanks. Some of the landslides eventually involved a subaqueous component that, in turn, generated a series of displacement waves tsunami- like waves produced by the fast entry of a ubaerial landmass into a water body within the fjord [Naranjo et al., 2009; Sepúlveda and Serey, 2009; Hermanns et al., 2013]. These waves, with run-ups several meters high along the shoreline, caused 10 fatalities. In addition, they severely damaged salmon farms, which constitute the main economic activity in the region, setting free millions of cultivated salmon with still unknown ecological consequences.

Pore-scale modeling of two-phase flow instabilities in porous media

Les problèmes d'écoulements multiphasiques en média poreux sont d'un grand intérêt pour de nombreuses applications scientifiques et techniques ; comme la séquestration de C02, l'extraction de pétrole et la dépollution des aquifères. La complexité intrinsèque des systèmes multiphasiques et l'hétérogénéité des formations géologiques sur des échelles multiples représentent un challenge majeur pour comprendre et modéliser les déplacements immiscibles dans les milieux poreux. Les descriptions à l'échelle supérieure basées sur la généralisation de l'équation de Darcy sont largement utilisées, mais ces méthodes sont sujettes à limitations pour les écoulements présentant de l'hystérèse. Les avancées récentes en terme de performances computationnelles et le développement de méthodes précises pour caractériser l'espace interstitiel ainsi que la distribution des phases ont favorisé l'utilisation de modèles qui permettent une résolution fine à l'échelle du pore. Ces modèles offrent un aperçu des caractéristiques de l'écoulement qui ne peuvent pas être facilement observées en laboratoire et peuvent être utilisé pour expliquer la différence entre les processus physiques et les modèles à l'échelle macroscopique existants. L'objet premier de la thèse se porte sur la simulation numérique directe : les équations de Navier-Stokes sont résolues dans l'espace interstitiel et la méthode du volume de fluide (VOF) est employée pour suivre l'évolution de l'interface. Dans VOF, la distribution des phases est décrite par une fonction fluide pour l'ensemble du domaine et des conditions aux bords particulières permettent la prise en compte des propriétés de mouillage du milieu poreux. Dans la première partie de la thèse, nous simulons le drainage dans une cellule Hele-Shaw 2D avec des obstacles cylindriques. Nous montrons que l'approche proposée est applicable même pour des ratios de densité et de viscosité très importants et permet de modéliser la transition entre déplacement stable et digitation visqueuse. Nous intéressons ensuite à l'interprétation de la pression capillaire à l'échelle macroscopique. Nous montrons que les techniques basées sur la moyenne spatiale de la pression présentent plusieurs limitations et sont imprécises en présence d'effets visqueux et de piégeage. Au contraire, une définition basée sur l'énergie permet de séparer les contributions capillaires des effets visqueux. La seconde partie de la thèse est consacrée à l'investigation des effets d'inertie associés aux reconfigurations irréversibles du ménisque causé par l'interface des instabilités. Comme prototype pour ces phénomènes, nous étudions d'abord la dynamique d'un ménisque dans un pore angulaire. Nous montrons que, dans un réseau de pores cubiques, les sauts et reconfigurations sont si fréquents que les effets d'inertie mènent à différentes configurations des fluides. A cause de la non-linéarité du problème, la distribution des fluides influence le travail des forces de pression, qui, à son tour, provoque une chute de pression dans la loi de Darcy. Cela suggère que ces phénomènes devraient être pris en compte lorsque que l'on décrit l'écoulement multiphasique en média poreux à l'échelle macroscopique. La dernière partie de la thèse s'attache à démontrer la validité de notre approche par une comparaison avec des expériences en laboratoire : un drainage instable dans un milieu poreux quasi 2D (une cellule Hele-Shaw avec des obstacles cylindriques). Plusieurs simulations sont tournées sous différentes conditions aux bords et en utilisant différents modèles (modèle intégré 2D et modèle 3D) afin de comparer certaines quantités macroscopiques avec les observations au laboratoire correspondantes. Malgré le challenge de modéliser des déplacements instables, où, par définition, de petites perturbations peuvent grandir sans fin, notre approche numérique apporte de résultats satisfaisants pour tous les cas étudiés. - Problems involving multiphase flow in porous media are of great interest in many scientific and engineering applications including Carbon Capture and Storage, oil recovery and groundwater remediation. The intrinsic complexity of multiphase systems and the multi scale heterogeneity of geological formations represent the major challenges to understand and model immiscible displacement in porous media. Upscaled descriptions based on generalization of Darcy's law are widely used, but they are subject to several limitations for flow that exhibit hysteric and history- dependent behaviors. Recent advances in high performance computing and the development of accurate methods to characterize pore space and phase distribution have fostered the use of models that allow sub-pore resolution. These models provide an insight on flow characteristics that cannot be easily achieved by laboratory experiments and can be used to explain the gap between physical processes and existing macro-scale models. We focus on direct numerical simulations: we solve the Navier-Stokes equations for mass and momentum conservation in the pore space and employ the Volume Of Fluid (VOF) method to track the evolution of the interface. In the VOF the distribution of the phases is described by a fluid function (whole-domain formulation) and special boundary conditions account for the wetting properties of the porous medium. In the first part of this thesis we simulate drainage in a 2-D Hele-Shaw cell filled with cylindrical obstacles. We show that the proposed approach can handle very large density and viscosity ratios and it is able to model the transition from stable displacement to viscous fingering. We then focus on the interpretation of the macroscopic capillary pressure showing that pressure average techniques are subject to several limitations and they are not accurate in presence of viscous effects and trapping. On the contrary an energy-based definition allows separating viscous and capillary contributions. In the second part of the thesis we investigate inertia effects associated with abrupt and irreversible reconfigurations of the menisci caused by interface instabilities. As a prototype of these phenomena we first consider the dynamics of a meniscus in an angular pore. We show that in a network of cubic pores, jumps and reconfigurations are so frequent that inertia effects lead to different fluid configurations. Due to the non-linearity of the problem, the distribution of the fluids influences the work done by pressure forces, which is in turn related to the pressure drop in Darcy's law. This suggests that these phenomena should be taken into account when upscaling multiphase flow in porous media. The last part of the thesis is devoted to proving the accuracy of the numerical approach by validation with experiments of unstable primary drainage in a quasi-2D porous medium (i.e., Hele-Shaw cell filled with cylindrical obstacles). We perform simulations under different boundary conditions and using different models (2-D integrated and full 3-D) and we compare several macroscopic quantities with the corresponding experiment. Despite the intrinsic challenges of modeling unstable displacement, where by definition small perturbations can grow without bounds, the numerical method gives satisfactory results for all the cases studied.

Changes in running mechanics and spring-mass behaviour during a 5-km time trial.

Research into the biomechanical manifestation of fatigue during exhaustive runs is increasingly popular but additional understanding of the adaptation of the spring-mass behaviour during the course of strenuous, self-paced exercises continues to be a challenge in order to develop optimized training and injury prevention programs. This study investigated continuous changes in running mechanics and spring-mass behaviour during a 5-km run. 12 competitive triathletes performed a 5-km running time trial (mean performance: ̴17 min 30 s) on a 200 m indoor track. Vertical and anterior-posterior ground reaction forces were measured every 200 m by a 5-m long force platform system, and used to determine spring-mass model characteristics. After a fast start, running velocity progressively decreased (- 11.6%; P<0.001) in the middle part of the race before an end spurt in the final 400-600 m. Stride length (- 7.4%; P<0.001) and frequency (- 4.1%; P=0.001) decreased over the 25 laps, while contact time (+ 8.9%; P<0.001) and total stride duration (+ 4.1%; P<0.001) progressively lengthened. Peak vertical forces (- 2.0%; P<0.01) and leg compression (- 4.3%; P<0.05), but not centre of mass vertical displacement (+ 3.2%; P>0.05), decreased with time. As a result, vertical stiffness decreased (- 6.0%; P<0.001) during the run, whereas leg stiffness changes were not significant (+ 1.3%; P>0.05). Spring-mass behaviour progressively changes during a 5-km time trial towards deteriorated vertical stiffness, which alters impact and force production characteristics.

FXYD Proteins Reverse Inhibition of the Na+-K+ Pump Mediated by Glutathionylation of Its {beta}1 Subunit.

The seven members of the FXYD protein family associate with the Na(+)-K(+) pump and modulate its activity. We investigated whether conserved cysteines in FXYD proteins are susceptible to glutathionylation and whether such reactivity affects Na(+)-K(+) pump function in cardiac myocytes and Xenopus oocytes. Glutathionylation was detected by immunoblotting streptavidin precipitate from biotin-GSH loaded cells or by a GSH antibody. Incubation of myocytes with recombinant FXYD proteins resulted in competitive displacement of native FXYD1. Myocyte and Xenopus oocyte pump currents were measured with whole-cell and two-electrode voltage clamp techniques, respectively. Native FXYD1 in myocytes and FXYD1 expressed in oocytes were susceptible to glutathionylation. Mutagenesis identified the specific cysteine in the cytoplasmic terminal that was reactive. Its reactivity was dependent on flanking basic amino acids. We have reported that Na(+)-K(+) pump β(1) subunit glutathionylation induced by oxidative signals causes pump inhibition in a previous study. In the present study, we found that β(1) subunit glutathionylation and pump inhibition could be reversed by exposing myocytes to exogenous wild-type FXYD3. A cysteine-free FXYD3 derivative had no effect. Similar results were obtained with wild-type and mutant FXYD proteins expressed in oocytes. Glutathionylation of the β(1) subunit was increased in myocardium from FXYD1(-/-) mice. In conclusion, there is a dependence of Na(+)-K(+) pump regulation on reactivity of two specifically identified cysteines on separate components of the multimeric Na(+)-K(+) pump complex. By facilitating deglutathionylation of the β(1) subunit, FXYD proteins reverse oxidative inhibition of the Na(+)-K(+) pump and play a dynamic role in its regulation.

Organ volume measurements: comparison between MRI and autopsy findings in infants following sudden unexpected death.

OBJECTIVE: To assess the accuracy of a semiautomated 3D volume reconstruction method for organ volume measurement by postmortem MRI. METHODS: This prospective study was approved by the institutional review board and the infants' parents gave their consent. Postmortem MRI was performed in 16 infants (1 month to 1 year of age) at 1.5 T within 48 h of their sudden death. Virtual organ volumes were estimated using the Myrian software. Real volumes were recorded at autopsy by water displacement. The agreement between virtual and real volumes was quantified following the Bland and Altman's method. RESULTS: There was a good agreement between virtual and real volumes for brain (mean difference: -0.03% (-13.6 to +7.1)), liver (+8.3% (-9.6 to +26.2)) and lungs (+5.5% (-26.6 to +37.6)). For kidneys, spleen and thymus, the MRI/autopsy volume ratio was close to 1 (kidney: 0.87±0.1; spleen: 0.99±0.17; thymus: 0.94±0.25), but with a less good agreement. For heart, the MRI/real volume ratio was 1.29±0.76, possibly due to the presence of residual blood within the heart. The virtual volumes of adrenal glands were significantly underestimated (p=0.04), possibly due to their very small size during the first year of life. The percentage of interobserver and intraobserver variation was lower or equal to 10%, but for thymus (15.9% and 12.6%, respectively) and adrenal glands (69% and 25.9%). CONCLUSIONS: Virtual volumetry may provide significant information concerning the macroscopic features of the main organs and help pathologists in sampling organs that are more likely to yield histological findings.

Measurement of the tree root growth using electrical impedance spectroscopy

Abstract

The Determinants of Exit in Argentina: Core and Peripheral Regions

This paper analyses the regional determinants of exit in Argentina. We find evidence of a dynamic revolving door by which past entrants increase current exits, particularly in the peripheral regions. In the central regions, current and past incumbents cause an analogous displacement effect. Also, exit shows a U-shaped relationship with respect to the informal economy, although the positive effect is weaker in the central regions. These findings point to the existence of a core-periphery structure in the spatial distribution of exits. Key words: firm exit, count data models, Argentina JEL: R12; R30; C33

Le Briançonnais, terrain exotique dans les Alpes ?

New reconstructions of the Western Alps from late Early Jurassic till early Tertiary are proposed. These reconstructions use deep lithospheric data gathered through recent seismic surveys and tomographic studies carried out in the Alps. The present day position, under the Po plain, of the southern limit of the European plate (fig. 1), allows to define the former geometry of the Brianconnais peninsula. The Brianconnais domain is regarded as an exotic terrane formerly belonging to the European margin until Late Jurassic, then transported eastward during the drift of Iberia (fig. 5). Therefore, on a present day Western Alps cross section, a duplication of the European continental margin can be recognized (fig. 10). Stratigraphic and sedimentological data along a zone linking the Pyrenean fracture zone to the Brianconnais, can be related to a rifting event starting in Oxfordian time. This event is responsible for the Late Jurassic till mid-Cretaceous drift of Iberia opening, first the northern Atlantic, then the Gulf of Biscay. Simultaneously, the drift of the Brianconnais will open the Valais ocean and close the Piemontese ocean. The resulting oblique collision zone between the Brianconnais and the Apulian margin generates HP/LT metamorphism starting in Early Cretaceous. The eastward drift of the Brianconnais peninsula will eventually bring it in front of a more northerly segment of the former European margin. The thrusting of the Brianconnais unto that margin takes place in early Tertiary (fig. 9), following the subduction of the Valais ocean. The present nappe pile results not only from continent/continent frontal collision, but also from important lateral displacement of terranes, the most important one being the Brianconnais. The dilemma of `'en echelon'' oceanic domains in the Alps is an outcome of these translations. A solution is found when considering the opening of a Cretaceous Valais ocean across the European margin, running out eastward into the Piemontese ocean, where the drift is taken up along a former transform fault and compensated by subduction under the Apulian margin (fig. 8). In the Western Alps we are then dealing with two oceans, the Piemontese and the Valaisan and a duplicated European margin. In the Eastern Alps the single Piemontese ocean is cut by newly created oceanic crust. All these elements will be incorporated into the Penninic structural domain which does not represent a former unique paleogeographic area, it is a composite accretionary domain squeezed between Europe and Apulia.

Dimorphic turning bias in spontaneous rotational movement

Turning biases have been associated with unbalanced hemispheric dopaminergic activity, and this activity has been correlated with cue-directed behaviors. Moreover, a sexual differentiation in hippocampal dopaminergic receptors following learning has been shown. In humans, pointing responses towards the starting point is commonly used to assess the accuracy of direction estimation after locomotion. Thus, it may be of interest for the field of spatial cognition to explore human sex differences in spontaneous turning bias when a body rotation is required. To this end, male and female blindfolded subjects were guided in a linear displacement and asked to rotate in order to point in the direction of the starting position. The main finding was a massive difference between men and women. 80% of women showed a turning bias to the right when 69% of men showed a bias to the left. Moreover, these preferences were not correlated with handedness. These results suggest basic preferences associated to sex might influence male and female performance in spatial cognition. They also suggest experimental procedures may be biased in favor of male or female strategies. Therefore, such preferences should be considered in order to gain further insight into the development of more balanced procedures.

Combining digital elevation model analysis and run-out modeling to characterize hazard posed by a potentially unstable rock slope at Turtle Mountain, Alberta, Canada

Turtle Mountain in Alberta, Canada has become an important field laboratory for testing different techniques related to the characterization and monitoring of large slope mass movements as the stability of large portions of the eastern face of the mountain is still questionable. In order to better quantify the volumes potentially unstable and the most probable failure mechanisms and potential consequences, structural analysis and runout modeling were preformed. The structural features of the eastern face were investigated using a high resolution digital elevation model (HRDEM). According to displacement datasets and structural observations, potential failure mechanisms affecting different portions of the mountain have been assessed. The volumes of the different potentially unstable blocks have been calculated using the Sloping Local Base Level (SLBL) method. Based on the volume estimation, two and three dimensional dynamic runout analyses have been performed. Calibration of this analysis is based on the experience from the adjacent Frank Slide and other similar rock avalanches. The results will be used to improve the contingency plans within the hazard area.

Field Testing of Integral Abutments, Interim Report, HR-399, 1998

Integral abutment bridges are constructed without an expansion joint in the superstructure of the bridge; therefore, the bridge girders, deck, abutment diaphragms, and abutments are monolithically constructed. The abutment piles in an integral abutment bridge are vertically orientated, and they are embedded into the pile cap. When this type of a bridge experiences thermal expansion or contraction, horizontal displacements are induced at the top of the abutment piles. The flexibility of the abutment piles eliminates the need to provide an expansion joint at the inside face to the abutments: Integral abutment bridge construction has been used in Iowa and other states for many years. This research is evaluating the performance of integral abutment bridges by investigating thermally induced displacements, strains, and temperatures in two Iowa bridges. Each bridge has a skewed alignment, contains five prestressed concrete girders that support a 30-ft wide roadway for three spans, and involves a water crossing. The bridges will be monitored for about two years. For each bridge, an instrumentation package includes measurement devices and hardware and software support systems. The measurement devices are displacement transducers, strain gages, and thermocouples. The hardware and software systems include a data-logger; multiplexers; directline telephone service and computer terminal modem; direct-line electrical power; lap-top computer; and an assortment of computer programs for monitoring, transmitting, and management of the data. Instrumentation has been installed on a bridge located in Guthrie County, and similar instrumentation is currently being installed on a bridge located in Story County. Preliminary test results for the bridge located in Guthrie County have revealed that temperature changes of the bridge deck and girders induce both longitudinal and transverse displacements of the abutments and significant flexural strains in the abutment piles. For an average temperature range of 73° F for the superstructure concrete in the bridge located in Guthrie County, the change in the bridge length was about 1 118 in. and the maximum, strong-axis, flexural-strain range for one of the abutment piles was about 400 micro-strains, which corresponds to a stress range of about 11,600 psi.

Approche minutieuse des fractures epiphysaires de l'humérus distal [Meticulous approach of the distal numeral epiphyseal fractures].

Epiphyseal fractures of the distal humerus and their treatment are common, source of many complications that may be immediate but also delayed. If the families accept the possibility of trauma to their child as part of his leisure activities or sports at times intense, they are much more intolerant if their child does not return after injury and repair, optimal function and physical performance unchanged. This intolerance, faced with the trauma but not its consequences, requires us to be particularly careful in the information given to the patient but also in monitoring the first few weeks. The treatment of epiphyseal fractures remains difficult and must be performed by trained teams, aware of the subtleties of pathology. It is essential to organize close monitoring and critical to be able to correct errors or secondary displacement in the first 15 days.

Influence of infusion line compliance on drug delivery rate during acute line loop formation.

OBJECTIVE: To determine whether infusion line compliance contributes to irregular drug delivery during vertical displacement of syringe pumps. DESIGN: Five different commercially available infusion lines were studied at infusion rates of 0.5, 1.0, and 1.5 ml/h. Zero drug delivery time was measured after acute line loop formation (70 cm) using an electronic balance. Compliance of each infusion line was calculated using a pressure transducer and measurement of the occlusion release bolus at 300 mmHg occlusion pressure. Finally, the influence of infusion line compliance on drug delivery during acute lowering of the syringe pump was studied using low- and high-compliance infusion lines. RESULTS: Acute line loop formation resulted in zero drug delivery time from 5.1 +/- 1.5 to 44.0 +/- 6.8 s at flow rates of 0.5 ml/h. Increased flow rates significantly reduced loop-induced flow variability. A close correlation was found between zero drug delivery time and calculated infusion line compliance at 0.5 ml/h (linear regression R2 = 0.79). Lowering of the syringe pump 50 cm prolonged zero drug delivery time from 295.8 +/- 20.7 s with the low-compliance tube to 463.3 +/- 24.0 s with the high-compliance infusion line. CONCLUSIONS: Infusion line compliance contributes to irregular drug delivery associated with vertical displacement of syringe pumps. Siphoning of the infusion line during patient care should be avoided, and flow rates of 1 ml/h or higher are recommended. Low-compliance infusion lines are indicated whenever highly short-acting vasoactive drugs at low delivery rates are administered.

A spatially explicit life cycle inventory of the global textile chain

Life cycle analyses (LCA) approaches require adaptation to reflect the increasing delocalization of production to emerging countries. This work addresses this challenge by establishing a country-level, spatially explicit life cycle inventory (LCI). This study comprises three separate dimensions. The first dimension is spatial: processes and emissions are allocated to the country in which they take place and modeled to take into account local factors. Emerging economies China and India are the location of production, the consumption occurs in Germany, an Organisation for Economic Cooperation and Development country. The second dimension is the product level: we consider two distinct textile garments, a cotton T-shirt and a polyester jacket, in order to highlight potential differences in the production and use phases. The third dimension is the inventory composition: we track CO2, SO2, NO (x), and particulates, four major atmospheric pollutants, as well as energy use. This third dimension enriches the analysis of the spatial differentiation (first dimension) and distinct products (second dimension). We describe the textile production and use processes and define a functional unit for a garment. We then model important processes using a hierarchy of preferential data sources. We place special emphasis on the modeling of the principal local energy processes: electricity and transport in emerging countries. The spatially explicit inventory is disaggregated by country of location of the emissions and analyzed according to the dimensions of the study: location, product, and pollutant. The inventory shows striking differences between the two products considered as well as between the different pollutants considered. For the T-shirt, over 70% of the energy use and CO2 emissions occur in the consuming country, whereas for the jacket, more than 70% occur in the producing country. This reversal of proportions is due to differences in the use phase of the garments. For SO2, in contrast, over two thirds of the emissions occur in the country of production for both T-shirt and jacket. The difference in emission patterns between CO2 and SO2 is due to local electricity processes, justifying our emphasis on local energy infrastructure. The complexity of considering differences in location, product, and pollutant is rewarded by a much richer understanding of a global production-consumption chain. The inclusion of two different products in the LCI highlights the importance of the definition of a product's functional unit in the analysis and implications of results. Several use-phase scenarios demonstrate the importance of consumer behavior over equipment efficiency. The spatial emission patterns of the different pollutants allow us to understand the role of various energy infrastructure elements. The emission patterns furthermore inform the debate on the Environmental Kuznets Curve, which applies only to pollutants which can be easily filtered and does not take into account the effects of production displacement. We also discuss the appropriateness and limitations of applying the LCA methodology in a global context, especially in developing countries. Our spatial LCI method yields important insights in the quantity and pattern of emissions due to different product life cycle stages, dependent on the local technology, emphasizing the importance of consumer behavior. From a life cycle perspective, consumer education promoting air-drying and cool washing is more important than efficient appliances. Spatial LCI with country-specific data is a promising method, necessary for the challenges of globalized production-consumption chains. We recommend inventory reporting of final energy forms, such as electricity, and modular LCA databases, which would allow the easy modification of underlying energy infrastructure.