Petrophysical modeling of high seismic velocity crust at the Namibian volcanic margin

Author(s):  
R.B. Trumbull ◽  
S.V. Sobolev ◽  
K. Bauer
2020 ◽  
Vol 25 (3) ◽  
pp. 415-423
Author(s):  
Ahmed Lachhab ◽  
El Mehdi Benyassine ◽  
Mohamed Rouai ◽  
Abdelilah Dekayir ◽  
Jean C. Parisot ◽  
...  

The tailings of Zeida's abandoned mine are found near the city of Midelt, in the middle of the high Moulouya watershed between the Middle and the High Atlas of Morocco. The tailings occupy an area of about 100 ha and are stored either in large mining pit lakes with clay-marl substratum or directly on a heavily fractured granite bedrock. The high contents of lead and arsenic in these tailings have transformed them into sources of pollution that disperse by wind, runoff, and seepage to the aquifer through faults and fractures. In this work, the main goal is to identify the pathways of contaminated water with heavy metals and arsenic to the local aquifers, water ponds, and Moulouya River. For this reason, geophysical surveys including electrical resistivity tomography (ERT), seismic refraction tomography (SRT) and very low-frequency electromagnetic (VLF-EM) methods were carried out over the tailings, and directly on the substratum outside the tailings. The result obtained from combining these methods has shown that pollutants were funneled through fractures, faults, and subsurface paleochannels and contaminated the hydrological system connecting groundwater, ponds, and the river. The ERT profiles have successfully shown the location of fractures, some of which extend throughout the upper formation to depths reaching the granite. The ERT was not successful in identifying fractures directly beneath the tailings due to their low resistivity which inhibits electrical current from propagating deeper. The seismic refraction surveys have provided valuable details on the local geology, and clearly identified the thickness of the tailings and explicitly marked the boundary between the Triassic formation and the granite. It also aided in the identification of paleochannels. The tailings materials were easily identified by both their low resistivity and low P-wave velocity values. Also, both resistivity and seismic velocity values rapidly increased beneath the tailings due to the compaction of the material and lack of moisture and have proven to be effective in identifying the upper limit of the granite. Faults were found to lie along the bottom of paleochannels, which suggest that the locations of these channels were caused by these same faults. The VLF-EM surveys have shown tilt angle anomalies over fractured areas which were also evinced by low resistivity area in ERT profiles. Finally, this study showed that the three geophysical methods were complementary and in good agreement in revealing the pathways of contamination from the tailings to the local aquifer, nearby ponds and Moulouya River.


2020 ◽  
Vol 72 (1) ◽  
Author(s):  
Yasuhira Aoyagi ◽  
Haruo Kimura ◽  
Kazuo Mizoguchi

Abstract The earthquake rupture termination mechanism and size of the ruptured area are crucial parameters for earthquake magnitude estimations and seismic hazard assessments. The 2016 Mw 7.0 Kumamoto Earthquake, central Kyushu, Japan, ruptured a 34-km-long area along previously recognized active faults, eastern part of the Futagawa fault zone and northernmost part of the Hinagu fault zone. Many researchers have suggested that a magma chamber under Aso Volcano terminated the eastward rupture. However, the termination mechanism of the southward rupture has remained unclear. Here, we conduct a local seismic tomographic inversion using a dense temporary seismic network to detail the seismic velocity structure around the southern termination of the rupture. The compressional-wave velocity (Vp) results and compressional- to shear-wave velocity (Vp/Vs) structure indicate several E–W- and ENE–WSW-trending zonal anomalies in the upper to middle crust. These zonal anomalies may reflect regional geological structures that follow the same trends as the Oita–Kumamoto Tectonic Line and Usuki–Yatsushiro Tectonic Line. While the 2016 Kumamoto Earthquake rupture mainly propagated through a low-Vp/Vs area (1.62–1.74) along the Hinagu fault zone, the southern termination of the earthquake at the focal depth of the mainshock is adjacent to a 3-km-diameter high-Vp/Vs body. There is a rapid 5-km step in the depth of the seismogenic layer across the E–W-trending velocity boundary between the low- and high-Vp/Vs areas that corresponds well with the Rokkoku Tectonic Line; this geological boundary is the likely cause of the dislocation of the seismogenic layer because it is intruded by serpentinite veins. A possible factor in the southern rupture termination of the 2016 Kumamoto Earthquake is the existence of a high-Vp/Vs body in the direction of southern rupture propagation. The provided details of this inhomogeneous barrier, which are inferred from the seismic velocity structures, may improve future seismic hazard assessments for a complex fault system composed of multiple segments.


2021 ◽  
Vol 13 (14) ◽  
pp. 2684
Author(s):  
Eldert Fokker ◽  
Elmer Ruigrok ◽  
Rhys Hawkins ◽  
Jeannot Trampert

Previous studies examining the relationship between the groundwater table and seismic velocities have been guided by empirical relationships only. Here, we develop a physics-based model relating fluctuations in groundwater table and pore pressure with seismic velocity variations through changes in effective stress. This model justifies the use of seismic velocity variations for monitoring of the pore pressure. Using a subset of the Groningen seismic network, near-surface velocity changes are estimated over a four-year period, using passive image interferometry. The same velocity changes are predicted by applying the newly derived theory to pressure-head recordings. It is demonstrated that the theory provides a close match of the observed seismic velocity changes.


Author(s):  
Fumiaki Nagashima ◽  
Hiroshi Kawase

Summary P-wave velocity (Vp) is an important parameter for constructing seismic velocity models of the subsurface structures by using microtremors and earthquake ground motions or any other geophysical exploration data. In order to reflect the ground survey information in Japan to the Vp structure, we investigated the relationships among Vs, Vp, and depth by using PS-logging data at all K-NET and KiK-net sites. Vp values are concentrated at around 500 m/s and 1,500 m/s when Vs is lower than 1,000 m/s, where these concentrated areas show two distinctive characteristics of unsaturated and saturated soil, respectively. Many Vp values in the layer shallower than 4 m are around 500 m/s, which suggests the dominance of unsaturated soil, while many Vp values in the layer deeper than 4 m are larger than 1,500 m/s, which suggests the dominance of saturated soil there. We also investigated those relationships for different soil types at K-NET sites. Although each soil type has its own depth range, all soil types show similar relationships among Vs, Vp, and depth. Then, considering the depth profile of Vp, we divided the dataset into two by the depth, which is shallower or deeper than 4 m, and calculated the geometrical mean of Vp and the geometrical standard deviation in every Vs bins of 200 m/s. Finally, we obtained the regression curves for the average and standard deviation of Vp estimated from Vs to get the Vp conversion functions from Vs, which can be applied to a wide Vs range. We also obtained the regression curves for two datasets with Vp lower and higher than 1,200 m/s. These regression curves can be applied when the groundwater level is known. In addition, we obtained the regression curves for density from Vs or Vp. An example of the application for those relationships in the velocity inversion is shown.


2019 ◽  
Vol 71 (1) ◽  
Author(s):  
Saeko Kita

AbstractI relocated the hypocenters of the 2018 M6.7 Hokkaido Eastern Iburi earthquake and its surrounding area, using a three-dimensional seismic structure, the double-difference relocation method, and the JMA earthquake catalog. After relocation, the focal depth of the mainshock became 35.4 km. As previous studies show, in south-central Hokkaido, the Hidaka collision zone is formed, and anomalous deep and thickened forearc crust material is subducting at depths of less than 70 km. The mainshock and its aftershocks are located at depths of approximately 10 to 40 km within the lower crust of the anomalous deep and thickened curst near the uppermost mantle material intrusions in the northwestern edge of this Hidaka collision zone. Like the two previous large events, the aftershocks of this event incline steeply eastward and appear to be distributed in the deeper extension of the Ishikari-teichi-toen fault zone. The highly inclined fault in the present study is consistent with a fault model by a geodetic analysis with InSAR. The aftershocks at depths of 10 to 20 km are located at the western edge of the high-attenuation (low-Qp) zone. These kinds of relationships between hypocenters and materials are the same as the 1970 and 1982 events in the Hidaka collision zone. The anomalous large focal depths of these large events compared with the average depth limit of inland earthquakes in Japan could be caused by the locally lower temperature in south-central Hokkaido. This event is one of the approximately M7 large inland earthquakes that occurred repeatedly at a recurrence interval of approximately 40 years and is important in the collision process in the Hidaka collision zone.


2019 ◽  
Vol 71 (1) ◽  
Author(s):  
Tohru Watanabe ◽  
Miho Makimura ◽  
Yohei Kaiwa ◽  
Guillaume Desbois ◽  
Kenta Yoshida ◽  
...  

AbstractElastic wave velocity and electrical conductivity in a brine-saturated granitic rock were measured under confining pressures of up to 150 MPa and microstructure of pores was examined with SEM on ion-milled surfaces to understand the pores that govern electrical conduction at high pressures. The closure of cracks under pressure causes the increase in velocity and decrease in conductivity. Conductivity decreases steeply below 10 MPa and then gradually at higher pressures. Though cracks are mostly closed at the confining pressure of 150 MPa, brine must be still interconnected to show observed conductivity. SEM observation shows that some cracks have remarkable variation in aperture. The aperture varies from ~ 100 nm to ~ 3 μm along a crack. FIB–SEM observation suggests that wide aperture parts are interconnected in a crack. Both wide and narrow aperture parts work parallel as conduction paths at low pressures. At high pressures, narrow aperture parts are closed but wide aperture parts are still open to maintain conduction paths. The closure of narrow aperture parts leads to a steep decrease in conductivity, since narrow aperture parts dominate cracks. There should be cracks in various sizes in the crust: from grain boundaries to large faults. A crack must have a variation in aperture, and wide aperture parts must govern the conduction paths at depths. A simple tube model was employed to estimate the fluid volume fraction. The fluid volume fraction of 10−4–10−3 is estimated for the conductivity of 10−2 S/m. Conduction paths composed of wide aperture parts are consistent with observed moderate fluctuations (< 10%) in seismic velocity in the crust.


Sign in / Sign up

Export Citation Format

Share Document