A matrix method for elastic wave problems

The solution to problems of elastic wave propagation in multilayered media, in which each layer is homogeneous and where the ensemble of layers has physical properties that vary only with one coordinate, may be given as the quotient of products of matrices. In the case of SH waves, the matrices are of order two; in the case of P-SV waves the matrices are of order four. The individual matrix elements are themselves determinants of order two or four in the two cases. The solution is obtained by means of Laplace's development by minors

Observations of multiple seismic events

Two or more dispersed wave trains each with constant amplitude will interfere giving a resultant wave train which is amplitude modulated, if the individual waves have their principal energies in a common frequency band and if the trains arrive with time separations small compared to their total length. The dispersive characteristics of the trains need not be the same. If the component trains are of comparable magnitude, the modulation due to interference becomes significant and a "beat" phenomenon occurs. Multiple trains of dispersed seismic surface waves may occur because of a temporal and/or spatial distribution at the source or because of multipath propagation. Each of these causal mechanisms influences the amplitude and phase spectra of the resultant wave train; derived properties such as phase velocities and amplitude ratios are also influenced. In the case of multipath propagation, wavelength dependent time delays may occur. Two cases of twin earthquakes are analyzed, and the significant features of interference are demonstrated. In one case, estimates are obtained for the amplitude ratio and time delay of the second shock with respect to the first. The interpretation of seismograms and spectra influenced by multiple events is discussed

Rank-Ordering Statistics of Extreme Events: Application to the Distribution of Large Earthquakes

Rank-ordering statistics provides a perspective on the rare, largest elements of a population, whereas the statistics of cumulative distributions are dominated by the more numerous small events. The exponent of a power law distribution can be determined with good accuracy by rank-ordering statistics from the observation of only a few tens of the largest events. Using analytical results and synthetic tests, we quantify the systematic and the random errors. We also study the case of a distribution defined by two branches, each having a power law distribution, one defined for the largest events and the other for smaller events, with application to the World-Wide (Harvard) and Southern California earthquake catalogs. In the case of the Harvard moment catalog, we make more precise earlier claims of the existence of a transition of the earthquake magnitude distribution between small and large earthquakes; the $b$-values are $b_2 = 2.3 \pm 0.3$ for large shallow earthquakes and $b_1 = 1.00 \pm 0.02$ for smaller shallow earthquakes. However, the cross-over magnitude between the two distributions is ill-defined. The data available at present do not provide a strong constraint on the cross-over which has a $50\%$ probability of being between magnitudes $7.1$ and $7.6$ for shallow earthquakes; this interval may be too conservatively estimated. Thus, any influence of a universal geometry of rupture on the distribution of earthquakes world-wide is ill-defined at best. We caution that there is no direct evidence to confirm the hypothesis that the large-moment branch is indeed a power law. In fact, a gamma distribution fits the entire suite of earthquake moments from the smallest to the largest satisfactorily. There is no evidence that the earthquakes of the Southern California catalog have a distribution with tw

Explaining coexistence of nitrogen fixing and non-fixing rhizobia in legume-rhizobia mutualism using mathematical modeling

In the mutualism established between legumes and soil bacteria known as rhizobia, bacteria from soil infect plants roots and reproduce inside root nodules where they fix atmospheric N2 for plant nutrition, receiving carbohydrates in exchange. Host-plant sanctions against non N2 fixing, cheating bacterial symbionts have been proposed to act in the legume-Rhizobium symbiosis, to preserve the mutualistic relationship. Sanctions include decreased rhizobial survival in nodules occupied by cheating rhizobia. Previously, a simple population model experimentally based showed that the coexistence of fixing and cheating rhizobia strains commonly found in field conditions is possible, and that the inclusion of sanctions leads to the extinction of cheating strains in soil. Here, we extend the previous model to include other factors that could complicate the sanction scenario, like horizontal transmission of symbiotic plasmids, turning non-nodulating strains into nodulating rhizobia, and competition between fixing and cheating strains for nodulation. In agreement with previous results, we show that plant populations persist even in the presence of cheating rhizobia without incorporating any sanction against the cheater populations in the model, under the realistic assumption that plants can at least get some amount of fixed N2 from the effectively mutualistic rhizobia occupying some nodules. Inclusion of plant sanctions leads to the unrealistic extinction of cheater strains in soil. Our results agree with increasing experimental evidence and theoretical work showing that mutualisms can persist in presence of cheating partners.Fil: Moyano, Gabriel Eduardo. Universidad Nacional de Córdoba. Facultad de Matemática, Astronomía y Física; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Marco, Diana Elizabeth. Universidad Nacional de Córdoba. Facultad de Cs.exactas Físicas y Naturales. Departamento de Matemáticas; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas; ArgentinaFil: Knopoff, Damián Alejandro. Universidad Nacional de Córdoba. Facultad de Matemática, Astronomia y Física. Sección Matemática; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Centro de Investigación y Estudios de Matemática. Universidad Nacional de Córdoba. Centro de Investigación y Estudios de Matemática; ArgentinaFil: Torres, German Ariel. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Nordeste. Instituto de Modelado e Innovación Tecnológica. Universidad Nacional del Nordeste. Facultad de Ciencias Exactas Naturales y Agrimensura. Instituto de Modelado e Innovación Tecnológica; ArgentinaFil: Turner, Cristina Vilma. Universidad Nacional de Córdoba. Facultad de Matemática, Astronomía y Física; Argentina. Consejo Nacional de Investigaciones Científicas y Técnicas. Centro Científico Tecnológico Conicet - Córdoba. Centro de Investigación y Estudios de Matemática. Universidad Nacional de Córdoba. Centro de Investigación y Estudios de Matemática; Argentin

First motions from seismic sources near a free surface

The radiation patterns of first motions are calculated for the sudden occurrence of an arbitrarily oriented fault (dislocation) at the surface of a half space; the dislocation in the fault plane is also arbitrarily oriented and is assumed to occur over a very small area of the fault plane. Initially the source is considered at a finite depth and the solution is obtained by allowing the depth to tend to zero. In general the results show a surprising directionality for the radiation of SV. In the focal plane projection the first motions of P and SH for a strike-slip fault show the familiar four-lobed radiation patterns. The first motions of SV show some reversals in polarity with angular distance from the source. The first motions for all components of motion for a dip-slip fault have characteristics governed strongly by the presence of the free surface, and hence differ markedly from the usual radiation patterns for a deeply imbedded source