Summary Three- and four-dimensional geological models provide important means for mineral exploration as they help to unravel the tectonic and depositional history of an area. They also serve to outline alteration zone patterns, establish ore-bearing fluid migration paths and identify potential areas for the deposition of minerals in economic concentrations. Geological models serve as constraints within the geophysical modelling process and contribute in the generation of new targets for exploration. In this work, we present a new 3D geological model of the Neves-Corvo mining area in the Iberian Pyrite Belt, Castro Verde, Portugal. The model provides higher resolution insight into a subsample of an area previously modelled using an updated (4D) geoscience data catalogue. The 3D model has been built using extensive ground-loop electromagnetic and seismic reflection data acquired by Lundin Mining, more than 1,000 surface drill-holes, 7,000 underground drill-holes, as well as updated geological cross-sections and surface geology data. Gravity and magnetic data were also used qualitatively to assess the model. The result is an explicit 3D model with a considerable level of detail honouring all available data which may be used for further mineral exploration.

The South Portuguese Terrane encompasses three major pre-orogenic mega-sequences: the Phyllite-Quartzite Group and the Tercenas Fm., the Volcanic-sedimentary Complex, and the Carrapateira Group, all three overlain by the syn-orogenic Baixo Alentejo Flysch Group, dealt with in Chap. 11. During the Mid-Late Devonian all the three mega-sequences belonged to an epicontinental sea. From the late Devonian to late Visean major changes took place: the region that became the IPB Domain was the locus of important lithospheric extension that led to the emplacement of important volumes of igneous rocks associated to a Volcanic-Sedimentary Complex, to which the VHMS ore deposits are associated; in Southwest Portugal the Tercenas Fm. shallow water sediments changed to a mixed-siliciclastic-carbonate shelf succession (Carrapateira Group); during the late Visean the extensional regime changed to crustal compression that caused the onset of huge deposits of flysch sediments (Baixo Alentejo Flysch Group) in successive depocenters migrating southward until the late Moscovian.

The knowledge of the anisotropic properties beneath the Iberian Peninsula and Northern Morocco has been dramatically improved since late 2007 with the analysis of the data provided by the dense TopoIberia broad-band seismic network, the increasing number of permanent stations operating in Morocco, Portugal and Spain, and the contribution of smaller scale/higher resolution experiments. Results from the two first TopoIberia deployments have evidenced a spectacular rotation of the fast polarization direction (FPD) along the Gibraltar Arc, interpreted as an evidence of mantle flow deflected around the high velocity slab beneath the Alboran Sea, and a rather uniform N100°E FPD beneath the central Iberian Variscan Massif, consistent with global mantle flow models taking into account contributions of surface plate motion, density variations and net lithosphere rotation. The results from the last Iberarray deployment presented here, covering the northern part of the Iberian Peninsula, also show a rather uniform FPD orientation close to N100°E, thus confirming the previous interpretation globally relating the anisotropic parameters to the LPO of mantle minerals generated by mantle flow at asthenospheric depths. However, the degree of anisotropy varies significantly, from delay time values of around 0.5s beneath NW Iberia to values reaching 2.0s in its NE corner. The anisotropic parameters retrieved from single events providing high quality data also show significant differences for stations located in the Variscan units of NW Iberia, suggesting that the region includes multiple anisotropic layers or complex anisotropy systems. These results allow to complete the map of the anisotropic properties of the westernmost Mediterranean region, which can now be considered as one of best constrained regions worldwide, with more than 300 sites investigated over an area extending from the Bay of Biscay to the Sahara platform.

Mainland Portugal, on the southwestern edge of the European continent, is located directly north of the boundary between the Eurasian and Nubian plates. It lies in a region of slow lithospheric deformation (< 5 mm yr(-1)), which has generated some of the largest earthquakes in Europe, both intraplate (mainland) and interplate (offshore). Some offshore earthquakes are nucleated on old and cold lithospheric mantle, at depths down to 60 km. The seismicity of mainland Portugal and its adjacent offshore has been repeatedly classified as diffuse. In this paper, we analyse the instrumental earthquake catalogue for western Iberia, which covers the period between 1961 and 2013. Between 2010 and 2012, the catalogue was enriched with data from dense broad-band deployments. We show that although the plate boundary south of Portugal is diffuse, in that deformation is accommodated along several distributed faults rather than along one long linear plate boundary, the seismicity itself is not diffuse. Rather, when located using high-quality data, earthquakes collapse into well-defined clusters and lineations. We identify and characterize the most outstanding clusters and lineations of epicentres and correlate them with geophysical and tectonic features (historical seismicity, topography, geologically mapped faults, Moho depth, free-air gravity, magnetic anomalies and geotectonic units). Both onshore and offshore, clusters and lineations of earthquakes are aligned preferentially NNE-SSW and WNW-ESE. Cumulative seismic moment and epicentre density decrease from south to north, with increasing distance from the plate boundary. Only few earthquake lineations coincide with geologically mapped faults. Clusters and lineations that do not match geologically mapped faults may correspond to previously unmapped faults (e.g. blind faults), rheological boundaries or distributed fracturing inside blocks that are more brittle and therefore break more easily than neighbour blocks. The seismicity map of western Iberia presented in this article opens important questions concerning the regional seismotectonics. This work shows that the study of low-magnitude earthquakes using dense seismic deployments is a powerful tool to study lithospheric deformation in slowly deforming regions, such as western Iberia, where high-magnitude earthquakes occur with long recurrence intervals.

P and S receiver functions (PRF and SRF) from 19 seismograph stations in the Gibraltar Arc and the Iberian Massif reveal new details of the regional deep structure. Within the high-velocity mantle body below southern Spain the 660-km discontinuity is depressed by at least 20 km. The Ps phase from the 410-km discontinuity is missing at most stations in the Gibraltar Arc. A thin (similar to 50 km) low-S-velocity layer atop the 410-km discontinuity is found under the Atlantic margin. At most stations the S410p phase in the SRFs arrives 1.0-2.5 s earlier than predicted by IASP91 model, but, for the propagation paths through the upper mantle below southern Spain, the arrivals of S410p are delayed by up to +1.5 s. The early arrivals can be explained by elevated Vp/Vs ratio in the upper mantle or by a depressed 410-km discontinuity. The positive residuals are indicative of a low (similar to 1.7 versus similar to 1.8 in IASP91) Vp/Vs ratio. Previously, the low ratio was found in depleted lithosphere of Precambrian cratons. From simultaneous inversion of the PRFs and SRFs we recognize two types of the mantle: 'continental' and 'oceanic'. In the 'continental' upper mantle the S-wave velocity in the high-velocity lid is 4.4-4.5 km s(-1), the S-velocity contrast between the lid and the underlying mantle is often near the limit of resolution (0.1 km s(-1)), and the bottom of the lid is at a depth reaching 90 100 km. In the 'oceanic' domain, the S-wave velocities in the lid and the underlying mantle are typically 4.2-4.3 and similar to 4.0 km s(-1), respectively. The bottom of the lid is at a shallow depth (around 50 km), and at some locations the lid is replaced by a low S-wave velocity layer. The narrow S-N-oriented band of earthquakes at depths from 70 to 120 km in the Alboran Sea is in the 'continental' domain, near the boundary between the 'continental' and 'oceanic' domains, and the intermediate seismicity may be an effect of ongoing destruction of the continental lithosphere.

Tese de doutoramento, Ciencias Geofisicas e da Geoinformacao (Geofisica), Universidade de Lisboa, Faculdade de Ciencias, 2012

Seismic recordings of IRIS/IDA/GSN station CMLA and of several temporary stations in the Azores archipelago are processed with P and S receiver function (PRF and SRF) techniques. Contrary to regional seismic tomography these methods provide estimates of the absolute velocities and of the Vp/Vs ratio up to a depth of ~300km. Joint inversion of PRFs and SRFs for a few data sets consistently reveals a division of the subsurface medium into four zones with a distinctly different Vp/Vs ratio: the crust ~20km thick with a ratio of ~1.9 in the lower crust, the high-Vs mantle lid with a strongly reduced Vp/Vs velocity ratio relative to the standard 1.8, the low-velocity zone (LVZ) with a velocity ratio of ~2.0, and the underlying upper-mantle layer with a standard velocity ratio. Our estimates of crustal thickness greatly exceed previous estimates (~10km). The base of the high-Vs lid (the Gutenberg discontinuity) is at a depth of ~80km. The LVZ with a reduction of S velocity of ~15% relative to the standard (IASP91) model is terminated at a depth of ~200km. The average thickness of the mantle transition zone (TZ) is evaluated from the time difference between the S410p and SKS660p, seismic phases that are robustly detected in the S and SKS receiver functions. This thickness is practically similar to the standard IASP91 value of 250km, and is characteristic of a large region of the North Atlantic outside the Azores plateau. Our data are indicative of a reduction of the S-wave velocity of several percent relative to the standard velocity in a depth interval from 460 to 500km. This reduction is found in the nearest vicinities of the Azores, in the region sampled by the PRFs, but, as evidenced by SRFs, it is missing at a distance of a few hundred kilometers from the islands. We speculate that this anomaly may correspond to the source of a plume which generated the Azores hotspot. Previously, a low S velocity in this depth range was found with SRF techniques beneath a few other hotspots.

The Faial earthquake (M L 5.8) that occurred on the 9th of July, 1998, in the Azores region (north Atlantic), caused nine casualties and severe destruction affecting more than 5,000 people. The main shock was located at sea, 10 km NE of the Faial Island, and triggered a seismic sequence that lasted for several weeks and was characterized by an unusual high p-value of 1.40 for the modified Omori law. We present here the results of a joint inversion of hypocenters and 1D velocity model performed on the data collected by the permanent network complemented with a temporary network installed shortly after the occurrence of the main event. The 1D velocity model shows a heterogeneous upper crust, testified by the observed differences in site effects at the stations, while the middle crust from ∼2.5 to 8 km in depth is quite homogeneous. The Moho is located at a depth of about 12–13 km and the Vp/Vs ratio is found to be around 1.78. The events at depth are mainly concentrated in the middle-lower crust (8–12 km), while their spatial distribution shows a main cluster, visible after relocation, SSE trending. This direction of elongation is consistent with one of the fault planes (N151°E) of the centroid moment tensor (CMT) solution for the main shock. The same plane is the preferred main shock fault plane inferred after a Coulomb failure function analysis on the aftershock distribution. The main event relocation points to a focal depth shallower than 5 km. The aftershocks pattern shows that several fault systems were reactivated by the stress perturbation induced by the main shock. Besides the two main tectonic directions, trending WNW–ESE and NNW–SSE, observed in the tectonics of Faial, Pico, and S. Jorge, there is also evidence of a new tectonic direction trending WSW–ENE.

—The project GEOALGAR, initiated in May 2000, is devoted to the geodynamic monitoring and seismic characterization of the Algarve region. A brief description of the project goals, as well as the first results concerning the analysis of the recent seismic digital data, from 1999 and 2000, are presented. After simultaneous inversion of the seismic data and the velocity model parameters, the relocation of the hypocenters was performed for two selected areas. Twenty-five earthquakes were used for Area 1 and 125 earthquakes for Area 2, selected from the period 1999–2000 and with magnitude ML ≥ 2.0. The results show that there are two main regions where there are more events: the Monchique region (inland) and another one in the area of Guadalquivir Bank (ranging from 36.4°N, 08°W to 36.8°N, 7.2°W); for Area 1, the hypocentral corrections are relatively small, with the focus slightly deeper than that in the old solution; for the more regional events, the hypocenters corrections are bigger, with the focus becoming more shallow. Fault-plane solutions for the recent events were also estimated, showing that best solutions are dominated by strike-slip movement consistent with a stress controlled by a horizontal compression in the NW-SE to NNW-SSE direction, with two exceptions showing reverse mechanisms, with maximum horizontal stress orientation slightly rotated to a N-S direction. These results are also in agreement with those presented in previous studies performed by different authors. The new epicentral locations show a more organized spatial distribution that could indicate a possible correlation with some known tectonic features. However the fault-plane solutions are considerably more difficult to correlate with the neotectonic features.