Ductile fabrics in the zone of active oblique convergence near the Alpine Fault, New Zealand: identifying the neotectonic overprint

The mid-crustal Alpine Schist in central Southern Alps, New Zealand has been exhumed during the past similar to3 m.y. on the hanging wall of the oblique-slip Alpine Fault. These rocks underwent ductile deformation during their passage through the similar to 150-km-wide Pacific-Australia plate boundary zone. Likely to be Cretaceous in age, peak metamorphism predates the largely Pliocene and younger oblique convergence that continues to uplift the Southern Alps today. Late Cenozoic ductile deformation constructively reinforced a pre-existing fabric that was well oriented to accommodate a dextral-transpressive overprint. Quartz microstructures below a recently exhumed brittle-ductile transition zone reflect a late Cenozoic increment of ductile strain that was distributed across deeper levels of the Pacific Plate. Deformation was transpressive, including a dextral-normal shear component that bends and rotates a delaminated panel of Pacific Plate crust onto the oblique footwall ramp of the Alpine Fault. Progressive ductile shear in mylonites at the base of the Pacific Plate overprints earlier fabrics in a dextral-reverse sense, a deformation that accompanies translation of the schists up the Alpine Fault. Ductile shear along that structure affects not only the 12-km-thick section of Alpine mylonites, but is distributed across several kilometres of overlying nonmylonitic rocks. (C) 2001 Elsevier Science Ltd. All rights reserved.

Folding processes and solitary waves in structural geology

Evolution of localized folding patterns in layered elastic and visco-elastic materials is reviewed in the context of compressed geological systems. The thin strut or plate embedded in a visco-elastic medium is used as an archetypal example to describe localized buckles which, in contrast to those from earlier formulations, appear in the absence of triggering imperfections. Structural and material effects are surveyed and important nonlinear characteristics are identified. A brief review of possible methods of analysis is conducted.

Analytical solutions of linear finite- and small-strain one-dimensional consolidation

Analytical solutions are presented for linear finite-strain one-dimensional consolidation of initially unconsolidated soil layers with surcharge loading for both one- and two-way drainage. These solutions complement earlier solutions for initially unconsolidated soil layers without surcharge and initially normally consolidated soil layers with surcharge. Small-strain solutions for the consolidation of initially unconsolidated soil layers with surcharge loading are also presented, and the relationship between the earlier solutions for initially unconsolidated soil without surcharge and the corresponding small-strain solutions, which was not addressed in the earlier work, is clarified. The new solutions for initially unconsolidated soil with surcharge loading can be applied to the analysis of low stress consolidation tests and to the partial validation of numerical solutions of non-linear finite-strain consolidation. They also clarify a formerly perplexing aspect of finite-strain solution charts first noted in numerical solutions. Copyright (C) 2004 John Wiley Sons, Ltd.

A sensitive vorticity gauge using rotated porphyroblasts, and its application to rocks adjacent to the Alpine Fault, New Zealand

Variable aspect ratio porphyroblasts deformed in non-coaxial flow. and internally containing rotated relicts of an external foliation, can be used to characterise plane strain flow regimes. The distribution obtained by plotting the orientation of the long axis of such grains, classified by aspect ratio, against the orientation of the internal foliation is potentially a sensitive gauge of both the bulk shear strain (as previously suggested) and kinematic vorticity number. We illustrate the method using rotated biotite porphyroblasts in the Alpine Schist: a sequence of mid-crustal rocks that have been ramped to the surface along the Alpine Fault. a major transpressional plate boundary. Results indicate that, at distances greater than or equal to similar to1 km from the fault, the rocks have undergone a combination of irrotational fattening and dextral-oblique, normal-sense shear, with a bulk shear strain of similar to0.6 and kinematic vorticity number of similar to0.2. The vorticity analysis is compatible with estimates of strongly oblate bulk strain of similar to 75% maximum shortening. Dextral-reverse transpressional flow characterises higher strain S-tectonite mylonite within similar to1 km of the Alpine Fault. These relationships provide insight into the kinematics of flow and distribution of strain in the hangingwall of the Alpine Fault and place constraints on numerical mechanical models for the exhumation of these mid-crustal rocks. (C) 2001 Elsevier Science Ltd. All rights reserved.

Three-dimensional orthotropic viscoelastic finite element model of a human ligament

Ligaments undergo finite strain displaying hyperelastic behaviour as the initially tangled fibrils present straighten out, combined with viscoelastic behaviour (strain rate sensitivity). In the present study the anterior cruciate ligament of the human knee joint is modelled in three dimensions to gain an understanding of the stress distribution over the ligament due to motion imposed on the ends, determined from experimental studies. A three dimensional, finite strain material model of ligaments has recently been proposed by Pioletti in Ref. [2]. It is attractive as it separates out elastic stress from that due to the present strain rate and that due to the past history of deformation. However, it treats the ligament as isotropic and incompressible. While the second assumption is reasonable, the first is clearly untrue. In the present study an alternative model of the elastic behaviour due to Bonet and Burton (Ref. [4]) is generalized. Bonet and Burton consider finite strain with constant modulii for the fibres and for the matrix of a transversely isotropic composite. In the present work, the fibre modulus is first made to increase exponentially from zero with an invariant that provides a measure of the stretch in the fibre direction. At 12% strain in the fibre direction, a new reference state is then adopted, after which the material modulus is made constant, as in Bonet and Burton's model. The strain rate dependence can be added, either using Pioletti's isotropic approximation, or by making the effect depend on the strain rate in the fibre direction only. A solid model of a ligament is constructed, based on experimentally measured sections, and the deformation predicted using explicit integration in time. This approach simplifies the coding of the material model, but has a limitation due to the detrimental effect on stability of integration of the substantial damping implied by the nonlinear dependence of stress on strain rate. At present, an artificially high density is being used to provide stability, while the dynamics are being removed from the solution using artificial viscosity. The result is a quasi-static solution incorporating the effect of strain rate. Alternate approaches to material modelling and integration are discussed, that may result in a better model.

Partitioning of deformation within a subduction channel during exhumation of high-pressure rocks: a case study from the Western Alps

The metamorphic belt of the Western Alps was subjected to widespread extensional tectonism at the end of the Eocene (ca. 45-35 Ma). Extension was accommodated by hinterland-directed movements along gently inclined extensional shear zones, which facilitated rapid exhumation of high-pressure and ultra-high-pressure rocks. This deformation resulted in a normal metamorphic sequence. Extension in the inner parts of the Western Alps was coeval with shortening at the front of the belt (foreland-directed thrusts), which took place during decompression, and emplaced higher grade metamorphic units over lower grade metamorphic rocks, thus forming an inverse metamorphic sequence. Two mechanisms for this extensional episode are discussed: (1) collapse of an overthickened lithosphere, and (2) internal readjustments within the orogenic wedge due to subduction channel dynamics. We favour the latter mechanism because it can account for the development of the observed inverse and normal metamorphic sequences along foreland-directed thrusts and hinterland-directed detachments, respectively. This hypothesis is supported by published structural, metamorphic and geochronological data from four geological transects through the Western Alps. This study also emphasizes the importance of post-shearing deformation (e.g. horizontal buckling versus vertical flattening), which can modify the distribution of hinterland- and foreland-directed shear zones in orogenic belts. (c) 2006 Elsevier Ltd. All rights reserved.

Kinematics of oblique collision and ramping inferred from microstructures and strain in middle crustal rocks, central Southern Alps, New Zealand

The hanging wall of the Alpine Fault near Franz Josef Glacier has been exhumed during the past similar to2-3 m.y. providing a sample of the ductilely deformed middle crust of a modem obliquely convergent orogen. Presently exposed rocks of the Pacific Plate are inferred to have undergone several phases of ductile deformation as they moved westward above a mid-crustal detachment. Initially they were transpressed across the outboard part of the orogen, resulting in oblate fabrics with a down-dip stretch. Later, they encountered the Alpine Fault, experiencing an oblique-slip backshearing on vertical planes. This escalator-like deformation tilted and thinned the incoming crust onto that crustal-scale oblique ramp. This style of hanging wall deformation may affect only the most rapidly uplifting, central part of the Southern Alps because of the low flexural rigidity of the crust in that region and its displacement over a relatively sharp ramp-angle at depth. A 3D transpressive flow affected mylonites locally near the fault, but their shear direction remained parallel to plate motion, ruling out ductile 'extrusion' as an important process in this orogen. Outside the mylonite zone, late Cenozoic shortening is inferred to be modest (30-40%), as measured from deformation of younger biotite grains. Oblique collision is dominated by translation on the Alpine Fault, and rocks migrate rapidly through the deforming zone, preventing the accumulation of large finite strains. Transpression may play a minor role in oblique collision. (C) 2001 Elsevier Science Ltd. All rights reserved.

Comparison of beam theory and finite-element analysis with in vivo bone strain data from the alligator cranium

The mechanical behavior of the vertebrate skull is often modeled using free-body analysis of simple geometric structures and, more recently, finite-element (FE) analysis. In this study, we compare experimentally collected in vivo bone strain orientations and magnitudes from the cranium of the American alligator with those extrapolated from a beam model and extracted from an FE model. The strain magnitudes predicted from beam and FE skull models bear little similarity to relative and absolute strain magnitudes recorded during in vivo biting experiments. However, quantitative differences between principal strain orientations extracted from the FE skull model and recorded during the in vivo experiments were smaller, and both generally matched expectations from the beam model. The differences in strain magnitude between the data sets may be attributable to the level of resolution of the models, the material properties used in the FE model, and the loading conditions (i.e., external forces and constraints). This study indicates that FE models and modeling of skulls as simple engineering structures may give a preliminary idea of how these structures are loaded, but whenever possible, modeling results should be verified with either in vitro or preferably in vivo testing, especially if precise knowledge of strain magnitudes is desired. (c) 2005 Wiley-Liss, Inc.

Strain energies due to nonplanar distortion of fullerenes and their dependence on structural motifs

Accurate strain energies due to nonplanar distortion of 114 isolated pentagon rule (IPR) fullerenes with 60-102 carbon atoms have been calculated based on B3LYP/6-31G(d) optimized structures. The calculated values of strain energy due to nonplanar distortion (E-np) are reproduced by three simple schemes based upon counts of 8, 16, and 30 distinct structural motifs composed of hexagons and pentagons. Using C-180 (I-h) and CN (I-h) (N is very large) as test molecules, the intrinsic limitations of the motif model based on six-membered rings (6-MRs) as the central unit have been discussed. On the basis of the relationship between the contributions of motifs to E-np and the number of five-membered rings (5-MRs) in motifs, we found that IPR fullerenes with dispersed 5-MRs present smaller nonplanar distortions.

Finite deformation model of simple shear of fault with microrotations: apparent strain localisation and en-echelon fracture pattern

Strain localisation is a widespread phenomenon often observed in shear and compressive loading of geomaterials, for example, the fault gouge. It is believed that the main mechanisms of strain localisation are strain softening and mismatch between dilatancy and pressure sensitivity. Observations show that gouge deformation is accompanied by considerable rotations of grains. In our previous work as a model for gouge material, we proposed a continuum description for an assembly of particles of equal radius in which the particle rotation is treated as an independent degree of freedom. We showed that there exist critical values of the model parameters for which the displacement gradient exhibits a pronounced localisation at the mid-surface layers of the fault, even in the absence of inelasticity. Here, we generalise the model to the case of finite deformations characteristic for the gouge deformation. We derive objective constitutive relationships relating the Jaumann rates of stress and moment stress to the relative strain and curvature rates, respectively. The model suggests that the pattern of localisation remains the same as in the linear case. However, the presence of the Jaumann terms leads to the emergence of non-zero normal stresses acting along and perpendicular to the shear layer (with zero hydrostatic pressure), and localised along the mid-line of the gouge; these stress components are absent in the linear model of simple shear. These additional normal stresses, albeit small, cause a change in the direction in which the maximal normal stresses act and in which en-echelon fracturing is formed.