The source of the 1722 Algarve earthquake: evidence from MCS and tsunami data

The 27 December 1722 Algarve earthquake destroyed a large area in southern Portugal generating a local tsunami that inundated the shallow areas of Tavira. It is unclear whether its source was located onshore or offshore and, in any case, what was the tectonic source responsible for the event. We analyze available historical information concerning macroseismicity and the tsunami to discuss the most probable location of the source. We also review available seismotectonic knowledge of the offshore region close to the probable epicenter, selecting a set of four candidate sources. We simulate tsunamis produced by these candidate sources assuming that the sea bottom displacement is caused by a compressive dislocation over a rectangular fault, as given by the half-space homogeneous elastic approach, and we use numerical modeling to study wave propagation and run-up. We conclude that the 27 December 1722 Tavira earthquake and tsunami was probably generated offshore, close to 37 degrees 01'N, 7 degrees 49'W.

The 1531 Lisbon Earthquake: A Tsunami In the Tagus Estuary?

On 26 January 1531, a strong-magnitude earthquake heavily impacted Lisbon downtown. Immediately after the earthquake, the eyewitnesses reported large waves in the Tagus estuary, mainly north of the city and along the northern bank of the river. Descriptions include large impacts on ships anchored in the estuary and even morphological changes in the riverbed. We present a synthesis of the available information concerning both the earthquake and the water disturbance as a basis for the discussion of the probable tectonic source and the magnitude of the associated river oscillations. We hypothesize that the initial disturbance of the water can be attributed to the coseismic deformation of the estuary riverbed, and we use a nonlinear shallow water model to simulate the tsunami propagation and inundation. We show that the Vila Franca de Xira fault is the most probable source of the 1531 event. The largest inundation effects of the model correlate well with the historical descriptions: the impact is relevant in the inner Tagus estuary, but inundation in downtown Lisbon is small.

The 1998 Faial earthquake, Azores: Evidence for a transform fault associated with the Nubia-Eurasia plate boundary?

With very few exceptions, M > 4 tectonic earthquakes in the Azores show normal fault solution and occur away from the islands. Exceptionally, the 1998 shock was pure strike-slip and occurred within the northern edge of the Pico-Faial Ridge. Fault plane solutions show two possible planes of rupture striking ENE-WSW (dextral) and NNW-SSE (sinistral). The former has not been recognised in the Azores, but is parallel to the transform direction related to the relative motion between the Eurasia and Nubia plates. Therefore, the main question we address in the present study is: do transform faults related to the Eurasia/Nubia plate boundary exist in the Azores? Knowing that the main source of strain is related to plate kinematics, we conclude that the sinistral strike-slip NNW-SSE fault plane solution is not consistent with either the fault dip (ca. 65, which is typical of a normal fault) or the ca. ENE-WSW direction of maximum extension; both are consistent with a normal fault, as observed in most major earthquakes on faults striking around NNW-SSE in the Azores. In contrast, the dextral strike-slip ENE-WSW fault plane solution is consistent with the transform direction related to the anticlockwise rotation of Nubia relative to Eurasia. Altogether, tectonic data, measured ground motion, observed destruction, and modelling are consistent with a dextral strike-slip source fault striking ENE-WSW. Furthermore, the bulk clockwise rotation measured by GPS is typical of bookshelf block rotations observed at the termination of such master strike-slip faults. Therefore, we suggest that the 1998 earthquake can be related to the WSW termination of a transform (ENE-WSW fault plane solution) associated with the Nubia-Eurasia diffuse plate boundary. (C) 2014 Elsevier B.V. All rights reserved.

Design optimization of cruciform specimens for biaxial fatigue loading

In order to correctly assess the biaxial fatigue material properties one must experimentally test different load conditions and stress levels. With the rise of new in-plane biaxial fatigue testing machines, using smaller and more efficient electrical motors, instead of the conventional hydraulic machines, it is necessary to reduce the specimen size and to ensure that the specimen geometry is appropriated for the load capacity installed. At the present time there are no standard specimen’s geometries and the indications on literature how to design an efficient test specimen are insufficient. The main goal of this paper is to present the methodology on how to obtain an optimal cruciform specimen geometry, with thickness reduction in the gauge area, appropriated for fatigue crack initiation, as a function of the base material sheet thickness used to build the specimen. The geometry is optimized for maximum stress using several parameters, ensuring that in the gauge area the stress is uniform and maximum with two limit phase shift loading conditions. Therefore the fatigue damage will always initiate on the center of the specimen, avoiding failure outside this region. Using the Renard Series of preferred numbers for the base material sheet thickness as a reference, the reaming geometry parameters are optimized using a derivative-free methodology, called direct multi search (DMS) method. The final optimal geometry as a function of the base material sheet thickness is proposed, as a guide line for cruciform specimens design, and as a possible contribution for a future standard on in-plane biaxial fatigue tests. © 2014, Gruppo Italiano Frattura. All rights reserved.

Optimization of cruciform specimens for biaxial fatigue loading with direct multi search

In order to correctly assess the biaxial fatigue material properties one must experimentally test different load conditions and stress levels. With the rise of new in-plane biaxial fatigue testing machines, using smaller and more efficient electrical motors, instead of the conventional hydraulic machines, it is necessary to reduce the specimen size and to ensure that the specimen geometry is appropriate for the load capacity installed. At the present time there are no standard specimen's geometries and the indications on literature how to design an efficient test specimen are insufficient. The main goal of this paper is to present the methodology on how to obtain an optimal cruciform specimen geometry, with thickness reduction in the gauge area, appropriate for fatigue crack initiation, as a function of the base material sheet thickness used to build the specimen. The geometry is optimized for maximum stress using several parameters, ensuring that in the gauge area the stress distributions on the loading directions are uniform and maximum with two limit phase shift loading conditions (delta = 0 degrees and (delta = 180 degrees). Therefore the fatigue damage will always initiate on the center of the specimen, avoiding failure outside this region. Using the Renard Series of preferred numbers for the base material sheet thickness as a reference, the reaming geometry parameters are optimized using a derivative-free methodology, called direct multi search (DMS) method. The final optimal geometry as a function of the base material sheet thickness is proposed, as a guide line for cruciform specimens design, and as a possible contribution for a future standard on in-plane biaxial fatigue tests

Educating for earthquake science and risk in a tectonically slowly deforming region

Over the past decade, scientists have been called to participate more actively in public education and outreach (E&O). This is particularly true in fields of significant societal impact, such as earthquake science. Local earthquake risk culture plays a role in the way that the public engages in educational efforts. In this article, we describe an adapted E&O program for earthquake science and risk. The program is tailored for a region of slow tectonic deformation, where large earthquakes are extreme events that occur with long return periods. The adapted program has two main goals: (1) to increase the awareness and preparedness of the population to earthquake and related risks (tsunami, liquefaction, fires, etc.), and (2) to increase the quality of earthquake science education, so as to attract talented students to geosciences. Our integrated program relies on activities tuned for different population groups who have different interests and abilities, namely young children, teenagers, young adults, and professionals.

Identification of tsunami-induced deposits using numerical modeling and rock magnetism techniques: A study case of the 1755 Lisbon tsunami in Algarve, Portugal

Storm- and tsunami-deposits are generated by similar depositional mechanisms making their discrimination hard to establish using classic sedimentologic methods. Here we propose an original approach to identify tsunami-induced deposits by combining numerical simulation and rock magnetism. To test our method, we investigate the tsunami deposit of the Boca do Rio estuary generated by the 1755 earthquake in Lisbon which is well described in the literature. We first test the 1755 tsunami scenario using a numerical inundation model to provide physical parameters for the tsunami wave. Then we use concentration (MS. SIRM) and grain size (chi(ARM), ARM, B1/2, ARM/SIRM) sensitive magnetic proxies coupled with SEM microscopy to unravel the magnetic mineralogy of the tsunami-induced deposit and its associated depositional mechanisms. In order to study the connection between the tsunami deposit and the different sedimentologic units present in the estuary, magnetic data were processed by multivariate statistical analyses. Our numerical simulation show a large inundation of the estuary with flow depths varying from 0.5 to 6 m and run up of similar to 7 m. Magnetic data show a dominance of paramagnetic minerals (quartz) mixed with lesser amount of ferromagnetic minerals, namely titanomagnetite and titanohematite both of a detrital origin and reworked from the underlying units. Multivariate statistical analyses indicate a better connection between the tsunami-induced deposit and a mixture of Units C and D. All these results point to a scenario where the energy released by the tsunami wave was strong enough to overtop and erode important amount of sand from the littoral dune and mixed it with reworked materials from underlying layers at least 1 m in depth. The method tested here represents an original and promising tool to identify tsunami-induced deposits in similar embayed beach environments.

Optimal design of piezolaminated structures using B-spline strip finite element models

A package of B-spline finite strip models is developed for the linear analysis of piezolaminated plates and shells. This package is associated to a global optimization technique in order to enhance the performance of these types of structures, subjected to various types of objective functions and/or constraints, with discrete and continuous design variables. The models considered are based on a higher-order displacement field and one can apply them to the static, free vibration and buckling analyses of laminated adaptive structures with arbitrary lay-ups, loading and boundary conditions. Genetic algorithms, with either binary or floating point encoding of design variables, were considered to find optimal locations of piezoelectric actuators as well as to determine the best voltages applied to them in order to obtain a desired structure shape. These models provide an overall economy of computing effort for static and vibration problems.

Impact of a 1755-like tsunami in Huelva, Spain

Coastal areas are highly exposed to natural hazards associated with the sea. In all cases where there is historical evidence for devastating tsunamis, as is the case of the southern coasts of the Iberian Peninsula, there is a need for quantitative hazard tsunami assessment to support spatial planning. Also, local authorities must be able to act towards the population protection in a preemptive way, to inform 'what to do' and 'where to go' and in an alarm, to make people aware of the incoming danger. With this in mind, we investigated the inundation extent, run-up and water depths, of a 1755-like event on the region of Huelva, located on the Spanish southwestern coast, one of the regions that was affected in the past by several high energy events, as proved by historical documents and sedimentological data. Modelling was made with a slightly modified version of the COMCOT (Cornell Multi-grid Coupled Tsunami Model) code. Sensitivity tests were performed for a single source in order to understand the relevance and influence of the source parameters in the inundation extent and the fundamental impact parameters. We show that a 1755-like event will have a dramatic impact in a large area close to Huelva inundating an area between 82 and 92 km(2) and reaching maximum run-up around 5 m. In this sense our results show that small variations on the characteristics of the tsunami source are not too significant for the impact assessment. We show that the maximum flow depth and the maximum run-up increase with the average slip on the source, while the strike of the fault is not a critical factor as Huelva is significantly far away from the potential sources identified up to now. We also show that the maximum flow depth within the inundated area is very dependent on the tidal level, while maximum run-up is less affected, as a consequence of the complex morphology of the area.

Projecto de execução de fundações e estruturas de uma ponte em betão armado pré-esforçado

Este Trabalho refere-se ao Projecto de Execução de Fundações e Estruturas de uma Ponte Rodoviária em betão armado pré-esforçado, realizado no âmbito do Trabalho Final de Mestrado em Engenharia Civil – Especialização em Estruturas, do Instituto Superior de Engenharia de Lisboa. O Projecto de Execução é composto de Peças Escritas e Peças Desenhadas. Nas Peças Escritas estão incluídos: Memória Justificativa e Descritiva; Cálculos Justificativos e Anexos. A ponte é composta por dois tabuleiros paralelos com 10,28m de largura cada um e afastados entre si de 0,10m. A obra é constituída de 8 tramos; os tramos correntes com 31m de comprimento e os tramos extremos com 25 e 20m de comprimento, perfazendo um comprimento total de 231m. A obra foi parcialmente isolada dos sismos pela introdução, em todos os pilares, de aparelhos de apoio de elevado amortecimento sísmico do tipo HDRB (High Damping Rubber Bearings). Encontram-se particularmente discriminadas e detalhadas neste projecto as seguintes situações: - Cálculo do Pré-esforço e respectivas perdas; - Acção das sobrecargas rodoviárias; - Diferença de comportamento da obra na entrada em serviço e no longo prazo; - Análise sísmica e do isolamento sísmico; - Estudo dos efeitos diferidos: retracção e fluência. Tendo as abordagens de cálculo e as verificações de segurança seguido a regulamentação nacional em vigor, nomeadamente RSA e REBAP, foi no entanto feita uma aproximação às regras do “Capacity Design” previstas no EC8, em que se privilegia a actuação do projectista sobre o comportamento da estrutura, procurando uma resposta não linear da mesma, visando garantir que: - A rotura não ocorrerá nos elementos de fundação; - Nos pilares a dissipação de energia se faz através de rótulas plásticas, evitando-se roturas associadas a esforços transversos. A aplicação destas regras neste Projecto demonstrou haver um agravamento substancial na definição dos esforços a que devem resistir alguns dos componentes da estrutura, designadamente os pilares e as fundações, originando soluções de secções de betão e armaduras bem mais exigentes do que aqueles que resultariam da simples verificação de segurança, pela comparação entre esforços actuante e esforços resistentes “secção a secção”, imposta pela actual regulamentação nacional.

Influência da modelação de núcleos de betão armado no cálculo automático de estruturas de edifícios

O presente trabalho aborda a questão do comportamento e influência estrutural da modelação de um núcleo de secção aberta de parede fina no modelo tridimensional de elementos finitos de um edifício de betão armado, composto por 5 pisos sujeito à acção sísmica regulamentar. Para o efeito, foram considerados os Eurocódigos estruturais, nomeadamente o Eurocódigo 0, Eurocódigo 1 e Eurocódigo 8. A fim de avaliar a importância do tipo de modelação destes elementos estruturais, elemento fundamental para a caracterização da resposta dinâmica, foi desenvolvido um estudo paramétrico, fazendo-se a análise de quatro formas distintas de modelação do mesmo núcleo. Para os mesmos casos de carregamento, é feita uma análise comparativa entre os resultados obtidos a nível de períodos, percentagem de armadura nos pilares e no próprio núcleo bem como os deslocamentos para os modelos em estudo, com recurso a um programa de cálculo automático. O edifício considerado neste trabalho foi modelado com base no programa comercial de cálculo automático SAP 2000 versão 14.4. Posteriormente apresentam-se e discutem-se os resultados obtidos nos vários modelos em causa e para as várias situações propostas ao longo do trabalho.

Caracterização mecânica de misturas betuminosas a aplicar em infra-estruturas de transportes

A maioria das infra-estruturas de transportes, nomeadamente os pavimentos rodoviários e aeroportuários, são constituídas por misturas betuminosas, o que permite um bom desempenho e uma adequada durabilidade, nas condições usuais de serviço. As misturas betuminosas são ainda amplamente utilizadas na construção de zonas de estacionamento de veículos, tendo-se verificado recentemente a sua aplicação também em infra-estruturas ferroviárias. Face à necessidade de melhorar o desempenho das vias-férreas, permitindo uma concepção mais durável de linhas de alta velocidade e uma redução dos custos da sua manutenção, tem-se vindo a desenvolver diversos estudos para promover a utilização de novos materiais, principalmente através da incorporação de misturas betuminosas. O presente trabalho tem como objectivo a caracterização do comportamento mecânico de misturas betuminosas a aplicar em infra-estruturas de transportes. Como metodologia para o estudo do comportamento mecânico das misturas betuminosas foram realizados em laboratório ensaios de cargas repetidas, nomeadamente, ensaios de flexão em quatro pontos para determinação da rigidez e da resistência à fadiga e ensaios de compressão triaxiais cíclicos para avaliação do comportamento à deformação permanente. A resistência à fadiga das misturas betuminosas em estudo foi avaliada através do ensaio de flexão em quatro pontos, com extensão controlada, e aplicação de um carregamento sinusoidal com diferentes frequências, de acordo com o procedimento de ensaio da norma europeia EN 12697-24 (2004 + A1: 2007). A resistência à deformação permanente das misturas betuminosas foi analisada através de ensaios de compressão triaxiais cíclicos, submetendo-as a uma tensão de confinamento estática pela aplicação parcial de vácuo e a uma pressão axial cíclica sob a forma rectangular, de acordo com a norma europeia EN 12697-25 (2004). O conhecimento destas propriedades mecânicas assume particular importância ao nível da formulação das misturas betuminosas, do dimensionamento de uma estrutura ou do estabelecimento de uma adequada solução para uma obra de reabilitação duma infra-estrutura de transportes. Para este estudo foi utilizado um modelo físico construído numa fossa no LNEC, com o propósito de serem testadas três substruturas ferroviárias não convencionais, utilizando sub-balastro betuminoso. A selecção das substruturas foi efectuada após uma análise de várias secções de estruturas já testadas e aplicadas noutros países, de forma a proporcionar comparações fiáveis entre elas. Os resultados obtidos mostraram que a mistura betuminosa AC20 base 50/70 (MB) aplicada na camada de sub-balastro é adequada para ser aplicada nas infra-estruturas de transportes pois apresenta um bom desempenho à fadiga e à deformação permanente. Através dos ensaios efectuados foi ainda possível entender a importante influência das características volumétricas, principalmente da porosidade para o bom comportamento da mistura betuminosa.

Design of a sea-level tsunami detection network for the Gulf of Cadiz

The devastating impact of the Sumatra tsunami of 26 December 2004, raised the question for scientists of how to forecast a tsunami threat. In 2005, the IOC-UNESCO XXIII assembly decided to implement a global tsunami warning system to cover the regions that were not yet protected, namely the Indian Ocean, the Caribbean and the North East Atlantic, the Mediterranean and connected seas (the NEAM region). Within NEAM, the Gulf of Cadiz is the more sensitive area, with an important record of devastating historical events. The objective of this paper is to present a preliminary design for a reliable tsunami detection network for the Gulf of Cadiz, based on a network of sea-level observatories. The tsunamigenic potential of this region has been revised in order to define the active tectonic structures. Tsunami hydrodynamic modeling and GIS technology have been used to identify the appropriate locations for the minimum number of sea-level stations. Results show that 3 tsunameters are required as the minimum number of stations necessary to assure an acceptable protection to the large coastal population in the Gulf of Cadiz. In addition, 29 tide gauge stations could be necessary to fully assess the effects of a tsunami along the affected coasts of Portugal, Spain and Morocco.

On the moroccan tsunami catalogue

A primary tool for regional tsunami hazard assessment is a reliable historical and instrumental catalogue of events. Morocco by its geographical situation, with two marine sides, stretching along the Atlantic coast to the west and along the Mediterranean coast to the north, is the country of Western Africa most exposed to the risk of tsunamis. Previous information on tsunami events affecting Morocco are included in the Iberian and/or the Mediterranean lists of tsunami events, as it is the case of the European GITEC Tsunami Catalogue, but there is a need to organize this information in a dataset and to assess the likelihood of claimed historical tsunamis in Morocco. Due to the fact that Moroccan sources are scarce, this compilation rely on historical documentation from neighbouring countries (Portugal and Spain) and so the compatibility between the new tsunami catalogue presented here and those that correspond to the same source areas is also discussed.

Response of a multi-domain continental margin to compression: Study from seismic reflection-refraction and numerical modelling in the Tagus Abyssal Plain

The effects of the Miocene through Present compression in the Tagus Abyssal Plain are mapped using the most up to date available to scientific community multi-channel seismic reflection and refraction data. Correlation of the rift basin fault pattern with the deep crustal structure is presented along seismic line IAM-5. Four structural domains were recognized. In the oceanic realm mild deformation concentrates in Domain I adjacent to the Tore-Madeira Rise. Domain 2 is characterized by the absence of shortening structures, except near the ocean-continent transition (OCT), implying that Miocene deformation did not propagate into the Abyssal Plain, In Domain 3 we distinguish three sub-domains: Sub-domain 3A which coincides with the OCT, Sub-domain 3B which is a highly deformed adjacent continental segment, and Sub-domain 3C. The Miocene tectonic inversion is mainly accommodated in Domain 3 by oceanwards directed thrusting at the ocean-continent transition and continentwards on the continental slope. Domain 4 corresponds to the non-rifted continental margin where only minor extensional and shortening deformation structures are observed. Finite element numerical models address the response of the various domains to the Miocene compression, emphasizing the long-wavelength differential vertical movements and the role of possible rheologic contrasts. The concentration of the Miocene deformation in the transitional zone (TC), which is the addition of Sub-domain 3A and part of 3B, is a result of two main factors: (1) focusing of compression in an already stressed region due to plate curvature and sediment loading; and (2) theological weakening. We estimate that the frictional strength in the TC is reduced in 30% relative to the surrounding regions. A model of compressive deformation propagation by means of horizontal impingement of the middle continental crust rift wedge and horizontal shearing on serpentinized mantle in the oceanic realm is presented. This model is consistent with both the geological interpretation of seismic data and the results of numerical modelling.