Mantle convection modeling with viscoelastic/brittle lithosphere: Numerical methodology and plate tectonic modeling

The earth's tectonic plates are strong, viscoelastic shells which make up the outermost part of a thermally convecting, predominantly viscous layer. Brittle failure of the lithosphere occurs when stresses are high. In order to build a realistic simulation of the planet's evolution, the complete viscoelastic/brittle convection system needs to be considered. A particle-in-cell finite element method is demonstrated which can simulate very large deformation viscoelasticity with a strain-dependent yield stress. This is applied to a plate-deformation problem. Numerical accuracy is demonstrated relative to analytic benchmarks, and the characteristics of the method are discussed.

Anisotropic convection model for the earth's mantle

The paper presents a theory for modeling flow in anisotropic, viscous rock. This theory has originally been developed for the simulation of large deformation processes including the folding and kinking of multi-layered visco-elastic rock (Muhlhaus et al. [1,2]). The orientation of slip planes in the context of crystallographic slip is determined by the normal vector - the director - of these surfaces. The model is applied to simulate anisotropic mantle convection. We compare the evolution of flow patterns, Nusselt number and director orientations for isotropic and anisotropic rheologies. In the simulations we utilize two different finite element methodologies: The Lagrangian Integration Point Method Moresi et al [8] and an Eulerian formulation, which we implemented into the finite element based pde solver Fastflo (www.cmis.csiro.au/Fastflo/). The reason for utilizing two different finite element codes was firstly to study the influence of an anisotropic power law rheology which currently is not implemented into the Lagrangian Integration point scheme [8] and secondly to study the numerical performance of Eulerian (Fastflo)- and Lagrangian integration schemes [8]. It turned out that whereas in the Lagrangian method the Nusselt number vs time plot reached only a quasi steady state where the Nusselt number oscillates around a steady state value the Eulerian scheme reaches exact steady states and produces a high degree of alignment (director orientation locally orthogonal to velocity vector almost everywhere in the computational domain). In the simulations emergent anisotropy was strongest in terms of modulus contrast in the up and down-welling plumes. Mechanisms for anisotropic material behavior in the mantle dynamics context are discussed by Christensen [3]. The dominant mineral phases in the mantle generally do not exhibit strong elastic anisotropy but they still may be oriented by the convective flow. Thus viscous anisotropy (the main focus of this paper) may or may not correlate with elastic or seismic anisotropy.

Apparent strain localization and shear wave dispersion in elastic fault gouge with microrotations

Shear deformation of fault gouge or other particulate materials often results in observed strain localization, or more precisely, the localization of measured deformation gradients. In conventional elastic materials the strain localization cannot take place therefore this phenomenon is attributed to special types of non-elastic constitutive behaviour. For particulate materials however the Cosserat continuum which takes care of microrotations independent of displacements is a more appropriate model. In elastic Cosserat continuum the localization in displacement gradients is possible under some combinations of the generalized Cosserat elastic moduli. The same combinations of parameters also correspond to a considerable dispersion in shear wave propagation which can be used for independent experimental verification of the proposed mechanism of apparent strain localization in fault gouge.

Surface characterization of 6H-SiC (0001) substrates in indentation and abrasive machining

Surface characterization of 6H-SiC (0001) substrates in indentation and abrasive machining was carried out to investigate microfracture, residual damage, and surface roughness associated with material removal and surface generation. Brittle versus plastic deformation was studied using Vickers indention and nano-indentation. To characterize the abrasive machining response, the 6H-SiC (0001) substrates were ground using diamond wheels with grit sizes of 25, 15 and 7 mum, and then polished with diamond suspensions of 3 and 0.05 mum. It is found that in indentation, there was a scale effect for brittle versus plastic deformation in 6H-SiC substrates. Also, in grinding, the scales of fracture and surface roughness of the substrates decreased with a decrease in diamond grit size. However, in polishing, a reduction in grit size of diamond suspensions gave no significant improvement in surface roughness. Furthermore, the results showed that fracture-free 6H-SiC (0001) surfaces were generated in polishing with the existence of the residual crystal defects, which were associated with the origin of defects in single crystal growth. (C) 2003 Elsevier Ltd. All rights reserved.

The Western Alps from the Jurassic to Oligocene: spatio-temporal constraints and evolutionary reconstructions

Despite extensive research in the last 150 years, the regional tectonic reconstruction of the Western Alps has remained controversial. The curved orogenic belt consists of several ribbon-like continental terranes (Sesia/Austroalpine, Internal Crystalline Massifs, Brianconnais), which are separated by two or more ophiolitic sutures (Piemonte, Valais, Antrona?, Lanzo/ Canavese?). High-pressure (HP) metamorphism of each terrane occurred during distinct orogenic episodes: at similar to65 Ma in the Sesia/Austroalpine, at similar to45 Ma in the Piemonte zone and at similar to35 Ma in the Internal Crystalline Massifs. It is suggested that these events reflect individual accretionary episodes, which together with kinematic indicators and the speed and direction of plate motions, provide constraints for the discussed reconstruction model. The model involves a prolonged orogenic history that took place during relative convergence of Europe and Adria (here considered as a promontory of the African plate). The first accretionary event involved the Sesia/Austroalpine terrane. Final closure of the Piemonte Ocean occurred during the Eocene (similar to45 Ma) and involved ultra-high-pressure (UHP) metamorphism of the Piemonte oceanic crust. Incorporation of the Brianconnais terrane in the accretionary wedge occurred thereafter, possibly during or after subduction of the Valais Ocean in the late Eocene (45-35 Ma). This subduction was terminated at ca. 35 Ma, when the Internal Crystalline Massifs (i.e. the assumed internal parts of the Brianconnais terrane) were buried into great depths and underwent HP and UHP metamorphism. (C) 2004 Elsevier B.V. All rights reserved.

Continental extension: From core complexes to rigid block faulting

Extension of overthickened continental crust is commonly characterized by an early core complex stage of extension followed by a later stage of crustal-scale rigid block faulting. These two stages are clearly recognized during the extensional destruction of the Alpine orogen in northeast Corsica, where rigid block faulting overprinting core complex formation eventually led to crustal separation and the formation of a new oceanic backarc basin (the Ligurian Sea). Here we investigate the geodynamic evolution of continental extension by using a novel, fully coupled thermomechanical numerical model of the continental crust. We consider that the dynamic evolution is governed by fault weakening, which is generated by the evolution of the natural-state variables (i.e., pressure, deviatoric stress, temperature, and strain rate) and their associated energy fluxes. Our results show the appearance of a detachment layer that controls the initial separation of the brittle crust on characteristic listric faults, and a core complex formation that is exhuming strongly deformed rocks of the detachment zone and relatively undeformed crustal cores. This process is followed by a transitional period, characterized by an apparent tectonic quiescence, in which deformation is not localized and energy stored in the upper crust is transferred downward and causes self-organized mobilization of the lower crust. Eventually, the entire crust ruptures on major crosscutting faults, shifting the tectonic regime from core complex formation to wholesale rigid block faulting.

Preliminary study of human breast tissue using synchrotron radiation combining WAXS and SAXS techniques

Using synchrotron radiation, we combined simultaneously wide angle X-ray scattering (WAXS) and small angle X-ray scattering (SAXS) techniques to obtain the scattering profiles of normal and neoplastic breast tissu-es samples at the momentum transfer range 6.28 nm(-1) <= Q(=4 pi.sin(theta/2)lambda) <= 50.26 nm(-1) and 0.15 nm(-1) <= Q <= 1.90 nm(-1), respectively. The results obtained show considerable differences between the scattering profiles of these tissues. We verified that the combination of some parameters (ratio between glandular and adipose peak intensity and third-order axial peak intensity) extracted from scattering profiles can be used for identifying breast cancer. (c) 2009 Elsevier Ltd. All rights reserved.

Mechanical loading modulates glutamate receptor subunit expression in bone

The cellular mechanisms coupling mechanical loading with bone remodeling remain unclear. In the CNS, the excitatory amino acid glutamate (Glu) serves as a potent neurotransmitter exerting its effects via various membrane Glu receptors (GluR). Nerves containing Glu exist close to bone cells expressing functional GluRs. Demonstration of a mechanically sensitive glutamate/aspartate transporter protein and the ability of glutamate to stimulate bone resorption in vitro suggest a role for glutamate linking mechanical load and bone remodeling. We used immunohistochemical techniques to identify the expression of N-methyl-D-aspartate acid (NMDA) and non-NMDA (AMPA or kainate) ionotropic GluR subunits on bone cells in vivo. In bone sections from young adult rats, osteoclasts expressed numerous GluR subunits including AMPA (GluR2/3 and GluR4), kainic acid (GluR567) and NMDA (NMDAR2A, NMDAR2B and NMDAR2C) receptor subtypes. Bone lining cells demonstrated immunoexpression for NMDAR2A, NMDAR2B, NMDAR2C, GluR567, GluR23, GuR2 and GluR4 subunits. Immunoexpression was not evident on osteocytes, chondrocytes or vascular channels. To investigate the effects of mechanical loading on GluR expression, we used a Materials Testing System (MTS) to apply 10 N sinusoidal axial compressive loads percutaneously to the right limbs (radius/ulna, tibia/fibula) of rats. Each limb underwent 300-load cycles/day (cycle rate, 1 Hz) for 4 consecutive days. Contralateral, non-loaded limbs served as controls. Mechanically loaded limbs revealed a load-induced loss of immunoexpression for GluR2/3, GluR4, GluR567 and NMDAR2A on osteoclasts and NMDAR2A, NMDAR2B, GluR2/3 and GluR4 on bone lining cells. Both neonatal rabbit and rat osteoclasts were cultured on bone slices to investigate the effect of the NMDA receptor antagonist, MK801, and the AMPA/kainic acid receptor antagonist, NBQX, on osteoclast resorptive activity in vitro. The inhibition of resorptive function seen suggested that both NMDAR and kainic acid receptor function are required for normal osteoclast function. While the exact role of ionotropic GluRs in skeletal tissue remains unclear, the modulation of GluR subunit expression by mechanical loading lends further support for participation of Glu as a mechanical loading effector. These ionotropic receptors appear to be functionally relevant to normal osteoclast resorptive activity. (C) 2005 Elsevier Inc. All rights reserved.

Effect of sodium dodecylbenzene sulfonate on the motion of three-phase contact lines on the Wilhelmy plate surface

The combined approach of the molecular-kinetic and hydrodynamic theories for description of the motion of three-phase gas-liquid-solid contact lines has been examined using the Wilhelmy plate method. The whole dynamic meniscus has been divided into molecular, hydrodynamic, and static-like regions. The Young-Laplace equation and the molecular-kinetic and hydrodynamic dewetting theories have been applied to describe the meniscus profiles and contact angle. The dissipative forces accompanying the dynamic dewetting have also been investigated. The experiments with a Wilhelmy plate made from an acrylic polymer sheet were carried out using a computerized apparatus for contact angle analysis (OCA 20, DataPhysics, Germany). The extrapolated dynamic contact angle versus velocity of the three-phase contact line for Milli-Q water and 5 x 10(-4) M SDBS solution was experimentally obtained and compared with the combined MHD models with low and moderate Reynolds numbers. The models predict similar results for the extrapolated contact angle. SDBS decreases the equilibrium contact angle and increases the molecular jumping length but does not affect the molecular frequency significantly. The hydrodynamic deformation of the meniscus, viscous dissipation, and friction were also influenced by the SDBS surfactant. (c) 2005 Elsevier Inc. All rights reserved.

Oculometric parameters of hyperopia in children with esotropic amblyopia

Purpose: To study the oculometric parameters of hyperopia in children with esotropic amblyopia, comparing amblyopic eyes with fellow eyes. Methods: Thirty-seven patients (5-8 years old) with bilateral hyperopia and esotropic amblyopia underwent a comprehensive ophthalmic examination, including cycloplegic refraction, keratometry and A-scan ultrasonography. Anterior chamber depth, lens thickness, vitreous chamber depth and total axial length were recorded. The refractive power of the crystalline lens was calculated using Bennett`s equations. Paired Student`s t-tests were used to compare ocular biometric measurements between amblyopic eyes and their fellow eyes. The associations of biometric parameters with refractive errors were assessed using Pearson correlation coefficients and linear regression. Multivariable models including axial length, corneal power and lens power were also constructed. Results: Amblyopic eyes were found to have significantly more hyperopic refraction, less corneal power, greater lens power, shorter vitreous chamber depth and shorter axial length, despite similar anterior chamber depth and lens thickness. The strongest correlation with refractive error was observed for the axial length/corneal radius ratio (r(36) = -0.92, p < 0.001 for amblyopic and r(36) = 0.87, p < 0.001 for fellow eyes). Axial length accounted for 39.2% (R(2)) of the refractive error variance in amblyopic eyes and 35.5% in fellow eyes. Adding corneal power to the model increased R(2) to 85.7% and 79.6%, respectively. A statistically significant correlation was found between axial length and corneal power, indicating decreasing corneal power with increasing axial length, and they were similar for amblyopic eyes (r(36) = 0.53,p < 0.001) and fellow eyes (r(36) = -0.57, p < 0.001). A statistically significant correlation was also found between axial length and lens power, indicating decreasing lens power with increasing axial length (r(36) = -0.72, p < 0.001 for amblyopic eyes and r(36) = -0.69, p < 0.001 for fellow eyes). Conclusions: We observed that the correlation among the major oculometric parameters and their individual contribution to hyperopia in esotropic children were similar in amblyopic and non-amblyopic eyes. This finding suggests that the counterbalancing effect of greater corneal and lens power associated with shorter axial length is similar in both eyes of patients with esotropic amblyopia.

Effects of Si content on defect band formation in hypoeutectic Al-Si die castings

Al-3-11% Si alloys have been high-pressure die-cast and characterized microstructurally. Alstruc was used to calculate the solidification characteristics and fraction of eutectic. Defect bands were observed at all Si contents, although their constitution, position and distinctiveness were a function of Si content. The defect bands contain a higher fraction Al-Si eutectic than the surroundings in all alloys, and porosity was additionally found in the band in AlSi3. With decreasing Si content, the defect bands formed closer to the casting surface, became more prevalent and also the width of the bands decreased. These differences are discussed by considering the effect of Si content on the distribution of solid in the mushy wall layers and on the feeding potentials of the alloys. The observations are consistent with the mechanism proposed by Gourlay et al. in which bands form due to deformation within the solidifying mushy wall layers. (c) 2005 Elsevier B.V. All rights reserved.

A Previously Undescribed Syndrome Combining Fibular Agenesis/Hypoplasia, Oligodactylous Clubfeet, Anonychia/Ungual Hypoplasia, and Other Defects

We describe an apparently new genetic syndrome in six members of a family living in a remote area in Northeastern Brazil. This syndrome comprises: short stature Clue to a marked decrease in the length of the lower limbs (predominantly mesomelic with fibular agenesis/marked hypoplasia), grossly malformed/deformed clubfeet with severe oligodactyly, tipper limbs with acromial dimples and variable motion limitation of the forearms and/or hands, severe nail hypoplasia/anonychia sometimes associated with mild brachydactyly and occasionally with pre-axial polydactyly. This syndrome is apparently distinct from the syndrome of brachydactyly-ectrodactyly with fibular aplasia or hypoplasia (OMIM 113310), the syndrome of fibular aplasia or hypoplasia, femoral bowing and poly-, syn-, and oligodactyly (OMIM 228930), and from other previously described conditions exhibiting fibular agenesis/hvpoplasia. (C) 2008 Wiley-Liss, Inc.

Dynamic strain measurement using piezoelectric polymer film

Piezoelectric polymers have been used to form the basis of dynamic strain gauges for the detection of stress waves. The linearity of response was tested using a split Hopkinson pressure bar arrangement. The results obtained illustrate the effectiveness of piezoelectric film strain gauges in the measurement of axial stress waves.

Determination of the force constant of a single-beam gradient trap by measurement of backscattered light

A single-beam gradient trap could potentially be used to hold a stylus for scanning force microscopy. With a view to development of this technique, we modeled the optical trap as a harmonic oscillator and therefore characterized it by its force constant. We measured force constants and resonant frequencies for 1-4-mu m-diameter polystyrene spheres in a single-beam gradient trap using measurements of back-scattered light. Force constants were determined with both Gaussian and doughnut laser modes, with powers of 3 and 1 mW, respectively. Typical values for spring constants were measured to be between 10(-6) and 4 x 10(-6) N/m. The resonant frequencies of trapped particles were measured to be between 1 and 10 kHz, and the rms amplitudes of oscillations were estimated to be around 40 nm. Our results confirm that the use of the doughnut mode for single-beam trapping is more efficient in the axial direction. (C) 1996 Optical Society of America.

Quantification of orientation of pore patterns in X-ray images of deformed clay

The use of X-ray imaging to assess variability in ceramic fabrication is common in archaeological studies aimed at examining ancient pottery technologies. In this paper, a method based on the measurement of individual pores orientation is presented. This method is successfully applied to ceramic specimens of known origin whose structure signified different deformation histories.