Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.

Fitting genetic models to twin data with binary and ordered categorical responses: A comparison of structural equation modelling and Bayesian hierarchical models

We compare Bayesian methodology utilizing free-ware BUGS (Bayesian Inference Using Gibbs Sampling) with the traditional structural equation modelling approach based on another free-ware package, Mx. Dichotomous and ordinal (three category) twin data were simulated according to different additive genetic and common environment models for phenotypic variation. Practical issues are discussed in using Gibbs sampling as implemented by BUGS to fit subject-specific Bayesian generalized linear models, where the components of variation may be estimated directly. The simulation study (based on 2000 twin pairs) indicated that there is a consistent advantage in using the Bayesian method to detect a correct model under certain specifications of additive genetics and common environmental effects. For binary data, both methods had difficulty in detecting the correct model when the additive genetic effect was low (between 10 and 20%) or of moderate range (between 20 and 40%). Furthermore, neither method could adequately detect a correct model that included a modest common environmental effect (20%) even when the additive genetic effect was large (50%). Power was significantly improved with ordinal data for most scenarios, except for the case of low heritability under a true ACE model. We illustrate and compare both methods using data from 1239 twin pairs over the age of 50 years, who were registered with the Australian National Health and Medical Research Council Twin Registry (ATR) and presented symptoms associated with osteoarthritis occurring in joints of the hand.

Exchange rate sensitivity of Australian international equity funds

This paper analyses the exchange rate exposure displayed by a sample of Australian international equity trusts (IET). Exchange rate exposure is also examined in the context of differing economic climates with particular emphasis on the Asian crisis in mid-1997. It is found that there is evidence of exchange rate exposure particularly in the context of a multiple exchange rate model. Exposure varies substantially between three alternative time periods with different exposure apparent subsequent to the Asian crisis than prior to this event.

The effect of heat transfer mode on heart rate responses and hysteresis during heating and cooling in the estuarine crocodile Crocodylus porosus

The effect of heating and cooling on heart rate in the estuarine crocodile Crocodylus porosus was studied in response to different heat transfer mechanisms and heat loads. Three heating treatments were investigated. C. porosus were: (1) exposed to a radiant heat source under dry conditions; (2) heated via radiant energy while half-submerged in flowing water at 23degreesC and (3) heated via convective transfer by increasing water temperature from 23degreesC to 35degreesC. Cooling was achieved in all treatments by removing the heat source and with C. porosus half-submerged in flowing water at 23degreesC. In all treatments, the heart rate of C. porosus increased markedly in response to heating and decreased rapidly with the removal of the heat source. Heart rate during heating was significantly faster than during cooling at any given body temperature, i.e. there was a significant heart rate hysteresis. There were two identifiable responses to heating and cooling. During the initial stages of applying or removing the heat source, there was a dramatic increase or decrease in heart rate ('rapid response'), respectively, indicating a possible cardiac reflex. This rapid change in heart rate with only a small change or no change in body temperature (

On the robustness of short-term interest rate models

This paper investigates the robustness of a range of short–term interest rate models. We examine the robustness of these models over different data sets, time periods, sampling frequencies, and estimation techniques. We examine a range of popular one–factor models that allow the conditional mean (drift) and conditional variance (diffusion) to be functions of the current short rate. We find that parameter estimates are highly sensitive to all of these factors in the eight countries that we examine. Since parameter estimates are not robust, these models should be used with caution in practice.

The rheology of flotation froths

A literature review has highlighted the need to measure flotation froth rheology in order to fully characterise the role of the froth in the flotation process. The initial investigation using a coaxial cylinder viscometer for froth rheology measurement led to the development of a new device employing a vane measuring head. The modified rheometer was used in industrial scale flotation tests at Mt. Isa Copper Concentrator. The measured froth rheograms show a non-Newtonian nature for the flotation froths (pseudoplastic flow). The evidence of the non-Newtonian flow has questioned the validity of application of the Laplace equation in froth motion modelling as used by a number of researchers, since the assumption of irrotational flow is violated. Correlations between the froth rheology and the froth retention time, water hold-up in the froth and concentrate grades have been found. These correlations are independent of air flow rate (test data at various air flow rates fall on one similar trend line). This implies that froth rheology may be used as a lumped parameter for other operating variables in flotation modelling and scale up. (C) 2003 Elsevier Science B.V. All rights reserved.

CRIM1 Regulates the Rate of Processing and Delivery of Bone Morphogenetic Proteins to the Cell Surface

The Crim1 gene is predicted to encode a transmembrane protein containing six von Willebrand-like cysteine-rich repeats (CRRs) similar to those in the BMP-binding antagonist Chordin (Chrd). In this study, we verify that CRIM1 is a glycosylated, Type I transmembrane protein and demonstrate that the extracellular CRR-containing domain can also be secreted, presumably via processing at the membrane. We have previously demonstrated Crim1 expression at sites consistent with an interaction with bone morphogenetic proteins (BMPs). Here we show that CRIM1 can interact with both BMP4 and BMP7 via the CRR-containing portion of the protein and in so doing acts as an antagonist in three ways. CRIM1 binding of BMP4 and -7 occurs when these proteins are co-expressed within the Golgi compartment of the cell and leads to (i) a reduction in the production and processing of preprotein to mature BMP, (ii) tethering of pre-BMP to the cell surface, and (iii) an effective reduction in the secretion of mature BMP. Functional antagonism was verified by examining the effect of coexpression of CRIM1 and BMP4 on metanephric explant culture. The presence of CRIM1 reduced the effective BMP4 concentration of the media, thereby acting as a BMP4 antagonist. Hence, CRIM1 modulates BMP activity by affecting its processing and delivery to the cell surface

Exceptionally high-rate nitrification in sequencing batch reactors treating high ammonia landfill leachate

The nitrogen removal capacity of a suspended culture system treating mature landfill leachate was investigated. Leachate containing high ammonium levels of 300-900 mg N/L was nitrified in a bench scale sequencing batch reactor. Leachate from four different landfills was treated over a two year period for the removal of nitrogen. In this time, a highly specific nitrifying culture was attained that delivered exceptionally high rates of ammonia removal. No sludge was wasted from the system to increase the throughput and up to 13 g/L of MLSS was obtained. Settleability of the purely nitrifying biomass was excellent with SVI less than 40 mL/g, even at the high sludge concentrations. Nitrification rates up to 246 mg NI(L h) (5.91 g N/(L d)) and specific nitrification rates of 36 mg N/(gVSS h) (880 mg N/(gVSS d)) were obtained. The loading to the system at this time allowed complete nitrification of the leachate with a hydraulic retention time of only 5 hours. Following these successful treatability studies, a full-scale plant was designed and built at one of the landfills investigated.

Fast, scalable master equation solution algorithms. IV. Lanczos iteration with diffusion approximation preconditioned iterative inversion

In this paper we propose a second linearly scalable method for solving large master equations arising in the context of gas-phase reactive systems. The new method is based on the well-known shift-invert Lanczos iteration using the GMRES iteration preconditioned using the diffusion approximation to the master equation to provide the inverse of the master equation matrix. In this way we avoid the cubic scaling of traditional master equation solution methods while maintaining the speed of a partial spectral decomposition. The method is tested using a master equation modeling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long-lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Fast, scalable master equation solution algorithms. III. Direct time propagation accelerated by a diffusion approximation preconditioned iterative solver

In this paper we propose a novel fast and linearly scalable method for solving master equations arising in the context of gas-phase reactive systems, based on an existent stiff ordinary differential equation integrator. The required solution of a linear system involving the Jacobian matrix is achieved using the GMRES iteration preconditioned using the diffusion approximation to the master equation. In this way we avoid the cubic scaling of traditional master equation solution methods and maintain the low temperature robustness of numerical integration. The method is tested using a master equation modelling the formation of propargyl from the reaction of singlet methylene with acetylene, proceeding through long lived isomerizing intermediates. (C) 2003 American Institute of Physics.

Calculation of bound and resonance states of HO2 for nonzero total angular momentum

Bound and resonance states of HO2 have been calculated quantum mechanically by the Lanczos homogeneous filter diagonalization method [Zhang and Smith, Phys. Chem. Chem. Phys. 3, 2282 (2001); J. Chem. Phys. 115, 5751 (2001)] for nonzero total angular momentum J = 1,2,3. For lower bound states, agreement between the results in this paper and previous work is quite satisfactory; while for high lying bound states and resonances these are the first reported results. A helicity quantum number V assignment (within the helicity conserving approximation) is performed and the results indicate that for lower bound states it is possible to assign the V quantum numbers unambiguously, but for resonances it is impossible to assign the V helicity quantum numbers due to strong mixing. In fact, for the high-lying bound states, the mixing has already appeared. These results indicate that the helicity conserving approximation is not good for the resonance state calculations and exact quantum calculations are needed to accurately describe the reaction dynamics for HO2 system. Analysis of the resonance widths shows that most of the resonances are overlapping and the interferences between them lead to large fluctuations from one resonance to another. In accord with the conclusions from earlier J = 0 calculations, this indicates that the dissociation of HO2 is essentially irregular. (C) 2003 American Institute of Physics.