Hybrid modelling of coupled pore fluid-solid deformation problems

A hybrid formulation for coupled pore fluid-solid deformation problems is proposed. The scheme is a hybrid in the sense that we use a vertex centered finite volume formulation for the analysis of the pore fluid and a particle method for the solid in our model. The pore fluid formally occupies the same space as the solid particles. The size of the particles is not necessarily equal to the physical size of materials. A finite volume mesh for the pore fluid flow is generated by Delaunay triangulation. Each triangle possesses an initial porosity. Changes of the porosity are specified by the translations of the mass centers of particles. Net pore pressure gradients are applied to the particle centers and are considered in the particle momentum balance. The potential of our model is illustrated by means of a simulation of coupled fracture and fluid flow developed in porous rock under biaxial compression condition.

Evolution of stress deficit and changing rates of seismicity in cellular automaton models of earthquake faults

We investigate the internal dynamics of two cellular automaton models with heterogeneous strength fields and differing nearest neighbour laws. One model is a crack-like automaton, transferring ail stress from a rupture zone to the surroundings. The other automaton is a partial stress drop automaton, transferring only a fraction of the stress within a rupture zone to the surroundings. To study evolution of stress, the mean spectral density. f(k(r)) of a stress deficit held is: examined prior to, and immediately following ruptures in both models. Both models display a power-law relationship between f(k(r)) and spatial wavenumber (k(r)) of the form f(k(r)) similar tok(r)(-beta). In the crack model, the evolution of stress deficit is consistent with cyclic approach to, and retreat from a critical state in which large events occur. The approach to criticality is driven by tectonic loading. Short-range stress transfer in the model does not affect the approach to criticality of broad regions in the model. The evolution of stress deficit in the partial stress drop model is consistent with small fluctuations about a mean state of high stress, behaviour indicative of a self-organised critical system. Despite statistics similar to natural earthquakes these simplified models lack a physical basis. physically motivated models of earthquakes also display dynamical complexity similar to that of a critical point system. Studies of dynamical complexity in physical models of earthquakes may lead to advancement towards a physical theory for earthquakes.

Accelerating seismic energy release and evolution of event time and size statistics: Results from two heterogeneous cellular automaton models

The evolution of event time and size statistics in two heterogeneous cellular automaton models of earthquake behavior are studied and compared to the evolution of these quantities during observed periods of accelerating seismic energy release Drier to large earthquakes. The two automata have different nearest neighbor laws, one of which produces self-organized critical (SOC) behavior (PSD model) and the other which produces quasi-periodic large events (crack model). In the PSD model periods of accelerating energy release before large events are rare. In the crack model, many large events are preceded by periods of accelerating energy release. When compared to randomized event catalogs, accelerating energy release before large events occurs more often than random in the crack model but less often than random in the PSD model; it is easier to tell the crack and PSD model results apart from each other than to tell either model apart from a random catalog. The evolution of event sizes during the accelerating energy release sequences in all models is compared to that of observed sequences. The accelerating energy release sequences in the crack model consist of an increase in the rate of events of all sizes, consistent with observations from a small number of natural cases, however inconsistent with a larger number of cases in which there is an increase in the rate of only moderate-sized events. On average, no increase in the rate of events of any size is seen before large events in the PSD model.

Energy and averages in the mechanics of granular materials

We derive a general thermo-mechanical theory for particulate materials consisting of granules of arbitrary whose material points possess three translational and three independent rotational degrees of freedom. Additional field variables are the translational and rotational granular temperatures, the kinetic energies shape and size. The kinematics of granulate is described within the framework of a polar continuum theory of the velocity and spin fluctuations respectively and the usual thermodynamic temperature. We distinguish between averages over particle categories (averages in mass/velocity and moment of inertia/spin space, respectively) and particle phases where the average extends over distinct subsets of particle categories (multi phase flows). The relationship between the thermal energy in the granular system and phonon energy in a molecular system is briefly discussed in the main body of the paper and discussed in detail in the Appendix A. (C) 2001 Elsevier Science B.V. All rights reserved.

Large amplitude folding in finely layered viscoelastic rock structures

We analyze folding phenomena in finely layered viscoelastic rock. Fine is meant in the sense that the thickness of each layer is considerably smaller than characteristic structural dimensions. For this purpose we derive constitutive relations and apply a computational simulation scheme (a finite-element based particle advection scheme; see MORESI et al., 2001) suitable for problems involving very large deformations of layered viscous and viscoelastic rocks. An algorithm for the time integration of the governing equations as well as details of the finite-element implementation is also given. We then consider buckling instabilities in a finite, rectangular domain. Embedded within this domain, parallel to the longer dimension we consider a stiff, layered plate. The domain is compressed along the layer axis by prescribing velocities along the sides. First, for the viscous limit we consider the response to a series of harmonic perturbations of the director orientation. The Fourier spectra of the initial folding velocity are compared for different viscosity ratios. Turning to the nonlinear regime we analyze viscoelastic folding histories up to 40% shortening. The effect of layering manifests itself in that appreciable buckling instabilities are obtained at much lower viscosity ratios (1:10) as is required for the buckling of isotropic plates (1:500). The wavelength induced by the initial harmonic perturbation of the director orientation seems to be persistent. In the section of the parameter space considered here elasticity seems to delay or inhibit the occurrence of a second, larger wavelength. Finally, in a linear instability analysis we undertake a brief excursion into the potential role of couple stresses on the folding process. The linear instability analysis also provides insight into the expected modes of deformation at the onset of instability, and the different regimes of behavior one might expect to observe.

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.

Relative motions of Africa, Iberia and Europe during Alpine orogeny

A revised kinematic model for the motions of Africa and Iberia relative to Europe since the Middle Jurassic is presented in order to provide boundary conditions for Alpine-Mediterranean reconstructions. These motions were calculated using up-to-date kinematic data predominantly based on magnetic isochrons in the Atlantic Ocean and published by various authors during the last 15 years. It is shown that convergence of Africa with respect to Europe commenced during the Cretaceous Normal Superchron (CNS), between chrons MO and 34 (120-83 Ma). This motion was subjected to fluctuations in convergence rates characterised by two periods of relatively rapid convergence (during Late Cretaceous and Eocene-Oligocene times) that alternated with periods of slower convergence (during the Paleocene and since the Early Miocene). Distinct changes in plate kinematics are recognised in the motion of Iberia with respect to Europe, indicated by: (1) a Late Jurassic-Early Cretaceous left-lateral strike-slip motion; (2) Late Cretaceous convergence; (3) Paleocene quiescence; (4) a short period of right-lateral strike-slip motion; and (5) final Eocene-Oligocene convergence. Based on these results, it is speculated that a collisional episode in the Alpine orogeny at ca. 65 Ma resulted in a dramatic decrease in the relative plate motions and that a slower motion since the Early Miocene promoted extension in the Mediterranean back-arc basins. (C) 2002 Elsevier Science B.V. All rights reserved.

Deconvolving the delta O-18 seawater component from subseasonal coral delta O-18 and Sr/Ca at Rarotonga in the southwestern subtropical Pacific for the period 1726 to 1997

To reconstruct oceanographic variations in the subtropical South Pacific, 271-year long subseasonal time series of Sr/Ca and delta(18)O were generated from a coral growing at Rarotonga (21.5degreesS, 159.5degreesW). In this case, coral Sr/Ca appears to be an excellent proxy for sea surface temperature (SST) and coral delta(18)O is a function of both SST and seawater delta(18)O composition (delta(18)O(sw)). Here, we focus on extracting the delta(18)O(sw) signal from these proxy records. A method is presented assuming that coral Sr/Ca is solely a function of SST and that coral delta(18)O is a function of both SST and delta(18)O(sw). This method separates the effects of delta(18)O(sw) from SST by breaking the instantaneous changes of coral delta(18)O into separate contributions by instantaneous SST and delta(18)O(sw) changes, respectively. The results show that on average delta(18)O(sw) at Rarotonga explains similar to39% of the variance in delta(18)O and that variations in SST explains the remaining similar to61% of delta(18)O variance. Reconstructed delta(18)O(sw) shows systematic increases in summer months (December-February) consistent with the regional pattern of variations in precipitation and evaporation. The delta(18)O(sw) also shows a positive linear correlation with satellite-derived estimated salinity for the period 1980 to 1997 (r = 0.72). This linear correlation between reconstructed delta(18)O(sw) and salinity makes it possible to use the reconstructed delta(18)O(sw) to estimate the past interannual and decadal salinity changes in this region. Comparisons of coral delta(18)O and delta(18)O(sw) at Rarotonga with the Pacific decadal oscillation index suggest that the decadal and interdecadal salinity and SST variability at Rarotonga appears to be related to basin-scale decadal variability in the Pacific. Copyright (C) 2002 Elsevier Science Ltd.

Longevity and stability of cratonic lithosphere: Insights from numerical simulations of coupled mantle convection and continental tectonics

[1] The physical conditions required to provide for the tectonic stability of cratonic crust and for the relative longevity of deep cratonic lithosphere within a dynamic, convecting mantle are explored through a suite of numerical simulations. The simulations allow chemically distinct continents to reside within the upper thermal boundary layer of a thermally convecting mantle layer. A rheologic formulation, which models both brittle and ductile behavior, is incorporated to allow for plate-like behavior and the associated subduction of oceanic lithosphere. Several mechanisms that may stabilize cratons are considered. The two most often invoked mechanisms, chemical buoyancy and/or high viscosity of cratonic root material, are found to be relatively ineffective if cratons come into contact with subduction zones. High root viscosity can provide for stability and longevity but only within a thick root limit in which the thickness of chemically distinct, high-viscosity cratonic lithosphere exceeds the thickness of old oceanic lithosphere by at least a factor of 2. This end-member implies a very thick mechanical lithosphere for cratons. A high brittle yield stress for cratonic lithosphere as a whole, relative to oceanic lithosphere, is found to be an effective and robust means for providing stability and lithospheric longevity. This mode does not require exceedingly deep strength within cratons. A high yield stress for only the crustal or mantle component of the cratonic lithosphere is found to be less effective as detachment zones can then form at the crust-mantle interface which decreases the longevity potential of cratonic roots. The degree of yield stress variations between cratonic and oceanic lithosphere required for stability and longevity can be decreased if cratons are bordered by continental lithosphere that has a relatively low yield stress, i.e., mobile belts. Simulations that combine all the mechanisms can lead to crustal stability and deep root longevity for model cratons over several mantle overturn times, but the dominant stabilizing factor remains a relatively high brittle yield stress for cratonic lithosphere.

Convective instability of 3-D fluid-saturated geological fault zones heated from below

We conduct a theoretical analysis to investigate the convective instability of 3-D fluid-saturated geological fault zones when they are heated uniformly from below. In particular, we have derived exact analytical solutions for the critical Rayleigh numbers of different convective flow structures. Using these critical Rayleigh numbers, three interesting convective flow structures have been identified in a geological fault zone system. It has been recognized that the critical Rayleigh numbers of the system have a minimum value only for the fault zone of infinite length, in which the corresponding convective flow structure is a 2-D slender-circle flow. However, if the length of the fault zone is finite, the convective flow in the system must be 3-D. Even if the length of the fault zone is infinite, since the minimum critical Rayleigh number for the 2-D slender-circle flow structure is so close to that for the 3-D convective flow structure, the system may have almost the same chance to pick up the 3-D convective flow structures. Also, because the convection modes are so close for the 3-D convective flow structures, the convective flow may evolve into the 3-D finger-like structures, especially for the case of the fault thickness to height ratio approaching zero. This understanding demonstrates the beautiful aspects of the present analytical solution for the convective instability of 3-D geological fault zones, because the present analytical solution is valid for any value of the ratio of the fault height to thickness. Using the present analytical solution, the conditions, under which different convective flow structures may take place, can be easily determined.

Neogene and Quaternary rollback evolution of the Tyrrhenian Sea, the Apennines, and the Sicilian Maghrebides

Reconstruction of the evolution of the Tyrrhenian Sea shows that the major stage of rifting associated with the opening of this basin began at similar to10 Ma. It involved two episodes of back arc extension, which were induced by the rollback of a west dipping subducting slab. The first period of extension (10-6 Ma) was prominent in the northern Tyrrhenian Sea and in the western part of the southern Tyrrhenian Sea. The second period of extension, mainly affected the southern Tyrrhenian Sea, began in the latest Messinian (6-5 Ma) and has been accompanied by subduction rollback at rates of 60-100 km Myr(-1). Slab reconstruction, combined with paleomagnetic and paleogeographic constraints, indicates that in the central Apennines, the latest Messinian (6-5 Ma) arrival of a carbonate platform at the subduction zone impeded subduction and initiated a slab tear and major strike-slip faults. These processes resulted in the formation of a narrow subducting slab beneath the Ionian Sea that has undergone faster subduction rollback and induced extreme rates of back arc extension.

Soil morphological control on saline and freshwater lake hydrogeochemistry in the Pantanal of Nhecolandia, Brazil

Joint pedological, geochemical, hydrological and geophysical investigations were performed to study the coexistence or saline and freshwater lakes in close proximity and similar climatic conditions in the Nhecolandia region, Pantanal wetlands in Brazil. The saline lakes are concentrically surrounded by green sandy loam horizons, which cause differential hydrological regimes. Mg-calcite, K-silicates, and amorphous silica precipitate in the soil cover, whereas Mg-silicates and more soluble Na-carbonates are concentrated in the topsoil along the shore of the saline lake. In saline solutions, some minor elements (As, Se) reach values above the water quality recommendations, whereas others are controlled and incorporated in solid phases (Ba, Sr). Locally, the destruction of the sandy loam horizons generates very acidic soil solution (pH similar to 3.5) through a process not yet understood. The soil distributions indicate that some freshwater lakes are former saline lakes. They are invaded by freshwater after destruction of the sandy loam green horizons, then the freshwater becomes enriched in K(+), SO(4)(2-), Fe, Al, and a stream of minor and trace elements. The formation of these green sandy loam horizons in the saline environment and their destruction in the non-saline one emphasizes the dynamic nature of this environment (C) 2008 Elsevier B.V. All rights reserved.

Subduction of the Nazca Ridge and the Inca Plateau: Insights into the formation of ore deposits in Peru

A large number of ore deposits that formed in the Peruvian Andes during the Miocene (15-5 Ma) are related to the subduction of the Nazea plate beneath the South American plate. Here we show that the spatial and temporal distribution of these deposits correspond with the arrival of relatively buoyant topographic anomalies, namely the Nazca Ridge in central Peru and the now-consumed Inca Plateau in northern Peru, at the subduction zone. Plate reconstruction shows a rapid metallogenic response to the arrival of the topographic anomalies at the subduction trench. This is indicated by clusters of ore deposits situated within the proximity of the laterally migrating zones of ridge subduction. It is accordingly suggested that tectonic changes associated with impingement of the aseismic ridge into the subduction zone may trigger the formation of ore deposits in metallogenically fertile suprasubduction environments. (c) 2005 Elsevier B.V All rights reserved.

Re-dating mid-Holocene betelnut (Areca catechu L.) and other plant use at Dongan, Papua New Guinea

Direct accelerator mass spectrometry (AMS) dating of anaerobically preserved plant remains from the Dongan site in New Guinea, combined with assessment of preservation condition, confirms earlier doubts about the antiquity of betel-nut (Areca catechu L.) found at the site. A possible sago leaf fragment is also identified as a modem contaminant. The mid-Holocene age of other fruit and nut remains is verified using these methods. The utility of AMS dating in combination with detailed archaeobotanical assessment is demonstrated, thus improving chronometric hygiene and with it knowledge of past plant use in Oceania.