Direct simulation of interface dynamics: linking capillary pressure, interfacial area and surface energy

We perform direct numerical simulations of drainage by solving Navier- Stokes equations in the pore space and employing the Volume Of Fluid (VOF) method to track the evolution of the fluid-fluid interface. After demonstrating that the method is able to deal with large viscosity contrasts and to model the transition from stable flow to viscous fingering, we focus on the definition of macroscopic capillary pressure. When the fluids are at rest, the difference between inlet and outlet pressures and the difference between the intrinsic phase average pressure coincide with the capillary pressure. However, when the fluids are in motion these quantities are dominated by viscous forces. In this case, only a definition based on the variation of the interfacial energy provides an accurate measure of the macroscopic capillary pressure and allows separating the viscous from the capillary pressure components.

Inheritance of ductile and brittle structures in the development of large rock slope instabilities: Examples from western Norway

The high density of slope failures in western Norway is due to the steep relief and to the concentration of various structures that followed protracted ductile and brittle tectonics. On the 72 investigated rock slope instabilities, 13 were developed in soft weathered mafic and phyllitic allochthons. Only the intrinsic weakness of such rocks increases the susceptibility to gravitational deformation. In contrast, the gravitational structures in the hard gneisses reactivate prominent ductile or/and brittle fabrics. At 30 rockslides along cataclinal slopes, weak mafic layers of foliation are reactivated as basal planes. Slope-parallel steep foliation forms back-cracks of unstable columns. Folds are specifically present in the Storfjord area, together with a clustering of potential slope failures. Folding increases the probability of having favourably orientated planes with respect to the gravitational forces and the slope. High water pressure is believed to seasonally build up along the shallow-dipping Caledonian detachments and may contribute to destabilization of the rock slope upwards. Regional cataclastic faults localized the gravitational structures at 45 sites. The volume of the slope instabilities tends to increase with the amount of reactivated prominent structures and the spacing of the latter controls the size of instabilities.

Structural analysis of Turtle mountain (Alberta) using digital elevation model: Toward a progressive failure

In 1903, the eastern slope of Turtle Mountain (Alberta) was affected by a 30 M m3-rockslide named Frank Slide that resulted in more than 70 casualties. Assuming that the main discontinuity sets, including bedding, control part of the slope morphology, the structural features of Turtle Mountain were investigated using a digital elevation model (DEM). Using new landscape analysis techniques, we have identified three main joint and fault sets. These results are in agreement with those sets identified through field observations. Landscape analysis techniques, using a DEM, confirm and refine the most recent geology model of the Frank Slide. The rockslide was initiated along bedding and a fault at the base of the slope and propagated up slope by a regressive process following a surface composed of pre-existing discontinuities. The DEM analysis also permits the identification of important geological structures along the 1903 slide scar. Based on the so called Sloping Local Base Level (SLBL) an estimation was made of the present unstable volumes in the main scar delimited by the cracks, and around the south area of the scar (South Peak). The SLBL is a method permitting a geometric interpretation of the failure surface based on a DEM. Finally we propose a failure mechanism permitting the progressive failure of the rock mass that considers gentle dipping wedges (30°). The prisms or wedges defined by two discontinuity sets permit the creation of a failure surface by progressive failure. Such structures are more commonly observed in recent rockslides. This method is efficient and is recommended as a preliminary analysis prior to field investigation.

Recrystallisation rims in zircon (Valle d'Arbedo, Switzerland): an integrated cathodoluminescence, LA-ICP-MS, SHRIMP, and TEM study

Recrystallization rims are a common feature of zircon crystals that underwent metamorphism. We present a microstructural and microchemical study of partially recrystallized zircon grains collected in polymetamorphic migmatites (Valle d'Arbedo, Ticino, Switzerland). The rims are bright in cathodo-luminescence (CL), with sharp and convex contacts characterized by inward-penetrating embayments transgressing igneous zircon cores. Laser ablation-inductively coupled plasma-mass spectrometry (LA-ICP-MS) data and transmission electron microscopy (TEM) imaging indicate that the rims are chemically and microstructurally different from the cores. The rims are strongly depleted in REE, with concentrations up to two orders of magnitude lower than in the cores, indicating a significant loss of REE during zircon recrystallization. Enrichment in non-formula elements, such as Ca, has not been observed in the rims. The microstructure of zircon cores shows a dappled intensity at and below the 100 nm scale, possibly due to radiation damage. Other defects such as pores and dislocations are absent in the core except at healed cracks. Zircon rims are mostly dapple-free, but contain nanoscale pores and strain centers, interpreted as fluid inclusions and chemical residues, respectively. Sensitive high-resolution ion microprobe (SHRIMP) U-Pb ages show that the recrystallization of the rims took place >200 Ma ago when the parent igneous zircon was not metamict. The chemical composition and the low-Ti content of the rims indicate that they form at sub-solidus temperatures (550-650 degrees C). Recrystallization rims in Valle d'Arbedo zircon are interpreted as the result of the migration of chemical reaction fronts in which fluid triggered in situ and contemporaneous interface-coupled dissolution-reprecipitation mechanisms. This study indicates that strong lattice strain resulting from the incorporation of a large amount of impurities and structural defects is not a necessary condition for zircon to recrystallize. Our observations suggest that the early formation of recrystallization rims played a major role in preserving zircon from the more recent Alpine metamorphic overprint.

A new technique to measure micromotion distribution around a cementless femoral stem.

The interfacial micromotion is closely associated to the long-term success of cementless hip prostheses. Various techniques have been proposed to measure them, but only a few number of points over the stem surface can be measured simultaneously. In this paper, we propose a new technique based on micro-Computer Tomography (μCT) to measure locally the relative interfacial micromotions between the metallic stem and the surrounding femoral bone. Tantalum beads were stuck at the stem surface and spread at the endosteal surface. Relative micromotions between the stem and the endosteal bone surfaces were measured at different loading amplitudes. The estimated error was 10μm and the maximal micromotion was 60μm, in the loading direction, at 1400N. This pilot study provided a local measurement of the micromotions in the 3 direction and at 8 locations on the stem surface simultaneously. This technique could be easily extended to higher loads and a much larger number of points, covering the entire stem surface and providing a quasi-continuous distribution of the 3D interfacial micromotions around the stem. The new measurement method would be very useful to compare the induced micromotions of different stem designs and to optimize the primary stability of cementless total hip arthroplasty.

Biomechanical characterization and in vitro mechanical injury of elderly human femoral head cartilage: comparison to adult bovine humeral head cartilage.

OBJECTIVES: In vitro mechanical injury of articular cartilage is useful to identify events associated with development of post-traumatic osteoarthritis (OA). To date, many in vitro injury models have used animal cartilage despite the greater clinical relevance of human cartilage. We aimed to characterize a new in vitro injury model using elderly human femoral head cartilage and compare its behavior to that of an existing model with adult bovine humeral head cartilage. DESIGN: Mechanical properties of human and bovine cartilage disks were characterized by elastic modulus and hydraulic permeability in radially confined axial compression, and by Young's modulus, Poisson's ratio, and direction-dependent radial strain in unconfined compression. Biochemical composition was assessed in terms of tissue water, solid, and glycosaminoglycan (GAG) contents. Responses to mechanical injury were assessed by observation of macroscopic superficial tissue cracks and histological measurements of cell viability following single injurious ramp loads at 7 or 70%/s strain rate to 3 or 14 MPa peak stress. RESULTS: Confined compression moduli and Young's moduli were greater in elderly human femoral cartilage vs adult bovine humeral cartilage whereas hydraulic permeability was less. Radial deformations of axially compressed explant disks were more anisotropic (direction-dependent) for the human cartilage. In both cartilage sources, tissue cracking and associated cell death during injurious loading was common for 14 MPa peak stress at both strain rates. CONCLUSION: Despite differences in mechanical properties, acute damage induced by injurious loading was similar in both elderly human femoral cartilage and adult bovine humeral cartilage, supporting the clinical relevance of animal-based cartilage injury models. However, inherent structural differences such as cell density may influence subsequent cell-mediated responses to injurious loading and affect the development of OA.

Seismic attenuation due to wave-induced fluid flow at microscopic and macroscopic scales

Wave-induced fluid flow at microscopic and mesoscopic scales arguably constitutes the major cause of intrinsic seismic attenuation throughout the exploration seismic and sonic frequency ranges. The quantitative analysis of these phenomena is, however, complicated by the fact that the governing physical processes may be dependent. The reason for this is that the presence of microscopic heterogeneities, such as micro-cracks or broken grain contacts, causes the stiffness of the so-called modified dry frame to be complex-valued and frequency-dependent, which in turn may affect the viscoelastic behaviour in response to fluid flow at mesoscopic scales. In this work, we propose a simple but effective procedure to estimate the seismic attenuation and velocity dispersion behaviour associated with wave-induced fluid flow due to both microscopic and mesoscopic heterogeneities and discuss the results obtained for a range of pertinent scenarios.

Post-Plutonic Magmatism: from the Source to the Emplacement of an Upper Crustal Dyke Swarm, Adamello Massif, Italy

Magmas of the arc-tholeiitic and calc-alkaline differentiation suites contribute substantially to the formation of continental crust in subduction zones. Different geochemical-petrological models have been put forward to achieve evolved magmas forming large volumes of tonalitic to granitic plutons, building an important part of the continental crust. Primary magmas produced in the mantle wedge overlying the subducted slab migrate through the mantle and the crust. During the transfer, magma can accumulate in intermediate reservoirs at different levels where crystallization leads to differentiation and the heat transfer from the magma, together with gained heat from solidification, lead to partial melting of the crust. Partial melts can be assimilated and mix with more primitive magma. Moreover, already formed crystal cumulates or crystal mushes can be recycled and reactivated to transfer to higher crustal levels. Magma transport in the crust involves fow through fractures within a brittle elastic rock. The solidified magma filled crack, a dyke, can crosscut previously formed geological structures and thus serves as a relative or absolute time marker. The study area is situated in the Adamello massif. The Adamello massif is a composite of plutons that were emplaced between 42 and 29 million years. A later dyke swarm intruded into the southern part of the Adamello Batholith. A fractionation model covering dyke compositions from picrobasalts to dacites results in the cummulative crystallization of 17% olivine, 2% Cr-rich spinel, 18% clinopyroxene, 41% amphibole, 4% plagioclase and 0.1% magnetite to achieve an andesitic composition out of a hydrous primitive picrobasalt. These rocks show a similar geochemical evolution as experimental data simulating fractional crystallization and associated magma differentiation at lower crustal depth (7-10 kbar). The peraluminous, corundum normative composition is one characteristic of more evolved dacitic magmas, which has been explained in a long lasting debate with two di_erent models. Melting of mafic crust or politic material provides one model, whereas an alternative is fractionation from primary mantle derived melts. Amphibole occurring in basaltic-andesitic and andesitic dyke rocks as fractionating cumulate phase extracted from lower crustal depth (6-7.5 kbar) is driving the magmas to peraluminous, corundum normative compositions, which are represented by tonalites forming most of the Adamello Batholith. Most primitive picrobasaltic dykes have a slightly steepened chondrite normalized rare earth elements (REE) pattern and the increased enrichment of light-REE (LREE) for andesites and dacites can be explained by the fractional crystallization model originating from a picrobasalt, taking the changing fractionating phase assemblage and temperature into account. The injection of hot basaltic magma (~1050°C) in a closely spaced dyke swarm increases the surface of the contact to the mainly tonalitic wallrock. Such a setting induces partial melting of the wall rock and selective assimilation. Partial melting of the tonalite host is further expressed through intrusion breccias from basaltic dykes. Heat conduction models with instantaneous magma injection for such a dyke swarm geometry can explain features of partial melting observed in the field. Geochemical data of minerals and bulk rock further underline the selective or bulk assimilation of the tonalite host rock at upper crustal levels (~2-3 kbar), in particular with regard to light ion lithophile elements (LILE) such as Sr, Ba and Rb. Primitive picrobasalts carry an immiscible felsic assimilant as enclaves that bring along refractory rutile and zircon with textures typically found in oceanic plagiogranites or high pressure/low-temperature metamorphic rocks in general. U-Pb data implies a lower Cretaceous age for zircon not yet described as assimilant in Eocene to Oligocene magmatic rocks of the Central Southern Alps. The distribution of post-plutonic dykes in large batholiths such as the Adamello is one of the key features for understanding the regional stress field during the post-batholith emplacement cooling history. The emplacement of the regional dyke swarm covering the southern part of the Adamello massif was associated with consistent left lateral strike-slip movement along magma dilatation planes, leading to en echelon segmentation of dykes. Through the dilation by magma of pre-existing weaknesses and cracks in an otherwise uniform host rock, the dyke propagation and according orientation in the horizontal plane adjusted continuously perpendicular to least compressive remote stress σ3, resulting in an inferred rotation of the remote principal stress field. Les magmas issus des zones de subduction contribuent substantiellement à la formation de la croûte continentale. Les plutons tonalitiques et granitiques représentent, en effet, une partie importante de la croûte continentale. Des magmas primaires produits dans le 'mantle wedge ', partie du manteau se trouvant au-dessus de la plaque plongeante dans des zones de subduction, migrent à travers le manteau puis la croûte. Pendant ce transfert, le magma peut s'accumuler dans des réservoirs intermédiaires à différentes profondeurs. Le stockage de magma dans ces réservoirs engendre, d'une part, la différentiation des magmas par cristallisation fractionnée et, d'autre part, une fusion partielle la croûte continentale préexistante associée au transfert de la chaleur des magmas vers l'encaissant. Ces liquides magmatiques issus de la croûte peuvent, ensuite, se mélanger avec des magmas primaires. Le transport du magma dans la croûte implique notamment un flux de magma à travers différentes fractures recoupant les roches encaissantes élastiques. Au cours de ce processus de migration, des cumulats de cristaux ou des agrégats de cristaux encore non-solidifiés, peuvent être recyclés et réactivés pour être transportés à des niveaux supérieures de la croûte. Le terrain d'étude est situé dans le massif d'Adamello. Celui-ci est composé de plusieurs plutons mis en place entre 42 et 29 millions d'années. Dans une phase tardive de l'activité magmatique liée à ce batholite, une série de filons de composition variable allant de picrobasalte à des compositions dacitiques s'est mise en place la partie sud du massif. Deux modèles sont proposés dans la littérature, pour expliquer la formation des magmas dacitiques caractérisés par des compositions peralumineux (i.e. à corindon normatif). Le premier modèle propose que ces magmas soient issus de la fusion de matériel mafique et pélitique présent dans la partie inférieur de la croûte, alors que le deuxième modèle suggère une évolution par cristallisation fractionnée à partir de liquides primaires issus du manteau. Un modèle de cristallisation fractionnée a pu être développé pour expliquer l'évolution des filons de l'Adamello. Ce modèle explique la formation des filons dacitiques par la cristallisation fractionnée de 17% olivine, 2% spinelle riche en Cr, 18% clinopyroxène, 41% amphibole, 4% plagioclase et 0.1% magnetite à partir de liquide de compositions picrobasaltiques. Ce modèle prend en considération les contraintes pétrologiques déduites de l'observation des différents filons ainsi que du champ de stabilité des différentes phases en fonction de la température. Ces roches montrent une évolution géochimique similaire aux données expérimentales simulant la cristallisation fractionnée de magmas évoluant à des niveaux inférieurs de la croûte (7-10 kbar). Le modèle montre, en particulier, le rôle prépondérant de l'amphibole, une phase qui contrôle en particulier le caractère peralumineux des magmas différentiés ainsi que leurs compositions en éléments en traces. Des phénomènes de fusion partielle de l'encaissant tonalitique lors de la mise en place de _lons mafiques sont observée sur le terrain. L'injection du magma basaltique chaud (~1050°C) sous forme de filons rapprochés augmente la surface du contact avec l'encaissante tonalitique. Une telle situation produit la fusion partielle des roches encaissantes nécessaire à l'incorporation d'enclaves mafiques observés au sein des tonalites. Pour comprendre les conditions nécessaires pour la fusion partielle des roches encaissantes, des modèles de conduction thermique pour une injection simultanée d'une série de filons ont été développées. Des données géochimiques sur les minéraux et sur les roches totales soulignent qu'au niveau supérieur de la croûte, l'assimilation sélective ou totale de l'encaissante tonalitique modifie la composition du liquide primaire pour les éléments lithophiles tel que le Sr, Ba et Rb. Un autre aspect important concernant la pétrologie des filons de l'Adamello est la présence d'enclaves felsiques dans les filons les plus primitifs. Ces enclaves montrent, en particulier, des textures proches de celles rencontrées dans des plagiogranites océaniques ou dans des roches métamorphiques de haute pression/basse température. Ces enclaves contiennent du zircon et du rutile. La datations de ces zircons à l'aide du géochronomètre U-Pb indique un âge Crétacé inférieur. Cet âge est important, car aucune roche de cet âge n'a été considérée comme un assimilant potentiel pour des roches magmatiques d'âge Eocène à Oligocène dans les Alpes Sud Centrales. La réparation spatiale des filons post-plutoniques dans des grands batholites tel que l'Adamello, est une caractéristique clé pour la compréhension des champs de contraintes lors du refroidissement du batholite. L'orientation des filons va, en particulier, indiqué la contrainte minimal au sein des roches encaissante. La mise en place de la série de filon recoupant la partie Sud du massif de l'Adamello est associée à un décrochement senestre, un décrochement que l'on peut lié aux contraintes tectoniques régionales auxquelles s'ajoutent l'effet de la dilatation produite par la mise en place du batholite lui-même. Ce décrochement senestre produit une segmentation en échelon des filons.

Origin of CO2 and carbonate veins in mantle-derived xenoliths in the Pannonian Basin

The origin and evolution of CO2 inclusions and calcite veins in peridotite xenoliths of the Pannonian Basin, Hungary, were investigated by means of petrographic investigation and stable isotope analyses. The fluid inclusions recovered in paragenetic olivine and clinopyroxene belong to distinct populations: type A (texturally early) inclusions with regular shapes (often with negative crystal forms) forming intragranular trails, type B (texturally late) inclusions defining randomly oriented trails that reach grain boundaries Type B inclusions are often associated with silicate melt (type C) inclusions Stable carbon isotope compositions in inclusion-hosted CO2 were obtained by vacuum crushing followed by conventional dual inlet as well as continuous flow mass spectrometry in order to eliminate possible lab artifacts. Olivines, clino- and orthopyroxenes of the host peridotite have oxygen isotope compositions from 5.3 to 6.0 parts per thousand (relative to V-SMOW). without any relationship with xenolith texture. Some of the xenoliths contained calcite in various forms veins and infillings in silicate globules in veins, secondary carbonate veins filling cracks and metasomatic veins with diffuse margins The former two carbonate types have delta C-13 values around -13 parts per thousand (relative to V-PDB) and low Sr contents (<05 wt %), whereas the third type,veins with high-temperature metasomatic features have a delta C-13 value of -5 0 parts per thousand and high Sr contents up to 34 wt.% In spite of the mantle-like delta C-13 value and the unusually high Sr content typical for mantle-derived carbonate, trace element compositions have proven a crustal origin. This observation supports the conclusions of earlier studies that the carbonate melt droplets found on peridotite xenoliths in the alkaline basalts represent mobilized sedimentary carbonate. The large delta C-13 range and the C-12-enrichment in the carbonates can be attributed to devolanlization of the migrating carbonate or infiltration of surficial fluids containing C-12-rich dissolved carbon Carbon isotope compositions of inclusion-hosted CO2 range from -17 8 to -4.8 parts per thousand (relative to V-PDB) with no relation to the amount of CO2 released by vacuum crushing. Low-delta C-13 values measured by stepwise heating under vacuum suggest that the carbon component is pristine and not related to surficial contamination, and that primary mantle fluids with delta C-13 values around -5 parts per thousand were at least partly preserved in the xenoliths Tectonic reworking and heating by the basaltic magma resulted in partial CO2 release and local C-13-depletion. (C) 2010 Elsevier B V All rights reserved

Probing functional groups at the gas-aerosol interface using heterogeneous titration reactions: a tool for predicting aerosol health effects?

The complex chemical and physical nature of combustion and secondary organic aerosols (SOAs) in general precludes the complete characterization of both bulk and interfacial components. The bulk composition reveals the history of the growth process and therefore the source region, whereas the interface controls--to a large extent--the interaction with gases, biological membranes, and solid supports. We summarize the development of a soft interrogation technique, using heterogeneous chemistry, for the interfacial functional groups of selected probe gases [N(CH(3))(3), NH(2)OH, CF(3)COOH, HCl, O(3), NO(2)] of different reactivity. The technique reveals the identity and density of surface functional groups. Examples include acidic and basic sites, olefinic and polycyclic aromatic hydrocarbon (PAH) sites, and partially and completely oxidized surface sites. We report on the surface composition and oxidation states of laboratory-generated aerosols and of aerosols sampled in several bus depots. In the latter case, the biomarker 8-hydroxy-2'-deoxyguanosine, signaling oxidative stress caused by aerosol exposure, was isolated. The increase in biomarker levels over a working day is correlated with the surface density N(i)(O3) of olefinic and/or PAH sites obtained from O(3) uptakes as well as with the initial uptake coefficient, γ(0), of five probe gases used in the field. This correlation with γ(0) suggests the idea of competing pathways occurring at the interface of the aerosol particles between the generation of reactive oxygen species (ROS) responsible for oxidative stress and cellular antioxidants.