scholarly journals Seismic Structure Beneath the Wellington Region from Receiver Functions

2021 ◽  
Author(s):  
◽  
Lucy Caroline Hall

<p>Seismic velocity structures, interpreted as being associated with the Hikurangi subduction system beneath the lower North Island of New Zealand, are imaged using stacked P wave receiver functions computed using teleseismic earthquakes. Receiver functions are a seismological technique that exploits the phenomenon of wave conversion. The upcoming P wave interacts with seismic velocity impedance contrasts below the receiving station to produce polarized P to SV converted phases. The time delay between the first arriving P wave and the SV converted phase is interpreted to infer the depth of interfaces and the velocity structure directly below the receiver, allowing estimates to be made of the physical properties of the interface. Passive seismic data were recorded at eighteen seismic stations deployed across a ~90km transect stretching across the breadth of lower North Island of New Zealand, from Kapiti Island, 5km off the west coast, to the eastern coast. The transect is oriented normal to the strike of the subducting Pacific Plate, as it dives beneath the overriding Australian Plate. Data were recorded at 10 broadband and 2 short period sensors, deployed as part of the Seismic Array Hikurangi Project (SAHKE 1 deployment), 3 Geonet (New Zealand Geonet Project) permanent short period stations, and 3 temporary stations from part of the 1991-1992 POMS project. Seismic data were recorded between November 2009 and March 2010 on the short period sensors and up to 18 months on the broadband sensor. Data recorded between November 2009 and November 2011 were utilised from the Geonet stations. P wave receiver functions are computed using the multi-taper correlation method using 389 > 6.0 Mw teleseismic earthquakes recorded at the individual seismic stations. A total of 1082 individual receiver functions from all the stations are stacked for both the individual stations and as a ‘super-stack’ across the complete transect, using the common conversion point (CCP) method. The CCP stack shows a distinct, thick low velocity layer (LVL), dipping to the west, from ~18km depth in the east to ~30km depth in the west. This is above a higher velocity layer, also dipping west, at depths of between ~22km and ~ 37km. The LVL is interpreted as being subducted sediments overlying the higher velocity plate interface. Structures towards the west indicate the presence of possibly imbricated features associated with the overriding plate. Deeper structures, down to a depth of 140km are evident, but have less clarity than the shallower features. Some of the deeper layers appear to be dipping towards the west, some to the east. The results of the CCP stack agree well with results from active source methods.</p>

2021 ◽  
Author(s):  
◽  
Lucy Caroline Hall

<p>Seismic velocity structures, interpreted as being associated with the Hikurangi subduction system beneath the lower North Island of New Zealand, are imaged using stacked P wave receiver functions computed using teleseismic earthquakes. Receiver functions are a seismological technique that exploits the phenomenon of wave conversion. The upcoming P wave interacts with seismic velocity impedance contrasts below the receiving station to produce polarized P to SV converted phases. The time delay between the first arriving P wave and the SV converted phase is interpreted to infer the depth of interfaces and the velocity structure directly below the receiver, allowing estimates to be made of the physical properties of the interface. Passive seismic data were recorded at eighteen seismic stations deployed across a ~90km transect stretching across the breadth of lower North Island of New Zealand, from Kapiti Island, 5km off the west coast, to the eastern coast. The transect is oriented normal to the strike of the subducting Pacific Plate, as it dives beneath the overriding Australian Plate. Data were recorded at 10 broadband and 2 short period sensors, deployed as part of the Seismic Array Hikurangi Project (SAHKE 1 deployment), 3 Geonet (New Zealand Geonet Project) permanent short period stations, and 3 temporary stations from part of the 1991-1992 POMS project. Seismic data were recorded between November 2009 and March 2010 on the short period sensors and up to 18 months on the broadband sensor. Data recorded between November 2009 and November 2011 were utilised from the Geonet stations. P wave receiver functions are computed using the multi-taper correlation method using 389 > 6.0 Mw teleseismic earthquakes recorded at the individual seismic stations. A total of 1082 individual receiver functions from all the stations are stacked for both the individual stations and as a ‘super-stack’ across the complete transect, using the common conversion point (CCP) method. The CCP stack shows a distinct, thick low velocity layer (LVL), dipping to the west, from ~18km depth in the east to ~30km depth in the west. This is above a higher velocity layer, also dipping west, at depths of between ~22km and ~ 37km. The LVL is interpreted as being subducted sediments overlying the higher velocity plate interface. Structures towards the west indicate the presence of possibly imbricated features associated with the overriding plate. Deeper structures, down to a depth of 140km are evident, but have less clarity than the shallower features. Some of the deeper layers appear to be dipping towards the west, some to the east. The results of the CCP stack agree well with results from active source methods.</p>


2021 ◽  
Author(s):  
V Lay ◽  
S Buske ◽  
SB Bodenburg ◽  
John Townend ◽  
R Kellett ◽  
...  

No description supplied


2020 ◽  
pp. 1-38
Author(s):  
Chao Xu ◽  
Pinbo Ding ◽  
Bangrang Di ◽  
Jianxin Wei

We investigated fluid effects on seismic responses using seismic data from a physical modeling experiment. Eight cubic samples with cavities quantitatively filled with air, oil, and water and sixteen non-fluid samples were set within a physical model. Both pre-stack and post-stack seismic responses of the samples were analyzed to quantitatively investigate the fluid effect on the seismic response. It was indicated that fluids could cause detectable changes in both pre-stack and post-stack seismic responses for tight rocks. At first, fluids filled within samples caused changes in pre-stack seismic responses. Visible differences could be detected between angle gathers of the samples filled with air, oil, and water. For the base reflections, the amplitudes at large angles of the air-filled and oil-filed samples are obviously stronger than those of the water-filled sample. In addition, the presence of fluids within samples led to significant changes in post-stack seismic reflections. For samples with similar P-wave impedances to the background, we found strong seismic reflections for the fluid samples and weak or even no reflections for the non-fluid samples. There was notable interference between the top and base reflections for the fluid samples while there was none for the non-fluid samples. Seismic velocities were estimated using the two-way travel times between the top and base reflections. The estimated seismic velocity gently declined with increasing water saturation until 90%. When the water saturation was more than 90%, the seismic velocity showed a steep increase.


2020 ◽  
Author(s):  
Vera Lay ◽  
Stefan Buske ◽  
Sascha Barbara Bodenburg ◽  
Franz Kleine ◽  
John Townend ◽  
...  

&lt;p&gt;The Alpine Fault along the West Coast of the South Island (New Zealand) is a major plate boundary that is expected to rupture in the next 50 years, likely as a magnitude 8 earthquake. The Deep Fault Drilling Project (DFDP) aims to deliver insight into the geological structure of this fault zone and its evolution by drilling and sampling the Alpine Fault at depth. &amp;#160;&lt;/p&gt;&lt;p&gt;Here we present results from a 3D seismic survey around the DFDP-2 drill site in the Whataroa Valley where the drillhole penetrated almost down to the fault surface. Within the glacial valley, we collected 3D seismic data to constrain valley structures that were obscured in previous 2D seismic data. The new data consist of a 3D extended vertical seismic profiling (VSP) survey using three-component receivers and a fibre optic cable in the DFDP-2B borehole as well as a variety of receivers at the surface.&lt;/p&gt;&lt;p&gt;The data set enables us to derive a reliable 3D P-wave velocity model by first-arrival travel time tomography. We identify a 100-460 m thick sediment layer (average velocity 2200&amp;#177;400 m/s) above the basement (average velocity 4200&amp;#177;500 m/s). Particularly on the western valley side, a region of high velocities steeply rises to the surface and mimics the topography. We interpret this to be the infilled flank of the glacial valley that has been eroded into the basement. In general, the 3D structures implied by the velocity model on the upthrown (Pacific Plate) side of the Alpine Fault correlate well with the surface topography and borehole findings.&lt;/p&gt;&lt;p&gt;A reliable velocity model is not only valuable by itself but it is also required as input for prestack depth migration (PSDM). We performed PSDM with a part of the 3D data set to derive a structural image of the subsurface within the Whataroa Valley. The top of the basement identified in the P-wave velocity model coincides well with reflectors in the migrated images so that we can analyse the geometry of the basement in detail.&lt;/p&gt;


2019 ◽  
Vol 219 (1) ◽  
pp. 313-327 ◽  
Author(s):  
Erin Cunningham ◽  
Vedran Lekic

SUMMARY Receiver functions are sensitive to sharp seismic velocity variations with depth and are commonly used to constrain crustal thickness. The H–κ stacking method of Zhu & Kanamori is often used to constrain both the crustal thickness (H) and ${V_P}$/${V_S}$ ratio ($\kappa $) beneath a seismic station using P-to-s converted waves (Ps). However, traditional H–κ stacks require an assumption of average crustal velocity (usually ${V_P}$). Additionally, large amplitude reverberations from low velocity shallow layers, such as sedimentary basins, can overprint sought-after crustal signals, rendering traditional H–$\ \kappa $ stacking uninterpretable. We overcome these difficulties in two ways. When S-wave reverberations from sediment are present, they are removed by applying a resonance removal filter allowing crustal signals to be clarified and interpreted. We also combine complementary Ps receiver functions, Sp receiver functions, and the post-critical P-wave reflection from the Moho (SPmp) to remove the dependence on an assumed average crustal ${V_P}$. By correcting for sediment and combining multiple data sets, the crustal thickness, average crustal P-wave velocity and crustal ${V_P}$/${V_S}$ ratio is constrained in geological regions where traditional H–$\ \kappa $ stacking fails, without making an initial P-wave velocity assumption or suffering from contamination by sedimentary reverberations.


Geophysics ◽  
2010 ◽  
Vol 75 (6) ◽  
pp. V77-V87 ◽  
Author(s):  
Rishi Bansal ◽  
Mike Matheney

Converted-wave (PS) data, when converted to PP time, develop time- and location-varying compression of the seismic wavelet due to a variable subsurface [Formula: see text] [Formula: see text]. The time-dependent compression distorts the wavelet in a seismic trace. The lack of a consistent seismic wavelet in a domain-converted PS volume can eventually lead to an erroneous joint PP/PS inversion result. Depth-converted seismic data also have wavelet distortion due to velocity-dependent wavelet stretch. A high value of seismic velocity produces more stretch in a seismic wavelet than a low value. Variable wavelet stretch renders the depth data unsuitable for attribute analysis. A filtering scheme is proposed that corrects for distortion in seismic wavelets due to domain conversions (PS to PP time and time-to-depth) of seismic data in an amplitude-preserving manner. The method uses a Fourier scaling theorem to predict the seismic wavelet in the converted domain and calculates a shaping filter for each time/depth sample that corrects for the distortion in the wavelet. The filter is applied to the domain-converted data using the method of nonstationary filtering. We provide analytical expressions for the squeeze factor [Formula: see text] that is used to predict the wavelet in the converted domain. The squeeze factor [Formula: see text] for PS to PP time conversion is a function of the subsurface [Formula: see text] whereas for PP time-to-depth conversion [Formula: see text] is dependent on subsurface P-wave velocity. After filtering, the squeezed wavelets in domain-converted PS data appear to have resulted from a constant subsurface [Formula: see text], which we denote as [Formula: see text]. Similarly, the filtered depth-converted data appear to have resulted from a constant subsurface P-wave velocity [Formula: see text].


2021 ◽  
Author(s):  
◽  
Pauline Maria Galea

<p>A shallow aftershock sequence in the Hawkes Bay region of the North Island, New Zealand (May 1990) was recorded with high quality on an L-shaped, 7-station array of 3-component, short-period seismographs at Wellington, such that the seismic waves travelled almost along strike of the subducted Pacific plate in this region. The arrival times at the stations of the Pn wave pulse from a number of aftershocks could be picked sufficiently accurately for a least-squares inversion to be carried out for wavefront speed, c, and incident azimuth, [theta]. The results show a high apparent velocity, 8.7 [plus or minus] 0.2 km/s, and an azimuth which is shifted by 6.0 [plus or minus] 2.5 degrees east of the true epicentre - station azimuth. The azimuthal anomaly, [delta][theta], has been interpreted as due to lateral refraction of Pn off the subducted slab. The effect of different geometries of the slab on the Pn wavefront characteristics (c and [delta][theta]) at Wellington have been explored through both simple geometrical considerations (in the case of a plane or cylindrical slab) as well as through 3-dimensional ray tracing (in the case of irregular curvature of the slab). It has been shown that a plane or cylindrical slab would require P-wave velocities of about 9.0 km/s to exist within it in order to fit both c and [delta][theta], whereas a model of the slab which departs from a regular cylinder and has a small updip component along strike can fit the observations with P-wave velocities of 8.75 km/s in the high velocity medium. This model has been proposed by Ansell and Bannister (1991) after detailed consideration of the shallow seismicity that defines the slab surface in the lower North Island. Information about the nature of the high velocity medium has been obtained by modeling the waveforms through generation of synthetic seismograms by the reflectivity technique of Kennett (1983). The large number of aftershocks within a small source region, and the sampling of much the same wavepath, meant that a sufficient number of seismograms had very similar and characteristic features that could be modelled. The typical seismogram of the data set had a simple Pn wavepulse, followed immediately by a complex. high frequency (up to 15 Hz) phase (here referred to as Phf) and a high amplitude, lower frequency phase that dominated the P-wavetrain (here referred to as P). A velocity profile that contained a layer of 8.75 km/s material at least 4 km thick, underlying "normal" mantle material of P-velocity 8.2 [plus or minus] 0.2 km/s. and whose surface lies approximately 18 km below the slab surface reproduced the observed seismogram features well. The presence of velocity gradients above and below the layer is not excluded. A gradual decrease in velocity below the layer in fact gives a better fit of the Pn pulse shape. By breaking down the synthetic seismogram into simpler versions. using Kennett's wavefield approximation technique, it has been shown that the Pn wave propagates through the high velocity layer, the Phf phase through the overlying layers as a sequence of reflections and refractions, and the P group as a reverberatory phase in a crustal waveguide, with its energy mostly in the form of free surface reflections and S to P conversion. These results have also been confirmed by ray tracing. Waveform modelling has also clearly shown that a low velocity layer (representing subducted sediment) on the top of the subducted slab produces a highly characteristic imprint on the synthetic seismogram, in the form of an energetic, reverberatory, lower frequency signal late in the P-wavetrain. Wavefield approximations show that this is also a crustal waveguide effect, with a strong component of mode conversion at the free surface, but P - S conversion appears to be the dominant mechanism. Seismograms very similar to such synthetic ones have been observed for the Weber aftershocks recorded at stations along the northern East Coast. The presence of such a low velocity layer in the East Coast region is thus implied, consistently with previous proposals. The petrological implications of the high velocity layer in the subducted Pacific plate are discussed. The most likely explanation is that it represents the maximum P velocity of an anisotropic layer within the Pacific upper mantle. It is proposed that the conditions of stress orientation, pressure and temperature at approximately 36 - 50 km depth in this region induces a strong realignment of olivine crystals with their fast direction along strike of the slab, normal to the maximum compressive stress axis. The upper mantle of the segment of the Pacific ocean just east of the Tonga - Kermadec trench and the North Island has been shown in this study to possess P-wave anisotropy, with the P-velocity reaching a maximum of 8.37 km/s in a direction N60 degrees E. This result was obtained by analysing a large set of ISC travel times from earthquakes along the Tonga - Kermadec - New Zealand subduction zone recorded at stations Niue, Rarotonga and the Chatham Islands. It is suggested that an enhancement of this anisotropy, accompanied by some re-orientation, takes place as the upper mantle medium is subjected to the new stress conditions in the initial stages of subduction.</p>


2021 ◽  
Author(s):  
◽  
Erik Ewig

<p>South Wanganui Basin (SWB), New Zealand, is located behind the southern end of the Hikurangi subduction system. One of the most marked geophysical characteristics of the basin is the -150 mGal Bouguer/isostatic gravity anomaly. Sediment fill can only partly explain this anomaly. 3-D gravity models show that the gravity anomaly associated with the basin is generally consistent with a downwarp model of the entire crust. However, the downwarp of the Moho has to be 3-4 times larger than the downwarp of the sediment-basement interface to fit the observed gravity anomaly. Hence a model of lithospheric shortening where ductile thickening of the crust increases with depth is proposed. Finite element modelling demonstrates that the crust, in order to produce the ductile downwarp, is best modelled with at least two distinct different layers. The model requires the top 15-20 km of the crust to behave purely elastic and the lower part (10 km thick) to be viscoelastic with a viscosity of 10[to the power of 21 pascal-seconds]. The existence of this ductile lower continental crust can be explained due to fluids released from the subducting slab accumulating in the lower crust. This is supported by receiver function analysis results. These results propose a 10+/-2 km thick low S-wave velocity layer in the lower crust. The vertical loading necessary to create the basin is high (up to 200MPa) and is difficult to explain by slab pull forces transmitted via a strongly coupled subduction interface alone. An additional driving mechanism proposed is a thickened mantle lithosphere inducing normal forces on the base of the crust. However, the exact origin of the basin remains a puzzling aspect. Receiver function analysis shows that the crust of the subducting Pacific plate underneath the mainland in the lower North Island is abnormally thick ([approximates]10 km) for oceanic crust. This matches with results from the 3-D gravity modelling. Further features discovered with the receiver function analysis are an up to 6 km thick low-velocity layer on top of the slab, which is interpreted as a zone of crushed crustal material with subducted sediments. Furthermore, a deep Moho (39.5+/-1.5 km) is proposed underneath the northern tip of theMarlborough sounds. Shallow seismic and gravity investigations of the southeastern corner of the SWB reveal a complex faulting regime with high-angle normal and reverse faults as well as a component of strike slip. The overall style of faulting in the SWB changes from the west to the east. There are the low-angle thrust faults of the Taranaki Fault zone in the west, the high-angle mostly reverse faults in the eastern part of the basin and the strike slip faults, with a component of vertical movement, at the eastern boundary within the Tararua Ranges.</p>


2021 ◽  
Author(s):  
◽  
Erik Ewig

<p>South Wanganui Basin (SWB), New Zealand, is located behind the southern end of the Hikurangi subduction system. One of the most marked geophysical characteristics of the basin is the -150 mGal Bouguer/isostatic gravity anomaly. Sediment fill can only partly explain this anomaly. 3-D gravity models show that the gravity anomaly associated with the basin is generally consistent with a downwarp model of the entire crust. However, the downwarp of the Moho has to be 3-4 times larger than the downwarp of the sediment-basement interface to fit the observed gravity anomaly. Hence a model of lithospheric shortening where ductile thickening of the crust increases with depth is proposed. Finite element modelling demonstrates that the crust, in order to produce the ductile downwarp, is best modelled with at least two distinct different layers. The model requires the top 15-20 km of the crust to behave purely elastic and the lower part (10 km thick) to be viscoelastic with a viscosity of 10[to the power of 21 pascal-seconds]. The existence of this ductile lower continental crust can be explained due to fluids released from the subducting slab accumulating in the lower crust. This is supported by receiver function analysis results. These results propose a 10+/-2 km thick low S-wave velocity layer in the lower crust. The vertical loading necessary to create the basin is high (up to 200MPa) and is difficult to explain by slab pull forces transmitted via a strongly coupled subduction interface alone. An additional driving mechanism proposed is a thickened mantle lithosphere inducing normal forces on the base of the crust. However, the exact origin of the basin remains a puzzling aspect. Receiver function analysis shows that the crust of the subducting Pacific plate underneath the mainland in the lower North Island is abnormally thick ([approximates]10 km) for oceanic crust. This matches with results from the 3-D gravity modelling. Further features discovered with the receiver function analysis are an up to 6 km thick low-velocity layer on top of the slab, which is interpreted as a zone of crushed crustal material with subducted sediments. Furthermore, a deep Moho (39.5+/-1.5 km) is proposed underneath the northern tip of theMarlborough sounds. Shallow seismic and gravity investigations of the southeastern corner of the SWB reveal a complex faulting regime with high-angle normal and reverse faults as well as a component of strike slip. The overall style of faulting in the SWB changes from the west to the east. There are the low-angle thrust faults of the Taranaki Fault zone in the west, the high-angle mostly reverse faults in the eastern part of the basin and the strike slip faults, with a component of vertical movement, at the eastern boundary within the Tararua Ranges.</p>


Sign in / Sign up

Export Citation Format

Share Document