Procedure for human saphenous veins ex vivo perfusion and external reinforcement.

The mainstay of contemporary therapies for extensive occlusive arterial disease is venous bypass graft. However, its durability is threatened by intimal hyperplasia (IH) that eventually leads to vessel occlusion and graft failure. Mechanical forces, particularly low shear stress and high wall tension, are thought to initiate and to sustain these cellular and molecular changes, but their exact contribution remains to be unraveled. To selectively evaluate the role of pressure and shear stress on the biology of IH, an ex vivo perfusion system (EVPS) was created to perfuse segments of human saphenous veins under arterial regimen (high shear stress and high pressure). Further technical innovations allowed the simultaneous perfusion of two segments from the same vein, one reinforced with an external mesh. Veins were harvested using a no-touch technique and immediately transferred to the laboratory for assembly in the EVPS. One segment of the freshly isolated vein was not perfused (control, day 0). The two others segments were perfused for up to 7 days, one being completely sheltered with a 4 mm (diameter) external mesh. The pressure, flow velocity, and pulse rate were continuously monitored and adjusted to mimic the hemodynamic conditions prevailing in the femoral artery. Upon completion of the perfusion, veins were dismounted and used for histological and molecular analysis. Under ex vivo conditions, high pressure perfusion (arterial, mean = 100 mm Hg) is sufficient to generate IH and remodeling of human veins. These alterations are reduced in the presence of an external polyester mesh.

A multilevel multiscale finite volume method

The Multiscale Finite Volume (MsFV) method has been developed to efficiently solve reservoir-scale problems while conserving fine-scale details. The method employs two grid levels: a fine grid and a coarse grid. The latter is used to calculate a coarse solution to the original problem, which is interpolated to the fine mesh. The coarse system is constructed from the fine-scale problem using restriction and prolongation operators that are obtained by introducing appropriate localization assumptions. Through a successive reconstruction step, the MsFV method is able to provide an approximate, but fully conservative fine-scale velocity field. For very large problems (e.g. one billion cell model), a two-level algorithm can remain computational expensive. Depending on the upscaling factor, the computational expense comes either from the costs associated with the solution of the coarse problem or from the construction of the local interpolators (basis functions). To ensure numerical efficiency in the former case, the MsFV concept can be reapplied to the coarse problem, leading to a new, coarser level of discretization. One challenge in the use of a multilevel MsFV technique is to find an efficient reconstruction step to obtain a conservative fine-scale velocity field. In this work, we introduce a three-level Multiscale Finite Volume method (MlMsFV) and give a detailed description of the reconstruction step. Complexity analyses of the original MsFV method and the new MlMsFV method are discussed, and their performances in terms of accuracy and efficiency are compared.

Rockfall characterisation and structural protection: A review

Rockfall is an extremely rapid process involving long travel distances. Due to these features, when an event occurs, the ability to take evasive action is practically zero and, thus, the risk of injury or loss of life is high. Damage to buildings and infrastructure is quite likely. In many cases, therefore, suitable protection measures are necessary. This contribution provides an overview of previous and current research on the main topics related to rockfall. It covers the onset of rockfall and runout modelling approaches, as well as hazard zoning and protection measures. It is the aim of this article to provide an in-depth knowledge base for researchers and practitioners involved in projects dealing with the rockfall protection of infrastructures, who may work in the fields of civil or environmental engineering, risk and safety, the earth and natural sciences.

The economic burden of alcohol dependence in europe.

AIMS: To determine the economic burden pertaining to alcohol dependence in Europe. METHODS: Database searching was combined with grey literature searching to identify costs and resource use in Europe relating to alcohol dependence as defined by the Diagnostic and Statistical Manual of Mental Disorders (DSM-IV) or the World Health Organisation's International Classification of Diseases (ICD-10). Searches combined MeSH headings for both economic terms and terms pertaining to alcohol dependence. Relevant outcomes included direct healthcare costs and indirect societal costs. Main resource use outcomes included hospitalization and drug costs. RESULTS: Compared with the number of studies of the burden of alcohol use disorders in general, relatively few focussed specifically on alcohol dependence. Twenty-two studies of variable quality were eligible for inclusion. The direct costs of alcohol dependence in Europe were substantial, the treatment costs for a single alcohol-dependent patient lying within the range euro1591-euro7702 per hospitalization and the annual total direct costs accounting for 0.04-0.31% of an individual country's gross domestic product (GDP). These costs were driven primarily by hospitalization; in contrast, the annual drug costs for alcohol dependence were low. The indirect costs were more substantial than the direct costs, accounting for up to 0.64% of GDP per country annually. Alcohol dependence may be more costly in terms of health costs per patient than alcohol abuse. CONCLUSIONS: This review confirms that alcohol dependence represents a significant burden for European healthcare systems and society. Difficulties in comparing across cost-of-illness studies in this disease area, however, prevent specific estimation of the economic burden.

Organisation and role of actin microfilaments during cellular growth and morphogenesis in the moss Physcomitrella patens

ABSTRACT The network of actin cytoskeleton is composed of actin filaments (F-actin) that are made by polymerisation of actin monomers and actin binding proteins. It is required for growth and morphogenesis of eukaryotic cells. The labelling of F-actin with constitutively expressed GFP-Talin (Kost et al., 1998) reveals the organisation of cellular actin networks in plants. Due to the lack of information on actin cytoskeleton through gametophytic development of the model moss plant Physcornitrella patens, stable transgenic lines overexpressing GFP-Talin were generated to detect F-actin structures. It is shown that the 35S promoter driven expression is not suitable for F-actin labelling in all cells. When it is replaced by the inducible heat-shock promoter Gmhsp17.3 from soybean, one hour mild heat stress at 37°C followed by recovery at 25°C is enough to induce efficient and transient labelling in all tissues without altering cellular morphology. The optimal observations of F-actin structures at different stages of moss development can be done between 12-18 hours after the induction. By using confocal microscopy, we demonstrate that stellated actin arrays were densely accumulated at the growing tip in regenerating protoplasts, apical protonemal cells and rhizoids and connected with a fine dispersed F-actin mesh. Following three-dimensional growth, the cortical star-like structures are widespread in the meristematic cells of developing bud and young gametophores. On the contrary, undulating networks of actin cables are found at the final stage of cell differentiation. During redifferentiation of mature leaf cells into protonemal filaments the rather stagnant web of actin cables is replaced by diffuse actin meshwork. In eukaryotes, nucleation of the actin monomers prior to their polymerization is driven by the seven-subunit ARP2/3 complex and formins. We cloned the gene encoding the ARP3 subunit of P. patens and generated arp3 mutants of the moss through gene disruption. The knockout of ARP3 affects the elongation of chloronemal cells and blocks further differentiation of caulonemal cells and rhizoids, and the gametophores are slightly stunted compared to wild-type. The arp mutants were created in the heat-shock inducible GFP-Talin strains allowing us to visualise a disorganised actin network and a lack of star-like actin cytoskeleton arrays. We conclude that ARP2/3 dependent nucleation of actin filaments is critical for the growth of filamentous cells, which in turn influences moss colonization. In complementation assays, the overexpression of Physcomitrella and Arab idopsis ARP3 genes in the moss arp3 mutant results in full recovery of wild type phenotype. In contrast the ARP3 subunit of fission yeast is not able to complement the moss arp3 mutant of moss indicating that regulation of the ARP2/3 dependent actin nucleation diverged in different kingdoms. RESUME Le réseau d'actine est composé de filaments de F-actine et d'un ensemble de protéines s'y attachant (Actin binding proteins). Le réseau d'actine est nécessaire à la croissance et à la morphogenèse de toutes les cellules eucaryotes. Chez les plantes, le marquage ainsi que l'étude de l'organisation du réseau d'actine ont été réalisés en utilisant une fusion GFP-Talin (Kost et al., 1998) exprimée sous le control d'un promoteur constitutif. Afin d'étudier les structures F-actine dans les cellules de Physcomitrella Patens et pour combler le manque d'information sur le développement des gamétophores, des lignées transgéniques stables surexprimant GFP-Talin ont été crées. Nous avons démontré que l'utilisation du promoteur 35S est inadéquate pour le marquage complet et homogène des filaments d'actine dans toutes les cellules de P. patens. Par contre, l'utilisation du promoteur inductible Gmhsp17.3 nous a permis de réaliser un marquage transitoire et général dans tous les tissus de la mousse. Une heure de choc thermique à 37°C suivis d'un temps de récupération de 12-18h à 25°C sont les conditions optimales (sans dommages cellulaires) pour l'observation des structures F-actine à différentes étapes de développement de la mousse. En utilisant la microscopie confocale, nous avons observé l'existence de structures F-actine accumulées en forme d'étoiles. Ces structures, qui sont liées au réseau de microfilaments d'actine, ont été observées dans les protoplastes en régénération, les cellules des protonema apicales ainsi que dans les rhizoïdes. En suivant la croissance tridimensionnelle, ces structures en étoiles ont été observées dans les cellules meristématiques des bourgeons et des jeunes gamétophores. Par contre, dans les cellules différentiées ces structures laissent place à des réseaux de câbles épais. Nous avons également remarqué que durant la redifferentiation des cellules foliaires le réseau de câbles de F-actine est remplacé par un réseau de F-actine diffus. Dans les cellules eucaryotes, la nucléation des filaments d'actirie précédant leur polymérisation est contrôlé par sept sous unités du complexe ARP2/3 et par des formines. Nous avons isolé le gène codant pour la sous unité ARP3 de P. patens et nous avons crée des mutants arp3 par intégration ciblée (Knockout). L'élongation des cellules chloronema est clairement affectée dans les mutants arp3. La différentiation des caulonemata et des rhizoïdes est bloquée et les gametophores sont légèrement plus courts comparé au type sauvage. A fin d'étudier l'organisation des filaments d'actines dans les mutants arp3, nous avons aussi réalisé un arp3-knockout dans la lignée Hsp-GFP-Talin. La nouvelle lignée générée nous a permis de visualiser une désorganisation du réseau d'actine et une absence complète de structures de F-actine accumulée en forme d'étoiles. Les résultats obtenus nous amènent à conclure que la nucléation (ARP2/3 dépendante) des filaments d'actine est indispensable à la croissance des cellules filamenteuses. Par conséquent, les filaments d'actine semblent avoir un rôle dans la colonisation des milieux par les mousses. Nous avons également procédé à des essais de complémentation du mutant arp3. La surexpression des gènes ARP3 de Physcomitrella et d'Arabidopsis dans les cellules du mutant arp3 rétabli complètement le phénotype WT. Par contre, le gène ARP3 des levures n'est pas suffisant pour complémenter la même mutation dans les cellules de mousses. Ce résultat démontre que les mécanismes de régulation de la nucléation des filaments d'actine (ARP2/3 dépendante) sont différents entre les différents groupes d'eucaryotes.

Algies après cure de hernie inguinale: que faire [Pain after inguinal hernia repair: what to do?]

Hernia repair one of the most frequently performed operations in general surgery. With the introduction of tension-free mesh repair, recurrence rates dropped well below 5% for open and laparoscopic procedures. However, chronic postoperative pain remains a widely neglected complication with a high socio-economic impact. It occurs in about 10-20% of patients after hernia repair. We review the different types of post-herniorrhaphy pain with the typical diagnostic features and we conclude with a pragmatic algorithm based on our clinical experience.

Constraining 3-D electrical resistance tomography with GPR reflection data forimproved aquifer characterization

Surface-based ground penetrating radar (GPR) and electrical resistance tomography (ERT) are common tools for aquifer characterization, because both methods provide data that are sensitive to hydrogeologically relevant quantities. To retrieve bulk subsurface properties at high resolution, we suggest incorporating structural information derived from GPR reflection data when inverting surface ERT data. This reduces resolution limitations, which might hinder quantitative interpretations. Surface-based GPR reflection and ERT data have been recorded on an exposed gravel bar within a restored section of a previously channelized river in northeastern Switzerland to characterize an underlying gravel aquifer. The GPR reflection data acquired over an area of 240×40 m map the aquifer's thickness and two internal sub-horizontal regions with different depositional patterns. The interface between these two regions and the boundary of the aquifer with then underlying clay are incorporated in an unstructured ERT mesh. Subsequent inversions are performed without applying smoothness constraints across these boundaries. Inversion models obtained by using these structural constraints contain subtle resistivity variations within the aquifer that are hardly visible in standard inversion models as a result of strong vertical smearing in the latter. In the upper aquifer region, with high GPR coherency and horizontal layering, the resistivity is moderately high (N300 Ωm). We suggest that this region consists of sediments that were rearranged during more than a century of channelized flow. In the lower low coherency region, the GPR image reveals fluvial features (e.g., foresets) and generally more heterogeneous deposits. In this region, the resistivity is lower (~200 Ωm), which we attribute to increased amounts of fines in some of the well-sorted fluvial deposits. We also find elongated conductive anomalies that correspond to the location of river embankments that were removed in 2002.

A numeric model to simulate solar individual ultraviolet exposure.

Exposure to solar ultraviolet (UV) light is the main causative factor for skin cancer. UV exposure depends on environmental and individual factors. Individual exposure data remain scarce and development of alternative assessment methods is greatly needed. We developed a model simulating human exposure to solar UV. The model predicts the dose and distribution of UV exposure received on the basis of ground irradiation and morphological data. Standard 3D computer graphics techniques were adapted to develop a rendering engine that estimates the solar exposure of a virtual manikin depicted as a triangle mesh surface. The amount of solar energy received by each triangle was calculated, taking into account reflected, direct and diffuse radiation, and shading from other body parts. Dosimetric measurements (n = 54) were conducted in field conditions using a foam manikin as surrogate for an exposed individual. Dosimetric results were compared to the model predictions. The model predicted exposure to solar UV adequately. The symmetric mean absolute percentage error was 13%. Half of the predictions were within 17% range of the measurements. This model provides a tool to assess outdoor occupational and recreational UV exposures, without necessitating time-consuming individual dosimetry, with numerous potential uses in skin cancer prevention and research.

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.

Inferring ultraviolet anatomical exposure patterns while distinguishing the relative contribution of radiation components

Exposure to solar ultraviolet (UV) radiation is the main causative factor for skin cancer. UV exposure depends on environmental and individual factors, but individual exposure data remain scarce. UV irradiance is monitored via different techniques including ground measurements and satellite observations. However it is difficult to translate such observations into human UV exposure or dose because of confounding factors (shape of the exposed surface, shading, behavior, etc.) A collaboration between public health institutions, a meteorological office and an institute specialized in computing techniques developed a model predicting the dose and distribution of UV exposure on the basis of ground irradiation and morphological data. Standard 3D computer graphics techniques were adapted to develop this tool, which estimates solar exposure of a virtual manikin depicted as a triangle mesh surface. The amount of solar energy received by various body locations is computed for direct, diffuse and reflected radiation separately. The radiation components are deduced from corresponding measurements of UV irradiance, and the related UV dose received by each triangle of the virtual manikin is computed accounting for shading by other body parts and eventual protection measures. The model was verified with dosimetric measurements (n=54) in field conditions using a foam manikin as surrogate for an exposed individual. Dosimetric results were compared to the model predictions. The model predicted exposure to solar UV adequately. The symmetric mean absolute percentage error was 13%. Half of the predictions were within 17% range of the measurements. This model allows assessing outdoor occupational and recreational UV exposures, without necessitating time-consuming individual dosimetry, with numerous potential uses in skin cancer prevention and research. Using this tool, we investigated solar UV exposure patterns with respect to the relative contribution of the direct, diffuse and reflected radiation. We assessed exposure doses for various body parts and exposure scenarios of a standing individual (static and dynamic postures). As input, the model used erythemally-weighted ground irradiance data measured in 2009 at Payerne, Switzerland. A year-round daily exposure (8 am to 5 pm) without protection was assumed. For most anatomical sites, mean daily doses were high (typically 6.2-14.6 SED) and exceeded recommended exposure values. Direct exposure was important during specific periods (e.g. midday during summer), but contributed moderately to the annual dose, ranging from 15 to 24% for vertical and horizontal body parts, respectively. Diffuse irradiation explained about 80% of the cumulative annual exposure dose. Acute diffuse exposures were also obtained for cloudy summer days. The importance of diffuse UV radiation should not be underestimated when advocating preventive measures. Messages focused on avoiding acute direct exposures may be of limited efficiency to prevent skin cancers associated with chronic exposure (e.g., squamous cell carcinomas).

Uroradiological screening for upper and lower urinary tract anomalies in patients with hypospadias: a systematic literature review.

INTRODUCTION: Hypospadias is associated with anomalies of the urinary tract, but the exact prevalence and significance of these anomalies are still controversial. OBJECTIVES: To assess the percentage of patients with hypospadias and associated urological anomalies, either requiring or not requiring medical or surgical attention. MATERIAL AND METHODS: We searched several databases using the following Mesh terms: hypospadias AND urination, ultrasonography, urinary tract/abnormalities, urinary bladder/radiography, ureteral obstruction, hydronephrosis or vesico-ureteral reflux. Type of uroradiological studies performed, type of urological anomalies, medical or surgical interventions, number of patients available, enrolled and undergoing uroradiological studies and number of patients with abnormal uroradiological exams were recorded. RESULTS: We found 24 studies. Four studies included 100% of available patients. In the other ones, the percentage of patients undergoing uroradiological screening varied from 12 to 82%. Frequency of anomalies varied from 0 to 56%. The most common anomalies were kidney position anomalies, vesico-ureteral reflux and hydronephrosis. CONCLUSIONS: The data published about screening patients with hypospadias for associated anomalies of their urinary tract are of poor quality. The clinical significance of the anomalies found is difficult to evaluate. We found no relationship between the severity of the hypospadias and associated anomalies of the upper or lower urinary tract.

A systematic review of the surgical treatment of large incisional hernia.

PURPOSE: Incisional hernia (IH) is one of the most frequent postoperative complications. Of all patients undergoing IH repair, a vast amount have a hernia which can be defined as a large incisional hernia (LIH). The aim of this study is to identify the preferred technique for LIH repair. METHODS: A systematic review of the literature was performed and studies describing patients with IH with a diameter of 10 cm or a surface of 100 cm2 or more were included. Recurrence hazards per year were calculated for all techniques using a generalized linear model. RESULTS: Fifty-five articles were included, containing 3,945 LIH repairs. Mesh reinforced techniques displayed better recurrence rates and hazards than techniques without mesh reinforcement. Of all the mesh techniques, sublay repair, sandwich technique with sublay mesh and aponeuroplasty with intraperitoneal mesh displayed the best results (recurrence rates of <3.6%, recurrence hazard <0.5% per year). Wound complications were frequent and most often seen after complex LIH repair. CONCLUSIONS: The use of mesh during LIH repair displayed the best recurrence rates and hazards. If possible mesh in sublay position should be used in cases of LIH repair.

A promising new device for the prevention of parastomal hernia.

Parastomal hernia (PSH) is the most frequent long-term stoma complication with serious negative effects on quality of life. Surgical revision is often required and has a substantial morbidity and recurrence rate. The development of PSH requires revisional surgery with a substantial perioperative morbidity and high failure rate in the long-term follow-up. Prophylactic parastomal mesh insertion during stoma creation has the potential to reduce the rate of PSH, but carries the risk of early and late mesh-related complications such as infection, fibrosis, mesh shrinkage, and/or bowel erosion. We developed a new stomaplasty ring (KORING), which is easy to implant, avoids potential mesh-related complications, and has a high potential of long-term prevention of PSH. Here we describe the technique and the first use.

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.

Rockfall susceptibility assessment and remote geological mapping with LiDAR point clouds

Characterizing the geological features and structures in three dimensions over inaccessible rock cliffs is needed to assess natural hazards such as rockfalls and rockslides and also to perform investigations aimed at mapping geological contacts and building stratigraphy and fold models. Indeed, the detailed 3D data, such as LiDAR point clouds, allow to study accurately the hazard processes and the structure of geologic features, in particular in vertical and overhanging rock slopes. Thus, 3D geological models have a great potential of being applied to a wide range of geological investigations both in research and applied geology projects, such as mines, tunnels and reservoirs. Recent development of ground-based remote sensing techniques (LiDAR, photogrammetry and multispectral / hyperspectral images) are revolutionizing the acquisition of morphological and geological information. As a consequence, there is a great potential for improving the modeling of geological bodies as well as failure mechanisms and stability conditions by integrating detailed remote data. During the past ten years several large rockfall events occurred along important transportation corridors where millions of people travel every year (Switzerland: Gotthard motorway and railway; Canada: Sea to sky highway between Vancouver and Whistler). These events show that there is still a lack of knowledge concerning the detection of potential rockfalls, making mountain residential settlements and roads highly risky. It is necessary to understand the main factors that destabilize rocky outcrops even if inventories are lacking and if no clear morphological evidences of rockfall activity are observed. In order to increase the possibilities of forecasting potential future landslides, it is crucial to understand the evolution of rock slope stability. Defining the areas theoretically most prone to rockfalls can be particularly useful to simulate trajectory profiles and to generate hazard maps, which are the basis for land use planning in mountainous regions. The most important questions to address in order to assess rockfall hazard are: Where are the most probable sources for future rockfalls located? What are the frequencies of occurrence of these rockfalls? I characterized the fracturing patterns in the field and with LiDAR point clouds. Afterwards, I developed a model to compute the failure mechanisms on terrestrial point clouds in order to assess the susceptibility to rockfalls at the cliff scale. Similar procedures were already available to evaluate the susceptibility to rockfalls based on aerial digital elevation models. This new model gives the possibility to detect the most susceptible rockfall sources with unprecented detail in the vertical and overhanging areas. The results of the computation of the most probable rockfall source areas in granitic cliffs of Yosemite Valley and Mont-Blanc massif were then compared to the inventoried rockfall events to validate the calculation methods. Yosemite Valley was chosen as a test area because it has a particularly strong rockfall activity (about one rockfall every week) which leads to a high rockfall hazard. The west face of the Dru was also chosen for the relevant rockfall activity and especially because it was affected by some of the largest rockfalls that occurred in the Alps during the last 10 years. Moreover, both areas were suitable because of their huge vertical and overhanging cliffs that are difficult to study with classical methods. Limit equilibrium models have been applied to several case studies to evaluate the effects of different parameters on the stability of rockslope areas. The impact of the degradation of rockbridges on the stability of large compartments in the west face of the Dru was assessed using finite element modeling. In particular I conducted a back-analysis of the large rockfall event of 2005 (265'000 m3) by integrating field observations of joint conditions, characteristics of fracturing pattern and results of geomechanical tests on the intact rock. These analyses improved our understanding of the factors that influence the stability of rock compartments and were used to define the most probable future rockfall volumes at the Dru. Terrestrial laser scanning point clouds were also successfully employed to perform geological mapping in 3D, using the intensity of the backscattered signal. Another technique to obtain vertical geological maps is combining triangulated TLS mesh with 2D geological maps. At El Capitan (Yosemite Valley) we built a georeferenced vertical map of the main plutonio rocks that was used to investigate the reasons for preferential rockwall retreat rate. Additional efforts to characterize the erosion rate were made at Monte Generoso (Ticino, southern Switzerland) where I attempted to improve the estimation of long term erosion by taking into account also the volumes of the unstable rock compartments. Eventually, the following points summarize the main out puts of my research: The new model to compute the failure mechanisms and the rockfall susceptibility with 3D point clouds allows to define accurately the most probable rockfall source areas at the cliff scale. The analysis of the rockbridges at the Dru shows the potential of integrating detailed measurements of the fractures in geomechanical models of rockmass stability. The correction of the LiDAR intensity signal gives the possibility to classify a point cloud according to the rock type and then use this information to model complex geologic structures. The integration of these results, on rockmass fracturing and composition, with existing methods can improve rockfall hazard assessments and enhance the interpretation of the evolution of steep rockslopes. -- La caractérisation de la géologie en 3D pour des parois rocheuses inaccessibles est une étape nécessaire pour évaluer les dangers naturels tels que chutes de blocs et glissements rocheux, mais aussi pour réaliser des modèles stratigraphiques ou de structures plissées. Les modèles géologiques 3D ont un grand potentiel pour être appliqués dans une vaste gamme de travaux géologiques dans le domaine de la recherche, mais aussi dans des projets appliqués comme les mines, les tunnels ou les réservoirs. Les développements récents des outils de télédétection terrestre (LiDAR, photogrammétrie et imagerie multispectrale / hyperspectrale) sont en train de révolutionner l'acquisition d'informations géomorphologiques et géologiques. Par conséquence, il y a un grand potentiel d'amélioration pour la modélisation d'objets géologiques, ainsi que des mécanismes de rupture et des conditions de stabilité, en intégrant des données détaillées acquises à distance. Pour augmenter les possibilités de prévoir les éboulements futurs, il est fondamental de comprendre l'évolution actuelle de la stabilité des parois rocheuses. Définir les zones qui sont théoriquement plus propices aux chutes de blocs peut être très utile pour simuler les trajectoires de propagation des blocs et pour réaliser des cartes de danger, qui constituent la base de l'aménagement du territoire dans les régions de montagne. Les questions plus importantes à résoudre pour estimer le danger de chutes de blocs sont : Où se situent les sources plus probables pour les chutes de blocs et éboulement futurs ? Avec quelle fréquence vont se produire ces événements ? Donc, j'ai caractérisé les réseaux de fractures sur le terrain et avec des nuages de points LiDAR. Ensuite, j'ai développé un modèle pour calculer les mécanismes de rupture directement sur les nuages de points pour pouvoir évaluer la susceptibilité au déclenchement de chutes de blocs à l'échelle de la paroi. Les zones sources de chutes de blocs les plus probables dans les parois granitiques de la vallée de Yosemite et du massif du Mont-Blanc ont été calculées et ensuite comparés aux inventaires des événements pour vérifier les méthodes. Des modèles d'équilibre limite ont été appliqués à plusieurs cas d'études pour évaluer les effets de différents paramètres sur la stabilité des parois. L'impact de la dégradation des ponts rocheux sur la stabilité de grands compartiments de roche dans la paroi ouest du Petit Dru a été évalué en utilisant la modélisation par éléments finis. En particulier j'ai analysé le grand éboulement de 2005 (265'000 m3), qui a emporté l'entier du pilier sud-ouest. Dans le modèle j'ai intégré des observations des conditions des joints, les caractéristiques du réseau de fractures et les résultats de tests géoméchaniques sur la roche intacte. Ces analyses ont amélioré l'estimation des paramètres qui influencent la stabilité des compartiments rocheux et ont servi pour définir des volumes probables pour des éboulements futurs. Les nuages de points obtenus avec le scanner laser terrestre ont été utilisés avec succès aussi pour produire des cartes géologiques en 3D, en utilisant l'intensité du signal réfléchi. Une autre technique pour obtenir des cartes géologiques des zones verticales consiste à combiner un maillage LiDAR avec une carte géologique en 2D. A El Capitan (Yosemite Valley) nous avons pu géoréferencer une carte verticale des principales roches plutoniques que j'ai utilisé ensuite pour étudier les raisons d'une érosion préférentielle de certaines zones de la paroi. D'autres efforts pour quantifier le taux d'érosion ont été effectués au Monte Generoso (Ticino, Suisse) où j'ai essayé d'améliorer l'estimation de l'érosion au long terme en prenant en compte les volumes des compartiments rocheux instables. L'intégration de ces résultats, sur la fracturation et la composition de l'amas rocheux, avec les méthodes existantes permet d'améliorer la prise en compte de l'aléa chute de pierres et éboulements et augmente les possibilités d'interprétation de l'évolution des parois rocheuses.