Investigation of young water inflow in karst aquifers using SF6–CFC–3H/He–85Kr–39Ar and stable isotope components

Karst aquifers are known for their wide distribution of water transfer velocities. From this observation, a multiple geochemical tracer approach seems to be particularly well suited to provide a significant assessment of groundwater flows, but the choice of adapted tracers is essential. In this study, several common tracers in karst aquifers such as physicochemical parameters, major ions, stable isotopes, and d13C to more specific tracers such as dating tracers – 14C, 3H, 3H–3He, CFC-12, SF6 and 85Kr, and 39Ar – were used, in a fractured karstic carbonated aquifer located in Burgundy (France). The information carried by each tracer and the best sampling strategy are compared on the basis of geochemical monitoring done during several recharge events and over longer time periods (months to years). This study’s results demonstrate that at the seasonal and recharge event time scale, the variability of concentrations is low for most tracers due to the broad spectrum of groundwater mixings. The tracers used traditionally for the study of karst aquifers, i.e., physicochemical parameters and major ions, efficiently describe hydrological processes such as the direct and differed recharge, but require being monitored at short time steps during recharge events to be maximized. From stable isotopes, tritium, and Cl� contents, the proportion of the fast direct recharge by the largest porosity was estimated using a binary mixing model. The use of tracers such as CFC-12, SF6, and 85Kr in karst aquifers provides additional information, notably an estimation of apparent age, but they require good preliminary knowledge of the karst system to interpret the results suitably. The CFC-12 and SF6 methods efficiently determine the apparent age of baseflow, but it is preferable to sample the groundwater during the recharge event. Furthermore, these methods are based on different assumptions such as regional enrichment in atmospheric SF6, excess air, and flow models among others. 85Kr and 39Ar concentrations can potentially provide a more direct estimation of groundwater residence time. Conversely, the 3H–3He method is inefficient in the karst aquifer for dating due to 3He degassing.

Systemic inflammation is higher in peripheral artery disease than in stable coronary artery disease.

OBJECTIVE The knowledge on the level of systemic inflammation in peripheral artery disease (PAD) is less well established than that in coronary artery disease (CAD). Systemic inflammation frequently coincides with atherosclerosis, but also with various traits of the metabolic syndrome (MetS). The individual contribution of CAD, PAD, and the MetS to inflammation is not known. METHODS We enrolled a total of 1396 patients, 460 patients with PAD Fontaine stages IIa-IV verified by duplex ultrasound (PAD group) and 936 patients free of limb claudication undergoing coronary angiography, of whom 507 had significant CAD with coronary stenoses ≥50% (CAD group), and 429 did not have significant CAD at angiography (control group). RESULTS C-reactive protein (CRP) was significantly higher in the PAD than in the CAD or in the control group (0.86 ± 1.85 mg/dl versus 0.44 ± 0.87 mg/dl and 0.39 ± 0.52 mg/dl, respectively, p < 0.001 for both comparisons). These significant differences were confirmed when patients with and subjects without the MetS were analyzed separately. In particular, within the PAD group, CRP was significantly higher in patients with the MetS than in subjects without the MetS (1.04 ± 2.01 vs. 0.67 ± 1.64 mg/dl; p = 0.001) and both, the presence of PAD and the MetS proved to be independently associated with CRP in analysis of covariance (F = 31.84; p < 0.001 and F = 10.52; p = 0.001, respectively). CONCLUSION Inflammatory activity in PAD patients is higher than in CAD patients and is particularly high in PAD patients affected by the MetS. Low grade systemic inflammation is independently associated with both the MetS and PAD.

Effect of large doses of parenteral vitamin D on glycaemic control and calcium/phosphate metabolism in patients with stable type 2 diabetes mellitus: a randomised, placebo-controlled, prospective pilot study.

OBJECTIVE Vitamin D (D₃) status is reported to correlate negatively with insulin production and insulin sensitivity in patients with type 2 diabetes mellitus (T2DM). However, few placebo-controlled intervention data are available. We aimed to assess the effect of large doses of parenteral D3 on glycosylated haemoglobin (HbA(₁c)) and estimates of insulin action (homeostasis model assessment insulin resistance: HOMA-IR) in patients with stable T2DM. MATERIALS AND METHODS We performed a prospective, randomised, double-blind, placebo-controlled pilot study at a single university care setting in Switzerland. Fifty-five patients of both genders with T2DM of more than 10 years were enrolled and randomised to either 300,000 IU D₃ or placebo, intramuscularly. The primary endpoint was the intergroup difference in HbA(₁c) levels. Secondary endpoints were: changes in insulin sensitivity, albuminuria, calcium/phosphate metabolism, activity of the renin-aldosterone axis and changes in 24-hour ambulatory blood pressure values. RESULTS After 6 months of D₃ supply, there was a significant intergroup difference in the change in HbA(₁c) levels (relative change [mean ± standard deviation] +2.9% ± 1.5% in the D₃ group vs +6.9% ± 2.1% the in placebo group, p = 0.041) as HOMA-IR decreased by 12.8% ± 5.6% in the D₃ group and increased by 10% ± 5.4% in the placebo group (intergroup difference, p = 0.032). Twenty-four-hour urinary albumin excretion decreased in the D₃ group from 200 ± 41 to 126 ± 39, p = 0.021). There was no significant intergroup difference for the other secondary endpoints. CONCLUSIONS D₃ improved insulin sensitivity (based on HOMA-IR) and affected the course of HbA(₁c) positively compared with placebo in patients with T2DM.

Lipid- and mechanosensitivities of sodium/hydrogen exchangers analyzed by electrical methods.

Sodium/hydrogen exchangers (NHEs) are ubiquitous ion transporters that serve multiple cell functions. We have studied two mammalian isoforms, NHE1 (ubiquitous) and NHE3 (epithelial-specific), by measuring extracellular proton (H+) gradients during whole-cell patch clamp with perfusion of the cell interior. Maximal Na(+)-dependent H+ fluxes (JH+) are equivalent to currents >20 pA for NHE1 in Chinese hamster ovary fibroblasts, >200 pA for NHE1 in guinea pig ventricular myocytes, and 5-10 pA for NHE3 in opossum kidney cells. The fluxes are blocked by an NHE inhibitor, ethylisopropylamiloride, and are absent in NHE-deficient AP-1 cells. NHE1 activity is stable with perfusion of nonhydrolyzable ATP [adenosine 5'-(beta,gamma-imido)triphosphate], is abolished by ATP depletion (2 deoxy-D-glucose with oligomycin or perfusion of apyrase), can be restored with phosphatidylinositol 4,5-bisphosphate, and is unaffected by actin cytoskeleton disruption (latrunculin or pipette perfusion of gelsolin). NHE3 (but not NHE1) is reversibly activated by phosphatidylinositol 3,4,5-trisphosphate. Both NHE1 and NHE3 activities are disrupted in giant patches during gigaohm seal formation. NHE1 (but not NHE3) is reversibly activated by cell shrinkage, even at neutral cytoplasmic pH without ATP, and inhibited by cell swelling. NHE1 in Chinese hamster ovary fibroblasts (but not NHE3 in opossum kidney cells) is inhibited by agents that thin the membrane (L-alpha-lysophosphatidylcholine and octyl-beta-D-glucopyranoside) and activated by cholesterol enrichment, which thickens membranes. Expressed in AP-1 cells, however, NHE1 is insensitive to these agents but remains sensitive to volume changes. Thus, changes of hydrophobic mismatch can modulate NHE1 but do not underlie its volume sensitivity.

Decreased phosphatidylcholine plasmalogens - A putative novel lipid signature in patients with stable coronary artery disease and acute myocardial infarction.

OBJECTIVE Glycerophospholipids and sphingolipids are structurally heterogeneous due to differences in the O- and N-linked fatty acids and head groups. Sphingolipids also show a heterogeneity in their sphingoid base composition which up to now has been little appreciated. The aim of this study was to investigate the association of certain glycerophospholipid and sphingolipid species with stable coronary artery disease (CAD) and acute myocardial infarction (AMI). METHODS The lipid profile in plasma from patients with stable CAD (n = 18) or AMI (n = 17) was compared to healthy subjects (n = 14). Sixty five glycerophospholipid and sphingolipid species were quantified by LC-MS. The relative distribution of these lipids into lipoprotein fractions was analyzed. RESULTS In the CAD cohort, 45 glycerophospholipid and sphingolipid species were significantly lower compared to healthy controls. In the AMI group, 42 glycerophospholipid and sphingolipid species were reduced. Four PC plasmalogens (PC33:1, PC33:2, PC33:3 and PC35:3) showed the most significant difference. Out of eleven analyzed sphingoid bases, four were lower in the CAD and six in the AMI group. Sphingosine-1-phosphate (S1P) levels were reduced in the AMI group whereas an atypical C16:1 S1P was lower in both groups. Phosphatidylcholine and sphingomyelin species were exclusively present in lipoprotein particles, whereas lysophosphatidylcholines were mainly found in the lipoprotein-free fraction. The observed differences were not explained by the use of statins as confirmed in a second, independent cohort. CONCLUSIONS Reduced levels of four PC plasmalogens (PC33:1, PC33:2, PC33:3 and PC35:3) were identified as a putatively novel lipid signature for CAD and AMI.

Preprocedural High-Sensitivity Cardiac Troponin T and Clinical Outcomes in Patients With Stable Coronary Artery Disease Undergoing Elective Percutaneous Coronary Intervention.

BACKGROUND Cardiac troponin detected by new-generation, highly sensitive assays predicts clinical outcomes among patients with stable coronary artery disease (SCAD) treated medically. The prognostic value of baseline high-sensitivity cardiac troponin T (hs-cTnT) elevation in SCAD patients undergoing elective percutaneous coronary interventions is not well established. This study assessed the association of preprocedural levels of hs-cTnT with 1-year clinical outcomes among SCAD patients undergoing percutaneous coronary intervention. METHODS AND RESULTS Between 2010 and 2014, 6974 consecutive patients were prospectively enrolled in the Bern Percutaneous Coronary Interventions Registry. Among patients with SCAD (n=2029), 527 (26%) had elevated preprocedural hs-cTnT above the upper reference limit of 14 ng/L. The primary end point, mortality within 1 year, occurred in 20 patients (1.4%) with normal hs-cTnT versus 39 patients (7.7%) with elevated baseline hs-cTnT (P<0.001). Patients with elevated hs-cTnT had increased risks of all-cause (hazard ratio 5.73; 95% confidence intervals 3.34-9.83; P<0.001) and cardiac mortality (hazard ratio 4.68; 95% confidence interval 2.12-10.31; P<0.001). Preprocedural hs-TnT elevation remained an independent predictor of 1-year mortality after adjustment for relevant risk factors, including age, sex, and renal failure (adjusted hazard ratio 2.08; 95% confidence interval 1.10-3.92; P=0.024). A graded mortality risk was observed across higher tertiles of elevated preprocedural hs-cTnT, but not among patients with hs-cTnT below the upper reference limit. CONCLUSIONS Preprocedural elevation of hs-cTnT is observed in one fourth of SCAD patients undergoing elective percutaneous coronary intervention. Increased levels of preprocedural hs-cTnT are proportionally related to the risk of death and emerged as independent predictors of all-cause mortality within 1 year. CLINICAL TRIAL REGISTRATION URL: http://www.clinicaltrials.gov. Unique identifier: NCT02241291.

A review of cataract in stable and conflict zones in Africa

Aim. To review and explore cataract prevalence in stable and unstable countries by examining published and unpublished ocular literature about Africa from 1980 onwards.^ Methods. Searches using OVID, Proquest Dissertations, WHO, and Ebsco Host were done. The review was restricted to articles utilizing WHO definitions of blindness and low vision. Random cluster sampling technique with a minimum sample size of 1,500, and reporting causes of blindness categorized by age and gender were inclusion considerations in the selected articles. ^ Results. Blindness and low vision increased with conflict. Women and the elderly were more likely to have vision impairing cataract. Cataract was the leading cause of blindness; the prevalence range was 22%–81% for the reviewed nations.^ Conclusion. Instability was connected to higher cataract prevalence and worse visual outcome across all characteristics examined except cataract surgical rates. ^

A comparative review of water disinfection methods appropriate for developing countries and their efficacy, cost-efficiency, and usability

Developing countries are heavily burdened by limited access to safe drinking water and subsequent water-related diseases. Numerous water treatment interventions combat this public health crisis, encompassing both traditional and less-common methods. Of these, water disinfection serves as an important means to provide safe drinking water. Existing literature discusses a wide range of traditional treatment options and encourages the use of multi-barrier approaches including coagulation-flocculation, filtration, and disinfection. Most sources do not delve into approaches specifically appropriate for developing countries, nor do they exclusively examine water disinfection methods.^ The objective of this review is to focus on an extensive range of chemical, physio-chemical, and physical water disinfection techniques to provide a compilation, description and evaluation of options available. Such an objective provides further understanding and knowledge to better inform water treatment interventions and explores alternate means of water disinfection appropriate for developing countries. Appropriateness for developing countries corresponds to the effectiveness of an available, easy to use disinfection technique at providing safe drinking water at a low cost.^ Among chemical disinfectants, SWS sodium hypochlorite solution is preferred over sodium hypochlorite bleach due to consistent concentrations. Tablet forms are highly recommended chemical disinfectants because they are effective and very easy to use, but also because they are stable. Examples include sodium dichloroisocyanurate, calcium hypochlorite, and chlorine dioxide, which vary in cost depending on location and availability. Among physio-chemical disinfection options, electrolysis which produces mixed oxidants (MIOX) provides a highly effective disinfection option with a higher upfront cost but very low cost over the long term. Among physical disinfection options, solar disinfection (SODIS) applications are effective, but they treat only a fixed volume of water at a time. They come with higher initial costs but very low on-going costs. Additional effective disinfection techniques may be suitable depending on the location, availability and cost.^

AN INVESTIGATION OF THE EFFECT OF FORMALDEHYDE SPECIES ON ENVIRONMENTAL MEASUREMENT METHODS UNDER AMBIENT CONDITIONS

Species variations in formaldehyde solutions and gases were investigated by means of infrared spectral analysis. Double beam infrared spectrometry in conjunction with sodium chloride wafer technique and solvent compensation technique were employed. Formaldehyde species in various solutions were investigated. Formalin 37% was stable for many months. Refrigeration had no effects on its stability. Spectral changes were detected in 1000 ppm formaldehyde solutions. The absorbances of very diluted solutions up to 100 ppm were lower than the detection limit of the instruments. Solvent compensation improved resolution, but was associated with an observed lack of repeatability. Formaldehyde species in animal chambers containing animals and in mobile home air were analyzed with the infrared spectrophotometer equipped with a 10 cm gas cell. Spectra were not different from the spectrum of clean air. A portable single beam infrared spectrometer with a 20 meter pathlength was used for reinvestigation. Indoor formaldehyde could not be detected in the spectral; conversely, an absorption peak at 3.58 microns was found in the spectra of 3 and 15 ppm formaldehyde gas in animal chambers. This peak did not appear in the spectrum of the control chamber. Because of concerns over measurement bias among various analytical methods for formaldehyde, side-by-side comparisons were conducted in both laboratory and field measurements. The chromotropic acid method with water and 1% sodium bisulfite as collection media, the pararosaniline method, and a single beam infrared spectrometer were compared. Measurement bias was elucidated and the extent of the effects of temperature and humidity was also determined. The problems associated with related methods were discussed. ^

Stable oxygen isotope record of Neogloboquadrina pachyderma of ODP Sites 177-1091 and 177-1094

Ocean Drilling Program (ODP) cores permit us to extend the study of millennial-scale climate variability beyond the time period that is generally accessible for piston cores (i.e., the last glacial cycle). ODP Leg 177 provided for the first time continuous high sedimentation rate cores along a north-south transect from 41°to 53°S across the main subdivisions of the Southern Ocean (Shipboard Scientific Party, 1999, doi:10.2973/odp.proc.ir.177.101.1999). The main purpose of this drilling was to investigate the Pleistocene and Holocene paleoceanographic history of this region, documented in the sedimentary records. ODP Sites 1094, 1093, 1091, and 1089 accumulated throughout the Pleistocene at rates >10 cm/k.y. and are the most detailed Pleistocene climatic records ever retrieved from the Southern Ocean. These sections provide a unique opportunity to fill an important gap in the knowledge of the paleoclimatic evolution of the high southern latitude regions. The composite sections at each site were generated shipboard using magnetic susceptibility, gamma ray attenuation (GRA) density, and reflectance data to correlate the drill holes and splice together an optimal (complete and undisturbed) record of the sedimentary sequence at each site. A preliminary magnetic polarity stratigraphy was generated on the 'archive' halves of the core sections from each hole, using the shipboard pass-through magnetometer after demagnetization at a single peak alternating field (Shipboard Scientific Party, 1999). During July 1998, we sampled core sections spanning the mid-Pleistocene interval (0.65-1.2 Ma) from Sites 1094, 1093, and 1091 at the ODP Bremen Core Repository and have since then analyzed the stable isotopic ratios of foraminifers in samples from Sites 1094 and 1091. Our goals for these studies are to establish detailed chronology for the mid-Pleistocene Southern Ocean records from Leg 177 using high-resolution stable isotope analyses, and furthermore, to trace the evolution of millennial-scale variability in proxy records from older glacial and interglacial periods characterized by higher-frequency variation. Here, we report on our stratigraphic results to date and describe the laboratory methods employed for sample preparation and stable isotope analysis. Furthermore, we provide tab-delimited text files of the age models.

Surface sediment samples from several fore-arc basins west and southwest of the Indonesian Archipelago, analyzed by planktonic foraminifera, stable oxygen and carbon isotopic signals and opal and CaCO3 contents in bulk sediment

A total of 69 surface sediment samples from several fore-arc basins located west and southwest of the Indonesian Archipelago was analyzed with respect to the faunal composition of planktonic foraminifera, the stable oxygen and carbon isotopic signal of a surface-dwelling (Globigerinoides ruber) and a thermocline-dwelling (Neogloboquadrina dutertrei) species, and the opal and CaCO3 contents in bulk sediment. Our results show that the distribution pattern of opal in surface sediments corresponds well to the upwelling-induced chlorophyll concentration in the upper water column and thus, represents a reliable proxy for marine productivity in the coastal upwelling area off S and SW Indonesia. Present-day oceanography and marine productivity are also reflected in the tropical to subtropical and upwelling assemblages of planktonic foraminifera in the surface sediments, which in part differ from previous studies in this region probably due to different coring methods and dissolution effects. The average stable oxygen isotopic values (d18O) of G. ruber in surface sediments vary between 2.9 per mill and 3.2 per mill from basin to basin and correspond to the oceanographic settings during the SE monsoon (July-October) off west Sumatra, whereas off southern Indonesia, they reflect the NW monsoon (December-March) or annual average conditions. The d18O values of N. dutertrei show a stronger interbasinal variation between 1.6 per mill and 2.2 per mill and correspond to the upper thermocline hydrology in July-October. In addition, the difference between the shell carbon isotopic values (d13C) of G. ruber and N. dutertrei (Delta d13C) appears to be an appropriate productivity recorder only in the non-upwelling areas off west Sumatra. Consequently, joint interpretation of the isotopic values of these species is distinctive for different fore-arc basins W and SW of Indonesia and should be considered in paleoceanographic studies.

Mineralogical and stable isotope analyses of carbonates from the eastern Black Sea

Authigenic carbonates associated with cold seeps provide valuable archives of changes in the long-term seepage activity. To investigate the role of shallow-buried hydrates on the seepage strength and fluid composition we analysed methane-derived carbonate precipitates from a high-flux hydrocarbon seepage area ("Batumi seep area") located on the south-eastern Black Sea slope in ca. 850 m. In a novel approach, we combined computerized X-ray tomography (CT) with mineralogical and isotope geochemical methods to get additional insights into the three-dimensional internal structure of the carbonate build-ups. X-ray diffractometry revealed the presence of two different authigenic carbonate phases, i.e. pure aragonitic rims associated with vital microbial mats and high-Mg calcite cementing the hemipelagic sediment. As indicated by the CT images, the initial sediment has been strongly deformed, first plastic then brittle, leading to brecciation of the progressively cemented sediment. The aragonitic rims on the other hand, represent a presumably recent carbonate growth phase since they cover the already deformed sediment. The stable oxygen isotope signature indicates that the high-Mg calcite cement incorporated pore water mixed with substantial hydrate water amounts. This points at a dominant role of high gas/fluid flux from decomposing gas hydrates leading to the deformation and cementation of the overlying sediment. In contrast, the aragonitic rims do not show an influence of 18O-enriched hydrate water. The differences in d18O between the presumably recent aragonite precipitates and the older high-Mg cements suggest that periods of hydrate dissociation and vigorous fluid discharge alternated with times of hydrate stability and moderate fluid flow. These results indicate that shallow-buried gas hydrates are prone to episodic decomposition with associated vigorous fluid flow. This might have a profound impact on the seafloor morphology resulting e.g. in the formation of carbonate pavements and pockmark-like structures but might also affect the local carbon cycle.

Multimodal non-rigid image registration methods for therapy and surgery planning

La planificación pre-operatoria se ha convertido en una tarea esencial en cirugías y terapias de marcada complejidad, especialmente aquellas relacionadas con órgano blando. Un ejemplo donde la planificación preoperatoria tiene gran interés es la cirugía hepática. Dicha planificación comprende la detección e identificación precisa de las lesiones individuales y vasos así como la correcta segmentación y estimación volumétrica del hígado funcional. Este proceso es muy importante porque determina tanto si el paciente es un candidato adecuado para terapia quirúrgica como la definición del abordaje a seguir en el procedimiento. La radioterapia de órgano blando es un segundo ejemplo donde la planificación se requiere tanto para la radioterapia externa convencional como para la radioterapia intraoperatoria. La planificación comprende la segmentación de tumor y órganos vulnerables y la estimación de la dosimetría. La segmentación de hígado funcional y la estimación volumétrica para planificación de la cirugía se estiman habitualmente a partir de imágenes de tomografía computarizada (TC). De igual modo, en la planificación de radioterapia, los objetivos de la radiación se delinean normalmente sobre TC. Sin embargo, los avances en las tecnologías de imagen de resonancia magnética (RM) están ofreciendo progresivamente ventajas adicionales. Por ejemplo, se ha visto que el ratio de detección de metástasis hepáticas es significativamente superior en RM con contraste Gd–EOB–DTPA que en TC. Por tanto, recientes estudios han destacado la importancia de combinar la información de TC y RM para conseguir el mayor nivel posible de precisión en radioterapia y para facilitar una descripción precisa de las lesiones del hígado. Con el objetivo de mejorar la planificación preoperatoria en ambos escenarios se precisa claramente de un algoritmo de registro no rígido de imagen. Sin embargo, la gran mayoría de sistemas comerciales solo proporcionan métodos de registro rígido. Las medidas de intensidad de voxel han demostrado ser criterios de similitud de imágenes robustos, y, entre ellas, la Información Mutua (IM) es siempre la primera elegida en registros multimodales. Sin embargo, uno de los principales problemas de la IM es la ausencia de información espacial y la asunción de que las relaciones estadísticas entre las imágenes son homogéneas a lo largo de su domino completo. La hipótesis de esta tesis es que la incorporación de información espacial de órganos al proceso de registro puede mejorar la robustez y calidad del mismo, beneficiándose de la disponibilidad de las segmentaciones clínicas. En este trabajo, se propone y valida un esquema de registro multimodal no rígido 3D usando una nueva métrica llamada Información Mutua Centrada en el Órgano (Organ-Focused Mutual Information metric (OF-MI)) y se compara con la formulación clásica de la Información Mutua. Esto permite mejorar los resultados del registro en áreas problemáticas incorporando información regional al criterio de similitud, beneficiándose de la disponibilidad real de segmentaciones en protocolos estándares clínicos, y permitiendo que la dependencia estadística entre las dos modalidades de imagen difiera entre órganos o regiones. El método propuesto se ha aplicado al registro de TC y RM con contraste Gd–EOB–DTPA así como al registro de imágenes de TC y MR para planificación de radioterapia intraoperatoria rectal. Adicionalmente, se ha desarrollado un algoritmo de apoyo de segmentación 3D basado en Level-Sets para la incorporación de la información de órgano en el registro. El algoritmo de segmentación se ha diseñado específicamente para la estimación volumétrica de hígado sano funcional y ha demostrado un buen funcionamiento en un conjunto de imágenes de TC abdominales. Los resultados muestran una mejora estadísticamente significativa de OF-MI comparada con la Información Mutua clásica en las medidas de calidad de los registros; tanto con datos simulados (p<0.001) como con datos reales en registro hepático de TC y RM con contraste Gd– EOB–DTPA y en registro para planificación de radioterapia rectal usando OF-MI multi-órgano (p<0.05). Adicionalmente, OF-MI presenta resultados más estables con menor dispersión que la Información Mutua y un comportamiento más robusto con respecto a cambios en la relación señal-ruido y a la variación de parámetros. La métrica OF-MI propuesta en esta tesis presenta siempre igual o mayor precisión que la clásica Información Mutua y consecuentemente puede ser una muy buena alternativa en aplicaciones donde la robustez del método y la facilidad en la elección de parámetros sean particularmente importantes. Abstract Pre-operative planning has become an essential task in complex surgeries and therapies, especially for those affecting soft tissue. One example where soft tissue preoperative planning is of high interest is liver surgery. It involves the accurate detection and identification of individual liver lesions and vessels as well as the proper functional liver segmentation and volume estimation. This process is very important because it determines whether the patient is a suitable candidate for surgical therapy and the type of procedure. Soft tissue radiation therapy is a second example where planning is required for both conventional external and intraoperative radiotherapy. It involves the segmentation of the tumor target and vulnerable organs and the estimation of the planned dose. Functional liver segmentations and volume estimations for surgery planning are commonly estimated from computed tomography (CT) images. Similarly, in radiation therapy planning, targets to be irradiated and healthy and vulnerable tissues to be protected from irradiation are commonly delineated on CT scans. However, developments in magnetic resonance imaging (MRI) technology are progressively offering advantages. For instance, the hepatic metastasis detection rate has been found to be significantly higher in Gd–EOB–DTPAenhanced MRI than in CT. Therefore, recent studies highlight the importance of combining the information from CT and MRI to achieve the highest level of accuracy in radiotherapy and to facilitate accurate liver lesion description. In order to improve those two soft tissue pre operative planning scenarios, an accurate nonrigid image registration algorithm is clearly required. However, the vast majority of commercial systems only provide rigid registration. Voxel intensity measures have been shown to be robust measures of image similarity, and among them, Mutual Information (MI) is always the first candidate in multimodal registrations. However, one of the main drawbacks of Mutual Information is the absence of spatial information and the assumption that statistical relationships between images are the same over the whole domain of the image. The hypothesis of the present thesis is that incorporating spatial organ information into the registration process may improve the registration robustness and quality, taking advantage of the clinical segmentations availability. In this work, a multimodal nonrigid 3D registration framework using a new Organ- Focused Mutual Information metric (OF-MI) is proposed, validated and compared to the classical formulation of the Mutual Information (MI). It allows improving registration results in problematic areas by adding regional information into the similitude criterion taking advantage of actual segmentations availability in standard clinical protocols and allowing the statistical dependence between the two modalities differ among organs or regions. The proposed method is applied to CT and T1 weighted delayed Gd–EOB–DTPA-enhanced MRI registration as well as to register CT and MRI images in rectal intraoperative radiotherapy planning. Additionally, a 3D support segmentation algorithm based on Level-Sets has been developed for the incorporation of the organ information into the registration. The segmentation algorithm has been specifically designed for the healthy and functional liver volume estimation demonstrating good performance in a set of abdominal CT studies. Results show a statistical significant improvement of registration quality measures with OF-MI compared to MI with both simulated data (p<0.001) and real data in liver applications registering CT and Gd–EOB–DTPA-enhanced MRI and in registration for rectal radiotherapy planning using multi-organ OF-MI (p<0.05). Additionally, OF-MI presents more stable results with smaller dispersion than MI and a more robust behavior with respect to SNR changes and parameters variation. The proposed OF-MI always presents equal or better accuracy than the classical MI and consequently can be a very convenient alternative within applications where the robustness of the method and the facility to choose the parameters are particularly important.

Computational Methods for Identification of Vibrating Structures

System identification deals with the problem of building mathematical models of dynamical systems based on observed data from the system" [1]. In the context of civil engineering, the system refers to a large scale structure such as a building, bridge, or an offshore structure, and identification mostly involves the determination of modal parameters (the natural frequencies, damping ratios, and mode shapes). This paper presents some modal identification results obtained using a state-of-the-art time domain system identification method (data-driven stochastic subspace algorithms [2]) applied to the output-only data measured in a steel arch bridge. First, a three dimensional finite element model was developed for the numerical analysis of the structure using ANSYS. Modal analysis was carried out and modal parameters were extracted in the frequency range of interest, 0-10 Hz. The results obtained from the finite element modal analysis were used to determine the location of the sensors. After that, ambient vibration tests were conducted during April 23-24, 2009. The response of the structure was measured using eight accelerometers. Two stations of three sensors were formed (triaxial stations). These sensors were held stationary for reference during the test. The two remaining sensors were placed at the different measurement points along the bridge deck, in which only vertical and transversal measurements were conducted (biaxial stations). Point estimate and interval estimate have been carried out in the state space model using these ambient vibration measurements. In the case of parametric models (like state space), the dynamic behaviour of a system is described using mathematical models. Then, mathematical relationships can be established between modal parameters and estimated point parameters (thus, it is common to use experimental modal analysis as a synonym for system identification). Stable modal parameters are found using a stabilization diagram. Furthermore, this paper proposes a method for assessing the precision of estimates of the parameters of state-space models (confidence interval). This approach employs the nonparametric bootstrap procedure [3] and is applied to subspace parameter estimation algorithm. Using bootstrap results, a plot similar to a stabilization diagram is developed. These graphics differentiate system modes from spurious noise modes for a given order system. Additionally, using the modal assurance criterion, the experimental modes obtained have been compared with those evaluated from a finite element analysis. A quite good agreement between numerical and experimental results is observed.

Analysis of Meshfree Methods for Lagrangian Fluid-Structure Interaction

In this dissertation a new numerical method for solving Fluid-Structure Interaction (FSI) problems in a Lagrangian framework is developed, where solids of different constitutive laws can suffer very large deformations and fluids are considered to be newtonian and incompressible. For that, we first introduce a meshless discretization based on local maximum-entropy interpolants. This allows to discretize a spatial domain with no need of tessellation, avoiding the mesh limitations. Later, the Stokes flow problem is studied. The Galerkin meshless method based on a max-ent scheme for this problem suffers from instabilities, and therefore stabilization techniques are discussed and analyzed. An unconditionally stable method is finally formulated based on a Douglas-Wang stabilization. Then, a Langrangian expression for fluid mechanics is derived. This allows us to establish a common framework for fluid and solid domains, such that interaction can be naturally accounted. The resulting equations are also in the need of stabilization, what is corrected with an analogous technique as for the Stokes problem. The fully Lagrangian framework for fluid/solid interaction is completed with simple point-to-point and point-to-surface contact algorithms. The method is finally validated, and some numerical examples show the potential scope of applications.