Local site effect microzonation of Lorca town (SE Spain)

Local site effect assessment based on subsurface ground conditions is often the key to evaluate urban seismic hazard. The site effect evaluation in Lorca town (south-eastern Spain) started with a classification of urban geology through the geological mapping at scale 1:10,000 and the use of geotechnical data and geophysical surveys. The 17 geological formations identified were classified into 5 geological/seismic formations according to their seismic amplification capacity obtained from ambient vibration measurements as well as from simultaneous strong motion records. The shear-wave velocity structure of each geological/seismic formation was evaluated by means of inversion of Rayleigh wave dispersion data obtained from vertical-component array records of ambient noise. Nakamura’s method was applied to determine a predominant period distribution map. The spectral amplification factors were fourfold the values recorded in a reference hard-rock site. Finally, the capability of this study for explaining the damage distribution caused by the May 11th, 2011 Lorca destructive earthquake (Mw \(=\) 5.2) was examined. The methods used in this work are of assistance to evaluate ground amplification phenomena in urban areas of complex geology as Lorca town due to future earthquakes with applicability on urban seismic risk management.     

Identification of potential subsidence related to pumping in the Almería basin (SE Spain)

Ground subsidence of detrital deposits in the Almería basin (SE Spain) was studied using the remote sensing technique of Differential Interferometry SAR (DInSAR). This basin is one of the most arid in Europe, receiving an average rainfall of 250 mm per year. Over the last 60 years the region has experienced an enormous agricultural and urban expansion, whose water demand has been largely supplied from groundwater, leading to the current situation of overexploitation of water resources. This paper outlines the likely relationship between groundwater abstraction and subsidence. To this end, 34 ERS and Envisat images, taken between 2003 and 2009, were analysed to estimate ground surface deformations, and hence, compared with water table variations measured in a number of piezometers in the basin. The analysis shows a clear parallelism between the variations in piezometric level and deformation of the ground surface. In addition, the zones of greatest subsidence coincide with those areas where groundwater abstractions are concentrated. Subsidence over the examined period varies from 10 to 30 mm, with extreme values as high as 50 mm, which translates to a rate of between 1·7 and 5 mm/year, reaching maximum rates of 8 mm/year at some points. Given such subsidence rates, damage to urban infrastructures are, for the moment, incipient. Copyright © 2011 John Wiley & Sons, Ltd.     

Physico-chemical evolution of groundwater in tectonically active areas. Application to the Leana hot spring (Murcia Region,SE Spain)

Seismic events can affect the physico-chemical characteristics of groundwater. These anomalies are of a pre-seismic, co-seismic and post-seismic nature and correspond to pulse variations, sudden increases and decreases without return to initial values and upward or downward changes in trend. Continuous and in situ conductivity and temperature monitoring and periodic water sampling at a hot spring associated with neotectonic activity are of great interest for establishing predictive methods. This method is limited to the seismic activity affecting the fracturing system with which the hot spring is associated. The Region of Murcia and surroundings (southeast Spain) was selected as the study area for exploring the nature of these influences on groundwater. A hot spring in the Leana spa (Murcia) was equipped and monitored during the period 2006–2008, allowing for the in situ determination of conductivity and temperature as well as of major and minor constituents at the laboratory. Due to its proximity and related with fault network, we suggest that 86 % of earthquakes located between 0 and 10 km may affect in situ parameters of groundwater, and 75 % may affect laboratory determinations. This percentage drops in more distant zones. Of all earthquakes that seem to influence groundwater, 55 % of the in situ parameter anomalies and 53 % of laboratory were of a pre-seismic nature.     

The Bullas (Murcia, SE Spain) earthquake, 29 January 2005

The intensity and focal mechanism of the Bullas earthquake of 29 January 2005 have been determined. The intensity assessment has been carried out using the EMS-98 scale and plotted on an intensity map. The focal mechanisms for the main shock and largest aftershock have been estimated from moment tensor inversion of body waves. The mechanisms obtained show left-lateral strike-slip faulting and foci at shallow depths, 3–4 km. These results are compared with those of the Bullas 2002 and Mula 1999 earthquakes in the same region.This paper has not been submitted elsewhere in identical or similar form, nor will it be during the first three months after its submission to Journal of Seismology.     

Impacts of human activities on recharge in a multilayered semiarid aquifer (Campo de Cartagena,SE Spain)

The development of intense agriculture in semiarid areas modifies intensity and spatial distribution of groundwater recharge by summing irrigation return flow to limited rainfall infiltration. Environmental tracers provide key information, but their interpretation is complicated by more complex groundwater flow patterns. In multilayered aquifers, the real origin of the groundwater samples is hard to assess because of local mixing processes occurring inside long‐screened boreholes. We use environmental tracers (¹⁴C, ¹³C, ²H, ¹⁸O, ³H) to investigate the long‐term evolution of recharge in the five‐layer Campo de Cartagena aquifer in South‐Eastern Spain, in addition to high‐resolution temperature loggings to identify the depth of origin of groundwater. Despite the complex background, this methodology allowed a reliable interpretation of the geochemistry and provided a better understanding of the groundwater flow patterns. The tritium method did not give good quantitative results because of the high variability of the recharge signal but remained an excellent indicator of recent recharge. Nonetheless, both pre‐anthropization and post‐anthropization recharge regime could be identified and quantified by radiocarbon. Before the development of agriculture, recharge varied from 17 mm.year^‐1 at the mountain ranges to 6 mm.year^‐1 in the plain, whereas the mean annual rainfall is about 300 mm. In response to the increase of agricultural activity, recharge fluxes to the plain were amplified and nowadays reach up to 210 mm.year^‐1 in irrigated areas. These values are strengthened by global water budget and local unsaturated zone studies. Copyright © 2013 John Wiley & Sons, Ltd.     

Indications of syngenetic origin of gold ore and ignimbrites near Rodalquilar (SE Spain)

Bulletin of Volcanology - The hydrothermal gold-alunite deposits of Rodalquilar are confined to a complex of very strongly altered rocks, the bulk of which is thought to be of ignimbritic origin,...     

Mathematical Modeling of Displacement Processes in Confined Flows in Coastal Aquifers

Mathematical models of fresh subsoil water flows in a confined aquifer to a sea (pit, pool, and the like) containing salt water are considered. To study these models mixed boundary value problems of the theory of analytical functions are formulated and solved with the use of the Polubarinova-Kochina method. The models were used to develop algorithms for calculating displacement in situations where subsoil water flows discharge into the sea at lateral inflow or inflow from below. The effect of the model physical parameters on the character and extent of displacement was analyzed with the use of obtained exact analytical relationships and numerical calculations. The hydrodynamic structure is described, and specific features of the flows being modeled are established.     

The Mar Menor lagoon (SE Spain): A singular natural ecosystem threatened by human activities

The Mar Menor lagoon is one of the most important ecological singularities in the Mediterranean area. At the same time, it is an area where many economic and industrial activities meet. The sum of the impacts of mining, agriculture and urban development in the surroundings to the lagoon during the last decades has affected its ecosystem. In this paper, we have reviewed the studies done by researchers over two decades regarding the impacts of human activities in the lagoon, reporting data of contamination levels and the effects generated in the ecosystem of the lagoon. The discharge of enriched nutrient waters seems to be the most important input in the lagoon. Changes to more friendly agricultural techniques in the surrounding areas are therefore necessary. After our review, we can conclude that it is necessary to have higher efforts from the public administrations and stricter environmental regulations in order to preserve, on one hand, the important ecological values of the lagoon, and on other hand, public health.     

Hydrogeochemical balance of forest umbrisol profiles (‘Sierra de Gata’, central western Spain)

Three techniques for obtaining soil water solutions (gravitational and matrical waters extracted using both in situ tension lysimeters and in vitro pressure chambers) and their later chemical analysis were performed in order to know the evolution of the soil‐solution composition when water moves down through the soil, from the Ah soil horizon to the BwC‐ or C‐horizons of forest soils located in western Spain. Additionally, ion concentrations and water volumes of input waters to soil (canopy washout) and exported waters (drainage solutions from C‐horizons) were determined to establish the net balance of solutes in order to determine the rates of leaching or retention of ions. A generalized process of sorption or retention of most components (even Cl^?) was observed, from the soil surface to the C‐horizon, in both gravitational and matrical waters, with H₄SiO₄, Mn²⁺, Na⁺, and SO₄^2? being the net exported components from the soil through the groundwater. These results enhance the role of the recycling effect in these forest soils. The net percentages of elements retained in these forest soils, considering the inputs and the outputs balance, were 68% K⁺, 85% Ca²⁺, 58% Mg²⁺, 7% Al³⁺, 5% Fe³⁺, 34% Zn²⁺, 57% Cl^?, and 20% NO₃^?, and about 75% of dissolved organic carbon was mineralized. Copyright © 2007 John Wiley & Sons, Ltd.     

Complex Isotopic and Hydrogeochemical Studies of Groundwater in the Kirov Region

Uranium-isotopic and multi-element hydrogeochemical methods are used to assess the natural pollution of fresh groundwater in horizons under development due to a modern intrusion of hypogene water with increased boron and fluorine contents in zones of old tectonic dislocations. The total index of the abovestandard (>MAC) pollution of groundwater in the region of the city of Kirov (the Kirov area) reaches 15. To reduce the effect of natural groundwater pollutants on the health of population, it is recommended that the groundwater intake should be regulated, purifiers should be used, and prospects revealed by the isotopic-and-hydrogeochemical data should be explored. The feasibility of modeling the groundwater formation and tracing the neotectonic dislocations in geological platforms is shown.     

Seismic hazard assessment of the Province of Murcia (SE Spain): analysis of source contribution to hazard

A probabilistic seismic hazard assessment of the Province of Murcia in terms of peak ground acceleration (PGA) and spectral accelerations [SA(T)] is presented in this paper. In contrast to most of the previous studies in the region, which were performed for PGA making use of intensity-to-PGA relationships, hazard is here calculated in terms of magnitude and using European spectral ground-motion models. Moreover, we have considered the most important faults in the region as specific seismic sources, and also comprehensively reviewed the earthquake catalogue. Hazard calculations are performed following the Probabilistic Seismic Hazard Assessment (PSHA) methodology using a logic tree, which accounts for three different seismic source zonings and three different ground-motion models. Hazard maps in terms of PGA and SA(0.1, 0.2, 0.5, 1.0 and 2.0 s) and coefficient of variation (COV) for the 475-year return period are shown. Subsequent analysis is focused on three sites of the province, namely, the cities of Murcia, Lorca and Cartagena, which are important industrial and tourism centres. Results at these sites have been analysed to evaluate the influence of the different input options. The most important factor affecting the results is the choice of the attenuation relationship, whereas the influence of the selected seismic source zonings appears strongly site dependant. Finally, we have performed an analysis of source contribution to hazard at each of these cities to provide preliminary guidance in devising specific risk scenarios. We have found that local source zones control the hazard for PGA and SA(T ≤ 1.0 s), although contribution from specific fault sources and long-distance north Algerian sources becomes significant from SA(0.5 s) onwards.     

Active faulting offshore SE Spain (Alboran Sea): Implications for earthquake hazard assessment in the Southern Iberian Margin

The southern margin of the Iberian Peninsula hosts the convergent boundary between the European and African Plates. The area is characterised by low to moderate magnitude shallow earthquakes, although large historical events have also occurred. In order to determine the possible sources of these events, we recently acquired swath-bathymetry, TOBI sidescan sonar and high-resolution seismic data on the Almería Margin (Eastern Alboran Sea). The new dataset reveals the offshore continuation of the NE–SW trending Carboneras Fault, a master fault in the Eastern Betic Shear Zone, and its associated structures (N150 and NS faults). These structures are active since they cut the Late Quaternary sedimentary units. The submarine Carboneras Fault zone is 100 km long, 5–10 km wide, and is divided into two N045 and N060 segments separated by an underlapping restraining stepover. Geomorphic features typically found in subaerial strike-slip faults, such as deflected drainage, water gaps, shutter ridges, pressure ridges and “en echelon” folds suggest a strike-slip motion combined with a vertical component along the submarine Carboneras Fault. Considering the NNW–SSE regional shortening axis, a left-lateral movement is deduced for the Carboneras Fault, whereas right-lateral and normal components are suggested for the associated N150 and NS faults, respectively. The offshore portion of this fault is at least twice as long as its onshore portion and together they constitute one of the longest structures in the southeastern Iberian Margin. Despite the fact that present day seismicity in the Almería margin seems to be associated with the N150 to NS faults, the Carboneras Fault is a potential source of large magnitude (M_w ∼7.2) events. Hence, the Carboneras Fault zone could pose a significant earthquake and tsunami hazard to the coasts of Spain and North Africa, and should therefore be considered in any hazard re-evaluation.     

Analysis of the weightiness of site effects on reinforced concrete (RC) building seismic behaviour: The Adra town example (SE Spain)

The damage distribution in Adra town (south‐eastern Spain) during the 1993 and 1994 Adra earthquakes (5.0 magnitude), that reached a maximum intensity degree of VII (European Macroseismic Scale (EMS scale)), was concentrated mainly in the south‐east zone of the town and the most relevant damage occurred in reinforced concrete (RC) buildings with four or five storeys. In order to evaluate the influence of ground condition on RC building behaviour, geological, geomorphological and geophysical surveys were carried out, and a detailed map of ground surface structure was obtained. Short‐period microtremor observations were performed in 160 sites on a 100m × 100m dimension grid and Nakamura's method was applied in order to determine a distribution map of soil predominant periods. Shorter predominant periods (0.1–0.3 s) were found in mountainous and neighbouring zones and larger periods (greater than 0.5 s) in thicker Holocene alluvial fans. A relationship T = (0.049 ± 0.001)N, where T is the natural period of swaying motion and N is the number of storeys, has been empirically obtained by using microtremor measurements at the top of 38 RC buildings (ranging from 2 to 9 storeys). 1‐D simulation of strong motion on different soil conditions and for several typical RC buildings were computed, using the acceleration record in Adra town of the 1993 earthquake. It is noteworthy that all the aforementioned results show the influence of site effects in the degree and distribution of observed building damage. Copyright © 2007 John Wiley & Sons, Ltd.     

Seismic properties of lower crustal xenoliths from El Hoyazo (SE Spain): Experimental evidence up to partial melting

Seismic techniques provide unique tools to investigate the structure and, in combination with petrological, geochemical and petrophysical study, the composition of the lower crust. Controversies can be solved with comparative study of metamorphic terrains or xenoliths that occur adjacent to areas where seismic refraction/reflection data are available. Xenoliths represent a direct sampling of the inaccessible lower crust at the time of the volcanism, whilst exposed crustal sections can only be used as analogue of present day lower crust.The present study is focused on the measurements of compressional wave velocities up to conditions exceeding the beginning of melting (950 °C at 500 MPa confining pressure) on three garnet–biotite–sillimanite metapelitic xenoliths recovered from the Neogene dacites of El Hoyazo (SE Spain). They preserve widespread interstitial rhyolitic glass as evidence of primary melt extraction and represent the best example of partially molten lower crust in the Alborán Domain. The influence of glass on Vp is primarily reflected by anomalous positive dVp/dT while heating with velocity increasing at 500 MPa from 4.98 to 5.50 km s^− 1 at room temperature to 5.85–6.79 km s^− 1 at 650–700 °C. This corresponds to the glass transition where all the grain boundaries and most of the pores within the glass are closed. After this point, the velocity decreases to 6.2–6.5 km s^− 1 at 950 °C where re-melting of the glass is achieved and additional partial melt produced. On cooling, the behavior is normal with negative dVp/dT. After the thermal treatment velocities are 30% higher (6.07–7.21 km s^− 1) and reveal that in the presence of intergranular melt velocity measurements at room temperature cannot be extrapolated to high temperatures.P-waves measured at melting conditions are in agreement with deep seismic refraction data and tomography in the area and corroborate the hypothesis that partial melts are actually present in Alborán lower crust.     

Paleoseismic indications along ‘aseismic’ fault segments in the guadalentín depression (SE Spain)

The Guadalentín Depression, located in SE Spain (Murcia Region), is bounded by two of the main NE-SW master faults of the Eastern Betics Cordilleras: The Lorca-Alhama and the Palomares left-lateral strike-slip faults. Available earthquake data indicate that, in the last 600 years, some sectors of the Lorca-Alhama Fault and the entire sector of the Palomares Fault have not been associated with significant historical seismicity. However, they show a wide range of diagnostic features of earthquake surface displacements on late Pleistocene and Holocene alluvial and colluvial surfaces. Aside from the left-lateral offsets recorded along 045–050 ° master fault strands of the Lorca-Alhama Fault, major paleoseismic surface displacements show different kinematics in relation to the broad orientation of the fault strands: (1) vertical normal displacements along 010–020 ° trending faults mainly preserved as degraded fault scarps of 2.5-1.8 m high (Aljibejo site); and (2) vertical reverse displacements, with average offsets of 0.2 – 1 m, along 065–080 ° subsidiary faults. In this last group, the younger one (Carraclaca Baths site) remains as a fault scarp of 0.8 m height affecting a cascade tufa which was active until the Spanish Roman Period (2nd Century B.C. to 6th Century A.D.). In other cases, reverse offsets resulted in smaller displacements (0.26 m) of paleosols, but show a recurrent behaviour (La Escarihuela site). The strongest earthquakes recorded in the study area did not exceed more than Mb 4.5 or MSK Intensity VIII (historical) with no evidence of coseismic rupture. Therefore, the preliminary data presented here seem to indicate that the paleoseismic activity on both faults is capable of producing coseismic surface displacements, probably reaching magnitudes of at least 6.5. These data show that paleoseismic studies based on geomorphological analyses are a useful tool in the assessment of the relative degree of activity of apparently ‘aseismic’ fault traces.