Long-range automaton models of earthquakes: Power-law accelerations, correlation evolution, and mode-switching

We introduce a conceptual model for the in-plane physics of an earthquake fault. The model employs cellular automaton techniques to simulate tectonic loading, earthquake rupture, and strain redistribution. The impact of a hypothetical crustal elastodynamic Green's function is approximated by a long-range strain redistribution law with a r(-p) dependance. We investigate the influence of the effective elastodynamic interaction range upon the dynamical behaviour of the model by conducting experiments with different values of the exponent (p). The results indicate that this model has two distinct, stable modes of behaviour. The first mode produces a characteristic earthquake distribution with moderate to large events preceeded by an interval of time in which the rate of energy release accelerates. A correlation function analysis reveals that accelerating sequences are associated with a systematic, global evolution of strain energy correlations within the system. The second stable mode produces Gutenberg-Richter statistics, with near-linear energy release and no significant global correlation evolution. A model with effectively short-range interactions preferentially displays Gutenberg-Richter behaviour. However, models with long-range interactions appear to switch between the characteristic and GR modes. As the range of elastodynamic interactions is increased, characteristic behaviour begins to dominate GR behaviour. These models demonstrate that evolution of strain energy correlations may occur within systems with a fixed elastodynamic interaction range. Supposing that similar mode-switching dynamical behaviour occurs within earthquake faults then intermediate-term forecasting of large earthquakes may be feasible for some earthquakes but not for others, in alignment with certain empirical seismological observations. Further numerical investigation of dynamical models of this type may lead to advances in earthquake forecasting research and theoretical seismology.

Load-unload response ratio and accelerating moment/energy release critical region scaling and earthquake prediction

The main idea of the Load-Unload Response Ratio (LURR) is that when a system is stable, its response to loading corresponds to its response to unloading, whereas when the system is approaching an unstable state, the response to loading and unloading becomes quite different. High LURR values and observations of Accelerating Moment/Energy Release (AMR/AER) prior to large earthquakes have led different research groups to suggest intermediate-term earthquake prediction is possible and imply that the LURR and AMR/AER observations may have a similar physical origin. To study this possibility, we conducted a retrospective examination of several Australian and Chinese earthquakes with magnitudes ranging from 5.0 to 7.9, including Australia's deadly Newcastle earthquake and the devastating Tangshan earthquake. Both LURR values and best-fit power-law time-to-failure functions were computed using data within a range of distances from the epicenter. Like the best-fit power-law fits in AMR/AER, the LURR value was optimal using data within a certain epicentral distance implying a critical region for LURR. Furthermore, LURR critical region size scales with mainshock magnitude and is similar to the AMR/AER critical region size. These results suggest a common physical origin for both the AMR/AER and LURR observations. Further research may provide clues that yield an understanding of this mechanism and help lead to a solid foundation for intermediate-term earthquake prediction.

The role of conserved CXXXRXG motif in the expression and function of the human norepinephrine transporter

Highly conserved motifs in the monoamine transporters, e.g. the human norepinephrine transporter (hNET) GXXXRXG motif which was the focus of the present study, are likely to be important structural features in determining function. This motif was investigated by mutating the glycines to glutamate (causing loss of function) and alanine, and the arginine to glycine. The effects of hG117A, hR121G and hG123A mutations on function were examined in COS-7 cells and compared to hNET. Substrate K-m values were decreased for hG117A and hG123A, and their K values for inhibition of [3 H]nisoxetine binding were decreased 3-4-fold and 4-6-fold, respectively. Transporter turnover was reduced to 65% of hNET for hG117A and hR121G and to 28% for hG123A, suggesting that substrate translocation is impaired. K values of nisoxetine and desipramine for inhibition of [H-3]norepinephrine uptake were increased by 5-fold for hG117A, with no change for cocaine. The K-i value of cocaine was increased by 3-fold for hG123A, with no change for nisoxetine and desipramine. However, there were no effects of the mutations on the K-d of [H-3]nisoxetine binding or K-i values of desipramine or cocaine for inhibition of [H-3]nisoxetine binding. Hence, glycine residues of the GXXXRXG motif are important determinants of NET expression and function, while the arginine residue does not have a major role. This study also showed that antidepressants and psychostimulants have different NET binding sites and provided the first evidence that different sites on the NET are involved in the binding of inhibitors and their competitive inhibition of substrate uptake. (C) 2002 Elsevier Science B.V. All rights reserved.

Composition and energy density of eggs from two species of freshwater turtle with twofold ranges in egg size

Lipid, protein, ash, carbohydrate and water content and energy density of eggs were measured from different clutches over a range of egg size in two species of freshwater turtle. Dry egg contents consisted of protein (54-60%), lipid (25-31%) and ash (5-6%) while carbohydrate was found to be negligible (

Morphology, ultrastructure, and implied function of ciliated sensory structures on the developmental stages of Merizocotyle icopae (Monogenea : Monocotylidae)

Experimental infections were used to track the fate of the dorsal sensilla of Merizocotyle icopae (Monogenea: Monocotylidae) from nasal tissue of the shovelnose ray, Rhinobatos typus (Rhinobatidae). Scanning and transmission electron microscopy revealed that 3 types of uniciliate dorsal sensilla exist at different times in the development of the monogenean. Type 1 sensilla have little or no invagination where the cilium exits the distal end of the dendrite and possess a ring of epidermis surrounding the cilium distal to the invagination. Type 2 sensilla have a deep invagination where the cilium exits the dendrite. Type 3 sensilla can be distinguished from the other types by the shape of the dendrite. The larvae have predominantly Type I dorsal sensilla, most of which are lost approximately 24 h after infection and a few Type 2 sensilla, which are retained. Additional Type 2 sensilla (termed Adult Type 2 sensilla), which are slightly different morphologically from the Type 2 sensilla of the larvae, form in later stages of development. Numerous Type 3 sensilla are unique to the dorsal surface of adults. Loss of all Type I sensilla upon attachment to the host, R. typus, suggests that these may be chemo- or mechanoreceptors responsible for host location by the swimming infective larvae. Type 2 sensilla appear to be important in the larvae, juveniles, and adults whereas the modality mediated by Type 3 is specific to adults. (C) 2003 Wiley-Liss, Inc.

Identification of residues and domains of Raf important for function in vivo and in vitro

Random mutagenesis and genetic screens for impaired Raf function in Caenorhabditis elegans were used to identify six loss-of-function alleles of lin-45 raf that result in a substitution of a single amino acid. The mutations were classified as weak, intermediate, and strong based on phenotypic severity. We engineered these mutations into the homologous residues of vertebrate Raf-1 and analyzed the mutant proteins for their underlying biochemical defects. Surprisingly, phenotype strength did not correlate with the catalytic activity of the mutant proteins. Amino acid substitutions Val-589 and Ser-619 severely compromised Raf kinase activity, yet these mutants displayed weak phenotypes in the genetic screen. Interestingly, this is because these mutant Raf proteins efficiently activate the MAPK (mitogen-activated protein kinase) cascade in living cells, a result that may inform the analysis of knockout mice. Equally intriguing was the observation that mutant proteins with non-functional Ras-binding domains, and thereby deficient in Ras-mediated membrane recruitment, displayed only intermediate strength phenotypes. This confirms that secondary mechanisms exist to couple Ras to Raf in vivo. The strongest phenotype in the genetic screens was displayed by a S508N mutation that again did not correlate with a significant loss of kinase activity or membrane recruitment by oncogenic Ras in biochemical assays. Ser-508 lies within the Raf-1 activation loop, and mutation of this residue in Raf-1 and the equivalent Ser-615 in B-Raf revealed that this residue regulates Raf binding to MEK. Further characterization revealed that in response to activation by epidermal growth factor, the Raf-S508N mutant protein displayed both reduced catalytic activity and aberrant activation kinetics: characteristics that may explain the C. elegans phenotype.

Impulse-response function of splanchnic circulation with model-independent constraints: theory and experimental validation

Modeling physiological processes using tracer kinetic methods requires knowledge of the time course of the tracer concentration in blood supplying the organ. For liver studies, however, inaccessibility of the portal vein makes direct measurement of the hepatic dual-input function impossible in humans. We want to develop a method to predict the portal venous time-activity curve from measurements of an arterial time-activity curve. An impulse-response function based on a continuous distribution of washout constants is developed and validated for the gut. Experiments with simultaneous blood sampling in aorta and portal vein were made in 13 anesthetized pigs following inhalation of intravascular [O-15] CO or injections of diffusible 3-O[ C-11] methylglucose (MG). The parameters of the impulse-response function have a physiological interpretation in terms of the distribution of washout constants and are mathematically equivalent to the mean transit time ( T) and standard deviation of transit times. The results include estimates of mean transit times from the aorta to the portal vein in pigs: (T) over bar = 0.35 +/- 0.05 min for CO and 1.7 +/- 0.1 min for MG. The prediction of the portal venous time-activity curve benefits from constraining the regression fits by parameters estimated independently. This is strong evidence for the physiological relevance of the impulse-response function, which includes asymptotically, and thereby justifies kinetically, a useful and simple power law. Similarity between our parameter estimates in pigs and parameter estimates in normal humans suggests that the proposed model can be adapted for use in humans.

Identification of intracellular calcium oxalate crystals in Chamelaucium uncinatum (Myrtaceae)

Intracellular inclusions in the pedicel and calyx-tube tissues of Chamelaucium uncinatum Schauer ( Myrtaceae) flowers are irregular in shape. They were shown, by polarised light and scanning electron microscopy, to be birefringent 8.9-29.5 mum druse (i.e. aggregate) crystals. Energy-dispersive X-ray spectroscopy showed that these crystals were predominantly composed of calcium. Histochemical and acid-solubility tests indicated that the crystals were calcium oxalate. Raman microprobe spectroscopy was used to confirm this chemical identity. The calcium oxalate crystals were located in xylem-vessel lumens and also in parenchyma cells adjacent to vascular tissues. Thus, the crystals may function to regulate soluble calcium concentrations in C. uncinatum tissues near sites where calcium is unloaded from the xylem.

Electron exchange coupling for single-donor solid-state spin qubits

Intervalley interference between degenerate conduction band minima has been shown to lead to oscillations in the exchange energy between neighboring phosphorus donor electron states in silicon [B. Koiller, X. Hu, and S. Das Sarma, Phys. Rev. Lett. 88, 027903 (2002); Phys. Rev. B 66, 115201 (2002)]. These same effects lead to an extreme sensitivity of the exchange energy on the relative orientation of the donor atoms, an issue of crucial importance in the construction of silicon-based spin quantum computers. In this article we calculate the donor electron exchange coupling as a function of donor position incorporating the full Bloch structure of the Kohn-Luttinger electron wave functions. It is found that due to the rapidly oscillating nature of the terms they produce, the periodic part of the Bloch functions can be safely ignored in the Heitler-London integrals as was done by Koiller, Hu, and Das Sarma, significantly reducing the complexity of calculations. We address issues of fabrication and calculate the expected exchange coupling between neighboring donors that have been implanted into the silicon substrate using an 15 keV ion beam in the so-called top down fabrication scheme for a Kane solid-state quantum computer. In addition, we calculate the exchange coupling as a function of the voltage bias on control gates used to manipulate the electron wave functions and implement quantum logic operations in the Kane proposal, and find that these gate biases can be used to both increase and decrease the magnitude of the exchange coupling between neighboring donor electrons. The zero-bias results reconfirm those previously obtained by Koiller, Hu, and Das Sarma.

An augmented Lanczos algorithm for the efficient computation of a dot-product of a function of a large sparse symmetric matrix

The Lanczos algorithm is appreciated in many situations due to its speed. and economy of storage. However, the advantage that the Lanczos basis vectors need not be kept is lost when the algorithm is used to compute the action of a matrix function on a vector. Either the basis vectors need to be kept, or the Lanczos process needs to be applied twice. In this study we describe an augmented Lanczos algorithm to compute a dot product relative to a function of a large sparse symmetric matrix, without keeping the basis vectors.

Analysis of the mouse transcriptome for genes involved in the function of the nervous system

We analyzed the mouse Representative Transcript and Protein Set for molecules involved in brain function. We found full-length cDNAs of many known brain genes and discovered new members of known brain gene families, including Family 3 G-protein coupled receptors, voltage-gated channels, and connexins. We also identified previously unknown candidates for secreted neuroactive molecules. The existence of a large number of unique brain ESTs suggests an additional molecular complexity that remains to be explored. A list of genes containing CAG stretches in the coding region represents a first step in the potential identification of candidates for hereditary neurological disorders.