Simulation of the influence of rate- and state-dependent friction on the macroscopic behavior of complex fault zones with the lattice solid model

In order to understand the earthquake nucleation process, we need to understand the effective frictional behavior of faults with complex geometry and fault gouge zones. One important aspect of this is the interaction between the friction law governing the behavior of the fault on the microscopic level and the resulting macroscopic behavior of the fault zone. Numerical simulations offer a possibility to investigate the behavior of faults on many different scales and thus provide a means to gain insight into fault zone dynamics on scales which are not accessible to laboratory experiments. Numerical experiments have been performed to investigate the influence of the geometric configuration of faults with a rate- and state-dependent friction at the particle contacts on the effective frictional behavior of these faults. The numerical experiments are designed to be similar to laboratory experiments by DIETERICH and KILGORE (1994) in which a slide-hold-slide cycle was performed between two blocks of material and the resulting peak friction was plotted vs. holding time. Simulations with a flat fault without a fault gouge have been performed to verify the implementation. These have shown close agreement with comparable laboratory experiments. The simulations performed with a fault containing fault gouge have demonstrated a strong dependence of the critical slip distance D-c on the roughness of the fault surfaces and are in qualitative agreement with laboratory experiments.

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.

A microstrip Yagi antenna using EBG structure

[1] In this paper a detailed design, development and performances of a 5 GHz microstrip Yagi antenna, which uses a two-dimensional (2-D) electromagnetic band gap (EBG) structure in the ground plane, are presented. The results indicate that the use of the EBG structure improves the radiation pattern of the antenna. The cross polarization is suppressed by properly choosing the period and dimensions of EBGs. Also, the broadside gain is improved in comparison with the analogous antenna without the EBGs.

Simultaneous observations of ionospheric quasiperiodic scintillations from short and long meridional baselines using VHF transmissions from transit satellites

For the first time it was possible to observe regular quasiperiodic scintillations (QPS) in VHF radio-satellite transmissions from orbiting satellites simultaneously at short (2.1 km) and long (121 km) meridional baselines in the vicinity of a typical mid-latitude station (Brisbane; 27.5degreesS and 152.9degreesE geog. and 35.6degrees invar.lat.), using three sites (St. Lucia-S, Taringa-T in Brisbane and Boreen Pt.-B, north of Brisbane). A few pronounced quasiperiodic (QP) events were recorded showing unambiguous regular structures at the sites which made it possible to deduce a time displacement of the regular fading minimum at S, T and B. The QP structure is highly dependent on the geometry of the ray-path from a satellite to the observer which is manifested as a change of a QP event from symmetrical to non-symmetrical for stations separated by 2.1 km, and to a radical change in the structure of the event over a distance of 121 km. It is suggested the short-duration intense QP events are due to a Fresnel diffraction (or a reflection mechanism) of radio-satellite signals by a single ionospheric irregularity in a form of an ellipsoid with a large ionization gradient along the major axis. The structure of a QP event depends on the angle of viewing of the irregular blob from a radio-satellite. In view of this it is suggested that the reported variety of the ionization formation, responsible for different types of QPS, is only apparent but not real. (C) 2003 Elsevier Science Ltd. All rights reserved.

In search of the ultimate ionosonde deduced equator

A further progress has been made in defining the ionosonde deduced equator (IDE) which characterises a latitudinal transition from the northern to southern hemisphere. It is now possible to define the global IDE location as the locus of the average position between geographic and geomagnetic equators. A more complete insight to the phenomenon of the third equator (i.e. after geographic and geomagnetic equators) was made possible due to availability of ionospheric height (h'F) data from three stations positioned close to the IDE in the American and the far-east sectors. The IDE ionospheric signature (or E-type signature), detected at these stations, consists of bi-annual h'F height increases. This signature however is not consistently observed during solar cycle and at times, particularly at sunspot minimum, a weak hemispheric signature is observed (i.e. the northern or southern hemisphere signature). In general, the height increase at the IDE are considerably smaller (by a factor of 4) than at other equatorial locations, indicating that the ionosphere at the IDE location becomes less disturbed. It is suggested that the equatorial longitudinal regions which can be associated with more consistent E-type signature are located in the central Pacific and at the east coast of America, close to the intersection points of the geographic and geomagnetic equators. (C) 2003 Elsevier Ltd. All rights reserved.

Speleothem U-series dating of semi-synchronous climate oscillations during the last deglaciation

Evidence for nearly synchronous climate oscillations during the last deglaciation has been found throughout the Northern Hemisphere but few records are based on independent time scales of calendar years. We present a rare uranium-series dated oxygen-carbon isotope record for a speleothem from Tangshan Cave, China, which demonstrates that abrupt deglacial climatic oscillations from 16 800 to 10 500 yr BP are semi-synchronous with those found in Greenland ice core records. Relatively rapid shifts in speleothem oxygen isotope ratios demonstrate that the intensity of the East Asian monsoon switched in parallel with the abrupt transitions separating the Bolling-Allerod, Younger Dryas, and pre-Boreal climatic reversals. However, the dated isotopic transitions appear to have lasted longer. Our results demonstrate the dominant role of atmospheric teleconnections in the rapid propagation of deglacial climatic signals on a hemispheric scale, and highlight the importance of U-series dated speleothems in the timing and characterization of abrupt climate change. (C) 2003 Elsevier B.V. All rights reserved.

Evaluating Tsunami Impact on the Gulf of Cadiz Coast (Northeast Atlantic).

The Gulf of Cadiz coasts are exposed to tsunamis. Emergency planning tools are now taking into account this fact, especially because a series of historical occurrences were strikingly significant, having left strong evidence behind, in the mareographic records, the geological evidence or simply the memory of the populations. The study area is a strip along the Algarve coast, south Portugal, an area known to have been heavily impacted by the 1 November 1755 event. In this study we use two different tsunami scenarios generated by the rupture of two thrust faults identified in the area, corresponding to 8.1-8.3 magnitude earthquakes. Tsunami propagation and inundation computation is performed using a non-linear shallow water code with bottom friction. Numerical modeling results are presented in terms of flow depth and current velocity with maximum values of 7 m and 8 m/s for inundation depth and flow speed, respectively. These results constitute a valuable tool for local authorities, emergency and decision planners to define the priority zones where tsunami mitigation measures must be implemented and to develop tsunami-resilient communities.

Interpretation of gravity data to delineate structural features connected to low-temperature geothermal resources at Northeastern Portugal

A great number of low-temperature geothermal fields occur in Northern-Portugal related to fractured rocks. The most important superficial manifestations of these hydrothermal systems appear in pull-apart tectonic basins and are strongly conditioned by the orientation of the main fault systems in the region. This work presents the interpretation of gravity gradient maps and 3D inversion model produced from a regional gravity survey. The horizontal gradients reveal a complex fault system. The obtained 3D model of density contrast puts into evidence the main fault zone in the region and the depth distribution of the granitic bodies. Their relationship with the hydrothermal systems supports the conceptual models elaborated from hydrochemical and isotopic water analyses. This work emphasizes the importance of the role of the gravity method and analysis to better understand the connection between hydrothermal systems and the fractured rock pattern and surrounding geology. (c) 2013 Elsevier B.V. All rights reserved.

Avaliação multitécnica em hidrogeologia ambiental e suas implicações na gestão de recursos hídricos subterrâneos

Mestrado em Engenharia Geotécnica e Geoambiente

Desenvolvimento de uma plataforma SIG para apoio ao estudo da sismicidade em Portugal

Trabalho Final de Mestrado elaborado no Laboratório Nacional de Engenharia Civil (LNEC) para a obtenção do grau de Mestre em Engenharia Civil pelo Instituto Superior de Engenharia de Lisboa no âmbito do protocolo de cooperação entre o ISEL e o LNEC

Prospecção de recursos hidrominerais e geotérmicos na Orla Ocidental Portuguesa: os casos das Termas dos Cucos e de Monte Real

A cobertura sedimentar da região Oeste portuguesa é constituída por uma série possante de sedimentos com uma variedade de fácies com idades compreendidas entre o Triásico Superior e o actual. Estes sedimentos foram depositados numa bacia alongada com direcção NNE‐SSW. A tectónica desta cobertura sedimentar é condicionada pelas falhas tardi‐Variscas que afectaram o substrato e pelo complexo evaporítico depositado na base das séries sedimentares. Séries evaporíticas espessas de idade Hetangiana formaram numerosas estruturas diapíricas. Na região Oeste de Portugal existem diversas nascentes minerais e termais usadas para hidroterapia. Esta tese tem por objectivo estudar as concessões hidrominerais das Termas dos Cucos e das Termas de Monte Real (Portugal Central), bem como as suas áreas envolventes. Estas actividades hidroterapêuticas são muito relevantes em termos sócio‐económicos para os concelhos de Torres Vedras e Leiria. Os estudos contemplados nesta tese (que incluíram trabalho de campo no domínio da geotectónica, geomorfologia e hidrogeologia) permitiram delinear modelos hidrogeológicos conceptuais, apoiados, ainda, pela re‐interpretação de estudos geofísicos e hidrogeológicos prévios. A caracterização destas áreas foi suportada por inventários hidrogeológicos, tendo sido determinante para o projecto dos furos de captação, incluindo localização e profundidade. Todos os dados compilados foram representados cartograficamente numa base de Sistemas de Informação Geográfica (SIG).

Caracterização de estruturas tectónicas activas da região sul de Portugal com recurso a ferramentas SIG: o caso da falha de São Marcos - Quarteira

Dissertação apresentada como requisito parcial para obtenção do grau de Mestre em Ciência e Sistemas de Informação Geográfica

Probabilistic tsunami hazard in the Northeast Atlantic from near- and far-field tectonic sources

In this article, we present the first study on probabilistic tsunami hazard assessment for the Northeast (NE) Atlantic region related to earthquake sources. The methodology combines the probabilistic seismic hazard assessment, tsunami numerical modeling, and statistical approaches. We consider three main tsunamigenic areas, namely the Southwest Iberian Margin, the Gloria, and the Caribbean. For each tsunamigenic zone, we derive the annual recurrence rate for each magnitude range, from Mw 8.0 up to Mw 9.0, with a regular interval, using the Bayesian method, which incorporates seismic information from historical and instrumental catalogs. A numerical code, solving the shallow water equations, is employed to simulate the tsunami propagation and compute near shore wave heights. The probability of exceeding a specific tsunami hazard level during a given time period is calculated using the Poisson distribution. The results are presented in terms of the probability of exceedance of a given tsunami amplitude for 100- and 500-year return periods. The hazard level varies along the NE Atlantic coast, being maximum along the northern segment of the Morocco Atlantic coast, the southern Portuguese coast, and the Spanish coast of the Gulf of Cadiz. We find that the probability that a maximum wave height exceeds 1 m somewhere in the NE Atlantic region reaches 60 and 100 % for 100- and 500-year return periods, respectively. These probability values decrease, respectively, to about 15 and 50 % when considering the exceedance threshold of 5 m for the same return periods of 100 and 500 years.

On the applicability of joint inversion of gravity and resistivity data to the study of a tectonic sedimentary basin in Northern Portugal

The Chaves basin is a pull-apart tectonic depression implanted on granites, schists, and graywackes, and filled with a sedimentary sequence of variable thickness. It is a rather complex structure, as it includes an intricate network of faults and hydrogeological systems. The topography of the basement of the Chaves basin still remains unclear, as no drill hole has ever intersected the bottom of the sediments, and resistivity surveys suffer from severe equivalence issues resulting from the geological setting. In this work, a joint inversion approach of 1D resistivity and gravity data designed for layered environments is used to combine the consistent spatial distribution of the gravity data with the depth sensitivity of the resistivity data. A comparison between the results from the inversion of each data set individually and the results from the joint inversion show that although the joint inversion has more difficulty adjusting to the observed data, it provides more realistic and geologically meaningful models than the ones calculated by the inversion of each data set individually. This work provides a contribution for a better understanding of the Chaves basin, while using the opportunity to study further both the advantages and difficulties comprising the application of the method of joint inversion of gravity and resistivity data.

Estimation of the Crustal Bulk Properties Beneath Mainland Portugal from P-Wave Teleseismic Receiver Functions

In this work, we present results from teleseismic P-wave receiver functions (PRFs) obtained in Portugal, Western Iberia. A dense seismic station deployment conducted between 2010 and 2012, in the scope of the WILAS project and covering the entire country, allowed the most spatially extensive probing on the bulk crustal seismic properties of Portugal up to date. The application of the H-κ stacking algorithm to the PRFs enabled us to estimate the crustal thickness (H) and the average crustal ratio of the P- and S-waves velocities V p/V s (κ) for the region. Observations of Moho conversions indicate that this interface is relatively smooth with the crustal thickness ranging between 24 and 34 km, with an average of 30 km. The highest V p/V s values are found on the Mesozoic-Cenozoic crust beneath the western and southern coastal domain of Portugal, whereas the lowest values correspond to Palaeozoic crust underlying the remaining part of the subject area. An average V p/V s is found to be 1.72, ranging 1.63-1.86 across the study area, indicating a predominantly felsic composition. Overall, we systematically observe a decrease of V p/V s with increasing crustal thickness. Taken as a whole, our results indicate a clear distinction between the geological zones of the Variscan Iberian Massif in Portugal, the overall shape of the anomalies conditioned by the shape of the Ibero-Armorican Arc, and associated Late Paleozoic suture zones, and the Meso-Cenozoic basin associated with Atlantic rifting stages. Thickened crust (30-34 km) across the studied region may be inherited from continental collision during the Paleozoic Variscan orogeny. An anomalous crustal thinning to around 28 km is observed beneath the central part of the Central Iberian Zone and the eastern part of South Portuguese Zone.