Surface Soil Effects Study Using Short-period Microtremor Observations in Almería City, Southern Spain

—?In Almería city large earthquakes occurred and many buildings were completely destroyed in these historical earthquakes. The actual population of Almería city is about 200,000 people. This population is rapidly increasing and new urbanizing areas are growing to the eastern part of the city where they are located in softer soil conditions. Consequently, the evaluation of surface soil conditions is very important from a standpoint of earthquake disaster mitigation. We have obtained a landform classification map developed by analysing aerial photos, large-scale topographic maps and 80 borehole data. Eleven unit areas, which have different soil conditions, were inferred from this research. Also, S-wave velocity prospecting tests were carried out at several sites within the city. The shear-velocity values of the ground vary from 1689?m/s in hard rock to 298?m/s in soft soil. These results are useful for understanding the uppermost soil characteristics and are used for soil classification. Finally, short-period microtremor observations were densely carried out in the research area and NAKAMURA's method (1989) was applied for determining predominant periods. Microtremors were observed at about 173 sites with mainly 400?m interval in rock sites and 200?m interval in relatively soft soil sites. From the result of these microtremor measurements, the predominant period determined at rock site, in the western part of the city and historic area, is very short, about 0.1?s, and very stable. However at soft soil sites, in the center of the city, near Zapillo Beach and in the newly developed urban area, the predominant period is about 1.0?s and even larger in concordance with the geological conditions. Finally, at medium soil sites, in the eastern part of the city, the predominant period is about 0.4?s and it appears very stable in the whole region. The difference of predominant periods between hard rock and soft soil sites is very clear and it has been observed that the distribution of predominant periods depends heavily on the surface soil conditions.     

Application of conditional geostatistical simulation to calculate the probability of occurrence of earthquakes belonging to a seismic series

Geostatistics offers various techniques of estimation and simulation that have been satisfactorily applied in solving geological problems. In this sense, conditional geostatistical simulation is applied to calculate the probability of occurrence of an earthquake with a lower than or equal magnitude to one determined during a seismic series. It is possible to calculate the energy of the next most probable earthquake from a specific time, given knowledge of the structure existing among earthquakes occurring prior to a specific moment.     

Simulating a seismic series using geostatistical and stochastic methods: application to the seismic series in the Alborán Sea (1997 June 24–?)

Seismic series can be taken as examples of correlated unstationary sets of time-stochastic sequences. We investigate the possibility of estimating what is most probable to occur subsequently, if we know the events that have occurred up to a given moment.
The stochastic methods can be used with data of the seismic series, irrespective of their genesis and origin. Using three stochastic methods, namely (1) simulating the likelihood of occurrence by conditional geostatistical simulation; (2) developing a stochastic analysis of the energy release by means of energy packages; and (3) calculating the occurrence time of the most probable next earthquake, we were able to simulate the occurrence of earthquakes that took place during the Alborán Sea seismic series (1997–1998).
We conclude that it is possible to set limits on the time of occurrence and energy release, understood as the magnitude of the most probable earthquake, during the development of a seismic sequence and prior to the actual occurrence of the earthquake.     

Simulation of the Seismic Response of Sedimentary Basins with Vertical Constant-Gradient Velocity for Incident <Emphasis Type="Italic">SH</Emphasis> Waves

The simulation of the seismic response of heterogeneous sedimentary basins under incident plane waves is computed using the Indirect Boundary Element Method (IBEM). To deal with these kinds of basins we used approximate analytical expressions for the two-dimensional Greens functions of a medium with constant-gradient wave propagation velocity. On the other hand, for the homogeneous half space underlying the sedimentary basin, the full space Greens functions were used. The response of semi-circular heterogeneous basins under incident SH waves is explored by means of the displacements in the frequency-space diagrams and synthetic seismograms. Moreover, we compared these results with those obtained for other homogeneous semi-circular models. The principal differences among them are pointed out. This simulation provided interesting results that displayed a complex amplification pattern in a rich spectrum of frequencies and locations. The maximum amplitudes levels were found around the edges of the heterogeneous sedimentary basins. In time domain some features characterize the seismic response of the basin which include enhanced generation and trapping of surface waves inside the sediments, and the reduced emission of seismic energy to the hard rock. In the heterogeneous models the lateral reflections of surface waves greatly influence the total displacements at the free surface in comparison with the homogeneous models where the displacements have a shared influence among both vertical and lateral reflections.     

Co-seismic landslide detection after M 7.4 earthquake on June 23, 2020, in Oaxaca,Mexico, based on rapid mapping method using high and medium resolution synthetic aperture radar (SAR) images

Landslides are the fourth most common natural disasters in the world, with Costa Rica and southern Mexico being the most affected regions of Central America (Froude and Petley, 2018). In this work, we propose a semi-automated method to detect earthquake-triggered landslides for rapid mapping after a disaster event using open Sentinel-1 data. We used high-resolution TerraSAR-X data and very high-resolution Spot-7 images to compare and evaluate the accuracy of landslide distribution maps generated from the semi-automated method, applied to the M 7.1 earthquake on June 23, 2017, in Oaxaca, Mexico. The outcomes showed better accuracy in descending orbits due to ‘windward-leeward’ physiographic conditions, with a 50.56% quality percentage. This shows a reasonably good capacity to detect co-seismic landslides. However, the breaching factor was also high because several features, such as bare soils and agricultural areas, were incorrectly identified as co-seismic landslides. Finally, this semi-automated method establishes a basis for future improvements in methodologies applied to construct rapid mapping inventories using medium SAR scales.

     

A Moore's cellular automaton model to get probabilistic seismic hazard maps for different magnitude releases: A case study for Greece

Cellular automata are simple mathematical idealizations of natural systems and they supply useful models for many investigations in natural science. Examples include sandpile models, forest fire models, and slider block models used in seismology. In the present paper, they have been used for establishing temporal relations between the energy releases of the seismic events that occurred in neighboring parts of the crust. The catalogue is divided into time intervals, and the region is divided into cells which are declared active or inactive by means of a threshold energy release criterion. Thus, a pattern of active and inactive cells which evolves over time is determined. A stochastic cellular automaton is constructed starting with these patterns, in order to simulate their spatio-temporal evolution, by supposing a Moore's neighborhood interaction between the cells. The best model is chosen by maximizing the mutual information between the past and the future states. Finally, a Probabilistic Seismic Hazard Map is given for the different energy releases considered. The method has been applied to the Greece catalogue from 1900 to 1999. The Probabilistic Seismic Hazard Maps for energies corresponding to m = 4 and m = 5 are close to the real seismicity after the data in that area, and they correspond to a background seismicity in the whole area. This background seismicity seems to cover the whole area in periods of around 25–50 years. The optimum cell size is in agreement with other studies; for m > 6 the optimum area increases according to the threshold of clear spatial resolution, and the active cells are not so clustered. The results are coherent with other hazard studies in the zone and with the seismicity recorded after the data set, as well as provide an interaction model which points out the large scale nature of the earthquake occurrence.     

Diffraction of P, S and Rayleigh waves by three-dimensional topographies

The diffraction of P, S and Rayleigh waves by 3-D topographies in an elastic half-space is studied using a simplified indirect boundary element method (IBEM). This technique is based on the integral representation of the diffracted elastic fields in terms of single-layer boundary sources. It can be seen as a numerical realization of Huygens principle because diffracted waves are constructed at the boundaries from where they are radiated by means of boundary sources. A Fredholm integral equation of the second kind for such sources is obtained from the stress-free boundary conditions. A simplified discretization scheme for the numerical and analytical integration of the exact Green's functions, which employs circles of various sizes to cover most of the boundary surface, is used.
The incidence of elastic waves on 3-D topographical profiles is studied. We analyse the displacement amplitudes in the frequency, space and time domains. The results show that the vertical walls of a cylindrical cavity are strong diffractors producing emission of energy in all directions. In the case of a mountain and incident P, SV and SH waves the results show a great variability of the surface ground motion. These spatial variations are due to the interference between locally generated diffracted waves. A polarization analysis of the surface displacement at different locations shows that the diffracted waves are mostly surface and creeping waves.     

Spatial and temporal analysis of a seismic series using a new version of the three point method: application to the 1989 Antequera (Spain) earthquakes

The three point method (TPM) has been successfully applied to several seismic series and has provided information about the spatial characteristics (azimuth and dip) of the fault planes activated in the rupture process. A new development of the TPM to determine temporal characteristics, is presented, to obtain the evolution of the fracturing process of an active fault system. For the analysis of the 158 microearthquakes and earthquakes that took place in the seismic series of Antequera in June 1989, the choice of a threshold magnitude (m_u = 2.5) has permitted the events related to the most relevant fractures to be distinguished. Secondly, only the events between two concentric spheres (here named Spatial Crown) with respect to a given earthquake, have been used in order to avoid taking into account earthquakes that are too close to each other, together with the very distant events that have little relation to the event analysed. The Spatial Crown has permitted some clear results in the Antequera series, where we have found that the fracturing process began fundamentally with N 80° E planes and evolved to N 65° W planes. Finally, an error analysis enables an estimation of the uncertainty in the results to be obtained from the errors in the data.