Imaging New Zealand's Crustal Structure Using Ambient Seismic Noise Recordings from Permanent and Temporary Instruments

<p>We use ambient seismic noise to image the crust and uppermost mantle, and to determine the spatiotemporal characteristics of the noise field itself, and examine the way in which those characteristics may influence imaging results. Surface wave information extracted from ambient seismic noise using cross-correlation methods significantly enhances our knowledge of the crustal and uppermost mantle shear-velocity structure of New Zealand. We assemble a large dataset of three-component broadband continuous seismic data from temporary and permanent seismic stations, increasing the achievable resolution of surface wave velocity maps in comparison to a previous study. Three-component data enables us to examine both Rayleigh and Love waves using noise cross-correlation functions. Employing a Monte Carlo inversion method, we invert Rayleigh and Love wave phase and group velocity dispersion curves separately for spatially averaged isotropic shear velocity models beneath the Northland Peninsula. The results yield first-order radial anisotropy estimates of 2% in the upper crust and up to 15% in the lower crust, and estimates of Moho depth and uppermost mantle velocity compatible with previous studies. We also construct a high-resolution, pseudo-3D image of the shear-velocity distribution in the crust and uppermost mantle beneath the central North Island using Rayleigh and Love waves. We document, for the first time, the lateral extent of low shear-velocity zones in the upper and mid-crust beneath the highly active Taupo Volcanic Zone, which have been reported previously based on spatially confined 1D shear-velocity profiles. Attributing these low shear-velocities to the presence of partial melt, we use an empirical relation to estimate an average percentage of partial melt of < 4:2% in the upper and middle crust. Analysis of the ambient seismic noise field in the North Island using plane wave beamforming and slant stacking indicates that higher mode Rayleigh waves can be detected, in addition to the fundamental mode. The azimuthal distributions of seismic noise sources inferred from beamforming are compatible with high near-coastal ocean wave heights in the period band of the secondary microseism (~7 s). Averaged over 130 days, the distribution of seismic noise sources is azimuthally homogeneous, indicating that the seismic noise field is well-suited to noise cross-correlation studies. This is underpinned by the good agreement of our results with those from previous studies. The effective homogeneity of the seismic noise field and the large dataset of noise cross-correlation functions we here compiled, provide the cornerstone for future studies of ambient seismic noise and crustal shear velocity structure in New Zealand.</p>

Imaging New Zealand's Crustal Structure Using Ambient Seismic Noise Recordings from Permanent and Temporary Instruments

<p>We use ambient seismic noise to image the crust and uppermost mantle, and to determine the spatiotemporal characteristics of the noise field itself, and examine the way in which those characteristics may influence imaging results. Surface wave information extracted from ambient seismic noise using cross-correlation methods significantly enhances our knowledge of the crustal and uppermost mantle shear-velocity structure of New Zealand. We assemble a large dataset of three-component broadband continuous seismic data from temporary and permanent seismic stations, increasing the achievable resolution of surface wave velocity maps in comparison to a previous study. Three-component data enables us to examine both Rayleigh and Love waves using noise cross-correlation functions. Employing a Monte Carlo inversion method, we invert Rayleigh and Love wave phase and group velocity dispersion curves separately for spatially averaged isotropic shear velocity models beneath the Northland Peninsula. The results yield first-order radial anisotropy estimates of 2% in the upper crust and up to 15% in the lower crust, and estimates of Moho depth and uppermost mantle velocity compatible with previous studies. We also construct a high-resolution, pseudo-3D image of the shear-velocity distribution in the crust and uppermost mantle beneath the central North Island using Rayleigh and Love waves. We document, for the first time, the lateral extent of low shear-velocity zones in the upper and mid-crust beneath the highly active Taupo Volcanic Zone, which have been reported previously based on spatially confined 1D shear-velocity profiles. Attributing these low shear-velocities to the presence of partial melt, we use an empirical relation to estimate an average percentage of partial melt of < 4:2% in the upper and middle crust. Analysis of the ambient seismic noise field in the North Island using plane wave beamforming and slant stacking indicates that higher mode Rayleigh waves can be detected, in addition to the fundamental mode. The azimuthal distributions of seismic noise sources inferred from beamforming are compatible with high near-coastal ocean wave heights in the period band of the secondary microseism (~7 s). Averaged over 130 days, the distribution of seismic noise sources is azimuthally homogeneous, indicating that the seismic noise field is well-suited to noise cross-correlation studies. This is underpinned by the good agreement of our results with those from previous studies. The effective homogeneity of the seismic noise field and the large dataset of noise cross-correlation functions we here compiled, provide the cornerstone for future studies of ambient seismic noise and crustal shear velocity structure in New Zealand.</p>

Microtremor surveys based on rotational seismology: theoretical analysis with focus on separation of Rayleigh and Love waves in general wavefield of microtremors

SUMMARY The applicability of rotational seismology to the general wavefield of microtremors is theoretically demonstrated based on a random process model of a 2-D wavefield. We show the effectiveness of taking the rotations (i.e. spatial differentiation) of microtremor waveforms in separating the Rayleigh and Love waves in a wavefield where waves are simultaneously arriving from various directions with different intensities. This means that a method based on rotational seismology (a rotational method) is capable of separating Rayleigh and Love waves without adopting a specific array geometry or imposing a specific assumption on the microtremor wavefield. This is an important feature of a rotational method because the spatial autocorrelation (SPAC) method, a conventional approach for determining phase velocities in microtremor array surveys, requires either the use of a circular array or the assumption of an isotropic wavefield (i.e. azimuthal averaging of correlations is required). Derivatives of the SPAC method additionally require the assumption that Rayleigh and Love waves are uncorrelated. We also show that it is possible to apply a rotational method to determine the characteristics of Love waves based on a simple three-point microtremor array that consists of translational (i.e. ordinary) three-component sensors. In later sections, we assume realistic data processing for microtremor arrays with translational sensors to construct a theoretical model to evaluate the effects of approximating spatial differentiation via finite differencing (i.e. array-derived rotation, ADR) and the effects of incoherent noise on analysis results. Using this model, it is shown that in a short-wavelength range compared to the distance for finite differencing (e.g. $\lambda &lt; 3h$, where $\lambda $ and $h$ are the wavelength and distance for finite differencing, respectively), the leakage of unwanted wave components can determine the analysis limit. It is also shown that in a long-wavelength range (e.g. $\lambda &gt; 3h$), the signal intensity gradually decreases, and thus the effects of incoherent noise increase (i.e. the signal-to-noise ratio, SNR decreases) and determine the analysis limit. We derive the relation between the SNR and wavelength. Although the analysis results quantitatively depend on the array geometry used for finite differencing, the qualitative understanding supported by mathematical expressions with a physically clear meaning can serve as a guideline for the treatment of data obtained from ADR.

Beamforming Reliability of DAS Ambient Noise Data and Wave Modes Identification

&lt;p&gt;Ambient noise generated by the anthropological activities in the urban environments may contain both Rayleigh and Love waves. Due to the differences in the physics of Rayleigh and Love waves, a pre-knowledge of the wave modes in the cross-correlogram is essential for an accurate inversion of the subsurface velocity model. Several studies (Martin and Biondi, 2017; Martin et al., 2017; Luo et al., 2020) demonstrated that only Rayleigh waves can be extracted by cross-correlation if the virtual source is colinear with the DAS array based on the assumption that the ambient noise sources are random and uniformly distributed. However, in realistic cases, ambient noise sources may come from a certain direction (e.g., Dou et al., 2017; Zhang et al., 2019). Moreover, the source propagation direction should be resolved and used to correct the apparent dispersion curves. Zhao et al. (2020) and van den Ende et al. (2020) proposed that beamforming results are not always reliable due to the measurements of DAS.&lt;/p&gt;&lt;p&gt;Based on the synthetic DAS ambient noise data recorded by a near &amp;#8220;L&amp;#8221; shape array (Source-West corner of the Stanford DAS-1 array), we prove that beamforming can resolve the source direction when the ambient sources are mainly coming from one direction. Two important processing procedures are that: check the polarity in the data and apply polarity flip on one part of the data; apply amplitude normalization on the data if strong amplitude difference exits in the data. Based on the source direction, the coordinate of the DAS array, and amplitude ratio of the data recorded by the two segments of the DAS array, we propose an inversion method to calculate the amplitude ratio of the Rayleigh and Love waves generated by the ambient sources.&lt;/p&gt;&lt;p&gt;We apply the method to two 100-second DAS ambient noise data recorded by the Stanford DAS-1 array. We first resolve the source propagation direction from the two data. The results indicate that the ambient noise in the data were mainly generated by the motor vehicles running on the Campus Drive in the northwest of the array. Then we invert for the Rayleigh and Love waves amplitude ratio using the proposed method. The ratios for the two data are 0.2 and 0.13, respectively. The results suggest that the ambient noise generated by motor vehicles running on the northwest corner of the Campus Drive mainly contain Love waves.&lt;/p&gt;

Source directionality of ambient seismic noise inferred from three-component beamforming

The increased use of ambient seismic noise for seismic imaging requires better understanding of the ambient seismic noise wavefield and its source locations and mechanisms. Although the source regions and mechanisms of Rayleigh waves have been studied extensively, characterization of Love wave source processes are sparse or absent. We present here the first systematic comparison of ambient seismic noise source directions within the primary (~10-20 s period) and secondary (~5-10 s period) microseism bands for both Rayleigh and Love waves in the Southern Hemisphere using vertical- and horizontal-component ambient seismic noise recordings from a dense temporary network of 68 broadband seismometers in New Zealand. Our analysis indicates that Rayleigh and Love waves within the primary microseism band appear to be mostly generated in different areas, whereas in the secondary microseism band they arrive from similar backazimuths. Furthermore, the source areas of surface waves within the secondary microseism band correlate well with modeled deep-water and near-coastal source regions. Key Points Rayleigh and Love wave source regions of the secondary microseism are co-located Rayleigh and Love wave source regions of the primary microseism differ strongly Observed and modeled source directions for the secondary microseism agree well ©2012. American Geophysical Union. All Rights Reserved.

Source directionality of ambient seismic noise inferred from three-component beamforming

The increased use of ambient seismic noise for seismic imaging requires better understanding of the ambient seismic noise wavefield and its source locations and mechanisms. Although the source regions and mechanisms of Rayleigh waves have been studied extensively, characterization of Love wave source processes are sparse or absent. We present here the first systematic comparison of ambient seismic noise source directions within the primary (~10-20 s period) and secondary (~5-10 s period) microseism bands for both Rayleigh and Love waves in the Southern Hemisphere using vertical- and horizontal-component ambient seismic noise recordings from a dense temporary network of 68 broadband seismometers in New Zealand. Our analysis indicates that Rayleigh and Love waves within the primary microseism band appear to be mostly generated in different areas, whereas in the secondary microseism band they arrive from similar backazimuths. Furthermore, the source areas of surface waves within the secondary microseism band correlate well with modeled deep-water and near-coastal source regions. Key Points Rayleigh and Love wave source regions of the secondary microseism are co-located Rayleigh and Love wave source regions of the primary microseism differ strongly Observed and modeled source directions for the secondary microseism agree well ©2012. American Geophysical Union. All Rights Reserved.

Local Coupling and Conversion of Surface Waves due to Earth’s Rotation. Part 2: Numerical Examples

Summary Rotation of the Earth affects the propagation of seismic waves. The global coupling of spheroidal and toroidal modes by the Coriolis force over time is described by normal-mode theory. The local action of the Coriolis force on the propagation of surface waves can be described by coefficients for the coupling between propagating Rayleigh and Love waves as derived by (Landau & Lifshitz 1959). Using global wavefield simulations we show how the Coriolis force leads to coupling and conversion between both surface wave types depending on latitude, propagation direction, frequency, and local velocity structure. Surface wave coupling is most efficient for periods where the modes have similar phase velocities, a condition that is equivalent to the selection rules of the angular degree in the normal-mode framework, a phenomenon that we refer to as resonant coupling. In the time-domain, resonant coupling gradually converts energy from one wave type–Rayleigh waves or Love wave–into the other, which then propagates independently. Due to the lateral heterogeneity, the condition of equal phase velocity renders the rotational coupling location-dependent. East-west oriented ray path segments and segments at high latitudes (across the Poles) only weakly couple the fundamental mode Rayleigh and Love waves while coupling is strongest for propagation along the meridians across the equator. At 250 s period, where Love and Rayleigh waves have similar phase velocities, the net energy transfer from Rayleigh to Love wave reaches about 10% for one orbit.