scholarly journals Determining seismic shear-velocity from ambient noise sources at regional and local scales

2021 ◽  
Author(s):  
◽  
Francesco Civilini
Keyword(s):  
Surface Wave ◽  
Ambient Noise ◽  
Shear Velocity ◽  
Volcanic Field ◽  
Geothermal Field ◽  
Near Surface ◽  
Velocity Models ◽  
Geothermal Fields ◽  
Cross Correlations ◽  
Velocity Changes

<p>We present three projects that use different bandwidths of the ambient noise spectrum to solve geophysical problems. Specifically, we use signals within the noise field to determine surface and shear wave velocities, image the shallow and deep crust, and monitor time-dependent deformation resulting from geothermal fluid injection and extraction.  Harrat Al-Madinah, a Cenozoic bimodal alkaline volcanic field in west-central Saudi Arabia, is imaged using shear-velocities obtained from natural ambient seismic noise. To our knowledge, this project is the first analysis of Saudi Arabia structure using ambient noise methods. Surface wave arrivals are extracted from a year's worth of station-pair cross-correlations, which are approximations of the empirical Green's function of the interstation path. We determine group and phase velocity surface wave dispersion maps with a 0.1 decimal degree resolution and resolve a zone of slow surface wave velocity south-east of the city of Medina, which is spatially correlated with the most recent historical eruption (the 1256 CE Medina eruption). Dispersion curves are calculated at each grid-point of the surface-wave velocity maps and inverted to obtain measurements of shear-velocity with depth. The 1D velocity models are then used to produce average shear-velocity models for the volcanic field. A shear-velocity increase ranging from 0.5 to 1.0 km/s, suggesting a layer interface, is detected at approximately 20 km depth and compared to P-wave measurement from a previous refraction study. We compute cross-section profiles by interpolating the inversions into a pseudo-3D model and resolve a zone of slow shear-velocity below the 1256 CE eruption location. These areas are also spatially correlated with low values of Bouguer gravity. We hypothesize that the low shear-velocity and gravity measurements are caused by fluids and fractures created from prior volcanic eruptions.   We use the coda of cross-correlations extracted from ambient noise to determine shear-velocity changes at Rotokawa and Ngatamariki, two electricity producing geothermal fields located in the North Island of New Zealand. Stacks of cross correlations between stations prior to the onset of production are compared to cross correlations of moving stacks in time periods of well stimulation and the onset of electricity production using the Moving Window Cross Spectral technique. An increase between 0.05% to 0.1% of shear-velocity is detected at Rotokawa coinciding with an increase of injection. The shear-velocity subsequently decreases by approximately 0.1% when the rate of production surpasses the rate of injection. A similar amplitude shear-velocity increase is detected at Ngatamariki during the beginning of injection. After the initial increase, the shear-velocity at Ngatamariki fluctuates in response to differences in injection and production rates. A straight-ray pseudo-tomography analysis is conducted at the geothermal fields, which reveals that localized positive velocity changes are co-located with injection wells.  Lastly, we use ambient noise and active sources at the Ngatamariki geothermal field to determine the structure of the top 200 meters using the Refraction Microtremor technique. We deployed a linear 72-channel array of vertical geophones with ten meter spacing at two locations of the geothermal field and determine average 1D and 2D shear-velocity profiles. We were able to image depths between 57 to 93 meters for 2D profiles and up to 165 meters for 1D profiles. A shear-velocity anomaly was detected across one of the lines that coincided with the inferred location of a fault determined from nearby well logs. This suggests that the method can be used to cheaply and quickly constrain near-surface geology at geothermal fields, where ambient noise is abundant and typical reflection and refraction surveys require large inputs of energy and are hindered by attenuation and scattering in near-surface layers.</p>

scholarly journals Determining seismic shear-velocity from ambient noise sources at regional and local scales

2021 ◽  
Author(s):  
◽  
Francesco Civilini
Keyword(s):  
Surface Wave ◽  
Ambient Noise ◽  
Shear Velocity ◽  
Volcanic Field ◽  
Geothermal Field ◽  
Near Surface ◽  
Velocity Models ◽  
Geothermal Fields ◽  
Cross Correlations ◽  
Velocity Changes

<p>We present three projects that use different bandwidths of the ambient noise spectrum to solve geophysical problems. Specifically, we use signals within the noise field to determine surface and shear wave velocities, image the shallow and deep crust, and monitor time-dependent deformation resulting from geothermal fluid injection and extraction.  Harrat Al-Madinah, a Cenozoic bimodal alkaline volcanic field in west-central Saudi Arabia, is imaged using shear-velocities obtained from natural ambient seismic noise. To our knowledge, this project is the first analysis of Saudi Arabia structure using ambient noise methods. Surface wave arrivals are extracted from a year's worth of station-pair cross-correlations, which are approximations of the empirical Green's function of the interstation path. We determine group and phase velocity surface wave dispersion maps with a 0.1 decimal degree resolution and resolve a zone of slow surface wave velocity south-east of the city of Medina, which is spatially correlated with the most recent historical eruption (the 1256 CE Medina eruption). Dispersion curves are calculated at each grid-point of the surface-wave velocity maps and inverted to obtain measurements of shear-velocity with depth. The 1D velocity models are then used to produce average shear-velocity models for the volcanic field. A shear-velocity increase ranging from 0.5 to 1.0 km/s, suggesting a layer interface, is detected at approximately 20 km depth and compared to P-wave measurement from a previous refraction study. We compute cross-section profiles by interpolating the inversions into a pseudo-3D model and resolve a zone of slow shear-velocity below the 1256 CE eruption location. These areas are also spatially correlated with low values of Bouguer gravity. We hypothesize that the low shear-velocity and gravity measurements are caused by fluids and fractures created from prior volcanic eruptions.   We use the coda of cross-correlations extracted from ambient noise to determine shear-velocity changes at Rotokawa and Ngatamariki, two electricity producing geothermal fields located in the North Island of New Zealand. Stacks of cross correlations between stations prior to the onset of production are compared to cross correlations of moving stacks in time periods of well stimulation and the onset of electricity production using the Moving Window Cross Spectral technique. An increase between 0.05% to 0.1% of shear-velocity is detected at Rotokawa coinciding with an increase of injection. The shear-velocity subsequently decreases by approximately 0.1% when the rate of production surpasses the rate of injection. A similar amplitude shear-velocity increase is detected at Ngatamariki during the beginning of injection. After the initial increase, the shear-velocity at Ngatamariki fluctuates in response to differences in injection and production rates. A straight-ray pseudo-tomography analysis is conducted at the geothermal fields, which reveals that localized positive velocity changes are co-located with injection wells.  Lastly, we use ambient noise and active sources at the Ngatamariki geothermal field to determine the structure of the top 200 meters using the Refraction Microtremor technique. We deployed a linear 72-channel array of vertical geophones with ten meter spacing at two locations of the geothermal field and determine average 1D and 2D shear-velocity profiles. We were able to image depths between 57 to 93 meters for 2D profiles and up to 165 meters for 1D profiles. A shear-velocity anomaly was detected across one of the lines that coincided with the inferred location of a fault determined from nearby well logs. This suggests that the method can be used to cheaply and quickly constrain near-surface geology at geothermal fields, where ambient noise is abundant and typical reflection and refraction surveys require large inputs of energy and are hindered by attenuation and scattering in near-surface layers.</p>

Near-surface Structure of Downtown Jinan, China: Application of Ambient Noise Tomography with a Dense Seismic Array

2019 ◽  
Vol 24 (4) ◽  
pp. 641-652
Author(s):  
Feng Liang ◽  
Zhihui Wang ◽  
Hailong Li ◽  
Kai Liu ◽  
Tao Wang
Keyword(s):  
Surface Wave ◽  
Human Activities ◽  
High Frequency ◽  
Urban Areas ◽  
Ambient Noise ◽  
Velocity Structure ◽  
Shear Velocity ◽  
Ambient Noise Tomography ◽  
Surface Wave Dispersion ◽  
Near Surface

Urban geophysics ups the ante in the world of applied geophysics, which requires innovative thinking and seemingly off-the-wall approaches, if for no other reason than the settings. Ambient-noise-tomography (ANT) can play a pivotal role in yielding subsurfa2ce information in urban areas, which is capable of dealing with challenges related to these scenarios ( e.g., human activities and low signal-to-noise ratio). In this study, the ANT was conducted to investigate the near-surface shear-velocity structure in the surrounding area of the Baotu Spring Park in downtown Jinan, Shandong Province, China. Quiet clear Rayleigh waves have been obtained by the cross-correlation, which indicates that strong human activities, such as moving vehicles and municipal engineering constructions, can produce approximately isotropic distribution of noise sources for high-frequency signals. The direct surface-wave tomographic method with period-dependent ray-tracing was used to invert all surface-wave dispersion data in the period band 0.2-1.5 s simultaneously for 3D variations of shear-velocity (Vs) structure. Our results show a good correspondence to the geological features with thinner Quaternary sediments, the geological structural characteristic of the limestone surrounded by the igneous which has the highest velocity than that of the limestone in the study area, and several concealed faults of which specific location has been detected at depth. The results demonstrate that it is possible to successfully use ANT with high-frequency signal in an urban environment provided a detailed planning and execution is implemented.

scholarly journals Using ambient seismic noise to study temporal and spatial surface wave velocity structures and ambient noise field characteristics of central South Island, New Zealand

2021 ◽  
Author(s):  
◽  
Rachel Heckels
Keyword(s):  
New Zealand ◽  
Surface Wave ◽  
Wave Velocity ◽  
Seismic Noise ◽  
Shear Velocity ◽  
Southern Alps ◽  
Ambient Seismic Noise ◽  
Near Surface ◽  
Surface Wave Velocity ◽  
Velocity Changes

<p>Ambient seismic noise is used to examine the spatial and temporal surface wave velocity structures and ambient seismic noise fields in the vicinity of different fault zone environments. This study focuses on two distinct regions of central South Island, New Zealand. The Canterbury Plains is a sedimentary basin with many minor faults, which was considered to have low seismic hazard prior to the 2010 – 2011 Canterbury earthquake sequence. We focus on the time period immediately following the 2010 Darfield earthquake, which ruptured the previously unmapped Greendale Fault. The second region of interest is the central Southern Alps. The locked portion of the Alpine Fault currently poses one of the largest seismic hazards for New Zealand. The wealth of data from both permanent and temporary seismic deployments in these regions make them ideal areas in which to assess the effectiveness of ambient noise for velocity modelling in regions surrounding faults at different stages of their seismic cycles.  Temporal velocity changes are measured following the Mw 7.1 Darfield earthquake of 4 September 2010 in the Canterbury Plains. Nine-component cross-correlations are computed from temporary and permanent seismic stations lying on and surrounding the Greendale Fault. Using the Moving-Window Cross-Spectral method, surface wave velocity changes are calculated for the four months immediately following the earthquake until 10 January 2011, for 0.1 — 1.0 Hz. An average increase in seismic velocity of 0.14 ± 0.04 % is determined throughout the region, providing the first such estimate of postseismic relaxation rates in Canterbury. Depth analyses further showed that velocity changes are confined to the uppermost 5 km of the subsurface and we attribute this to postseismic relaxation via crack-healing of the Greendale Fault and throughout the surrounding region.  Rayleigh and Love wave dispersion is examined throughout the Canterbury region. Multi-component cross-correlation functions are analysed for group and phase dispersion curves. These are inverted using frequency-time analysis for 2-D phase and group velocity maps of Rayleigh and Love waves. A high-velocity zone to the southeast of the region coincides with volcanic rocks of Banks Peninsula. Dispersion curves generated from the surface wave tomography are further inverted for one-dimensional shear velocity profiles. These models show a thin, low-velocity near surface layer consistent with the basin sediments, which thins towards the foothills of the Southern Alps. A near-surface damage zone is identified along the length of the Greendale Fault, with consistent reduced Vs velocities to depth of up to 5 km.  Surface and shear wave velocity maps are computed for the central Southern Alps to image the seismic structure of the region. Tomographic surface maps at periods of 5 – 12 s are produced from dispersion measurements of three-component cross-correlation functions. At periods of 5 – 8 s a strong NE-SW trending velocity contrast highlights the Alpine Fault. One-dimensional shear velocity models, computed from the surface wave maps, are in agreement with previous models produced by other conventional methods. An analysis of surface wave amplitudes through signal-to-noise ratios of cross-correlations reveals strong directional effects. Calculated signal-to-noise ratios are up to eight times higher for surface waves travelling north-west than for waves travelling to the south or east. We attribute this to a combination of more energetic ocean wave signals from the Southern Ocean compared to the Tasman Sea.</p>

Shear wave velocity changes induced by earthquakes and rainfall at the Rotokawa and Ngatamariki geothermal fields, Taupō Volcanic Zone, New Zealand

2019 ◽  
Vol 221 (1) ◽  
pp. 97-114
Author(s):  
F Civilini ◽  
M K Savage ◽  
J Townend
Keyword(s):  
Shear Wave ◽  
Wave Velocity ◽  
Shear Wave Velocity ◽  
Ambient Noise ◽  
Shear Velocity ◽  
Taupo Volcanic Zone ◽  
Fluid Injection ◽  
Volcanic Zone ◽  
Geothermal Fields ◽  
Velocity Changes

SUMMARY Fluid injection for geothermal production has the potential to produce subsidence and microseismicity that can incur heavy financial cost or hazard. Due to this, novel ways to monitor subsurface deformation to supplement existing methods are highly sought after. We use seismic ambient noise to obtain time-dependent measurements of shear velocity within the geothermal reservoirs of Rotokawa and Ngatamariki, two producing geothermal fields in the Taupō Volcanic Zone located in the central North Island of New Zealand and operated by Mercury Energy. We investigate the relationship between shear wave velocity changes and geothermal injection by selecting time periods at the fields when injection and production volumes were significantly altered: 2009–2010 at Rotokawa, when geothermal injection was quadrupled due to the start-up of a new power station, and 2012–2013 at Ngatamariki, the beginning of geothermal injection for electricity production at that field. Shear wave velocity changes are computed from the ambient noise cross-correlation coda using the Moving Window Cross-Spectral (MWCS) technique, with a reference stack encompassing all data prior to the change in injection rate and moving stacks of 10–50 d. Gradual positive and negative shear velocity changes with a periodicity of approximately 12 months were observed at both sites, with maximum amplitude of 0.06 ± 0.04 and –0.08 ± 0.03 per cent at Rotokawa and 0.07 ± 0.03 and –0.06 ± 0.02 per cent at Ngatamariki. We hypothesize that these changes are due to seasonal rainfall, as seismic velocities computed by ambient noise are sensitive to the filling and emptying of near-surface pore space. In addition to these gradual responses, we found several sharp negative changes in velocity that reach minimum values over a few days and then gradually equilibrate to prior values over a few weeks. The amplitude of these responses is between –0.03 and –0.07 per cent and coincides with regional and local earthquakes. We hypothesize that these responses are primarily produced by the creation of new fractures, the same mechanism that produces gradual groundwater level decreases at regional distances from earthquake epicentres. We analyse a periodic signal within the time-delay measurements and determine that it is at least in part caused by the MWCS window size smoothing the cross-coherence of the ambient seismic signal. We do not observe shear wave velocity changes coinciding with geothermal injection, which may suggest that the signal has lower amplitude compared to the seasonal and seismic responses. We use bandstop filters and polynomial curve fitting to remove the contribution of the seasonal signal, but see no evidence of a shear wave velocity response due to geothermal fluid injection.

Inferring shear-velocity structure of the upper 200 m using cultural ambient noise at the Ngatamariki geothermal field, Central North Island, New Zealand

2016 ◽  
Vol 4 (3) ◽  
pp. SJ87-SJ101 ◽  
Author(s):  
Francesco Civilini ◽  
Aasha Pancha ◽  
Martha Kane Savage ◽  
Steven Sewell ◽  
John Townend
Keyword(s):  
New Zealand ◽  
Ambient Noise ◽  
Velocity Structure ◽  
Shear Velocity ◽  
Geothermal Field ◽  
S Wave ◽  
Near Surface ◽  
Shallow Subsurface ◽  
Spatially Correlated ◽  
Geothermal Wells

We have determined subsurface structure using the refraction microtremor (ReMi) method at the Ngatamariki geothermal field, Central North Island, New Zealand. The local geology is such that refraction and reflection studies are hindered by energy scattering and attenuation in the near-surface layers. The ReMi method uses surface waves from ambient noise and active sources to determine S-wave velocities in the shallow subsurface. We have deployed two lines of 72-channel, 10 Hz vertical geophones with 10 m spacing, and we were able to model near-surface S-wave velocity to depths of 57–93 m for 2D profiles and as much as 165 m for 1D profiles. Shear-velocity anomalies were detected on one line that were spatially correlated with a fault. The location of the fault was previously inferred from stratigraphic offset in the geothermal wells, suggesting that the ReMi method can provide important constraints on near-surface geology in noisy geothermal settings.

scholarly journals Using ambient seismic noise to study temporal and spatial surface wave velocity structures and ambient noise field characteristics of central South Island, New Zealand

2021 ◽  
Author(s):  
◽  
Rachel Heckels
Keyword(s):  
New Zealand ◽  
Surface Wave ◽  
Wave Velocity ◽  
Seismic Noise ◽  
Shear Velocity ◽  
Southern Alps ◽  
Ambient Seismic Noise ◽  
Near Surface ◽  
Surface Wave Velocity ◽  
Velocity Changes

<p>Ambient seismic noise is used to examine the spatial and temporal surface wave velocity structures and ambient seismic noise fields in the vicinity of different fault zone environments. This study focuses on two distinct regions of central South Island, New Zealand. The Canterbury Plains is a sedimentary basin with many minor faults, which was considered to have low seismic hazard prior to the 2010 – 2011 Canterbury earthquake sequence. We focus on the time period immediately following the 2010 Darfield earthquake, which ruptured the previously unmapped Greendale Fault. The second region of interest is the central Southern Alps. The locked portion of the Alpine Fault currently poses one of the largest seismic hazards for New Zealand. The wealth of data from both permanent and temporary seismic deployments in these regions make them ideal areas in which to assess the effectiveness of ambient noise for velocity modelling in regions surrounding faults at different stages of their seismic cycles.  Temporal velocity changes are measured following the Mw 7.1 Darfield earthquake of 4 September 2010 in the Canterbury Plains. Nine-component cross-correlations are computed from temporary and permanent seismic stations lying on and surrounding the Greendale Fault. Using the Moving-Window Cross-Spectral method, surface wave velocity changes are calculated for the four months immediately following the earthquake until 10 January 2011, for 0.1 — 1.0 Hz. An average increase in seismic velocity of 0.14 ± 0.04 % is determined throughout the region, providing the first such estimate of postseismic relaxation rates in Canterbury. Depth analyses further showed that velocity changes are confined to the uppermost 5 km of the subsurface and we attribute this to postseismic relaxation via crack-healing of the Greendale Fault and throughout the surrounding region.  Rayleigh and Love wave dispersion is examined throughout the Canterbury region. Multi-component cross-correlation functions are analysed for group and phase dispersion curves. These are inverted using frequency-time analysis for 2-D phase and group velocity maps of Rayleigh and Love waves. A high-velocity zone to the southeast of the region coincides with volcanic rocks of Banks Peninsula. Dispersion curves generated from the surface wave tomography are further inverted for one-dimensional shear velocity profiles. These models show a thin, low-velocity near surface layer consistent with the basin sediments, which thins towards the foothills of the Southern Alps. A near-surface damage zone is identified along the length of the Greendale Fault, with consistent reduced Vs velocities to depth of up to 5 km.  Surface and shear wave velocity maps are computed for the central Southern Alps to image the seismic structure of the region. Tomographic surface maps at periods of 5 – 12 s are produced from dispersion measurements of three-component cross-correlation functions. At periods of 5 – 8 s a strong NE-SW trending velocity contrast highlights the Alpine Fault. One-dimensional shear velocity models, computed from the surface wave maps, are in agreement with previous models produced by other conventional methods. An analysis of surface wave amplitudes through signal-to-noise ratios of cross-correlations reveals strong directional effects. Calculated signal-to-noise ratios are up to eight times higher for surface waves travelling north-west than for waves travelling to the south or east. We attribute this to a combination of more energetic ocean wave signals from the Southern Ocean compared to the Tasman Sea.</p>

scholarly journals Crustal structures beneath the Eastern and Southern Alps from ambient noise tomography

2020 ◽  
Author(s):  
Ehsan Qorbani ◽  
Dimitri Zigone ◽  
Mark R. Handy ◽  
Götz Bokelmann ◽  
Keyword(s):  
High Resolution ◽  
Crustal Structure ◽  
Ambient Noise ◽  
Eastern Alps ◽  
Southern Alps ◽  
Ambient Noise Tomography ◽  
Tauern Window ◽  
Velocity Models ◽  
Cross Correlations ◽  
Velocity Anomalies

Abstract. We study the crustal structure under the Eastern and Southern Alps using ambient noise tomography. We use cross-correlations of ambient seismic noise between pairs of 71 permanent stations and 19 stations of the EASI profile to derive new high-resolution 3-D shear-velocity models for the crust. Continuous records from 2014 and 2015 are cross-correlated to estimate Green's functions of Rayleigh and Love waves propagating between the station pairs. Group velocities extracted from the cross-correlations are inverted to obtain isotropic 3-D Rayleigh and Love-wave shear-wave velocity models. Our high resolution models image several velocity anomalies and contrasts and reveal details of the crustal structure. Velocity variations at short periods correlate very closely with the lithologies of tectonic units at the surface and projected to depth. Low-velocity zones, associated with the Po and Molasse sedimentary basins, are imaged well to the south and north of the Alps, respectively. We find large high-velocity zones associated with the crystalline basement that forms the core of the Tauern Window. Small-scale velocity anomalies are also aligned with geological units such as the Ötztal and the Gurktal nappes of the Austroalpine nappes. Clear velocity contrasts in the Tauern Window along vertical cross-sections of the velocity model show the depth extent of the tectonic units and their bounding faults. A mid-crustal velocity contrast is interpreted as a manifestation of intracrustal decoupling in the Eastern Alps and decoupling between the Southern and Eastern Alps.

Sign in / Sign up

Export Citation Format

Share Document