Implications of fault constitutive properties for earthquake prediction.

The rate- and state-dependent constitutive formulation for fault slip characterizes an exceptional variety of materials over a wide range of sliding conditions. This formulation provides a unified representation of diverse sliding phenomena including slip weakening over a characteristic sliding distance Dc, apparent fracture energy at a rupture front, time-dependent healing after rapid slip, and various other transient and slip rate effects. Laboratory observations and theoretical models both indicate that earthquake nucleation is accompanied by long intervals of accelerating slip. Strains from the nucleation process on buried faults generally could not be detected if laboratory values of Dc apply to faults in nature. However, scaling of Dc is presently an open question and the possibility exists that measurable premonitory creep may precede some earthquakes. Earthquake activity is modeled as a sequence of earthquake nucleation events. In this model, earthquake clustering arises from sensitivity of nucleation times to the stress changes induced by prior earthquakes. The model gives the characteristic Omori aftershock decay law and assigns physical interpretation to aftershock parameters. The seismicity formulation predicts large changes of earthquake probabilities result from stress changes. Two mechanisms for foreshocks are proposed that describe observed frequency of occurrence of foreshock-mainshock pairs by time and magnitude. With the first mechanism, foreshocks represent a manifestation of earthquake clustering in which the stress change at the time of the foreshock increases the probability of earthquakes at all magnitudes including the eventual mainshock. With the second model, accelerating fault slip on the mainshock nucleation zone triggers foreshocks.

Slip complexity in earthquake fault models.

We summarize studies of earthquake fault models that give rise to slip complexities like those in natural earthquakes. For models of smooth faults between elastically deformable continua, it is critical that the friction laws involve a characteristic distance for slip weakening or evolution of surface state. That results in a finite nucleation size, or coherent slip patch size, h*. Models of smooth faults, using numerical cell size properly small compared to h*, show periodic response or complex and apparently chaotic histories of large events but have not been found to show small event complexity like the self-similar (power law) Gutenberg-Richter frequency-size statistics. This conclusion is supported in the present paper by fully inertial elastodynamic modeling of earthquake sequences. In contrast, some models of locally heterogeneous faults with quasi-independent fault segments, represented approximately by simulations with cell size larger than h* so that the model becomes "inherently discrete," do show small event complexity of the Gutenberg-Richter type. Models based on classical friction laws without a weakening length scale or for which the numerical procedure imposes an abrupt strength drop at the onset of slip have h* = 0 and hence always fall into the inherently discrete class. We suggest that the small-event complexity that some such models show will not survive regularization of the constitutive description, by inclusion of an appropriate length scale leading to a finite h*, and a corresponding reduction of numerical grid size.

The organization of seismicity on fault networks.

Although models of homogeneous faults develop seismicity that has a Gutenberg-Richter distribution, this is only a transient state that is followed by events that are strongly influenced by the nature of the boundaries. Models with geometrical inhomogeneities of fracture thresholds can limit the sizes of earthquakes but now favor the characteristic earthquake model for large earthquakes. The character of the seismicity is extremely sensitive to distributions of inhomogeneities, suggesting that statistical rules for large earthquakes in one region may not be applicable to large earthquakes in another region. Model simulations on simple networks of faults with inhomogeneities of threshold develop episodes of lacunarity on all members of the network. There is no validity to the popular assumption that the average rate of slip on individual faults is a constant. Intermediate term precursory activity such as local quiescence and increases in intermediate-magnitude activity at long range are simulated well by the assumption that strong weakening of faults by injection of fluids and weakening of asperities on inhomogeneous models of fault networks is the dominant process; the heat flow paradox, the orientation of the stress field, and the low average stress drop in some earthquakes are understood in terms of the asperity model of inhomogeneous faulting.

Geometric incompatibility in a fault system.

Interdependence between geometry of a fault system, its kinematics, and seismicity is investigated. Quantitative measure is introduced for inconsistency between a fixed configuration of faults and the slip rates on each fault. This measure, named geometric incompatibility (G), depicts summarily the instability near the fault junctions: their divergence or convergence ("unlocking" or "locking up") and accumulation of stress and deformations. Accordingly, the changes in G are connected with dynamics of seismicity. Apart from geometric incompatibility, we consider deviation K from well-known Saint Venant condition of kinematic compatibility. This deviation depicts summarily unaccounted stress and strain accumulation in the region and/or internal inconsistencies in a reconstruction of block- and fault system (its geometry and movements). The estimates of G and K provide a useful tool for bringing together the data on different types of movement in a fault system. An analog of Stokes formula is found that allows determination of the total values of G and K in a region from the data on its boundary. The phenomenon of geometric incompatibility implies that nucleation of strong earthquakes is to large extent controlled by processes near fault junctions. The junctions that have been locked up may act as transient asperities, and unlocked junctions may act as transient weakest links. Tentative estimates of K and G are made for each end of the Big Bend of the San Andreas fault system in Southern California. Recent strong earthquakes Landers (1992, M = 7.3) and Northridge (1994, M = 6.7) both reduced K but had opposite impact on G: Landers unlocked the area, whereas Northridge locked it up again.

Development of a bovine X chromosome linkage group and painting probes to assess cattle, sheep, and goat X chromosome segment homologies.

The X chromosome linkage group is conserved in placental mammals. However, X chromosome morphological differences, due to internal chromosome rearrangements, exist among mammalian species. We have developed bovine chromosome painting probes for Xp and Xq to assess segment homologies between the submetacentric bovine X chromosome and the acrocentric sheep and goat X chromosomes. These painting probes and their corresponding DNA libraries were developed by chromosome micromanipulation, DNA micropurification, microcloning, and PCR amplification. The bovine Xp painting probe identified an interstitially located homologous segment in the sheep and goat Xq region, most probably resulting from chromosome inversion. Ten type II (microsatellite) markers obtained from the bovine Xq library and five other X chromosome assigned, but unlinked, markers were used to generate a linkage map for Xq spanning 89.4 centimorgans. The chromosome painting probes and molecular markers generated in this study would be useful for comparative mapping and tracing of internal X chromosome rearrangements in all ruminant species and would contribute to the understanding of mammalian sex chromosome evolution.

A de novo missense mutation of the beta subunit of the epithelial sodium channel causes hypertension and Liddle syndrome, identifying a proline-rich segment critical for regulation of channel activity.

Liddle syndrome is a mendelian form of hypertension characterized by constitutively elevated renal Na reabsorption that can result from activating mutations in the beta or gamma subunit of the epithelial Na channel. All reported mutations have deleted the last 45-76 normal amino acids from the cytoplasmic C terminus of one of these channel subunits. While these findings implicate these terminal segments in the normal negative regulation of channel activity, they do not identify the amino acid residues that are critical targets for these mutations. Potential targets include the short highly conserved Pro-rich segments present in the C terminus of beta and gamma subunits; these segments are similar to SH3-binding domains that mediate protein-protein interaction. We now report a kindred with Liddle syndrome in which affected patients have a mutation in codon 616 of the beta subunit resulting in substitution of a Leu for one of these highly conserved Pro residues. The functional significance of this mutation is demonstrated both by the finding that this is a de novo mutation appearing concordantly with the appearance of Liddle syndrome in the kindred and also by the marked activation of amiloride-sensitive Na channel activity seen in Xenopus oocytes expressing channels containing this mutant subunit (8.8-fold increase compared with control oocytes expressing normal channel subunits; P = 0.003). These findings demonstrate a de novo missense mutation causing Liddle syndrome and identify a critical channel residue important for the normal regulation of Na reabsorption in humans.

Invariant chain made in Escherichia coli has an exposed N-terminal segment that blocks antigen binding to HLA-DR1 and a trimeric C-terminal segment that binds empty HLA-DR1.

Invariant chain (Ii), a membrane glycoprotein, binds class II major histocompatibility complex (MHC) glycoproteins, probably via its class II-associated Ii peptide (CLIP) segment, and escorts them toward antigen-containing endosomal compartments. We find that a soluble, trimeric ectodomain of Ii expressed and purified from Escherichia coli blocks peptide binding to soluble HLA-DR1. Proteolysis indicates that Ii contains two structural domains. The C-terminal two-thirds forms an alpha-helical domain that trimerizes and interacts with empty HLA-DR1 molecules, augmenting rather than blocking peptide binding. The N-terminal one-third, which inhibits peptide binding, is proteolytically susceptible over its entire length. In the trimer, the N-terminal domains act independently with each CLIP segment exposed and free to bind an MHC class II molecule, while the C-terminal domains act as a trimeric unit.

Expression of PTPH1, a rat protein tyrosine phosphatase, is restricted to the derivatives of a specific diencephalic segment.

Studies to date have identified only a few proteins that are expressed in a segment-specific manner within the mammalian brain. Here we report that a nonreceptor protein tyrosine phosphatase, PTPH1, is selectively expressed in the adult thalamus. Expression of PTPH1 mRNA is detected in most, but not all, thalamic nuclei. Nuclei that are derived embryonically from the dorsal thalamus and project to the neocortex express this gene, whereas those derived from the ventral thalamus do not. PTPH1 mRNA expression is also restricted to the dorsal thalamus during development and, thus, can serve as a specific marker for the dorsal thalamic nuclei. Since the subcellular localization of PTPH1 protein is not known, its functional role is not clear. However, the restriction of its expression to the thalamic nuclei that have thalamocortical connections suggests that PTPH1 may play a role in the maintenance of these connections or in determining the physiological properties of thalamic relay nuclei.

Binding of the transcription activator NRI (NTRC) to a supercoiled DNA segment imitates association with the natural enhancer: an electron microscopic investigation.

Electron microscopic visualization indicates that the transcription activator NRI (NTRC) binds with exceptional selectivity and efficiency to a sequence-induced superhelical (spiral) segment inserted upstream of the glnA promoter, accounting for its observed ability to substitute for the natural glnA enhancer. The cooperative binding of NRI to the spiral insert leads to protein oligomerization which, at higher concentration, promotes selective coating of the entire superhelical segment with protein. Localization of NRI at apical loops is observed with negatively supercoiled plasmid DNA. With a linear plasmid, bending of DNA is observed. We confirm that NRI is a DNA-bending protein, consistent with its high affinity for spiral DNA. These results prove that spiral DNA without any homology to the NRI-binding sequence site can substitute for the glnA enhancer by promoting cooperative activator binding to DNA and facilitating protein oligomerization. Similar mechanisms might apply to other prokaryotic and eukaryotic activator proteins that share the ability to bend DNA and act efficiently as multimers.