Ocean bottom seismic noise : applications for the crust knowledge, interaction ocean-atmosphere and instrumental behaviour

Tese de doutoramento, Ciências Geofísicas e da Geoinformação (Geofisíca), Universidade de Lisboa, Faculdade de Ciências, 2014

Neotectonics of the southwest Portugal mainland : implications on the regional seismic hazard

Tese de doutoramento, Geologia (Geodinâmica Interna), Universidade de Lisboa, Faculdade de Ciências, 2015

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.

Conditional probability applied to the definition of limestone geochemical facies: A new Methodological Approach

A new method, based on linear correlation and phase diagrams was successfully developed for processes like the sedimentary process, where the deposition phase can have different time duration - represented by repeated values in a series - and where the erosion can play an important rule deleting values of a series. The sampling process itself can be the cause of repeated values - large strata twice sampled - or deleted values: tiny strata fitted between two consecutive samples. What we developed was a mathematical procedure which, based upon the depth chemical composition evolution, allows the establishment of frontiers as well as the periodicity of different sedimentary environments. The basic tool isn't more than a linear correlation analysis which allow us to detect the existence of eventual evolution rules, connected with cyclical phenomena within time series (considering the space assimilated to time), with the final objective of prevision. A very interesting discovery was the phenomenon of repeated sliding windows that represent quasi-cycles of a series of quasi-periods. An accurate forecast can be obtained if we are inside a quasi-cycle (it is possible to predict the other elements of the cycle with the probability related with the number of repeated and deleted points). We deal with an innovator methodology, reason why it's efficiency is being tested in some case studies, with remarkable results that shows it's efficacy. Keywords: sedimentary environments, sequence stratigraphy, data analysis, time-series, conditional probability.

High resolution stratigraphy and miocene facies correlation in Lisbon and Setúbal Peninsula (Lower Tagus basin, Portugal)

Correlation between facies associations (marine, estuarine and distal fluviatile environments) and disconformities, observed between Foz da Fonte (SW of Setúbal Peninsula) and Santa Iria da Azóia (NE of Lisbon) are presented. The precise definition of the marine-continental facies relationships improved very much the chronology of the depositional sequence boundaries. Tectonic and eustatic controls are discussed on the basis of subsidence rates variation.

Facies, depositional environment, and palaeoecology of the Middle Triassic Cassina Beds (Meride Limestone, Monte San Giorgio, Switzerland).

The Ladinian Cassina beds belong to the fossiliferous levels of the world-famous Middle Triassic Monte San Giorgio Lagerstatte (UNESCO World Heritage List, Canton Ticino, Southern Alps). Although they are a rich archive for the depositional environment of an important thanatocoenosis, previous excavations focused on vertebrates and particularly on marine reptiles. In 2006, the Museo Cantonale di Storia Naturale (Lugano) started a new research project focusing for the first time on microfacies, micropalaeontological, palaeoecological and taphonomic analyses. So far, the upper third of the sequence has been excavated on a surface of around 40 m(2), and these new data complete those derived from new vertebrate finds (mainly fishes belonging to Saurichthys, Archaeosemionotus, Eosemionotus and Peltopleurus), allowing a better characterization of the basin. Background sedimentation on an anoxic to episodically suboxic seafloor resulted in a finely laminated succession of black shales and limestones, bearing a quasi-anaerobic biofacies, which is characterized by a monotypic benthic foraminiferal meiofauna and has been documented for the first time from the whole Monte San Giorgio sequence. Event deposition, testified by turbidites and volcaniclastic layers, is related to sediment input from basin margins and to distant volcanic eruptions, respectively. Fossil nekton points to an environment with only limited connection to the open sea. Terrestrial macroflora remains document the presence of emerged areas covered with vegetation and probably located relatively far away. Proliferation of benthic microbial mats is inferred on the basis of microfabrics, ecological considerations and taphonomic (both biostratinomic and diagenetic) features of the new vertebrate finds, whose excellent preservation is ascribed to sealing by biofilms. The occurrence of allochthonous elements allows an insight into the shallow-waters of the adjoining time-equivalent Salvatore platform. Finally, the available biostratigraphic data are critically reviewed.

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.

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.

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.

The application of ichnofabrics towards bridging the dichotomy between siliciclastic and carbonate shelf facies: examples from the Upper Jurassic Fulmar Formation (UK) and the Jubaila Formation (Saudi Arabia)

An initial study of the ichnofabrics of the Upper Jurassic (Kimmeridgian) Jubaila Formation of Saudi Arabia shows that the ichnofabrics are closely matched to the relatively well-described ichnofabrics of the contemporary Fulmar Formation of the UK Continental Shelf (North Sea), in respect of the lower shoreface/offshore transition facies to offshore facies. The ichnology and ichnofabrics of the Lower Jubaila Formation show that deposition took place on an open-marine platform on the Arabian craton subject to periodic storm activity, but under a persisting equilibrium between sediment accumulation and subsidence. This is consistent with the moderately deep-marine foraminiferal assemblages and the presence of calcareous nannofossils. Cyclicity is absent, though storm beds may be grouped, in contrast with the genetic sequences present in the rift and halokinetic scenario of the North Sea. In contrast with the siliciclastic setting hardgrotinds (with Gastrochaenolites), more common firmground omission surfaces, and micritic mudstones with Chondrites and Zoophycos are notable features of the carbonate facies. In siliciclastic successions (parasequences) the latter ichnotaxa are generally regarded as having been deposited in rather deeper water, but in the carbonate Jubaila Formation are interpreted as being associated with local areas of lower turbulence. Likewise, the hardgrounds and firmgrounds, which have not been traced laterally, are tentatively regarded to be of local significance.

A comparison of cross-hole electrical and seismic data in fractured rock

Cross-hole anisotropic electrical and seismic tomograms of fractured metamorphic rock have been obtained at a test site where extensive hydrological data were available. A strong correlation between electrical resistivity anisotropy and seismic compressional-wave velocity anisotropy has been observed. Analysis of core samples from the site reveal that the shale-rich rocks have fabric-related average velocity anisotropy of between 10% and 30%. The cross-hole seismic data are consistent with these values, indicating that observed anisotropy might be principally due to the inherent rock fabric rather than to the aligned sets of open fractures. One region with velocity anisotropy greater than 30% has been modelled as aligned open fractures within an anisotropic rock matrix and this model is consistent with available fracture density and hydraulic transmissivity data from the boreholes and the cross-hole resistivity tomography data. However, in general the study highlights the uncertainties that can arise, due to the relative influence of rock fabric and fluid-filled fractures, when using geophysical techniques for hydrological investigations.

A fuzzy multi-criteria decision approach for enhanced auto-tracking of seismic events

The main activity carried out by the geophysicist when interpreting seismic data, in terms of both importance and time spent is tracking (or picking) seismic events. in practice, this activity turns out to be rather challenging, particularly when the targeted event is interrupted by discontinuities such as geological faults or exhibits lateral changes in seismic character. In recent years, several automated schemes, known as auto-trackers, have been developed to assist the interpreter in this tedious and time-consuming task. The automatic tracking tool available in modem interpretation software packages often employs artificial neural networks (ANN's) to identify seismic picks belonging to target events through a pattern recognition process. The ability of ANNs to track horizons across discontinuities largely depends on how reliably data patterns characterise these horizons. While seismic attributes are commonly used to characterise amplitude peaks forming a seismic horizon, some researchers in the field claim that inherent seismic information is lost in the attribute extraction process and advocate instead the use of raw data (amplitude samples). This paper investigates the performance of ANNs using either characterisation methods, and demonstrates how the complementarity of both seismic attributes and raw data can be exploited in conjunction with other geological information in a fuzzy inference system (FIS) to achieve an enhanced auto-tracking performance.