Analysis and visualization of seismic data using mutual information

Seismic data is difficult to analyze and classical mathematical tools reveal strong limitations in exposing hidden relationships between earthquakes. In this paper, we study earthquake phenomena in the perspective of complex systems. Global seismic data, covering the period from 1962 up to 2011 is analyzed. The events, characterized by their magnitude, geographic location and time of occurrence, are divided into groups, either according to the Flinn-Engdahl (F-E) seismic regions of Earth or using a rectangular grid based in latitude and longitude coordinates. Two methods of analysis are considered and compared in this study. In a first method, the distributions of magnitudes are approximated by Gutenberg-Richter (G-R) distributions and the parameters used to reveal the relationships among regions. In the second method, the mutual information is calculated and adopted as a measure of similarity between regions. In both cases, using clustering analysis, visualization maps are generated, providing an intuitive and useful representation of the complex relationships that are present among seismic data. Such relationships might not be perceived on classical geographic maps. Therefore, the generated charts are a valid alternative to other visualization tools, for understanding the global behavior of earthquakes.

An interpretation of the deep seismic lines from the Penninic Alps of Valais (Switzerland)

The upper part of three deep seismic lines running across the Penninic Swiss Alps of Valais have been studied. Numerous reflectors illustrate the nappe structure of this internal part of the orogen. These reflectors, even at great depths (20-25 km), can be correlated with outcropping geological features and are most likely produced by lithological boundaries rather than by mylonites zones, which are hardly reflective in such an environment. Our interpretations, largely constrained by projections of the outcropping geology, have improved our knowledge of the deep structure of this segment of the Alpine belt, enhancing the importance of the backfolding and the crustal scale deformation phase which produced the Rawil-Valpelline depression and the Aar-Toce culmination. Furthermore we have here the possibility of correlating seismic patterns produced by ductile folds with the outcropping structures.

Dealing with Moral Hazard and Adverse Selection Simultaneously

Although Insurers Face Adverse Selection and Moral Hazard When They Set Insurance Contracts, These Two Types of Asymmetrical Information Have Been Given Separate Treatments Sofar in the Economic Literature. This Paper Is a First Attempt to Integrate Both Problems Into a Single Model. We Show How It Is Possible to Use Time in Order to Achieve a First-Best Allocation of Risks When Both Problems Are Present Simultaneously.

Moral Hazard and Marshallian Inefficiency: Evidence from Tunisia.

In the presence of moral hazard, received agency theory predicts the Marshallian inefficiency of agricultural tenancy contracts, meaning that inputs per hectare on sharecropped land will differ from that on owned land. in this paper, we test for the presence of Marshallian inefficiency using a unique data set collected in the Tunisian village of El Oulja in 1993.

Modelling and Analysis of Bivariate Lifetime Data using Reversed Hazard Rates

Department of Statistics, Cochin University of Science and Technology

Seismic behavior of capacity designed masonry walls in low seismicity regions

Eurocode 8 representing a new generation of structural design codes in Europe defines ‎requirements for the design of buildings against earthquake action. In Central and ‎Western Europe, the newly defined earthquake zones and corresponding design ground ‎acceleration values, will lead in many cases to earthquake actions which are remarkably ‎higher than those defined so far by the design codes used until now in Central Europe. ‎ In many cases, the weak points of masonry structures during an earthquake are the corner ‎regions of the walls. Loading of masonry walls by earthquake action leads in most cases ‎to high shear forces. The corresponding bending moment in such a wall typically causes a ‎significant increase of the eccentricity of the normal force in the critical wall cross ‎section. This in turn leads ultimately to a reduction of the size of the compression zone in ‎unreinforced walls and a high concentration of normal stresses and shear stresses in the ‎corner regions. ‎ Corner-Gap-Elements, consisting of a bearing beam located underneath the wall and ‎made of a sufficiently strong material (such as reinforced concrete), reduce the effect of ‎the eccentricity of the normal force and thus restricts the pinching effect of the ‎compression zone. In fact, the deformation can be concentrated in the joint below the ‎bearing beam. According to the principles of the Capacity Design philosophy, the ‎masonry itself is protected from high stresses as a potential cause of brittle failure. ‎ Shaking table tests at the NTU Athens Earthquake Engineering Laboratory have proven ‎the effectiveness of the Corner-Gap-Element. The following presentation will cover the ‎evaluation of various experimental results as well as a numerical modeling of the ‎observed phenomena.‎

3D Simultaneous Inversion for Seismic Structure and Local Hypocenters in Germany under Consideration of Seismic Anisotropy in the Upper Mantel

The main task of this work has been to investigate the effects of anisotropy onto the propagation of seismic waves along the Upper Mantle below Germany and adjacent areas. Refraction- and reflexion seismic experiments proved the existence of Upper Mantle anisotropy and its influence onto the propagation of Pn-waves. By the 3D tomographic investigations that have been done here for the crust and the upper mantle, considering the influence of anisotropy, a gap for the investigations in Europe has been closed. These investigations have been done with the SSH-Inversionprogram of Prof. Dr. M. Koch, which is able to compute simultaneously the seismic structure and hypocenters. For the investigation, a dataset has been available with recordings between the years 1975 to 2003 with a total of 60249 P- and 54212 S-phase records of 10028 seismic events. At the beginning, a precise analysis of the residuals (RES, the difference between calculated and observed arrivaltime) has been done which confirmed the existence of anisotropy for Pn-phases. The recognized sinusoidal distribution has been compensated by an extension of the SSH-program by an ellipse with a slow and rectangular fast axis with azimuth to correct the Pn-velocities. The azimuth of the fast axis has been fixed by the application of the simultaneous inversion at 25° - 27° with a variation of the velocities at +- 2.5 about an average value at 8 km/s. This new value differs from the old one at 35°, recognized in the initial residual analysis. This depends on the new computed hypocenters together with the structure. The application of the elliptical correction has resulted in a better fit of the vertical layered 1D-Model, compared to the results of preceding seismological experiments and 1D and 2D investigations. The optimal result of the 1D-inversion has been used as initial starting model for the 3D-inversions to compute the three dimensional picture of the seismic structure of the Crust and Upper Mantle. The simultaneous inversion has showed an optimization of the relocalization of the hypocenters and the reconstruction of the seismic structure in comparison to the geology and tectonic, as described by other investigations. The investigations for the seismic structure and the relocalization have been confirmed by several different tests. First, synthetic traveltime data are computed with an anisotropic variation and inverted with and without anisotropic correction. Further, tests with randomly disturbed hypocenters and traveltime data have been proceeded to verify the influence of the initial values onto the relocalization accuracy and onto the seismic structure and to test for a further improvement by the application of the anisotropic correction. Finally, the results of the work have been applied onto the Waldkirch earthquake in 2004 to compare the isotropic and the anisotropic relocalization with the initial optimal one to verify whether there is some improvement.

Genetic testing with primary prevention and moral hazard

We develop a model where a free genetic test reveals whether the individual tested has a low or high probability of developing a disease. A costly prevention effort allows high-risk agents to decrease the probability of developing the disease. Agents are not obliged to take the test, but must disclose its results to insurers. Insurers offer separating contracts which take into account the individual risk, so that taking the test is associated to a discrimination risk. We study the individual decisions to take the test and to undertake the prevention effort as a function of the effort cost and of its e¢ ciency. We obtain that, if effort is observable by insurers, agents undertake the test only if the effort cost is neither too large nor too low. If the effort cost is not observable by insurers, they face a moral hazard problem which induces them to under-provide insurance. We obtain the counterintuitive result that moral hazard increases the value of the test if the effort cost is low enough. Also, agents may perform the test for lower levels of prevention e¢ ciency when effort is not observable

Reconocimiento geofísico del deslizamiento de Rosiana (depresión de Tirajana, Gran Canaria)

The landslide of Rosiana is considered the largest slope movement amongst those known in historical times in Gran Canana, Canary Islands. It has been activated at least 4 times in the last century, and in the movement of 1956, when about 3.106 m3 of materials were involved, 250 people had to be evacuated and many buildings were destroyed. The present geological hazard has lead to specific studies of the phenomenon which, once characterised, can be used as a guide for the scientific and technical works that are to be made in this or similar areas. This paper wants to increase the knowledge about the unstable mass of Rosiana by using geophysical techniques based on the method of seismic by refraction. The geophysical measues have been interpreted with the aid of the available geomorphologic data, thus obtaining a first approximation to the geometry of the slope movements

Geostatistics for constrained variables: positive data, compositions and probabilities. Applications to environmental hazard monitoring

Aquesta tesi estudia com estimar la distribució de les variables regionalitzades l'espai mostral i l'escala de les quals admeten una estructura d'espai Euclidià. Apliquem el principi del treball en coordenades: triem una base ortonormal, fem estadística sobre les coordenades de les dades, i apliquem els output a la base per tal de recuperar un resultat en el mateix espai original. Aplicant-ho a les variables regionalitzades, obtenim una aproximació única consistent, que generalitza les conegudes propietats de les tècniques de kriging a diversos espais mostrals: dades reals, positives o composicionals (vectors de components positives amb suma constant) són tractades com casos particulars. D'aquesta manera, es generalitza la geostadística lineal, i s'ofereix solucions a coneguts problemes de la no-lineal, tot adaptant la mesura i els criteris de representativitat (i.e., mitjanes) a les dades tractades. L'estimador per a dades positives coincideix amb una mitjana geomètrica ponderada, equivalent a l'estimació de la mediana, sense cap dels problemes del clàssic kriging lognormal. El cas composicional ofereix solucions equivalents, però a més permet estimar vectors de probabilitat multinomial. Amb una aproximació bayesiana preliminar, el kriging de composicions esdevé també una alternativa consistent al kriging indicador. Aquesta tècnica s'empra per estimar funcions de probabilitat de variables qualsevol, malgrat que sovint ofereix estimacions negatives, cosa que s'evita amb l'alternativa proposada. La utilitat d'aquest conjunt de tècniques es comprova estudiant la contaminació per amoníac a una estació de control automàtic de la qualitat de l'aigua de la conca de la Tordera, i es conclou que només fent servir les tècniques proposades hom pot detectar en quins instants l'amoni es transforma en amoníac en una concentració superior a la legalment permesa.

Ovenbird (Seiurus aurocapilla) Territory Placement Near Seismic Lines is Influenced by Forest Regeneration and Conspecific Density

The boreal forest of western Canada is being dissected by seismic lines used for oil and gas exploration. The vast amount of edge being created is leading to concerns that core habitat will be reduced for forest interior species for extended periods of time. The Ovenbird (Seiurus aurocapilla) is a boreal songbird known to be sensitive to newly created seismic lines because it does not include newly cut lines within its territory. We examined multiple hypotheses to explain potential mechanisms causing this behavior by mapping Ovenbird territories near lines with varying states of vegetation regeneration. The best model to explain line exclusion behavior included the number of neighboring conspecifics, the amount of bare ground, leaf-litter depth, and canopy closure. Ovenbirds exclude recently cut seismic lines from their territories because of lack of protective cover (lower tree and shrub cover) and because of reduced food resources due to large areas of bare ground. Food reduction and perceived predation risk effects seem to be mitigated once leaf litter (depth and extent of cover) and woody vegetation cover are restored to forest interior levels. However, as conspecific density increases, lines are more likely to be used as landmarks to demarcate territorial boundaries, even when woody vegetation cover and leaf litter are restored. This behavior can reduce territory density near seismic lines by changing the spatial distribution of territories. Landmark effects are longer lasting than the effects from reduced food or perceived predation risk because canopy height and tree density take >40 years to recover to forest interior levels. Mitigation of seismic line impacts on Ovenbirds should focus on restoring forest cover as quickly as possible after line cutting.