Validity of the acoustic approximation in full-waveform seismic crosshole tomography

Acoustic waveform inversions are an increasingly popular tool for extracting subsurface information from seismic data. They are computationally much more efficient than elastic inversions. Naturally, an inherent disadvantage is that any elastic effects present in the recorded data are ignored in acoustic inversions. We investigate the extent to which elastic effects influence seismic crosshole data. Our numerical modeling studies reveal that in the presence of high contrast interfaces, at which P-to-S conversions occur, elastic effects can dominate the seismic sections, even for experiments involving pressure sources and pressure receivers. Comparisons of waveform inversion results using a purely acoustic algorithm on synthetic data that is either acoustic or elastic, show that subsurface models comprising small low-to-medium contrast (?30%) structures can be successfully resolved in the acoustic approximation. However, in the presence of extended high-contrast anomalous bodies, P-to-S-conversions may substantially degrade the quality of the tomographic images. In particular, extended low-velocity zones are difficult to image. Likewise, relatively small low-velocity features are unresolved, even when advanced a priori information is included. One option for mitigating elastic effects is data windowing, which suppresses later arriving seismic arrivals, such as shear waves. Our tests of this approach found it to be inappropriate because elastic effects are also included in earlier arriving wavetrains. Furthermore, data windowing removes later arriving P-wave phases that may provide critical constraints on the tomograms. Finally, we investigated the extent to which acoustic inversions of elastic data are useful for time-lapse analyses of high contrast engineered structures, for which accurate reconstruction of the subsurface structure is not as critical as imaging differential changes between sequential experiments. Based on a realistic scenario for monitoring a radioactive waste repository, we demonstrated that acoustic inversions of elastic data yield substantial distortions of the tomograms and also unreliable information on trends in the velocity changes.

Stiffness tomography by atomic force microscopy.

The atomic force microscope is a convenient tool to probe living samples at the nanometric scale. Among its numerous capabilities, the instrument can be operated as a nano-indenter to gather information about the mechanical properties of the sample. In this operating mode, the deformation of the cantilever is displayed as a function of the indentation depth of the tip into the sample. Fitting this curve with different theoretical models permits us to estimate the Young's modulus of the sample at the indentation spot. We describe what to our knowledge is a new technique to process these curves to distinguish structures of different stiffness buried into the bulk of the sample. The working principle of this new imaging technique has been verified by finite element models and successfully applied to living cells.

Magmatic fluids in the Breccia-hosted epithermal Au-Ag deposit of Rosia Montana, Romania

The breccia-hosted epithermal Au-Ag deposit of Rosia Montana is located 7 kin northeast of Abrud, in the northern part of the South Apuseni Mountains, Romania. Estimated total reserves of 214.91 million metric toils (Mt) of ore at 1.46 g/t An and 6.9 g/t Ag (10.1 Moz of An and 47.6 Moz of Ag) make Rosia Montana one of the largest gold deposits in Europe. At this location, Miocene calc-alkaline magmatic and hydrothermal activity was associated with local extensional tectonics within a strike-slip regime related to the indentation of the Adriatic microplate into the European plate during the Carpathian orogenesis. The host rocks of the magmatic complex consist of pre-Mesozoic metamorphosed continental crust covered by Cretaceous turbiditic sediment (flysch). Magmatic activity at Rosia Montana and its surroundings occurred in several pulses and lasted about 7 m.y, Rosia Montana is a breccia-hosted epithermal system related to strong phreatomagmatic activity due to the shallow emplacement of the Montana dacite. The Montana dacite intruded Miocene volcaniclastic material (volcaniclastic breccias) and crops out at Cetate and Carnic Hills. Current mining is focused primarily on the Cetate open pit, which was mapped in detail, leading to the recognition of three distinct breccia bodies: the dacite breccia with a dominantly hydrothermal matrix, the gray polymict breccia with a greater proportion of sand-sized matrix support, and the black polymict breccia, which reached to the surface, contains carbonized tree trunks and has a dominantly barren elastic matrix. The hydrothermal alteration is pervasive. Adularia alteration with a phyllic overprint is ubiquitous; silicification and argillic alteration occur locally. Mineralization consists of quartz, adularia, carbonates (commonly Mn-rich), pyrite, Fe-poor sphalerite, galena, chalcopyrite, tetrahedrite, and native gold and occurs as disseminations, as well as in veins and filling vugs within the Montana dacite and the different breccias. The age of mineralization (12.85 +/- 0.07 Ma) was determined by Ar-40- Ar-39 dating on hydrothermal adularia crystals from vugs in the dacite breccia in the Cetate open pit. Microthermometric measurements of fluid inclusions in quartz phenocrysts from the Montana dacite revealed two fluid types that are absent from the hydrothermal breccia and must have been trapped at depth prior to dacite dome emplacement: brine inclusions (32-55 -wt % NaCl equiv, homogenizing at T-h > 460 degrees C) and intermediate density fluids (4.9-15.6 wt % NaCl equiv, T, between 345 degrees-430 degrees C). Secondary aqueous fluid inclusion assemblages in the phenocrysts have salinities of 0.2 to 2.2 wt percent NaCl equiv and T-h of 200 degrees to 280 degrees C. Fluid inclusion assemblages in hydrothermal quartz from breccias and veins have salinities of 0.2 to 3.4 wt percent NaCl equiv and T-h, from 200 degrees to 270 degrees C. The oxygen isotope composition of several zones of an ore-related epithermal quartz crystal indicate a very constant delta O-18 of 4.5 to 5.0 per mil for the mineralizing fluid, despite significant salinity and temperature variation over time. Following microthermometry, selected fluid inclusion assemblages were analyzed by laser ablation-inductively coupled-plasma mass spectrometry (LA-ICMS). Despite systematic differences in salinity between phenocryst-hosted fluids trapped at depth and fluids from quartz in the epithermal breccias, all fluids have overlapping major and trace cation ratios, including identical Na/K/Rb/Sr/Cs/Ba. Consistent with the constant near-magmatic oxygen isotope composition of the hydrothermal fluids, these data strongly indicate a common magmatic component of these chemically conservative solutes in all fluids. Cu, Pb, Zn, and Mn show variations in concentration relative to the relatively non-reactive alkalis, reflecting the precipitation of sulfide minerals together with An in the epithermal breccia, and possibly of Cu in an inferred subjacent porphyry environment. The magmatic-hydrothermal processes responsible for epithermal Au-Ag mineralization at Rosia Montana are, however, not directly related to the formation of the spatially associated porphyry Cu-Au deposit of Rosia Poieni, which occurred lout 3 m.y. later.

Simulation of surface waves in porous media

We present a novel numerical algorithm for the simulation of seismic wave propagation in porous media, which is particularly suitable for the accurate modelling of surface wave-type phenomena. The differential equations of motion are based on Biot's theory of poro-elasticity and solved with a pseudospectral approach using Fourier and Chebyshev methods to compute the spatial derivatives along the horizontal and vertical directions, respectively. The time solver is a splitting algorithm that accounts for the stiffness of the differential equations. Due to the Chebyshev operator the grid spacing in the vertical direction is non-uniform and characterized by a denser spatial sampling in the vicinity of interfaces, which allows for a numerically stable and accurate evaluation of higher order surface wave modes. We stretch the grid in the vertical direction to increase the minimum grid spacing and reduce the computational cost. The free-surface boundary conditions are implemented with a characteristics approach, where the characteristic variables are evaluated at zero viscosity. The same procedure is used to model seismic wave propagation at the interface between a fluid and porous medium. In this case, each medium is represented by a different grid and the two grids are combined through a domain-decomposition method. This wavefield decomposition method accounts for the discontinuity of variables and is crucial for an accurate interface treatment. We simulate seismic wave propagation with open-pore and sealed-pore boundary conditions and verify the validity and accuracy of the algorithm by comparing the numerical simulations to analytical solutions based on zero viscosity obtained with the Cagniard-de Hoop method. Finally, we illustrate the suitability of our algorithm for more complex models of porous media involving viscous pore fluids and strongly heterogeneous distributions of the elastic and hydraulic material properties.

Effect of intradialytic resistance band exercise on physical function in patients on maintenance hemodialysis: a pilot study.

Although physical activity is recommended in patients on maintenance hemodialysis (MHD), randomized controlled trials testing the effects of exercise in this population have given conflicting results. In general, aerobic exercises mostly failed to produce improvements in physical function, whereas resistance exercises, although less studied, appeared to be more promising. The use of sophisticated materials such as leg press and free weights may preclude widespread application of resistance training in patients on MHD. Simple and cheap elastic bands may thus be an attractive alternative. We tested the feasibility of a supervised intradialytic resistance band exercise training program, and its effects on physical function, in patients on MHD. A total of 11 unselected adult patients on MHD from our center, aged 70 ± 10.7 (mean ± standard deviation) years, including 8 men and 3 women, accepted to follow the program under the supervision of qualified physiotherapists. Thirty-six exercise sessions of moderate intensity (twice a week, mean duration 40 minutes each, during 4.5 to 6 months), mainly involving leg muscles against an elastic resistance, were performed. The exercise program was well tolerated and all patients completed it. Statistically significant improvements were observed in the following tests: Tinetti test, 23.9 ± 3.9 points before versus 25.7 ± 3.5 points after the program (P = .022); the Timed Up and Go test, 12.1 ± 6.6 versus 10 ± 5.8 seconds (P = .0156). Improvements in the 6-minute walk distance and in the one-leg balance tests just failed to reach statistical significance. In this single-center pilot study, an intradialytic resistance band exercise program was feasible, well tolerated, and showed encouraging results on physical function.

The clinical use of quantitative ultrasound (QUS) in the detection and management of osteoporosis

For the detection and management of osteoporosis and osteoporosis-related fractures, quantitative ultrasound (QUS) is emerging as a relatively low-cost and readily accessible alternative to dual-energy X-ray absorptiometry (DXA) measurement of bone mineral density (BMD) in certain circumstances. The following is a brief, but thorough review of the existing literature with respect to the use of QUS in 6 settings: 1) assessing fragility fracture risk; 2) diagnosing osteoporosis; 3) initiating osteoporosis treatment; 4) monitoring osteoporosis treatment; 5) osteoporosis case finding; and 6) quality assurance and control. Many QUS devices exist that are quite different with respect to the parameters they measure and the strength of empirical evidence supporting their use. In general, heel QUS appears to be most tested and most effective. Overall, some, but not all, heel QUS devices are effective assessing fracture risk in some, but not all, populations, the evidence being strongest for Caucasian females over 55 years old. Otherwise, the evidence is fair with respect to certain devices allowing for the accurate diagnosis of likelihood of osteoporosis, and generally fair to poor in terms of QUS use when initiating or monitoring osteoporosis treatment. A reasonable protocol is proposed herein for case-finding purposes, which relies on a combined assessment of clinical risk factors (CR.F) and heel QUS. Finally, several recommendations are made for quality assurance and control.

Comparison of aortic elasticity determined by cardiovascular magnetic resonance imaging in obese versus lean adults.

The vascular properties of large vessels in the obese have not been adequately studied. We used cardiovascular magnetic resonance imaging to quantify the cross-sectional area and elastic properties of the ascending thoracic and abdominal aorta in 21 clinically healthy obese young adult men and 25 men who were age-matched lean controls. Obese subjects had greater maximal cross-sectional area of the ascending thoracic aorta (984 +/- 252 vs 786 +/- 109 mm(2), p <0.01) and of the abdominal aorta (415 +/- 71 vs 374 +/- 51 mm(2), p <0.05). When indexed for height the differences persisted, but when indexed for body surface area, a significant difference between groups was found only for the maximal abdominal aortic cross-sectional area. The obese subjects also had decreased abdominal aortic elasticity, characterized by 24% lower compliance (0.0017 +/- 0.0004 vs 0.0021 +/- 0.0005 mm(2)/kPa/mm, p <0.01), 22% higher stiffness index beta (6.0 +/- 1.5 vs 4.9 +/- 0.7, p <0.005), and 41% greater pressure-strain elastic modulus (72 +/- 25 vs 51 +/- 9, p <0.005). At the ascending thoracic aorta, only the pressure-strain elastic modulus was different between obese and lean subjects (85 +/- 42 vs 65 +/- 26 kPa, respectively; p <0.05), corresponding to a 31% difference-but arterial compliance and stiffness index were not significantly different between groups. In clinically healthy young adult obese men, obesity is associated with increased cross-sectional aortic area and decreased aortic elasticity.

High-precision dating and origin of synsedimentary volcanism in the Late Carboniferous Salvan-Dorenaz basin (Aiguilles-Rouges Massif, Western Alps)

Late Variscan volcanic activity is documented in the Late Carboniferous Salvan-Dorenaz sedimentary basin and in the neighboring basement units of the Aiguilles-Rouges and Mont-Blanc crystalline massifs (Western Alps). Precise U/Pb isotopic dating, zircon morphology and geochemical analyses indicate that volcanism occurred during short-lived pulses and that coexisting crustal and mantle sources were involved in the production of melts. Volcanic and subvolcanic products were emplaced along major N-S to NNE-SSW transtensional fracture zones, similar to the ones that governed intense basement exhumation and that favored the formation and filling of the Late Carboniferous Salvan-Dorenaz continental basin. In the Aiguilles-Rouges massif, dacitic flows outcropping at the base of the Salvan-Dorenaz basin erupted at 308 +/- 3 Ma; they represent the surface equivalent of the nearby Vallorcine peraluminous granite and associated rhyolitic dykes (311 +/- 17 Ma). In the Mont Blanc massif, calc-alkaline rhyolitic dykes were emplaced simultaneously (307 +/- 2 Ma) at shallow crustal levels, but they derive from deeper magma sources denoting enhanced mantellic activity. Recently identified tuffs and volcaniclastic layers embedded at different levels of the Salvan-Dorenaz stratigraphic record testify a 295 +3/-4 Ma old episode of highly explosive volcanism from distant volcanic centers, possibly located in the Aar-Gotthard massifs (Central Alps). Their zircon typology is highly heterogeneous. documenting wall-rock contamination of the melts and/or admixture of crustal sediments, whereas consistent subpopulations point to high-temperature magmas of deep-seated origin and alkaline affinity. The dated volcanic layers from the Salvan-Dorenaz basin set the beginning of the detrital sedimentation at 308 +/- 3 Ma and constrain the deposition of 1.5-1.7 km thick of elastic sediments within a time span of 10-15 Ma. These results infer minimum, long-term subsidence rates during basin evolution in the order of >0.1 mm/a, while in the surrounding basement units estimated exhumation rates are in the range of 1 mm/a. All dated rocks contain inherited zircon populations about 350, 450 or 600 Ma old.

Enterocystoplasty using modified pedicled, detubularized, de-epithelialized sigmoid patches in the mini-pig model.

Enterocystoplasty and gastrocystoplasty have been developed to restore adequate bladder capacity but they still result in the undesired harmful contact of urine with the intestinal mucosa. In an attempt to prevent this nonphysiological interface we performed several modified enterocystoplasties in the mini-pig model using a pedicled, detubularized, de-epithelialized sigmoid patch. Five techniques of patch coverage were used to evaluate urothelial growth or survival on the sigmoid patch. All but 1 patch had intestinal mucosa remnants with mucocele formation. When no coverage was applied or a biodegradable polyglactin mesh was temporarily covering the pedicled, detubularized, de-epithelialized sigmoid patch, severe shrinkage occurred. Partial cover with autologous urothelium islets seemed to decrease shrinkage. Adequate urothelium survival and satisfactory elastic properties of the patch were observed when total cover of the sigmoid patch was achieved with a sheet of homologous urothelium recovering autologous urothelial islets or when the patch was applied to protruding urothelium obtained following sagittal posterior detrusor myotomy.

Incidence and management of pulmonary embolism following spinal surgery occurring while under chemical thromboprophylaxis.

Patients undergoing spinal surgery are at risk of developing thromboembolic complications even though lower incidences have been reported as compared to joint arthroplasty surgery. Deep vein thrombosis (DVT) has been studied extensively in the context of spinal surgery but symptomatic pulmonary embolism (PE) has engaged less attention. We prospectively followed a consecutive cohort of 270 patients undergoing spinal surgery at a single institution. From these patients, only 26 were simple discectomies, while the largest proportion (226) was fusions. All patients received both low molecular weight heparin (LMWH) initiated after surgery and compressive stockings. PE was diagnosed with spiral chest CT. Six patients developed symptomatic PE, five during their hospital stay. In three of the six patients the embolic event occurred during the first 3 postoperative days. They were managed by the temporary insertion of an inferior vena cava (IVC) filter thus allowing for a delay in full-dose anticoagulation until removal of the filter. None of the PE patients suffered any bleeding complication as a result of the introduction of full anticoagulation. Two patients suffered postoperative haematomas, without development of neurological symptoms or signs, requiring emergency evacuation. The overall incidence of PE was 2.2% rising to 2.5% after exclusion of microdiscectomy cases. The incidence of PE was highest in anterior or combined thoracolumbar/lumbar procedures (4.2%). There is a large variation in the reported incidence of PE in the spinal literature. Results from the only study found in the literature specifically monitoring PE suggest an incidence of PE as high as 2.5%. Our study shows a similar incidence despite the use of LMWH. In the absence of randomized controlled trials (RCT) it is uncertain if this type of prophylaxis lowers the incidence of PE. However, other studies show that the morbidity of LMWH is very low. Since PE can be a life-threatening complication, LMWH may be a worthwhile option to consider for prophylaxis. RCTs are necessary in assessing the efficacy of DVT and PE prophylaxis in spinal patients.

Assessing dystrophies and other muscle diseases at the nanometer scale by atomic force microscopy.

AIM: Atomic force microscopy nanoindentation of myofibers was used to assess and quantitatively diagnose muscular dystrophies from human patients. MATERIALS & METHODS: Myofibers were probed from fresh or frozen muscle biopsies from human dystrophic patients and healthy volunteers, as well as mice models, and Young's modulus stiffness values were determined. RESULTS: Fibers displaying abnormally low mechanical stability were detected in biopsies from patients affected by 11 distinct muscle diseases, and Young's modulus values were commensurate to the severity of the disease. Abnormal myofiber resistance was also observed from consulting patients whose muscle condition could not be detected or unambiguously diagnosed otherwise. DISCUSSION & CONCLUSION: This study provides a proof-of-concept that atomic force microscopy yields a quantitative read-out of human muscle function from clinical biopsies, and that it may thereby complement current muscular dystrophy diagnosis.

Thoracic aortic dissection associated with cocaine abuse.

Cardiovascular complications of cocaine abuse include myocardial ischemia and infarction, dysrhythmias, cardiomyopathies and aortic dissection. The case in point pertains to a 26-year-old, Caucasian male, substance abuser who suffered a thoracic aortic dissection following the use of crack cocaine. The autopsy and histological findings showed a connective tissue abnormality including a focal microcystic medial necrosis and a fragmentation of the elastic fibers in the arterial walls. Blood concentrations of cocaine and benzoylecgonine, taken individually, were considered to be within a potentially toxic range. Blood concentrations of methadone also indicated use of this drug at the same time. The small amounts of morphine found in the blood and urine were compatible with heroine or morphine use more than 24 h before death.

Three-dimensional hydrodynamic numerical modelling of fold nappe formation, basement-cover deformation and slab detachment with applications to the helvetic nappe system (W Switzerland)

Many three-dimensional (3-D) structures in rock, which formed during the deformation of the Earth's crust and lithosphere, are controlled by a difference in mechanical strength between rock units and are often the result of a geometrical instability. Such structures are, for example, folds, pinch-and-swell structures (due to necking) or cuspate-lobate structures (mullions). These struc-tures occur from the centimeter to the kilometer scale and the related deformation processes con-trol the formation of, for example, fold-and-thrust belts and extensional sedimentary basins or the deformation of the basement-cover interface. The 2-D deformation processes causing these structures are relatively well studied, however, several processes during large-strain 3-D defor-mation are still incompletely understood. One of these 3-D processes is the lateral propagation of these structures, such as fold and cusp propagation in a direction orthogonal to the shortening direction or neck propagation in direction orthogonal to the extension direction. Especially, we are interested in fold nappes which are recumbent folds with amplitudes usually exceeding 10 km and they have been presumably formed by ductile shearing. They often exhibit a constant sense of shearing and a non-linear increase of shear strain towards their overturned limb. The fold axes of the Morcles fold nappe in western Switzerland plunges to the ENE whereas the fold axes in the more eastern Doldenhorn nappe plunges to the WSW. These opposite plunge direc-tions characterize the Rawil depression (Wildstrubel depression). The Morcles nappe is mainly the result of layer parallel contraction and shearing. During the compression the massive lime-stones were more competent than the surrounding marls and shales, which led to the buckling characteristics of the Morcles nappe, especially in the north-dipping normal limb. The Dolden-horn nappe exhibits only a minor overturned fold limb. There are still no 3-D numerical studies which investigate the fundamental dynamics of the formation of the large-scale 3-D structure including the Morcles and Doldenhorn nappes and the related Rawil depression. We study the 3-D evolution of geometrical instabilities and fold nappe formation with numerical simulations based on the finite element method (FEM). Simulating geometrical instabilities caused by sharp variations of mechanical strength between rock units requires a numerical algorithm that can accurately resolve material interfaces for large differences in material properties (e.g. between limestone and shale) and for large deformations. Therefore, our FE algorithm combines a nu-merical contour-line technique and a deformable Lagrangian mesh with re-meshing. With this combined method it is possible to accurately follow the initial material contours with the FE mesh and to accurately resolve the geometrical instabilities. The algorithm can simulate 3-D de-formation for a visco-elastic rheology. The viscous rheology is described by a power-law flow law. The code is used to study the 3-D fold nappe formation, the lateral propagation of folding and also the lateral propagation of cusps due to initial half graben geometry. Thereby, the small initial geometrical perturbations for folding and necking are exactly followed by the FE mesh, whereas the initial large perturbation describing a half graben is defined by a contour line inter-secting the finite elements. Further, the 3-D algorithm is applied to 3-D viscous nacking during slab detachment. The results from various simulations are compared with 2-D resulats and a 1-D analytical solution. -- On retrouve beaucoup de structures en 3 dimensions (3-D) dans les roches qui ont pour origines une déformation de la lithosphère terrestre. Ces structures sont par exemple des plis, des boudins (pinch-and-swell) ou des mullions (cuspate-lobate) et sont présentés de l'échelle centimétrique à kilométrique. Mécaniquement, ces structures peuvent être expliquées par une différence de résistance entre les différentes unités de roches et sont généralement le fruit d'une instabilité géométrique. Ces différences mécaniques entre les unités contrôlent non seulement les types de structures rencontrées, mais également le type de déformation (thick skin, thin skin) et le style tectonique (bassin d'avant pays, chaîne d'avant pays). Les processus de la déformation en deux dimensions (2-D) formant ces structures sont relativement bien compris. Cependant, lorsque l'on ajoute la troisiéme dimension, plusieurs processus ne sont pas complètement compris lors de la déformation à large échelle. L'un de ces processus est la propagation latérale des structures, par exemple la propagation de plis ou de mullions dans la direction perpendiculaire à l'axe de com-pression, ou la propagation des zones d'amincissement des boudins perpendiculairement à la direction d'extension. Nous sommes particulièrement intéressés les nappes de plis qui sont des nappes de charriage en forme de plis couché d'une amplitude plurikilométrique et étant formées par cisaillement ductile. La plupart du temps, elles exposent un sens de cisaillement constant et une augmentation non linéaire de la déformation vers la base du flanc inverse. Un exemple connu de nappes de plis est le domaine Helvétique dans les Alpes de l'ouest. Une de ces nap-pes est la Nappe de Morcles dont l'axe de pli plonge E-NE tandis que de l'autre côté de la dépression du Rawil (ou dépression du Wildstrubel), la nappe du Doldenhorn (équivalent de la nappe de Morcles) possède un axe de pli plongeant O-SO. La forme particulière de ces nappes est due à l'alternance de couches calcaires mécaniquement résistantes et de couches mécanique-ment faibles constituées de schistes et de marnes. Ces différences mécaniques dans les couches permettent d'expliquer les plissements internes à la nappe, particulièrement dans le flanc inver-se de la nappe de Morcles. Il faut également noter que le développement du flanc inverse des nappes n'est pas le même des deux côtés de la dépression de Rawil. Ainsi la nappe de Morcles possède un important flanc inverse alors que la nappe du Doldenhorn en est presque dépour-vue. A l'heure actuelle, aucune étude numérique en 3-D n'a été menée afin de comprendre la dynamique fondamentale de la formation des nappes de Morcles et du Doldenhorn ainsi que la formation de la dépression de Rawil. Ce travail propose la première analyse de l'évolution 3-D des instabilités géométriques et de la formation des nappes de plis en utilisant des simulations numériques. Notre modèle est basé sur la méthode des éléments finis (FEM) qui permet de ré-soudre avec précision les interfaces entre deux matériaux ayant des propriétés mécaniques très différentes (par exemple entre les couches calcaires et les couches marneuses). De plus nous utilisons un maillage lagrangien déformable avec une fonction de re-meshing (production d'un nouveau maillage). Grâce à cette méthode combinée il nous est possible de suivre avec précisi-on les interfaces matérielles et de résoudre avec précision les instabilités géométriques lors de la déformation de matériaux visco-élastiques décrit par une rhéologie non linéaire (n>1). Nous uti-lisons cet algorithme afin de comprendre la formation des nappes de plis, la propagation latérale du plissement ainsi que la propagation latérale des structures de type mullions causé par une va-riation latérale de la géométrie (p.ex graben). De plus l'algorithme est utilisé pour comprendre la dynamique 3-D de l'amincissement visqueux et de la rupture de la plaque descendante en zone de subduction. Les résultats obtenus sont comparés à des modèles 2-D et à la solution analytique 1-D. -- Viele drei dimensionale (3-D) Strukturen, die in Gesteinen vorkommen und durch die Verfor-mung der Erdkruste und Litosphäre entstanden sind werden von den unterschiedlichen mechani-schen Eigenschaften der Gesteinseinheiten kontrolliert und sind häufig das Resulat von geome-trischen Istabilitäten. Zu diesen strukturen zählen zum Beispiel Falten, Pich-and-swell Struktu-ren oder sogenannte Cusbate-Lobate Strukturen (auch Mullions). Diese Strukturen kommen in verschiedenen Grössenordungen vor und können Masse von einigen Zentimeter bis zu einigen Kilometer aufweisen. Die mit der Entstehung dieser Strukturen verbundenen Prozesse kontrol-lieren die Entstehung von Gerbirgen und Sediment-Becken sowie die Verformung des Kontaktes zwischen Grundgebirge und Stedimenten. Die zwei dimensionalen (2-D) Verformungs-Prozesse die zu den genannten Strukturen führen sind bereits sehr gut untersucht. Einige Prozesse wäh-rend starker 3-D Verformung sind hingegen noch unvollständig verstanden. Einer dieser 3-D Prozesse ist die seitliche Fortpflanzung der beschriebenen Strukturen, so wie die seitliche Fort-pflanzung von Falten und Cusbate-Lobate Strukturen senkrecht zur Verkürzungsrichtung und die seitliche Fortpflanzung von Pinch-and-Swell Strukturen othogonal zur Streckungsrichtung. Insbesondere interessieren wir uns für Faltendecken, liegende Falten mit Amplituden von mehr als 10 km. Faltendecken entstehen vermutlich durch duktile Verscherung. Sie zeigen oft einen konstanten Scherungssinn und eine nicht-lineare zunahme der Scherverformung am überkipp-ten Schenkel. Die Faltenachsen der Morcles Decke in der Westschweiz fallen Richtung ONO während die Faltenachsen der östicher gelegenen Doldenhorn Decke gegen WSW einfallen. Diese entgegengesetzten Einfallrichtungen charakterisieren die Rawil Depression (Wildstrubel Depression). Die Morcles Decke ist überwiegend das Resultat von Verkürzung und Scherung parallel zu den Sedimentlagen. Während der Verkürzung verhielt sich der massive Kalkstein kompetenter als der Umliegende Mergel und Schiefer, was zur Verfaltetung Morcles Decke führ-te, vorallem in gegen Norden eifallenden überkippten Schenkel. Die Doldenhorn Decke weist dagegen einen viel kleineren überkippten Schenkel und eine stärkere Lokalisierung der Verfor-mung auf. Bis heute gibt es keine 3-D numerischen Studien, die die fundamentale Dynamik der Entstehung von grossen stark verformten 3-D Strukturen wie den Morcles und Doldenhorn Decken sowie der damit verbudenen Rawil Depression untersuchen. Wir betrachten die 3-D Ent-wicklung von geometrischen Instabilitäten sowie die Entstehung fon Faltendecken mit Hilfe von numerischen Simulationen basiert auf der Finite Elemente Methode (FEM). Die Simulation von geometrischen Instabilitäten, die aufgrund von Änderungen der Materialeigenschaften zwischen verschiedenen Gesteinseinheiten entstehen, erfortert einen numerischen Algorithmus, der in der Lage ist die Materialgrenzen mit starkem Kontrast der Materialeigenschaften (zum Beispiel zwi-schen Kalksteineinheiten und Mergel) für starke Verfomung genau aufzulösen. Um dem gerecht zu werden kombiniert unser FE Algorithmus eine numerische Contour-Linien-Technik und ein deformierbares Lagranges Netz mit Re-meshing. Mit dieser kombinierten Methode ist es mög-lich den anfänglichen Materialgrenzen mit dem FE Netz genau zu folgen und die geometrischen Instabilitäten genügend aufzulösen. Der Algorithmus ist in der Lage visko-elastische 3-D Ver-formung zu rechnen, wobei die viskose Rheologie mit Hilfe eines power-law Fliessgesetzes beschrieben wird. Mit dem numerischen Algorithmus untersuchen wir die Entstehung von 3-D Faltendecken, die seitliche Fortpflanzung der Faltung sowie der Cusbate-Lobate Strukturen die sich durch die Verkürzung eines mit Sediment gefüllten Halbgraben bilden. Dabei werden die anfänglichen geometrischen Instabilitäten der Faltung exakt mit dem FE Netz aufgelöst wäh-rend die Materialgranzen des Halbgrabens die Finiten Elemente durchschneidet. Desweiteren wird der 3-D Algorithmus auf die Einschnürung während der 3-D viskosen Plattenablösung und Subduktion angewandt. Die 3-D Resultate werden mit 2-D Ergebnissen und einer 1-D analyti-schen Lösung verglichen.

Continuing Medical Education Self-Assessment Questions Vol. 50, No. 1 - January 2014. Actual evidence in the medical approach to adolescents with idiopathic scoliosis

Idiopathic scoliosis (IS) is a three-dimensional deformity of the spine and trunk. The most common form involve ado- lescents (AIS). The prevalence for AIS is 2-3% of the population, with 1 out of 6 patients requiring treatment of which 25% progress to surgery. Physical and rehabilitation medicine (PRM) plays a primary role in the so-called conservative treatment of adolescents with AIS, since all the therapeutic tools used (exercises and braces) fall into the PRM domain. According to a Cochrane systematic review there is evidence in favor of bracing, even if it is of low quality. Another shows that there is evidence in favor of exercises as an adjunctive treatment, but of low quality. Three meta-analysis have been published on bracing: one shows that bracing does not reduce surgery rates, but studies with bracing plus exercises were not included and had the highest effectiveness; another shows that full time is better than part-time bracing; the last focuses on observational studies following the SRS criteria and shows that not all full time rigid bracing are the same: some have the highest effectiveness, others have less than elastic and nighttime bracing. Two very important RCTs failed in recruitment, showing that in the field of bracing for scoliosis RCTs are not accepted by the patients. Consensuses by the international Society on Scoliosis Orthopedic and Rehabilitation Treatment (SOSORT) show that there is no agree- ment among experts either on the best braces or on their biomechanical action, and that compliance is a matter of clinical more than patients' behavior (there is strong agreement on the management criteria to achieve best results with bracing). A systematic review of all the existing studies shows effectiveness of exercises, and that auto-correction is the main goal of exercises. A systematic review shows that there are no studies on manual treatment. Research on conservative treat- ment of AIS has continuously decreased since the 1980s, but this trend changed only recently. The SOSORT Guidelines offers the actual standard of conservative care.

Seismic characterization of heterogeneous sedimentary and fractured rocks based on Biot's theory of poroelasticity

Understanding and quantifying seismic energy dissipation, which manifests itself in terms of velocity dispersion and attenuation, in fluid-saturated porous rocks is of considerable interest, since it offers the perspective of extracting information with regard to the elastic and hydraulic rock properties. There is increasing evidence to suggest that wave-induced fluid flow, or simply WIFF, is the dominant underlying physical mechanism governing these phenomena throughout the seismic, sonic, and ultrasonic frequency ranges. This mechanism, which can prevail at the microscopic, mesoscopic, and macroscopic scale ranges, operates through viscous energy dissipation in response to fluid pressure gradients and inertial effects induced by the passing wavefield. In the first part of this thesis, we present an analysis of broad-band multi-frequency sonic log data from a borehole penetrating water-saturated unconsolidated glacio-fluvial sediments. An inherent complication arising in the interpretation of the observed P-wave attenuation and velocity dispersion is, however, that the relative importance of WIFF at the various scales is unknown and difficult to unravel. An important generic result of our work is that the levels of attenuation and velocity dispersion due to the presence of mesoscopic heterogeneities in water-saturated unconsolidated clastic sediments are expected to be largely negligible. Conversely, WIFF at the macroscopic scale allows for explaining most of the considered data while refinements provided by including WIFF at the microscopic scale in the analysis are locally meaningful. Using a Monte-Carlo-type inversion approach, we compare the capability of the different models describing WIFF at the macroscopic and microscopic scales with regard to their ability to constrain the dry frame elastic moduli and the permeability as well as their local probability distribution. In the second part of this thesis, we explore the issue of determining the size of a representative elementary volume (REV) arising in the numerical upscaling procedures of effective seismic velocity dispersion and attenuation of heterogeneous media. To this end, we focus on a set of idealized synthetic rock samples characterized by the presence of layers, fractures or patchy saturation in the mesocopic scale range. These scenarios are highly pertinent because they tend to be associated with very high levels of velocity dispersion and attenuation caused by WIFF in the mesoscopic scale range. The problem of determining the REV size for generic heterogeneous rocks is extremely complex and entirely unexplored in the given context. In this pilot study, we have therefore focused on periodic media, which assures the inherent self- similarity of the considered samples regardless of their size and thus simplifies the problem to a systematic analysis of the dependence of the REV size on the applied boundary conditions in the numerical simulations. Our results demonstrate that boundary condition effects are absent for layered media and negligible in the presence of patchy saturation, thus resulting in minimum REV sizes. Conversely, strong boundary condition effects arise in the presence of a periodic distribution of finite-length fractures, thus leading to large REV sizes. In the third part of the thesis, we propose a novel effective poroelastic model for periodic media characterized by mesoscopic layering, which accounts for WIFF at both the macroscopic and mesoscopic scales as well as for the anisotropy associated with the layering. Correspondingly, this model correctly predicts the existence of the fast and slow P-waves as well as quasi and pure S-waves for any direction of wave propagation as long as the corresponding wavelengths are much larger than the layer thicknesses. The primary motivation for this work is that, for formations of intermediate to high permeability, such as, for example, unconsolidated sediments, clean sandstones, or fractured rocks, these two WIFF mechanisms may prevail at similar frequencies. This scenario, which can be expected rather common, cannot be accounted for by existing models for layered porous media. Comparisons of analytical solutions of the P- and S-wave phase velocities and inverse quality factors for wave propagation perpendicular to the layering with those obtained from numerical simulations based on a ID finite-element solution of the poroelastic equations of motion show very good agreement as long as the assumption of long wavelengths remains valid. A limitation of the proposed model is its inability to account for inertial effects in mesoscopic WIFF when both WIFF mechanisms prevail at similar frequencies. Our results do, however, also indicate that the associated error is likely to be relatively small, as, even at frequencies at which both inertial and scattering effects are expected to be at play, the proposed model provides a solution that is remarkably close to its numerical benchmark. -- Comprendre et pouvoir quantifier la dissipation d'énergie sismique qui se traduit par la dispersion et l'atténuation des vitesses dans les roches poreuses et saturées en fluide est un intérêt primordial pour obtenir des informations à propos des propriétés élastique et hydraulique des roches en question. De plus en plus d'études montrent que le déplacement relatif du fluide par rapport au solide induit par le passage de l'onde (wave induced fluid flow en anglais, dont on gardera ici l'abréviation largement utilisée, WIFF), représente le principal mécanisme physique qui régit ces phénomènes, pour la gamme des fréquences sismiques, sonique et jusqu'à l'ultrasonique. Ce mécanisme, qui prédomine aux échelles microscopique, mésoscopique et macroscopique, est lié à la dissipation d'énergie visqueuse résultant des gradients de pression de fluide et des effets inertiels induits par le passage du champ d'onde. Dans la première partie de cette thèse, nous présentons une analyse de données de diagraphie acoustique à large bande et multifréquences, issues d'un forage réalisé dans des sédiments glaciaux-fluviaux, non-consolidés et saturés en eau. La difficulté inhérente à l'interprétation de l'atténuation et de la dispersion des vitesses des ondes P observées, est que l'importance des WIFF aux différentes échelles est inconnue et difficile à quantifier. Notre étude montre que l'on peut négliger le taux d'atténuation et de dispersion des vitesses dû à la présence d'hétérogénéités à l'échelle mésoscopique dans des sédiments clastiques, non- consolidés et saturés en eau. A l'inverse, les WIFF à l'échelle macroscopique expliquent la plupart des données, tandis que les précisions apportées par les WIFF à l'échelle microscopique sont localement significatives. En utilisant une méthode d'inversion du type Monte-Carlo, nous avons comparé, pour les deux modèles WIFF aux échelles macroscopique et microscopique, leur capacité à contraindre les modules élastiques de la matrice sèche et la perméabilité ainsi que leur distribution de probabilité locale. Dans une seconde partie de cette thèse, nous cherchons une solution pour déterminer la dimension d'un volume élémentaire représentatif (noté VER). Cette problématique se pose dans les procédures numériques de changement d'échelle pour déterminer l'atténuation effective et la dispersion effective de la vitesse sismique dans un milieu hétérogène. Pour ce faire, nous nous concentrons sur un ensemble d'échantillons de roches synthétiques idéalisés incluant des strates, des fissures, ou une saturation partielle à l'échelle mésoscopique. Ces scénarios sont hautement pertinents, car ils sont associés à un taux très élevé d'atténuation et de dispersion des vitesses causé par les WIFF à l'échelle mésoscopique. L'enjeu de déterminer la dimension d'un VER pour une roche hétérogène est très complexe et encore inexploré dans le contexte actuel. Dans cette étude-pilote, nous nous focalisons sur des milieux périodiques, qui assurent l'autosimilarité des échantillons considérés indépendamment de leur taille. Ainsi, nous simplifions le problème à une analyse systématique de la dépendance de la dimension des VER aux conditions aux limites appliquées. Nos résultats indiquent que les effets des conditions aux limites sont absents pour un milieu stratifié, et négligeables pour un milieu à saturation partielle : cela résultant à des dimensions petites des VER. Au contraire, de forts effets des conditions aux limites apparaissent dans les milieux présentant une distribution périodique de fissures de taille finie : cela conduisant à de grandes dimensions des VER. Dans la troisième partie de cette thèse, nous proposons un nouveau modèle poro- élastique effectif, pour les milieux périodiques caractérisés par une stratification mésoscopique, qui prendra en compte les WIFF à la fois aux échelles mésoscopique et macroscopique, ainsi que l'anisotropie associée à ces strates. Ce modèle prédit alors avec exactitude l'existence des ondes P rapides et lentes ainsi que les quasis et pures ondes S, pour toutes les directions de propagation de l'onde, tant que la longueur d'onde correspondante est bien plus grande que l'épaisseur de la strate. L'intérêt principal de ce travail est que, pour les formations à perméabilité moyenne à élevée, comme, par exemple, les sédiments non- consolidés, les grès ou encore les roches fissurées, ces deux mécanismes d'WIFF peuvent avoir lieu à des fréquences similaires. Or, ce scénario, qui est assez commun, n'est pas décrit par les modèles existants pour les milieux poreux stratifiés. Les comparaisons des solutions analytiques des vitesses des ondes P et S et de l'atténuation de la propagation des ondes perpendiculaires à la stratification, avec les solutions obtenues à partir de simulations numériques en éléments finis, fondées sur une solution obtenue en 1D des équations poro- élastiques, montrent un très bon accord, tant que l'hypothèse des grandes longueurs d'onde reste valable. Il y a cependant une limitation de ce modèle qui est liée à son incapacité à prendre en compte les effets inertiels dans les WIFF mésoscopiques quand les deux mécanismes d'WIFF prédominent à des fréquences similaires. Néanmoins, nos résultats montrent aussi que l'erreur associée est relativement faible, même à des fréquences à laquelle sont attendus les deux effets d'inertie et de diffusion, indiquant que le modèle proposé fournit une solution qui est remarquablement proche de sa référence numérique.