To see the other types of publications on this topic, follow the link: 3-D S-wave velocity.
Journal articles on the topic '3-D S-wave velocity'
Create a spot-on reference in APA, MLA, Chicago, Harvard, and other styles
Consult the top 50 journal articles for your research on the topic '3-D S-wave velocity.'
Next to every source in the list of references, there is an 'Add to bibliography' button. Press on it, and we will generate automatically the bibliographic reference to the chosen work in the citation style you need: APA, MLA, Harvard, Chicago, Vancouver, etc.
You can also download the full text of the academic publication as pdf and read online its abstract whenever available in the metadata.
Browse journal articles on a wide variety of disciplines and organise your bibliography correctly.
1
Agudo, Òscar Calderón, Nuno Vieira da Silva, George Stronge, and Michael Warner. "Mitigating elastic effects in marine 3-D full-waveform inversion." Geophysical Journal International 220, no. 3 (December 18, 2019): 2089–104. http://dx.doi.org/10.1093/gji/ggz569.
Full textAbstract:
SUMMARY The potential of full-waveform inversion (FWI) to recover high-resolution velocity models of the subsurface has been demonstrated in the last decades with its application to field data. But in certain geological scenarios, conventional FWI using the acoustic wave equation fails in recovering accurate models due to the presence of strong elastic effects, as the acoustic wave equation only accounts for compressional waves. This becomes more critical when dealing with land data sets, in which elastic effects are generated at the source and recorded directly by the receivers. In marine settings, in which sources and receivers are typically within the water layer, elastic effects are weaker but can be observed most easily as double mode conversions and through their effect on P-wave amplitudes. Ignoring these elastic effects can have a detrimental impact on the accuracy of the recovered velocity models, even in marine data sets. Ideally, the elastic wave equation should be used to model wave propagation, and FWI should aim to recover anisotropic models of velocity for P waves (vp) and S waves (vs). However, routine three-dimensional elastic FWI is still commercially impractical due to the elevated computational cost of modelling elastic wave propagation in regions with low S-wave velocity near the seabed. Moreover, elastic FWI using local optimization methods suffers from cross-talk between different inverted parameters. This generally leads to incorrect estimation of subsurface models, requiring an estimate of vp/vs that is rarely known beforehand. Here we illustrate how neglecting elasticity during FWI for a marine field data set that contains especially strong elastic heterogeneities can lead to an incorrect estimation of the P-wave velocity model. We then demonstrate a practical approach to mitigate elastic effects in 3-D yielding improved estimates, consisting of using a global inversion algorithm to estimate a model of vp/vs, employing matching filters to remove elastic effects from the field data, and performing acoustic FWI of the resulting data set. The quality of the recovered models is assessed by exploring the continuity of the events in the migrated sections and the fit of the latter with the recovered velocity model.
2
Zhang, Yan, Ping Wang, and Chongbin Guo. "Low-velocity impact failure mechanism analysis of 3-D braided composites with Hilbert–Huang transform." Journal of Industrial Textiles 46, no. 5 (July 28, 2016): 1241–56. http://dx.doi.org/10.1177/1528083715619955.
Full textAbstract:
This paper aimed to reveal the low-velocity impact responses characteristics and failure mechanism of 3-D braided composites with experimental and frequency domain analysis method, respectively. The low-velocity impact tests were carried out by Instron® 9250 drop-weight instrument with five different impact velocities from 1 m/s to 6 m/s. The results showed that the peak load and absorbed energy increased with the increase of impact velocity. The load–time curves which were in time domain were transformed into frequency domain with Hilbert–Huang Transform (HHT) method. Combined the failure morphologies of 3-D braided composites with frequency domain analysis results, it could be precisely found out the failure mechanism of 3-D braided composites. At the impact velocity of 1 m/s, the 3-D braided composites only had elastic deformations. With the increase of impact velocity, resin crack was the main failure mode of 3-D braided composites. The frequency of impact stress waves which caused the elastic deformation and resin crack mainly located at 0–10 kHz and 60 kHz. When the impact velocity increased to 6 m/s, fiber tows breakage was the main failure mode, and the frequency of impact stress wave located at 15–20 kHz.
3
Mathisen, Mark E., Paul Cunningham, Jesse Shaw, Anthony A. Vasiliou, J. H. Justice, and N. J. Guinzy. "Crosswell seismic radial survey tomograms and the 3-D interpretation of a heavy oil steamflood." GEOPHYSICS 60, no. 3 (May 1995): 651–59. http://dx.doi.org/10.1190/1.1443804.
Full textAbstract:
S‐wave, P‐wave, and Poisson’s ratio tomograms have been used to interpret the 3-D distribution of rock and fluid properties during an early phase of a California heavy oil sand steamflood. Four lines of good quality crosswell seismic data, with source to receiver offsets ranging from 287 to 551 ft (87 to 168 m), were acquired in a radial pattern around a high temperature cemented receiver cable in four days. Processing, first‐arrival picking, and good quality tomographic reconstructions were completed despite offset‐related variations in data quality between the long and short lines. Interpretation was based on correlations with reservoir models, log, core, temperature, and steam injection data. S‐wave tomograms define the 3-D distribution of the “high flow” turbidite channel facies, the “moderate‐low flow” levee facies, porosity, and structural dip. The S‐wave tomograms also define an area with anomalously low S‐wave velocity, which correlates with low shear log velocities and suggests that pressure‐related dilation and compaction may be imageable. P‐wave tomograms define the same reservoir lithology and structure as the S‐wave tomograms and the 3-D distribution of low compressional velocity zones formed by previous steam‐heat injection and the formation of gas. The low P‐wave velocity zones, which are laterally continuous in the “high flow” channel facies near the top of most zones, indicate that the steam‐heat‐gas distribution is controlled by stratification. The stratigraphic control of gas‐bearing zones inferred from P‐wave tomograms is confirmed by Poisson’s ratio tomograms which display low Poisson’s ratios indicative of gas (<0.35) in the same zones as the low P‐wave velocities. The interpretation results indicate that radial survey tomograms can be tied at a central well and used to develop an integrated 3-D geoscience‐engineering reservoir model despite offset‐related variations in data quality. The laterally continuous, stratification‐controlled, low P‐wave velocity zones, which extend up‐dip, suggest that significant amounts of steam‐heat are not heating the surrounding reservoir volume but are flowing updip along “high flow” channels.
4
Movaghari, R., and G. Javan Doloei. "3-D crustal structure of the Iran plateau using phase velocity ambient noise tomography." Geophysical Journal International 220, no. 3 (December 17, 2019): 1555–68. http://dx.doi.org/10.1093/gji/ggz537.
Full textAbstract:
SUMMARY More accurate crustal structure models will help us to better understand the tectonic convergence between Arabian and Eurasian plates in the Iran plateau. In this study, the crustal and uppermost mantle velocity structure of the Iran plateau is investigated using ambient noise tomography. Three years of continuous data are correlated to retrieve Rayleigh wave empirical Green's functions, and phase velocity dispersion curves are extracted using the spectral method. High-resolution Rayleigh wave phase velocity maps are presented at periods of 8–60 s. The tomographic maps show a clear consistency with geological structures such as sedimentary basins and seismotectonic zones, especially at short periods. A quasi-3-D shear wave velocity model is determined from the surface down to 100 km beneath the Iran plateau. A transect of the shear wave velocity model has been considered along with a profile extending across the southern Zagros, the Sanandaj-Sirjan Zone (SSZ), the Urumieh-Dokhtar Magmatic Arc (UDMA) and Central Iran and Kopeh-Dagh (KD). Obvious crustal thinning and thickening are observable along the transect of the shear wave velocity model beneath Central Iran and the SSZ, respectively. The observed shear wave velocities beneath the Iran plateau, specifically Central Iran, support the slab break-off idea in which low density asthenospheric materials drive towards the upper layers, replacing materials in the subcrustal lithosphere.
5
Rowbotham, Peter S., and Neil R. Goulty. "Wavefield separation by 3-D filtering in crosshole seismic reflection processing." GEOPHYSICS 59, no. 7 (July 1994): 1065–71. http://dx.doi.org/10.1190/1.1443662.
Full textAbstract:
In processing crosshole seismic reflection data, it is necessary to separate the upgoing and downgoing primary reflections from each other, from the direct waves, and from other wave types in the data. We have implemented a 3-D f-k-k filter for wavefield separation that is applied in a single pass. The complete data set is treated as a data volume, with each sample defined by the three coordinates of source depth, receiver depth, and time. The filter works well because upgoing primaries, downgoing primaries, and direct waves lie in different quadrants in f-k-k space. The strongest multiples, including mode‐converted multiples, lie in the same quadrants in f-k-k space as the direct waves, so they are readily rejected together. Tube waves and mode‐converted primaries are also suppressed as most of the energy in these wave types lies outside the pass volume for P‐wave primaries. Some head wave and S‐wave primary energy will be passed by the filter; however, these waves tend to have low amplitudes and late arrival times, respectively, and will be smeared out on imaging with the P‐wave velocity field. We have processed a real crosshole data set using two different methods of wavefield separation: applying 2-D f-k filtering to common source gathers and applying a 3-D f-k-k filter to the whole data set. The migrated image produced after 3-D f-k-k filtering contains less coherent noise and consequently shows improved continuity of reflectors.
6
Chen, Kai-Xun, Po-Fei Chen, Li-Wei Chen, Huajian Yao, Hongjian Fang, and Po-Li Su. "South Ilan Plain High-Resolution 3-D S-Wave Velocity from Ambient Noise Tomography." Terrestrial, Atmospheric and Oceanic Sciences 27, no. 3 (2016): 375. http://dx.doi.org/10.3319/tao.2016.01.29.02(tem).
Full text
7
Earp, S., A. Curtis, X. Zhang, and F. Hansteen. "Probabilistic neural network tomography across Grane field (North Sea) from surface wave dispersion data." Geophysical Journal International 223, no. 3 (August 8, 2020): 1741–57. http://dx.doi.org/10.1093/gji/ggaa328.
Full textAbstract:
SUMMARY Surface wave tomography uses measured dispersion properties of surface waves to infer the spatial distribution of subsurface properties such as shear wave velocities. These properties can be estimated vertically below any geographical location at which surface wave dispersion data are available. As the inversion is significantly non-linear, Monte Carlo methods are often used to invert dispersion curves for shear wave velocity profiles with depth to give a probabilistic solution. Such methods provide uncertainty information but are computationally expensive. Neural network (NN) based inversion provides a more efficient way to obtain probabilistic solutions when those solutions are required beneath many geographical locations. Unlike Monte Carlo methods, once a network has been trained it can be applied rapidly to perform any number of inversions. We train a class of NNs called mixture density networks (MDNs), to invert dispersion curves for shear wave velocity models and their non-linearized uncertainty. MDNs are able to produce fully probabilistic solutions in the form of weighted sums of multivariate analytic kernels such as Gaussians, and we show that including data uncertainties as additional inputs to the MDN gives substantially more reliable velocity estimates when data contains significant noise. The networks were applied to data from the Grane field in the Norwegian North sea to produce shear wave velocity maps at several depth levels. Post-training we obtained probabilistic velocity profiles with depth beneath 26 772 locations to produce a 3-D velocity model in 21 s on a standard desktop computer. This method is therefore ideally suited for rapid, repeated 3-D subsurface imaging and monitoring.
8
Frankel, Arthur. "Three-dimensional simulations of ground motions in the San Bernardino Valley, California, for hypothetical earthquakes on the San Andreas fault." Bulletin of the Seismological Society of America 83, no. 4 (August 1, 1993): 1020–41. http://dx.doi.org/10.1785/bssa0830041020.
Full textAbstract:
Abstract Three-dimensional finite difference simulations of elastic waves in the San Bernardino Valley were performed for two hypothetical earthquakes on the San Andreas fault: a point source with moment magnitude M5 and an extended rupture with M6.5. A method is presented for incorporating a source with arbitrary focal mechanism in the grid. Synthetics from the 3-D simulations are compared with those derived from 2-D (vertical cross section) and 1-D (flat-layered) models. The synthetic seismograms from the 3-D and 2-D simulations exhibit large surface waves produced by conversion of incident S waves at the edge of the basin. Seismograms from the flat-layered model do not contain these converted surface waves and underestimate the duration of shaking. The seismograms from the 3-D simulations have larger amplitude coda than do the seismograms from the 2-D case because of the presence of off-azimuth surface wave arrivals in the 3-D simulations that are not included in the 2-D simulations. Snapshots of the wavefield of the 3-D simulation show that these off-azimuth arrivals represent surface waves reflected from the edges of the basin. The anelastic attenuation of the sediments is a key parameter controlling the overall duration of motion. Some of the coda energy at rock sites near the basin edges represents leakage of surface wave energy out of the basin. For the M6.5 earthquake simulation, the largest ground velocities occur where surface waves reflected from the edge of the basin interfere constructively with the trapped waves that follow the direct S-wave. Maps of maximum ground velocity are produced for two directions of rupture propagation. The largest velocities occur in localized portions of the basin. The location of the largest velocities changes with the rupture propagation direction. Contours of maximum shaking are also dependent on asperity positions and radiation pattern.
9
Kato, Kenichi, Keiiti Aki, and Ta-Liang Teng. "3-D simulations of surface wave propagation in the Kanto sedimentary basin, Japan—part 1: Application of the surface wave Gaussian beam method." Bulletin of the Seismological Society of America 83, no. 6 (December 1, 1993): 1676–99. http://dx.doi.org/10.1785/bssa0830061676.
Full textAbstract:
Abstract The goal of this study is to simulate the displacement waveforms of 8 sec-period Love waves observed at Tokyo and Yokohama stations in the Kanto sedimentary basin during the Izu-hanto-toho-oki Earthquake of 29 June 1980 (M0 = 7 × 1025 dyne · cm). The surface wave Gaussian beam method is applied for this purpose. On the basis of the 3-D seismic velocity and density structure of the Kanto basin and assuming that the earth medium is laterally homogeneous outside the Kanto basin, waveforms of Love waves are synthesized. The synthesized seismograms underestimate the observed peak amplitudes at Yokohama station. This is primarily because the station is located in the direction of the nodal plane of the Love-wave radiation. As indicated by Yamanaka et al. (1992), a Quaternary basin exists in the Sagami Bay between the source location and the Kanto basin. We include the Sagami basin in our model by the approximation of a circular low-velocity region. Excellent agreement between observed and synthesized waveforms was achieved not only for amplitude but also for phase for the early parts of the wave trains at both stations. We conclude that the low velocity Quaternary basin in the Sagami Bay acts like a lens to focus surface wave energy resulting in high amplitudes. The later arriving waves, in particular the long duration observed at Tokyo station, however, cannot be adequately explained by this method. One possible reason for the failure of simulating the later phases is that this method disregards the secondary Love waves converted from S-waves and/or surface waves at a laterally discontinuous boundary. Although the surface wave Gaussian beam method cannot adequately predict the duration of observed seismograms, it provides us with a satisfactory prediction of amplitudes and phases for early arrivals in laterally slowly-varying media at drastically lower computation costs and less memory requirements than does other methods.
10
Pujol, Jose, and Richard Aster. "Joint hypocentral determination and the detection of low-velocity anomalies. An example from the Phlegraean Fields earthquakes." Bulletin of the Seismological Society of America 80, no. 1 (February 1, 1990): 129–39. http://dx.doi.org/10.1785/bssa0800010129.
Full textAbstract:
Abstract Arrival time data from the Phlegraean Fields (Italy) earthquake swarm recorded by the University of Wisconsin array in 1983 to 1984 were reanalyzed using a joint hypocentral determination (JHD) technique. The P- and S-wave station corrections computed as part of the JHD analysis show a circular pattern of central positive values surrounded by negative values whose magnitudes increase with distance from the center of the pattern. This center roughly coincides with the point of the maximum uplift (almost 2 m) associated with the swarm. Corrections range from −0.85 to 0.10 sec for P-wave arrivals and from −1.09 to 0.70 sec for S-wave arrivals. We interpret these patterns of corrections as caused by a localized low-velocity anomaly in the epicentral area, which agrees with the results of a previous 3-D velocity inversion of the same data set. The relocated (JHD) epicenters show less scatter than the epicenters obtained in the velocity inversion, and move more of the seismic activity to the vicinity of the only presently active fumarolic feature. The capability of the JHD technique to detect low-velocity anomalies and at the same time to give reliable locations, particularly epicenters, was verified using synthetic data generated for a 3-D velocity model roughly resembling the model obtained by velocity inversion.
11
Kaviani, Ayoub, Anne Paul, Ali Moradi, Paul Martin Mai, Simone Pilia, Lapo Boschi, Georg Rümpker, Yang Lu, Zheng Tang, and Eric Sandvol. "Crustal and uppermost mantle shear wave velocity structure beneath the Middle East from surface wave tomography." Geophysical Journal International 221, no. 2 (February 21, 2020): 1349–65. http://dx.doi.org/10.1093/gji/ggaa075.
Full textAbstract:
SUMMARY We have constructed a 3-D shear wave velocity (Vs) model for the crust and uppermost mantle beneath the Middle East using Rayleigh wave records obtained from ambient-noise cross-correlations and regional earthquakes. We combined one decade of data collected from 852 permanent and temporary broad-band stations in the region to calculate group-velocity dispersion curves. A compilation of >54 000 ray paths provides reliable group-velocity measurements for periods between 2 and 150 s. Path-averaged group velocities calculated at different periods were inverted for 2-D group-velocity maps. To overcome the problem of heterogeneous ray coverage, we used an adaptive grid parametrization for the group-velocity tomographic inversion. We then sample the period-dependent group-velocity field at each cell of a predefined grid to generate 1-D group-velocity dispersion curves, which are subsequently inverted for 1-D Vs models beneath each cell and combined to approximate the 3-D Vs structure of the area. The Vs model shows low velocities at shallow depths (5–10 km) beneath the Mesopotamian foredeep, South Caspian Basin, eastern Mediterranean and the Black Sea, in coincidence with deep sedimentary basins. Shallow high-velocity anomalies are observed in regions such as the Arabian Shield, Anatolian Plateau and Central Iran, which are dominated by widespread magmatic exposures. In the 10–20 km depth range, we find evidence for a band of high velocities (>4.0 km s–1) along the southern Red Sea and Arabian Shield, indicating the presence of upper mantle rocks. Our 3-D velocity model exhibits high velocities in the depth range of 30–50 km beneath western Arabia, eastern Mediterranean, Central Iranian Block, South Caspian Basin and the Black Sea, possibly indicating a relatively thin crust. In contrast, the Zagros mountain range, the Sanandaj-Sirjan metamorphic zone in western central Iran, the easternmost Anatolian plateau and Lesser Caucasus are characterized by low velocities at these depths. Some of these anomalies may be related to thick crustal roots that support the high topography of these regions. In the upper mantle depth range, high-velocity anomalies are obtained beneath the Arabian Platform, southern Zagros, Persian Gulf and the eastern Mediterranean, in contrast to low velocities beneath the Red Sea, Arabian Shield, Afar depression, eastern Turkey and Lut Block in eastern Iran. Our Vs model may be used as a new reference crustal model for the Middle East in a broad range of future studies.
12
Deschenes, Michael R., Clinton M. Wood, Liam M. Wotherspoon, Brendon A. Bradley, and Ethan Thomson. "Development of Deep Shear Wave Velocity Profiles in the Canterbury Plains, New Zealand." Earthquake Spectra 34, no. 3 (August 2018): 1065–89. http://dx.doi.org/10.1193/122717eqs267m.
Full textAbstract:
Deep (typically > 1,000 m) shear wave velocity ( V S) profiles were developed across the Canterbury region of New Zealand at nine strong-motion stations using a combination of active and passive surface wave methods. A multimode, multimethod joint inversion process, which included Rayleigh and Love wave dispersion and horizontal-to-vertical spectral ratio data, was used to develop the V S profiles at each site. A priori geologic information was used in defining preliminary constraints on the complex geologic layering of the deep basin underlying the region, including velocity reversals in locations where interbedded terrestrial gravels and marine sediments are present. Shear wave profiles developed as part of this study had characteristics comparable to the profiles from 14 Christchurch sites detailed in a separate study. The profiles developed in the two studies were combined to form region-specific V S profiles for typical deposits, which can be used to improve the accuracy of current three-dimensional (3-D) crustal velocity models of the region.
13
Block, Lisa V., C. H. Cheng, Michael C. Fehler, and W. Scott Phillips. "Seismic imaging using microearthquakes induced by hydraulic fracturing." GEOPHYSICS 59, no. 1 (January 1, 1994): 102–12. http://dx.doi.org/10.1190/geo1992-0156.
Full textAbstract:
Seismic imaging using microearthquakes induced by hydraulic fracturing produces a three-dimensional (3-D), S-wave velocity model of the fractured zone, improves the calculated locations of the microearthquakes, and may lead to better estimates of fractureplane orientations, fracture density, and water flow paths. Such information is important for predicting the amount of heat energy that may be extracted from geothermal reservoir. A fractured zone was created at the Los Alamos Hot Dry Rock Reservoir in north-central New Mexico within otherwise impermeable basement rock by injecting [Formula: see text] of water into a borehole under high pressure at a depth of 3.5 km. Induced microearthquakes were observed using four borehole seismometers. The P-wave and S-wave arrival times have been inverted to find the 3-D velocity structures and the microearthquake locations and origin times. The inversion was implemented using the separation of parameters technique, and constraints were incorporated to require smooth velocity structures and to restrict the velocities within the fractured region to be less than or equal to the velocities of the unfractured basement rock. The rms amval time residuals decrease by 11–15 percent during the joint hypocenter-velocity inversion. The average change in the microearthquake locations is 20–27 m, depending on the smoothing parameter used. Tests with synthetic data imply that the absolute locations may improve by as much as 35 percent, while the relative locations may improve by 40 percent. The general S-wave velocity patterns are reliable, but the absolute velocity values are not uniquely determined. However, studies of inversions using various degrees of smoothing suggest that the S-wave velocities decrease by at least 13 percent in the most intensely fractured regions of the reservoir. The P-wave velocities are poorly constrained because the P-wave traveltime perturbations caused by the fluid-filled fractures are small compared to the amval time noise level. The significant difference in the relative signal-to-noise levels of the P-wave and S-wave arrival time data, coupled with the limited ray coverage, can produce a bias in the computed [Formula: see text] ratios, and corresponding systematic rotation of the microearthquake cluster. These adverse effects were greatly reduced by applying a [Formula: see text] lower bound based on the [Formula: see text] ratio of the unfractured basement rock.
14
Schwartz, Susan Y., and Glenn D. Nelson. "Loma Prieta aftershock relocation with S-P travel times: Effects of 3-D structure and true error estimates." Bulletin of the Seismological Society of America 81, no. 5 (October 1, 1991): 1705–25. http://dx.doi.org/10.1785/bssa0810051705.
Full textAbstract:
Abstract Aftershocks of the 18 October 1989 Loma Prieta, California, earthquake are located using S-P arrival-time measurements from stations of the PASSCAL aftershock deployment. We demonstrate the effectiveness of using S-P arrival-time data in locating earthquakes recorded by a sparse three-component network. Events are located using the program QUAKE3D (Nelson and Vidale, 1990) with both 2-D and 3-D velocity models that have been developed independently for this region. The dense coverage of the area around the Loma Prieta rupture zone by instruments of the California Network (CALNET) has allowed the U.S. Geological Survey (USGS) to find P-wave earthquake locations for both velocity models, which we compare with our solutions. We also perform synthetic calculations to estimate realistic location errors resulting from uncertainties in both the 3-D velocity structure and the timing of arrivals. These calculations provide a comparison of location accuracies obtained using S-P arrival times, S and P arrival times, and P times alone. We estimate average absolute errors in epicentral location and in depth for the Loma Prieta aftershocks to be just under 2 km and 1 km, respectively, using S-P phase data and the sparse PASSCAL instrument coverage. The synthetic tests show that these errors are much smaller than those predicted using P-wave data alone and are nearly the same as those predicted using S- and P-phase data separately. This suggests that future aftershock recording deployments with sparse networks of three-component data can retrieve accurate event locations even if absolute timing is problematic. We find moderate differences between our locations and those determined by the USGS from a larger network of stations; however, common characteristics in both seismicity patterns are apparent. Neither set of locations yields earthquake patterns that can be easily interpreted in terms of simple faulting geometries. The absence of a simple pattern in both sets of earthquake locations indicates that this complexity is not the result of earthquake mislocation but is a genuine feature of the seismicity. A deep southwesterly dipping plane and a near-vertical fault extending from the surface to at least 7-km depth beneath the surface trace of the San Andreas Fault are imaged by both sets of earthquake locations. Although earthquake locations indicate the existence of many more fault segments, the complexity of this region requires that a definitive picture of the faulting geometry will have to await improvement in our knowledge of the P- and S-wave velocity structures.
15
Guo, B., C. H. Thurber, S. W. Roecker, J. Townend, C. Rawles, C. J. Chamberlain, C. M. Boese, S. Bannister, J. Feenstra, and J. D. Eccles. "3-D P- and S-wave velocity structure along the central Alpine Fault, South Island, New Zealand." Geophysical Journal International 209, no. 2 (February 16, 2017): 935–47. http://dx.doi.org/10.1093/gji/ggx059.
Full text
16
Quispe, Selene, Kosuke Chimoto, Hiroaki Yamanaka, Hernando Tavera, Fernando Lazares, and Zenon Aguilar. "Estimation of S-Wave Velocity Profiles at Lima City, Peru Using Microtremor Arrays." Journal of Disaster Research 9, no. 6 (December 1, 2014): 931–38. http://dx.doi.org/10.20965/jdr.2014.p0931.
Full textAbstract:
Microtremor exploration was performed around seismic recording stations at five sites in Lima city, Peru in order to know the site amplification at these sites. The Spatial Autocorrelation (SPAC) method was applied to determine the observed phase velocity dispersion curve, which was subsequently inverted in order to estimate the 1-D S-wave velocity structure. From these results, the theoretical amplification factor was calculated to evaluate the site effect at each site. S-wave velocity profiles at alluvial gravel sites have S-wave velocities ranging from ∼500 to ∼1500 m/s which gradually increase with depth, while Vs profiles at sites located on fine alluvial material such as sand and silt have Swave velocities that vary between ∼200 and ∼500 m/s. The site responses of all Vs profiles show relatively high amplification levels at frequencies larger than 3 Hz. The average transfer function was calculated to make a comparison with values within the existing amplification map of Lima city. These calculations agreed with the proposed site amplification ranges.
17
van Stiphout, A. M., S. Cottaar, and A. Deuss. "Receiver function mapping of mantle transition zone discontinuities beneath Alaska using scaled 3-D velocity corrections." Geophysical Journal International 219, no. 2 (August 2, 2019): 1432–46. http://dx.doi.org/10.1093/gji/ggz360.
Full textAbstract:
SUMMARY The mantle transition zone is the region between the globally observed major seismic velocity discontinuities around depths of 410 and 660 km and is important for determining the style of convection and mixing between the upper and the lower mantle. In this study, P-to-S converted waves, or receiver functions, are used to study these discontinuities beneath the Alaskan subduction zone, where the Pacific Plate subducts underneath the North American Plate. Previous tomographic models do not agree on the depth extent of the subducting slab, therefore improved imaging of the Earth structure underneath Alaska is required. We use 27 800 high quality radial receiver functions to make common conversion point stacks. Upper mantle velocity anomalies are accounted for by two recently published regional tomographic S-wave velocity models. Using these two tomographic models, we show that the discontinuity depths within our CCP stacks are highly dependent on the choice of velocity model, between which velocity anomaly magnitudes vary greatly. We design a quantitative test to show whether the anomalies in the velocity models are too strong or too weak, leading to over- or undercorrected discontinuity depths. We also show how this test can be used to rescale the 3-D velocity corrections in order to improve the discontinuity topography maps. After applying the appropriate corrections, we find a localized thicker mantle transition zone and an uplifted 410 discontinuity, which show that the slab has clearly penetrated into the mantle transition zone. Little topography is seen on the 660 discontinuity, indicating that the slab has probably not reached the lower mantle. In the southwest, P410s arrivals have very small amplitudes or no significant arrival at all. This could be caused by water or basalt in the subducting slab, reducing the strength at the 410, or by topography on the 410 discontinuity, preventing coherent stacking. In the southeast of Alaska, a thinner mantle transition zone is observed. This area corresponds to the location of a slab window, and thinning of the mantle transition zone may be caused by hot mantle upwellings.
18
Grechka, Vladimir, and Ilya Tsvankin. "3-D moveout inversion in azimuthally anisotropic media with lateral velocity variation: Theory and a case study." GEOPHYSICS 64, no. 4 (July 1999): 1202–18. http://dx.doi.org/10.1190/1.1444627.
Full textAbstract:
Reflection moveout recorded over an azimuthally anisotropic medium (e.g., caused by vertical or dipping fractures) varies with the azimuth of the source‐receiver line. Normal‐moveout (NMO) velocity, responsible for the reflection traveltimes on conventional‐length spreads, forms an elliptical curve in the horizontal plane. While this result remains valid in the presence of arbitrary anisotropy and heterogeneity, the inversion of the NMO ellipse for the medium parameters has been discussed so far only for horizontally homogeneous models above a horizontal or dipping reflector. Here, we develop an analytic moveout correction for weak lateral velocity variation in horizontally layered azimuthally anisotropic media. The correction term is proportional to the curvature of the zero‐offset traveltime surface at the common midpoint and, therefore, can be estimated from surface seismic data. After the influence of lateral velocity variation on the effective NMO ellipses has been stripped, the generalized Dix equation can be used to compute the interval ellipses and evaluate the magnitude of azimuthal anisotropy (measured by P-wave NMO velocity) within the layer of interest. This methodology was applied to a 3-D “wide‐azimuth” data set acquired over a fractured reservoir in the Powder River Basin, Wyoming. The processing sequence included 3-D semblance analysis (based on the elliptical NMO equation) for a grid of common‐midpoint “supergathers,” spatial smoothing of the effective NMO ellipses and zero‐offset traveltimes, correction for lateral velocity variation, and generalized Dix differentiation. Our estimates of depth‐varying fracture trends in the survey area, based on the interval P-wave NMO ellipses, are in good agreement with the results of outcrop and borehole measurements and the rotational analysis of four‐ component S-wave data.
19
Cruz-Hernández, Favio, Luis A. Gallardo, Marco Calò, Raúl R. Castro, and José M. Romo-Jones. "Ambient noise tomography in the Cerro Prieto Basin, Baja California, Mexico from laterally constrained surface wave inversion." Geophysical Journal International 229, no. 3 (January 19, 2022): 1586–603. http://dx.doi.org/10.1093/gji/ggac017.
Full textAbstract:
SUMMARY We apply a new technique for a laterally constrained surface wave inversion (LCSWI) to estimate the 3-D sedimentary structure of the Cerro Prieto Basin, Baja California, Mexico. The basin contains the Cerro Prieto geothermal field, which is considered one of the most productive in the world. The data used consist of group velocity measurements of Rayleigh waves extracted from cross-correlations of ambient noise registered at 12 stations distributed in an 18 × 12 km area. We estimated an S-wave velocity model that clearly shows three relatively homogeneous geoseismic units that correlate with the stratigraphic column reported in previous studies. The deepest geoseismic unit is the most heterogeneous and shows low-velocity zones likely associated with fluids. The resulting velocity model shows similarity with the conceptual geological model of the geothermal field reported in the literature and recent geophysical studies that suggest the potential existence of another deeper reservoir west of the current exploitation area.
20
Duan, Chenglong, David Lumley, and Hejun Zhu. "Estimation of micro-earthquake source locations based on full adjoint P and S wavefield imaging." Geophysical Journal International 226, no. 3 (June 3, 2021): 2116–44. http://dx.doi.org/10.1093/gji/ggab203.
Full textAbstract:
SUMMARY Locating micro-earthquakes with high resolution and accuracy is a challenge for traveltime inversion, which has uncertainty on the order of a Fresnel zone (many wavelengths). We develop a wave-equation imaging method to increase resolution and reduce location errors to less than a wavelength, but requires very densely deployed receiver arrays with wide aperture and considerable computational cost. Instead of using acoustic data or direct P wave arrivals only, we use elastic multicomponent data and present a new method that uses the full P and S adjoint wavefields to image the microseismic source locations. We separate the P and S waves from the data, and extrapolate the P and S wavefields of each receiver subarray by solving the P and S adjoint wave equations in parallel. We formulate three source imaging conditions by multiplying over subarrays the adjoint P wavefield (IP), S wavefield (IS) and cross-correlated P and S wavefields (IPS). We perform numerical experiments on the highly realistic SEG SEAM4D reservoir model using surface acquisition array geometries. Results for 2-D and 3-D microseismic source estimations show clean images without noisy artefacts at shallow depths. In particular, IPS provides the highest resolution source location image, while IP is limited by the P wavelength and IS is influenced by small coda artefacts. The major-axis alignment and resolution of the source location image are determined by the hypocentral location with respect to the receiver array and illumination-angle coverage, respectively. We discuss the impacts of S-wave attenuation and frequency bandwidth on the source location images. Noise tests indicate that the imaging results are relatively insensitive to ambient noise, as is observed for the surface monitoring data. Using smoothed velocity models, the imaging results are similar to the results using the true realistically heterogeneous velocity model. The 90 per cent confidence ellipse of the source location due to Gaussian-distributed velocity errors shows a larger depth error as the source becomes deeper, while the horizontal error does not change as much.
21
Tsuno, Seiji, Andi Muhamad Pramatadie, Yadab P. Dhakal, Kosuke Chimoto, Wakana Tsutsumi, and Hiroaki Yamanaka. "Long-Period Ground Motions Observed in the Northern Part of Kanto Basin, During the 2011 off the Pacific Coast of Tohoku Earthquake, Japan." Journal of Disaster Research 8, sp (September 1, 2013): 781–91. http://dx.doi.org/10.20965/jdr.2013.p0781.
Full textAbstract:
During the 2011 off the Pacific coast of Tohoku earthquake (Mw 9.0), strong ground motions were observed at many seismic stations in the Tokyo Metropolitan Area located about 200 km away from the southern edge of the earthquake source fault. Large earthquake responses in high-rise buildings having long natural periods of several seconds were also observed. The largest ground responses for a period of 4 to 5 seconds were observed locally in Oyama (K-NET TCG012) and Koga (K-NET IBR009) on the border between Tochigi and Ibaraki Prefectures in the northern part of Kanto basin. Geophysical information in these areas was not accurate enough, however, to evaluate these ground motions. To understand S-wave velocity structures, we performed array microtremors observations at TCG012 seismic station in Oyama. We applied the Spatial Autocorrelation (SPAC) method to array microtremors data for vertical components. Rayleigh wave phase velocity from 0.3 to 1.6 km/s was obtained for a period of 0.25 to 3 seconds. We inverted phase velocity to a S-wave velocity structure reaching to bedrock at a depth of 1.6 km, using a Genetic Algorithm. The estimated structure explained the first peak of the H/V spectral ratio of microtremors well by the ellipticity of fundamentalmode Rayleigh wave. To evaluate long-period ground motions observed around Oyama during the main shock, we estimated earthquake ground motions by 1-D analysis, showing agreements with and the differences from those observed. As a result, velocity calculated at IBR008 located midway between the Tsukuba Mountains and Oyama, explained that observed for main phases and later phases. However, velocity calculated at TCG012 did not explain that observed for later phases. According to the emphasis of airy phases for group velocity of Love wave using the estimated S-wave velocity structure and the principal axis for later phases obtained by PCA corresponding to the vibration direction of Love wave propagating from the earthquake source fault and through the Tsukuba Mountains, long-period ground motions of a period of 3 to 5 seconds observed at TCG012 lasting for 200 seconds after the arrival of main phases, consist of Love wave.
22
Hosseini, Kasra, Karin Sigloch, Maria Tsekhmistrenko, Afsaneh Zaheri, Tarje Nissen-Meyer, and Heiner Igel. "Global mantle structure from multifrequency tomography using P, PP and P-diffracted waves." Geophysical Journal International 220, no. 1 (September 17, 2019): 96–141. http://dx.doi.org/10.1093/gji/ggz394.
Full textAbstract:
SUMMARY In global-scale seismic tomography, teleseismic P and PP waves mainly constrain structures in the upper two thirds of the mantle, whereas core-diffracted waves (Pdiff) constrain the lower third. This study is the first to invert a very large data set of Pdiff waves, up to the highest possible frequencies. This results in tomographic resolution matching and exceeding that of global S-wave tomographies, which have long been the models of choice for interpreting lowermost mantle structure. We present three new global tomography models of 3-D isotropic P-wave velocity in the earth’s mantle. Multifrequency cross-correlation traveltimes are measured on all phases in passbands from 30 s dominant period to the highest frequencies that produce satisfactory fits (≈3 s). Model DETOX-P1 fits ≈2.5 M traveltimes from teleseismic P waves. DETOX-P2 fits the same data, plus novel measurements of ≈1.4 M traveltimes of Pdiff waves. DETOX-P3 fits the same data as DETOX-P2, plus ≈ 1.2 M PP traveltimes. Synthetics up to 1 s dominant period are computed by full wave propagation in a spherically symmetric earth using the spectral-element method AxiSEM. Traveltimes are linked to 3-D velocity perturbations (dVP/VP) by finite-frequency Fréchet kernels, parametrized on an adaptive tetrahedral grid of ≈400 000 vertices spaced by ≈80 km in the best-sampled regions. To complete spatial coverage, the waveform cross-correlation measurements are augmented by ≈5.7 million analyst-picked, teleseismic P arrival times. P, Pdiff and PP traveltimes are jointly inverted for 3-D isotropic P-velocity anomalies in the mantle and for events corrections, by least squares solution of an explicit matrix–vector equation. Inclusion of Pdiff traveltimes (in DETOX-P2, -P3) improves the spatial sampling of the lowermost mantle 100- to 1000-fold compared to teleseismic P waves (DETOX-P1). Below ≈2400 km depth, seismically slow anomalies are clustered at southern and equatorial latitudes, in a dozen or more intensely slow patches of 600–1400 km diameter. These features had long been classed into two large low shear velocity provinces (LLVP), which now appears questionable. Instead, patches of intensely slow anomalies in the lowermost mantle seem to form a nearly continuous, globe-spanning chain beneath the southern hemisphere, according to our increased resolution of LLVP-internal subdivisions and newly imaged patches beneath South America. Our tomography also supports the existence of whole-mantle plumes beneath Iceland, Ascension, Afar, Kerguelen, Canary, Azores, Easter, Galapagos, Hawaii, French Polynesia and the Marquesas. Seismically fast structure in the lowermost mantle is imaged as narrowly elongated belts under Eastern Asia and the Americas, presumably reflecting the palaeo-trench geometries of subduction zones and arcs that assembled Eastern Asia and the American Cordilleras in Palaeozoic and early Mesozoic times. Mid-mantle structure is primarily constrained by teleseismic P waves, but Pdiff data have a stabilizing effect, for example, sharpening the geometries of subducted slabs under the Americas, Eurasia and the Northern Pacific in the upper 2000 km. PP traveltimes contribute complementary constraints in the upper and mid mantle, but they also introduce low-velocity artefacts beneath the oceans, through downward smearing of lithospheric structure. Our three new global P-wave models can be accessed and interactively visualized through the SubMachine web portal (http://submachine.earth.ox.ac.uk/).
23
Civilini, F., W. D. Mooney, M. K. Savage, J. Townend, and H. Zahran. "Crustal imaging of northern Harrat Rahat, Saudi Arabia, from ambient noise tomography." Geophysical Journal International 219, no. 3 (August 19, 2019): 1532–49. http://dx.doi.org/10.1093/gji/ggz380.
Full textAbstract:
SUMMARY Harrat Rahat is a volcanic field located in west-central Saudi Arabia and is the site of the most recent eruption in the country (1256 CE). An earthquake swarm at a nearby volcanic field in 2009 prompted the need for new hazard models for this region, which includes the holy city of Medina. Tomography studies can be used to infer material properties of the subsurface such as partial melt, and are instrumental for volcanic hazard assessment. Regional earthquakes have been used to determine mantle structure, but such crustal models are often hindered by an insufficient number of earthquakes in the plate interior. We use ambient seismic noise to compute Rayleigh and Love surface-wave dispersion maps between 5 and 12 s for northern Harrat Rahat. The surface-wave maps are inverted to produce shear-wave velocities using a neighbourhood algorithm and interpolated into a pseudo-3-D model. The distributions of surface-wave and shear-wave velocities are heterogenous, varying between ±3 and 8 per cent. However, low velocities are not restricted to the Harrat. We observed a difference between Rayleigh- and Love-wave velocities that extends north from the site of the 1256 CE eruption and coincides with a low gravity anomaly. We obtain a shear-wave velocity increase of 10–15 per cent between 15 and 25 km depth consistent with the Conrad discontinuity, the interface between andesitic upper crust and the mafic lower crust of the Arabian Shield. The average velocities of the upper and lower crust are estimated to be 3.64 and 3.95 km s–1 using Rayleigh waves and 3.53 and 4.16 km s–1 using Love waves, which are in good agreement with the results of other geophysical studies of this area. The magnitude of the low-velocity anomalies, their location away from the Harrat, and the lack of reversals in the shear-velocity inversions suggest that the presence of a crustal magma chamber is not likely. If a magma chamber exists, it is smaller than can be imaged with a secondary microseism source (approximately 15 km wavelength), deeper than 30 km, or shallower than 5 km with a small velocity contrast.
24
Tsvankin, Ilya. "Reflection moveout and parameter estimation for horizontal transverse isotropy." GEOPHYSICS 62, no. 2 (March 1997): 614–29. http://dx.doi.org/10.1190/1.1444170.
Full textAbstract:
Transverse isotropy with a horizontal axis of symmetry (HTI) is the simplest azimuthally anisotropic model used to describe fractured reservoirs that contain parallel vertical cracks. Here, I present an exact equation for normal‐moveout (NMO) velocities from horizontal reflectors valid for pure modes in HTI media with any strength of anisotropy. The azimuthally dependent P‐wave NMO velocity, which can be obtained from 3-D surveys, is controlled by the principal direction of the anisotropy (crack orientation), the P‐wave vertical velocity, and an effective anisotropic parameter equivalent to Thomsen's coefficient δ. An important parameter of fracture systems that can be constrained by seismic data is the crack density, which is usually estimated through the shear‐wave splitting coefficient γ. The formalism developed here makes it possible to obtain the shear‐wave splitting parameter using the NMO velocities of P and shear waves from horizontal reflectors. Furthermore, γ can be estimated just from the P‐wave NMO velocity in the special case of the vanishing parameter ε, corresponding to thin cracks and negligible equant porosity. Also, P‐wave moveout alone is sufficient to constrain γ if either dipping events are available or the velocity in the symmetry direction is known. Determination of the splitting parameter from P‐wave data requires, however, an estimate of the ratio of the P‐to‐S vertical velocities (either of the split shear waves can be used). Velocities and polarizations in the vertical symmetry plane of HTI media, that contains the symmetry axis, are described by the known equations for vertical transverse isotropy (VTI). Time‐related 2-D P‐wave processing (NMO, DMO, time migration) in this plane is governed by the same two parameters (the NMO velocity from a horizontal reflector and coefficient ε) as in media with a vertical symmetry axis. The analogy between vertical and horizontal transverse isotropy makes it possible to introduce Thomsen parameters of the “equivalent” VTI model, which not only control the azimuthally dependent NMO velocity, but also can be used to reconstruct phase velocity and carry out seismic processing in off‐symmetry planes.
25
Corchete, V. "Crustal and upper mantle structure beneath the South China Sea and Indonesia." GSA Bulletin 133, no. 1-2 (May 28, 2020): 177–84. http://dx.doi.org/10.1130/b35641.1.
Full textAbstract:
Abstract A three-dimensional (3-D) S-velocity model for the crust and upper mantle beneath the South China Sea and Indonesia is presented, determined by means of Rayleigh wave analysis, in the depth range from 0 km to 400 km. The crustal and lithospheric mantle structure of this study area was previously investigated using several methods and databases. Due to their low resolution, a 3-D structure for this area has not been previously determined. The determination of such a 3-D S-velocity model is the goal of the present study. The most conspicuous features of the crust and upper mantle structure include the S-velocity difference between the Java Sea and the Banda Sea regions and a transitional boundary between these two regions. This model confirms the principal structural features revealed in previous studies: an oceanic crust structure in the center of the South China Sea, crustal thinning from the northern continental margin of the South China Sea to this oceanic crust, and the existence of a high-velocity layer in the lower crust of the northern continental margin. This study concludes that the north of the South China Sea is a nonvolcanic-type continental margin, solving the open question of whether the continental margin of the northern South China Sea is volcanic or nonvolcanic. A new map of the asthenosphere’s base is also presented.
26
Kassaras, I., V. Kapetanidis, A. Karakonstantis, and P. Papadimitriou. "Deep structure of the Hellenic lithosphere from teleseismic Rayleigh-wave tomography." Geophysical Journal International 221, no. 1 (January 8, 2020): 205–30. http://dx.doi.org/10.1093/gji/ggz579.
Full textAbstract:
SUMMARY This research provides new constraints on the intermediate depth upper-mantle structure of the Hellenic lithosphere using a three-step Rayleigh-wave tomography. Broadband waveforms of about 1000 teleseismic events, recorded by ∼200 permanent broad-band stations between 2010 and 2018 were acquired and processed. Through a multichannel cross-correlation technique, the fundamental mode Rayleigh-wave phase-velocity dispersion curves in the period range 30–90 s were derived. The phase-velocities were inverted and a 3-D shear velocity model was obtained down to the depth of 140 km. The applied method has provided 3-D constraints on large-scale characteristics of the lithosphere and the upper mantle of the Hellenic region. Highlighted resolved features include the continental and oceanic subducting slabs in the region, the result of convergence between Adria and Africa plates with the Aegean. The boundary between the oceanic and continental subduction is suggested to exist along a trench-perpendicular line that connects NW Peloponnese with N. Euboea, bridging the Hellenic Trench with the North Aegean Trough. No clear evidence for trench-perpendicular vertical slab tearing was resolved along the western part of Hellenic Subduction Zone; however, subcrustal seismicity observed along the inferred continental–oceanic subduction boundary indicates that such an implication should not be excluded. The 3-D shear velocity model supports an N–S vertical slab tear beneath SW Anatolia that justifies deepening, increase of dip and change of dip direction of the Wadati-Benioff Zone. Low velocities found at depths <50 km beneath the island and the backarc, interrelated with recent/remnant volcanism in the Aegean and W. Anatolia, are explained by convection from a shallow asthenosphere.
27
De Martin, F., E. Chaljub, P. Thierry, P. Sochala, F. Dupros, E. Maufroy, B. Hadri, A. Benaichouche, and F. Hollender. "Influential parameters on 3-D synthetic ground motions in a sedimentary basin derived from global sensitivity analysis." Geophysical Journal International 227, no. 3 (August 3, 2021): 1795–817. http://dx.doi.org/10.1093/gji/ggab304.
Full textAbstract:
SUMMARY Which physical parameters are the most influential when predicting earthquake ground motions in a 3-D sedimentary basin? We answer quantitatively by doing a global sensitivity analysis of two quantities of interest: the peak ground motions (PGMs) and a time–frequency representation (the S transform) of ground motions resulting from the synthetic anelastic responses of the EUROSEISTEST. This domain of interest is modeled by two layers with uncertain depth-dependent mechanical properties and is illuminated by a plane S-wave propagating vertically upward in an uncertain homogeneous elastic bedrock. The global sensitivity analysis is conducted on 800+ physics-based simulations of the EUROSEISTEST requiring 8+ million core-hours (i.e. ≈ 900 yr of mono-core computation). The analysis of the PGMs at the free surface displays the spatial influence of the uncertain input parameters over the entire basin scale, while the analysis of the time–frequency representation shows their influence at a specific location inside the basin. The global sensitivity analysis done on the PGMs points out that their most influential parameter in the middle of the basin is the quality factor QS (it controls up to 80 per cent of the PGMs in certain locations where the sediments thickness is larger than 200 m). On the other hand, the geological layering configuration (here represented by the depth of a geological interface controlling the geological layering) strongly influences the PGMs close to the basin edges, up to 90 per cent. We also found that the shear wave velocity at the free surface of the basin and the one of the bedrock underlying the basin are to be considered on an equal footing, both influencing the PGMs in the middle of the basin and close to its edges. We highlight that the bedrock to basin amplification of the PGMs shows a clear increase with respect to the thickness of the sediments, but this amplification saturates from 200 m of sediments around the value of three and is frequency dependent. This PGMs amplification starts from about one tenth of the mean S-wavelength propagating in the basin. The global sensitivity analysis done on the S transform of the ground motions shows that (i) the own effect of the parameters fully controls the first S-wave train and mostly controls the direct arrival of the basin-induced surfaces waves, (ii) the quality factor QS controls 40–60 per cent of the decay of amplitude of coda waves, the remaining part being mainly controlled by interaction effects due to the coupling effect of several parameters and (iii) the interaction effects between the parameters increases with time, suggesting under the hypotheses of our study that the own effects control the ballistic wave propagation while the interaction effects control the diffusive wave propagation.
28
Asgharzadeh, Mehdi, Pouya Ahmadi, Andrej Bóna, and Milovan Urosevic. "Orthotropic anisotropy analysis and parameter estimation from 3-D vertical seismic profile data." Geophysical Journal International 229, no. 2 (December 11, 2021): 1338–56. http://dx.doi.org/10.1093/gji/ggab502.
Full textAbstract:
SUMMARY Orthotropic velocity models are commonly considered as to have the lowest complexity (or highest symmetry) that can realistically describe the combined effect of Earth’s velocity heterogeneity in the vertical and in the horizontal (radial) directions. Velocity heterogeneity in the vertical direction is commonly caused by horizontal stratification of subsurface rock layers with different elastic properties or caused intrinsically by microlayered rocks such as shales. In the horizontal direction, azimuthal variation of the Earth’s stress field can initially lead to preferential deformation of pore space in the rocks and eventual development of aligned (vertical) fractures that are commonly regarded as the main cause for seismic azimuthal anisotropy. Building an appropriate orthotropic velocity model that can explain seismic velocity variations in vertical planes and along any azimuthal direction can therefore enhance the quality of seismic data processing and imaging. This, however, requires knowledge of several anisotropy parameters defined in the three symmetry planes of the orthotropic model and also the orientation of these planes. A combination of P- and S-wave slowness and polarization data derived from vertical seismic profile (VSP) measurements has been routinely used to estimate parameters of less complex anisotropic media such as transverse isotropy with vertical axes of symmetry (VTI). The accuracy of anisotropy parameter estimation by these methods depends primarily on the availability of data that are used in the inversion and the size of anisotropy parameters. In this study, we develop a benchmark for the accuracy of orthotropic parameters that are estimated from 3-D VSP measurements using P-wave slowness only method. Through numerical analysis, we model a wide range of fracture-induced orthotropic symmetries by adding vertical fractures to a background VTI media and show that the accuracy of our orthotropic parameter estimation method is dependent on anisotropy in the host media, intensity of fracturing, noise level in slowness data and slowness data coverage in both the dip and azimuthal planes. We then expand our accuracy analysis work to multilayer media by examining the accuracy of P-wave slowness method for orthotropic anisotropy parameter estimation by inverting finite-difference data generated in a pseudo 3-D VSP experiment. We show that, for a typical 3-D VSP experiment, reasonable estimates of orthotropic parameters in the model can only be obtained with slowness data taken from Xmax/Z ≥ 1 and line azimuth interval ϕ ≤ 10° where Xmax is the maximum source to receiver offset and Z is the receiver depth. We finally describe how field VSP measurements can be inverted for orthotropic anisotropy parameters.
29
Wright, Corwin J., Neil P. Hindley, Lars Hoffmann, M. Joan Alexander, and Nicholas J. Mitchell. "Exploring gravity wave characteristics in 3-D using a novel S-transform technique: AIRS/Aqua measurements over the Southern Andes and Drake Passage." Atmospheric Chemistry and Physics 17, no. 13 (July 13, 2017): 8553–75. http://dx.doi.org/10.5194/acp-17-8553-2017.
Full textAbstract:
Abstract. Gravity waves (GWs) transport momentum and energy in the atmosphere, exerting a profound influence on the global circulation. Accurately measuring them is thus vital both for understanding the atmosphere and for developing the next generation of weather forecasting and climate prediction models. However, it has proven very difficult to measure the full set of GW parameters from satellite measurements, which are the only suitable observations with global coverage. This is particularly critical at latitudes close to 60° S, where climate models significantly under-represent wave momentum fluxes. Here, we present a novel fully 3-D method for detecting and characterising GWs in the stratosphere. This method is based around a 3-D Stockwell transform, and can be applied retrospectively to existing observed data. This is the first scientific use of this spectral analysis technique. We apply our method to high-resolution 3-D atmospheric temperature data from AIRS/Aqua over the altitude range 20–60 km. Our method allows us to determine a wide range of parameters for each wave detected. These include amplitude, propagation direction, horizontal/vertical wavelength, height/direction-resolved momentum fluxes (MFs), and phase and group velocity vectors. The latter three have not previously been measured from an individual satellite instrument. We demonstrate this method over the region around the Southern Andes and Antarctic Peninsula, the largest known sources of GW MFs near the 60° S belt. Our analyses reveal the presence of strongly intermittent highly directionally focused GWs with very high momentum fluxes (∼ 80–100 mPa or more at 30 km altitude). These waves are closely associated with the mountains rather than the open ocean of the Drake Passage. Measured fluxes are directed orthogonal to both mountain ranges, consistent with an orographic source mechanism, and are largest in winter. Further, our measurements of wave group velocity vectors show clear observational evidence that these waves are strongly focused into the polar night wind jet, and thus may contribute significantly to the missing momentum at these latitudes. These results demonstrate the capabilities of our new method, which provides a powerful tool for delivering the observations required for the next generation of weather and climate models.
30
Hable, Sarah, Karin Sigloch, Eléonore Stutzmann, Sergey Kiselev, and Guilhem Barruol. "Tomography of crust and lithosphere in the western Indian Ocean from noise cross-correlations of land and ocean bottom seismometers." Geophysical Journal International 219, no. 2 (July 26, 2019): 924–44. http://dx.doi.org/10.1093/gji/ggz333.
Full textAbstract:
SUMMARY We use seismic noise cross-correlations to obtain a 3-D tomography model of SV-wave velocities beneath the western Indian Ocean, in the depth range of the oceanic crust and uppermost mantle. The study area covers 2000 × 2000 km2 between Madagascar and the three spreading ridges of the Indian Ocean, centred on the volcanic hotspot of La Réunion. We use seismograms from 38 ocean bottom seismometers (OBSs) deployed by the RHUM-RUM project and 10 island stations on La Réunion, Madagascar, Mauritius, Rodrigues, and Tromelin. Phase cross-correlations are calculated for 1119 OBS-to-OBS, land-to-OBS, and land-to-land station pairs, and a phase-weighted stacking algorithm yields robust group velocity measurements in the period range of 3–50 s. We demonstrate that OBS correlations across large interstation distances of >2000 km are of sufficiently high quality for large-scale tomography of ocean basins. Many OBSs yielded similarly good group velocity measurements as land stations. Besides Rayleigh waves, the noise correlations contain a low-velocity wave type propagating at 0.8–1.5 km s−1 over distances exceeding 1000 km, presumably Scholte waves travelling through seafloor sediments. The 100 highest-quality group velocity curves are selected for tomographic inversion at crustal and lithospheric depths. The inversion is executed jointly with a data set of longer-period, Rayleigh-wave phase and group velocity measurements from earthquakes, which had previously yielded a 3-D model of Indian Ocean lithosphere and asthenosphere. Robust resolution tests and plausible structural findings in the upper 30 km validate the use of noise-derived OBS correlations for adding crustal structure to earthquake-derived tomography of the oceanic mantle. Relative to crustal reference model CRUST1.0, our new shear-velocity model tends to enhance both slow and fast anomalies. It reveals slow anomalies at 20 km depth beneath La Réunion, Mauritius, Rodrigues Ridge, Madagascar Rise, and beneath the Central Indian spreading ridge. These structures can clearly be associated with increased crustal thickness and/or volcanic activity. Locally thickened crust beneath La Réunion and Mauritius is probably related to magmatic underplating by the hotspot. In addition, these islands are characterized by a thickened lithosphere that may reflect the depleted, dehydrated mantle regions from which the crustal melts where sourced. Our tomography model is available as electronic supplement.
31
Magrini, Fabrizio, Giovanni Diaferia, Islam Fadel, Fabio Cammarano, Mark van der Meijde, and Lapo Boschi. "3-D shear wave velocity model of the lithosphere below the Sardinia–Corsica continental block based on Rayleigh-wave phase velocities." Geophysical Journal International 220, no. 3 (December 24, 2019): 2119–30. http://dx.doi.org/10.1093/gji/ggz555.
Full textAbstract:
SUMMARY Rayleigh-wave dispersion curves from both ambient noise and teleseismic events allow us to provide the first high-resolution 3-D shear wave velocity (VS) model of the crust and upper mantle below the Sardinia–Corsica microplate, an important continental block for understanding the evolution of the central-western Mediterranean. For a wide range of periods (from 3 to ∼30 s), the phase velocities of the study area are systematically higher than those measured within the Italian peninsula, in agreement with a colder geotherm. Relative and absolute variations in the VS allow us to detect a very heterogeneous upper crust down to 8 km, as opposed to a relatively homogeneous middle and lower crust. The isosurface at 4.1 km s−1 is consistent with a rather flat Moho at a depth of 28.0 ± 1.8 km (2σ). The lithospheric mantle is relatively cold, and we constrain the thermal lithosphere–asthenosphere boundary at ∼100 km. We find our estimate consistent with a continental geotherm based on a surface heat flow of 60 mW m−2. Our results suggest that most of the lithosphere endured the complex history of deformation experienced by the study area and imply, in general, that deep tectonic processes do not easily destabilize the deeper portion of the continental lithosphere, despite leaving a clear surface signature.
32
Buland, Arild, and Henning Omre. "Bayesian linearized AVO inversion." GEOPHYSICS 68, no. 1 (January 2003): 185–98. http://dx.doi.org/10.1190/1.1543206.
Full textAbstract:
A new linearized AVO inversion technique is developed in a Bayesian framework. The objective is to obtain posterior distributions for P‐wave velocity, S‐wave velocity, and density. Distributions for other elastic parameters can also be assessed—for example, acoustic impedance, shear impedance, and P‐wave to S‐wave velocity ratio. The inversion algorithm is based on the convolutional model and a linearized weak contrast approximation of the Zoeppritz equation. The solution is represented by a Gaussian posterior distribution with explicit expressions for the posterior expectation and covariance; hence, exact prediction intervals for the inverted parameters can be computed under the specified model. The explicit analytical form of the posterior distribution provides a computationally fast inversion method. Tests on synthetic data show that all inverted parameters were almost perfectly retrieved when the noise approached zero. With realistic noise levels, acoustic impedance was the best determined parameter, while the inversion provided practically no information about the density. The inversion algorithm has also been tested on a real 3‐D data set from the Sleipner field. The results show good agreement with well logs, but the uncertainty is high.
33
Furumura, Takashi, and Brian L. N. Kennett. "Propagation of distinct Love-wave pulses from regional to teleseismic distances in continental and oceanic environments." Geophysical Journal International 221, no. 1 (January 16, 2020): 665–82. http://dx.doi.org/10.1093/gji/ggaa028.
Full textAbstract:
SUMMARY Surface waves are usually dispersive with long wave trains and steady decay of amplitude with distance. However, if the group velocity is nearly constant for a span of periods a strong pulse is produced that retains its amplitude for large distances. This situation arises for the fundamental mode of Love waves in the period band 40–500 s for crust and mantle structures with a positive gradient of S wave speed in the uppermost mantle. Such a distinct Love-wave pulse with limited dispersion observed at teleseismic distance is termed the G wave in honour of Gutenberg. The long-period G-wave pulse caused by large earthquakes carries a large amount of energy to substantial distances, with significant effects across the globe, for example event triggering. A similar G-type Love-wave pulse with a much shorter-period of 10–20 s is generated for crustal structures without thick sediment. Such pulses produce anomalously large ground displacement at near-regional distances with, for example an overestimate of surface wave magnitude. We investigate the generation and propagation mechanism of the G-type Love-wave pulses in the crust and upper-mantle with the analysis of observed strong motion records from the Mw 6.2 2016 Central Tottori earthquake and the Mw 9.0 2011 Off Tohoku earthquake in Japan, in conjunction with 3-D finite-difference simulation of seismic wave propagation and analysis of dispersion curves.
34
Kennett, B. L. N., M. G. Bostock, and J. K. Xie. "Guided-wave tracking in 3-D: A tool for interpreting complex regional seismograms." Bulletin of the Seismological Society of America 80, no. 3 (June 1, 1990): 633–42. http://dx.doi.org/10.1785/bssa0800030633.
Full textAbstract:
Abstract The tracking of guided-wave trains by raytracing in 3-D structures can provide a means of interpreting complex seismograms at regional distances. The method relies on the interpretation of Lg as the constructive interference of multiple S reflections within the crust, and strong scattering can be simulated by the inclusion of secondary sources. This approach has been applied to Californian events observed in the southwestern United States whose records at regional distances frequently exhibit an extended and complex Lg coda with significant late energy arriving with group velocities of 2.5 km/sec or less. The general character of this energy precludes an interpretation as the result of stochastic scattering processes or dispersion in low velocity surface sediments; therefore some alternative explanation needs to be sought. The observed wavetrains of events from two sites in southwestern California (Coalinga and North Baja) at stations LAC, MNV, and ELK have been compared with the predictions from tracking guided-wave patterns. Multipathing of energy by major changes in the thickness of the crustal wave guide is consistent with the extended and often “pulse-like” nature of the Lg coda. Examples of plausible multipathing mechanisms include reflection from the ocean-continent transition, scattering associated with topographic features such as the Sierra Nevada, and resonance within the narrow corridor of Baja California.
35
Wu, Tengfei, Shuangxi Zhang, Zijun Cao, Mengkui Li, Yujin Hua, Xiaoying Fu, and Yu Wei. "Lithospheric structure of Hubei Province, central China, from Rayleigh wave tomography: insight into the spatial contact relationship between the Yangtze Plate and the eastern Qinling-Dabie orogenic belt." Geophysical Journal International 221, no. 3 (March 3, 2020): 1669–83. http://dx.doi.org/10.1093/gji/ggaa102.
Full textAbstract:
SUMMARY Knowledge about the spatial contact relationship between the Yangtze Plate and the eastern Qinling-Dabie orogenic belt can not only provide a scientific basis for the exploration of mineral resources, disaster prevention and earthquake prediction in the region, but also play an important role in reconstructing the geological process of the central orogenic belt. Hence, high-resolution lithospheric velocity model is essential to address these geological problems. In this study, using waveform data from 48 stations in Hubei Province and adjacent regions, central China, we invert for a 3-D S-wave velocity structure model of the crust and upper mantle from Rayleigh wave tomography. Our model reveals the complex subduction pattern of the Yangtze Plate to the north and the thrust-nappe tectonics of the Qinling-Dabie orogenic belt along the Mianlue suture with different scales and different deformation strengths. In addition, in the central part of Hubei Province, the local Yangtze slab has been broken into several pieces, among which the upwelling low-velocity anomalies appear. Moreover, the southern margin of the Dabie orogenic belt has undergone thrusting-nappe movement, and a series of associated structures are formed in the northern margin of the middle Yangtze platform. The contact zone between the two blocks in this area is composed of a series of thrust faults with dextrorotation slip component. Finally, based on the 3-D S-wave velocity image of Hubei Province and its vertical cross-section profiles along three different directions, three dynamic models are proposed to explain the spatial contact relationship between the Yangtze Plate and the eastern Qinling-Dabie orogenic belt in different regions.
36
Liu, Xin, and Dapeng Zhao. "Seismic evidence for a plume-modified oceanic lithosphere–asthenosphere system beneath Cape Verde." Geophysical Journal International 225, no. 2 (January 11, 2021): 872–86. http://dx.doi.org/10.1093/gji/ggab012.
Full textAbstract:
SUMMARY We determine a new 3-D shear wave velocity (Vs) model down to 400 km depth beneath the Cape Verde hotspot that is far from plate boundaries. This Vs model is obtained by using a new method of jointly inverting P- and S-wave receiver functions, Rayleigh-wave phase-velocity data and S-wave arrival times of teleseismic events. Two Vs discontinuities at ∼15 and ∼60 km depths are revealed beneath volcanic islands, which are interpreted as the Moho discontinuity and the Gutenberg (G) discontinuity. Between the north and south islands, obvious high-Vs anomalies exist in the uppermost mantle down to a depth of ∼100–150 km beneath the Atlantic Ocean, whereas obvious low-Vs anomalies exist in the uppermost mantle beneath the volcanic islands including the active Fogo volcano. These low-Vs anomalies merge into a significant column-like low-Vs zone at depths of ∼150–400 km beneath the Cape Verde swell. We propose that these features in the upper mantle reflect a plume-modified oceanic lithosphere–asthenosphere system beneath the Cape Verde hotspot.
37
Lindh, Per, and Polina Lemenkova. "Shear bond and compressive strength of clay stabilised with lime/cement jet grouting and deep mixing: A case of Norvik, Nynäshamn." Nonlinear Engineering 11, no. 1 (January 1, 2022): 693–710. http://dx.doi.org/10.1515/nleng-2022-0269.
Full textAbstract:
Abstract The strength of soil can significantly increase by stabilisation with binders. Adding binders in correct proportions to improve soil parameters is of paramount importance for earthworks. In this article, we presented a framework to explore strength characteristics of soil stabilised by several binders and evaluated using applied geophysical methods by estimated P-wave velocities. The core of our work is a systematic assessment of the effects on clay stabilisation from various binders on shear and compressive strength. The binders were combined from four stabilising agents: (i) CEM II/A, a Portland limestone cement; (ii) burnt lime; (iii) lime kiln dust (LKD) limited up to 50%; and (iv) cement kiln dust (CKD). Shear strength has shown a nonlinear dependence as an exponential curve with P-waves. Natural frequency analysis was modelled to simulate resonant frequencies as eigen values. Variations in strength proved that CEM II/A-M (Recipe A, 100% CEM II) has the best performance for weak soil stabilisation followed by the combinations: Recipe B (70% CEM II/A-M, 30% LKD), Recipe C with added 80% CEM II/A-M and 20% CKD, and Recipe D (70% CEM II/A-M 30% CKD). Recipe B has shown high values with maximum uniaxial compressive strength (UCS) at 13.8 MPa. The Recipe C was less effective with the highest value of UCS as 8.8 MPa. The least strength was shown in Recipe D, where UCS has maximal values of 3.7 MPa. The specimens stabilised by Recipe B demonstrated the highest P-wave velocity at 2,350 m/s, while Recipe C and Recipe D showed the highest P-wave velocity at 1,900 and 1,550 m/s. All specimens shown a gain of UCS with sharply increased P-wave speed during the 3 days of curing. The study contributes to the development of methods of soil testing in civil engineering.
38
Cheng, Ningya, and Chuen Hon Cheng. "Estimations of formation velocity, permeability, and shear‐wave anisotropy using acoustic logs." GEOPHYSICS 61, no. 2 (March 1996): 437–43. http://dx.doi.org/10.1190/1.1443971.
Full textAbstract:
Field data sets collected by an array monopole acoustic logging tool and a shear wave logging tool are processed and interpreted. The P‐ and S‐wave velocities of the formation are determined by threshold detection with cross‐correlation correction from the full waveform and the shear‐wave log, respectively. The array monopole acoustic logging data are also processed using the extended Prony’s method to estimate the borehole Stoneley wave phase velocity and attenuation as a function of frequency. The well formation between depths of 2950 and 3150 ft (899 and 960 m) can be described as an isotropic elastic medium. The inverted [Formula: see text] from the Stoneley wave phase velocity is in excellent agreement with the shear‐wave log results in this section. The well formation between the depths of 3715 and 3780 ft (1132 and 1152 m) can be described as a porous medium with shear‐wave velocity anisotropy about 10% to 20% and with the symmetry axis perpendicular to the borehole axis. The disagreement between the shear‐wave velocity from the Stoneley wave inversion and the direct shear‐wave log velocity in this section is beyond the errors in the measurements. Estimated permeabilities from low‐frequency Stoneley wave velocity and attenuation data are in good agreement with the core measurements. Also it is proven that the formation permeability is not the cause of the discrepancy. From the estimated “shear/pseudo‐Rayleigh” phase velocities in the array monopole log and the 3-D finite‐difference synthetics in the anisotropic formation, the discrepancy can best be explained as shear‐wave anisotropy.
39
Lu, Y., and Y. Ben-Zion. "Validation of seismic velocity models in southern California with full-waveform simulations." Geophysical Journal International 229, no. 2 (January 11, 2022): 1232–54. http://dx.doi.org/10.1093/gji/ggab534.
Full textAbstract:
SUMMARY Crustal seismic velocity models provide essential information for many applications including earthquake source properties, simulations of ground motion and related derivative products. We present a systematic workflow for assessing the accuracy of velocity models with full-waveform simulations. The framework is applied to four regional seismic velocity models for southern California: CVM-H15.11, CVM-S4.26, CVM-S4.26.M01 that includes a shallow geotechnical layer, and the model of Berg et al. For each model, we perform 3-D viscoelastic wave propagation simulations for 48 virtual seismic noise sources (down to 2 s) and 44 moderate-magnitude earthquakes (down to 2 s generally and 0.5 s for some cases) assuming a minimum shear wave velocity of 200 m s–1. The synthetic waveforms are compared with observations associated with both earthquake records and noise cross-correlation data sets. We measure, at multiple period bands for well-isolated seismic phases, traveltime delays and normalized zero-lag cross-correlation coefficients between the synthetic and observed data. The obtained measurements are summarized using the mean absolute derivation of time delay and the mean correlation coefficient. These two metrics provide reliable statistical representations of model quality with consistent results in all data sets. In addition to assessing the overall (average) performance of different models in the entire study area, we examine spatial variations of the models’ quality. All examined models show good phase and waveform agreements for surface waves at periods longer than 5 s, and discrepancies at shorter periods reflecting small-scale heterogeneities and near-surface structures. The model performing best overall is CVM-S4.26.M01. The largest misfits for both body and surface waves are in basin structures and around large fault zones. Inaccuracies generated in these areas may affect tomography and model simulation results at other regions. The seismic velocity models for southern California can be improved by adding better resolved structural representations of the shallow crust and volumes around the main faults.
40
Wei, Yu, Shuangxi Zhang, Mengkui Li, Tengfei Wu, Yujin Hua, Yu Zhang, and Jianfeng Cai. "Regional lithospheric deformation beneath the East Qinling-Dabie orogenic belt based on ambient noise tomography." Geophysical Journal International 228, no. 2 (October 1, 2021): 1294–312. http://dx.doi.org/10.1093/gji/ggab393.
Full textAbstract:
SUMMARY The Qinling–Dabie orogenic belt, which contain the arc-shaped Dabbashan orocline and is the world's largest belt of HP/UHP metamorphic rocks, formed by a long-term complex amalgamation process between the North China Block and the Yangtze Block. To understand the collision processes and tectonic evolution, we constructed a 3-D S-wave velocity model from the surface to a depth of ∼120 km in the eastern Qinling-Dabie orogenic belt and its adjacent region by inverting 5–70 s phase velocity dispersion data of Rayleigh waves extracted from ambient noise data. Our 3-D model reveals low velocities in the middle–lower crust and high velocities in the upper mantle beneath the orogenic belt, suggesting the delamination of the lower crust. Our results support a two-stage exhumation model for the HP/UHP rocks in the study area. First-stage exhumation was caused by the slab breaking away from the subducted Yangtze Block during the Early–Middle Triassic. Partial melting of the lithospheric mantle caused by slab breakoff-related asthenospheric upwelling weakened the lithospheric mantle beneath the orogenic belt, and continued convergence of the two continental blocks led to further thickening of the lower crust. Such processes promoted lower-crust delamination, which triggered the second-stage exhumation of the HP/UHP rocks. In the Dabbashan orocline, two deep-rooted high-velocity domes, that is, Hannan–Micang and Shennong–Huangling domes, acted as a pair of indenters during the formation stage. High-velocity lower crust was observed beneath the Dabbashan orocline. In addition, our 3-D model reveals that high-velocity lithospheric mantle extends from the Sichuan Basin to the Dabbashan orocline, with a subhorizontal distribution, providing strong support for the high-velocity lower crust. We also observed the destruction of lithospheric mantle beneath the Yangtze Block; the destruction area is bounded by the North–South Gravity Lineament, suggesting that the destruction mechanism of the Yangtze Block may be similar to the North China Block.
41
Guerin, Gauthier, Diane Rivet, Anne Deschamps, Christophe Larroque, Aurélien Mordret, Jean-Xavier Dessa, and Xavier Martin. "High resolution ambient noise tomography of the Southwestern Alps and the Ligurian margin." Geophysical Journal International 220, no. 2 (October 21, 2019): 806–20. http://dx.doi.org/10.1093/gji/ggz477.
Full textAbstract:
SUMMARY The Southwestern Alps and the Ligurian margin is a region of moderate seismicity with a high rate of small to moderate events. Identifying the active faults in this very densely populated region is critical to better assess the hazard and mitigate the risk. An accurate 3-D velocity model of the shallow to middle crust is a fundamental step to better locate the seismicity, and hence, the faults from which it originates. We performed ambient noise surface-wave tomography based on all available continuous seismological data from the French and Italian permanent networks (RESIF, INGV, RSNI), and current and past temporary experiments (AlpArray, CASSAT, SISVAR, RISVAL). In addition to these available data, we deployed three more stations to improve the spatial resolution in a region with sparse seismic station coverage. Overall, we used 55 inland seismic stations, 5 oceans bottom seismometers and 2 offshore cabled site/sensors. Data span the 2014–2018 time period. Time series from all available components were cross-correlated to reconstruct both Rayleigh and Love-wave Green's functions. For each station-pair Rayleigh and Love group velocity dispersion curves were semi-automatically picked using a frequency–time analysis. Then we regionalize these group velocities to build 2-D Rayleigh and Love velocity-maps between 1.5 and 9 s period. Using a two-step inversion, we estimate the best 3-D shear wave velocity model. The first step is based on a Neighbourhood Algorithm to recover the best three layers’ velocity model at each cell of the model. We then use this three-layer model as a starting model in a perturbational method based on finite elements. At periods up to 5 s, the spatial variation of the velocity is well correlated with the effective geology of the area. Lower velocities are observed in areas where the sedimentary cover is thicker, such as the Var and Paillon valley near Nice, or in the subalpine domain in the northwestern part of the region. Higher velocities are retrieved in areas where massifs are present, such as the Argentera-Mercantour massifs in the northeastern, or the Esterel massif in the southwestern part of the region.
42
Lu, Y., L. Stehly, R. Brossier, and A. Paul. "Imaging Alpine crust using ambient noise wave-equation tomography." Geophysical Journal International 222, no. 1 (March 24, 2020): 69–85. http://dx.doi.org/10.1093/gji/ggaa145.
Full textAbstract:
SUMMARY We present an improved crustal Vs model and Moho depth map using ambient noise wave-equation tomography. The so-called ‘ambient noise wave-equation tomography’ is a method to invert seismic ambient noise phase dispersion data based on elastic waveform simulation, which accounts for 3-D and finite-frequency effects. We use cross-correlations of up to 4 yr of continuous vertical-component ambient seismic noise recordings from 304 high-quality broad-band stations in the Alpine region. We use model LSP_Eucrust1.0 obtained from traditional ambient noise tomography as initial model, and we iteratively improve the initial model by minimizing frequency-dependent phase traveltime differences between the observed and synthetic waveforms of Rayleigh waves in the period range 10–50 s. We obtain the final model after 15 iterations with ∼65 per cent total misfit reduction compared to the initial model. At crustal depth, the final model significantly enhances the amplitudes and adjusts the shapes of velocity anomalies. At Moho and upper-mantle depth, the final model corrects an obvious systematic velocity shift of the initial model. The resulting isovelocity Moho map confirms a Moho step along the external side of the external crystalline massifs of the northwestern Alps and reveals underplated gabbroic plutons in the lower most crust of the central and eastern Alps. Ambient noise wave-equation tomography turns out to be a useful tool to refine shear wave velocity models obtained by traditional ambient noise tomography based on ray theory.
43
Rosenberger, L. J. "Pectoral fin locomotion in batoid fishes: undulation versus oscillation." Journal of Experimental Biology 204, no. 2 (January 15, 2001): 379–94. http://dx.doi.org/10.1242/jeb.204.2.379.
Full textAbstract:
This study explores the dichotomy between undulatory (passing multiple waves down the fin or body) and oscillatory (flapping) locomotion by comparing the kinematics of pectoral fin locomotion in eight species of batoids (Dasyatis americana, D. sabina, D. say, D. violacea, Gymnura micrura, Raja eglanteria, Rhinobatos lentiginosus and Rhinoptera bonasus) that differ in their swimming behavior, phylogenetic position and lifestyle. The goals of this study are to describe and compare the pectoral fin locomotor behavior of the eight batoid species, to clarify how fin movements change with swimming speed for each species and to analyze critically the undulation/oscillation continuum proposed by Breder using batoids as an example. Kinematic data were recorded for each species over a range of swimming velocities (1–3 disc lengths s(−1)). The eight species in this study vary greatly in their swimming modes. Rhinobatos lentiginosus uses a combination of axial-based and pectoral-fin-based undulation to move forward through the water, with primary thrust generated by the tail. The pectoral fins are activated in short undulatory bursts for increasing swimming speed and for maneuvering. Raja eglanteria uses a combination of pectoral and pelvic locomotion, although only pectoral locomotion is analyzed here. The other six species use pectoral locomotion exclusively to propel themselves through the water. Dasyatis sabina and D. say have the most undulatory fins with an average of 1.3 waves per fin length, whereas Rhinoptera bonasus has the most oscillatory fin behavior with 0.4 waves per fin length. The remaining species range between these two extremes in the degree of undulation present on their fins. There is an apparent trade-off between fin-beat frequency and amplitude. Rhinoptera bonasus has the lowest frequency and the highest fin amplitude, whereas Rhinobatos lentiginosus has the highest frequency and the lowest amplitude among the eight species examined. The kinematic variables that batoids modify to change swimming velocity vary among different species. Rhinobatos lentiginosus increases its tail-beat frequency to increase swimming speed. In contrast, the four Dasyatis species increase swimming speed by increasing frequency and wavespeed, although D. americana also changes wave number. Raja eglanteria modifies its swimming velocity by changing wavespeed and wave number. Rhinoptera bonasus increases wavespeed, Gymnura micrura decreases wave number, and both Rhinoptera bonasus and Gymnura micrura increase fin-tip velocity to increase swimming velocity. Batoid species fall onto a continuum between undulation and oscillation on the basis of the number of waves present on the fins.
44
Papadopoulos, Ilias, Kafele Reddock, Jevan Manzano, and Joan L. Latchman. "The Trinidad and Tobago Microzonation Project: Port of Spain." Geophysical Journal International 222, no. 3 (June 3, 2020): 1936–51. http://dx.doi.org/10.1093/gji/ggaa275.
Full textAbstract:
SUMMARY In this study, we present the results from the microzonation study conducted in Port of Spain (PoS), capital of the Republic of Trinidad & Tobago. A dense grid of single-site recordings was used to determine the fundamental frequency of soil above bedrock, while a grid of 26 array recordings comprised the database for finding the 1-D shear wave velocity, with depth. The resonant frequency was found to range from <1.0 Hz, for the deeper sediments to the south, near the coast, to above 4.0 Hz, on the northern outskirts of the city, closer to the rock formations. The array data processing revealed a shear wave velocity less than 360 m s–1, for the alluvial deposits, whilst for the harder formations, the velocity was at least 1000 m s–1. To validate the results, a parametric investigation, using synthetic seismograms of ambient noise for simplified 1-D models of the PoS basin sediments, was conducted. A 3-D geological model of the basin was developed, by integrating the experimental results with the simulated data. The model suggests a gradual increase, from north to south, in sediment depth down to ∼160 m. In order to understand and explain the variation of the resonance frequency, a review of the historical development of the area, for the past 250 yr, revealed large-scale, non-engineered land reclamation in the 19th and 20th centuries, resulting in areas with anomalously high amplification of seismic motion.
45
Tran, Khiem T., Justin Sperry, Michael McVay, Scott J. Wasman, and David Horhota. "Shear Wave Velocity Profiles of Roadway Substructures from Multichannel Analysis of Surface Waves and Waveform Tomography." Transportation Research Record: Journal of the Transportation Research Board 2655, no. 1 (January 2017): 36–44. http://dx.doi.org/10.3141/2655-06.
Full textAbstract:
Assessment of roadway subsidence caused by embedded low-velocity anomalies is critical to the health and safety of the traveling public. Surface-based seismic techniques are often used to assess roadways because of data acquisition convenience and large depths of characterization. To mitigate the negative impact of closing a traffic lane under traditional seismic testing, a new test system that uses a land streamer is presented. The main advantages of the system are the elimination of the need to couple the geophones to the roadway, the use of only one source at the end of the geophone array, and the movement of the whole test system along the roadway quickly. For demonstration, experimental data were collected on asphalt pavement overlying a backfilled sinkhole that was experiencing further subsidence. For the study, a 24-channel land streamer and a propelled energy generator to generate seismic energy were used. The test system was pulled by a pickup truck along the roadway and the data were collected with 81 shots at every 3 m for a road segment of 277.5 m, with a total data acquisition time of about 1 h. The measured seismic data set was analyzed by the standard multichannel analysis of surface waves (MASW) and advanced two-dimensional (2-D) waveform tomography methods. Eighty-one one-dimensional shear wave velocity (VS) profiles from the MASW were combined to obtain a single 2-D profile. The waveform tomography method was able to characterize subsurface structures at a high resolution (1.5- × 1.5-m cells) along the test length to a depth of 22.5 m. Very low S-wave velocity was obtained at the repaired sinkhole location. The 2-D VS profiles from the MASW and waveform tomography methods are consistent. Both methods were able to delineate high- and low-velocity soil layers and variable bedrock.
46
Tong, P., D. Zhao, D. Yang, X. Yang, J. Chen, and Q. Liu. "Wave-equation based traveltime seismic tomography – Part 2: Application to the 1992 Landers earthquake (<i>M</i><sub>w</sub> 7.3) area." Solid Earth Discussions 6, no. 2 (August 25, 2014): 2567–613. http://dx.doi.org/10.5194/sed-6-2567-2014.
Full textAbstract:
Abstract. High-resolution 3-D P and S wave crustal velocity and Poisson's ratio models of the 1992 Landers earthquake (Mw 7.3) area are determined iteratively by a wave-equation based traveltime seismic tomography (WETST) technique as developed in the first paper. The details of data selection, synthetic arrival-time determination, and trade-off analysis of damping and smoothing parameters are presented to show the performance of this new tomographic inversion method. A total of 78 523 P wave and 46 999 S wave high-quality arrival-time data from 2041 local earthquakes recorded by 275 stations during the period of 1992–2013 is used to obtain the final tomographic models which costs around 10 000 CPU h. Checkerboard resolution tests are conducted to verify the reliability of inversion results for the chosen seismic data and the wave-equation based traveltime seismic tomography method. Significant structural heterogeneities are revealed in the crust of the 1992 Lander earthquake area which may be closely related to the local seismic activities. Strong variations of velocity and Poisson's ratio exist in the source regions of the Landers and three other strong earthquakes in this area. Most seismicity occurs in areas with high-velocity and low Poisson's ratio, which may be associated with the seismogenic layer. Pronounced low-velocity anomalies revealed in the lower crust along the Elsinore, the San Jacinto and the San Andreas faults may reflect the existence of fluids in the lower crust. The recovery of these strong heterogeneous structures are facilitated by the use of full wave equation solvers and WETST and verifies their ability in generating high-resolution tomographic models.
47
Tong, P., D. Zhao, D. Yang, X. Yang, J. Chen, and Q. Liu. "Wave-equation-based travel-time seismic tomography – Part 2: Application to the 1992 Landers earthquake (<i>M</i><sub>w</sub> 7.3) area." Solid Earth 5, no. 2 (November 26, 2014): 1169–88. http://dx.doi.org/10.5194/se-5-1169-2014.
Full textAbstract:
Abstract. High-resolution 3-D P and S wave crustal velocity and Poisson's ratio models of the 1992 Landers earthquake (Mw 7.3) area are determined iteratively by a wave-equation-based travel-time seismic tomography (WETST) technique. The details of data selection, synthetic arrival-time determination, and trade-off analysis of damping and smoothing parameters are presented to show the performance of this new tomographic inversion method. A total of 78 523 P wave and 46 999 S wave high-quality arrival-time data from 2041 local earthquakes recorded by 275 stations during the period of 1992–2013 are used to obtain the final tomographic models, which cost around 10 000 CPU hours. Checkerboard resolution tests are conducted to verify the reliability of inversion results for the chosen seismic data and the wave-equation-based travel-time seismic tomography method. Significant structural heterogeneities are revealed in the crust of the 1992 Landers earthquake area which may be closely related to the local seismic activities. Strong variations of velocity and Poisson's ratio exist in the source regions of the Landers and three other nearby strong earthquakes. Most seismicity occurs in areas with high-velocity and low Poisson's ratio, which may be associated with the seismogenic layer. Pronounced low-velocity anomalies revealed in the lower crust along the Elsinore, the San Jacinto, and the San Andreas faults may reflect the existence of fluids in the lower crust. The recovery of these strong heterogeneous structures is facilitated by the use of full wave equation solvers and WETST and verifies their ability in generating high-resolution tomographic models.
48
Gupta, Sandeep, P. Mahesh, Nagaraju Kanna, K. Sivaram, and Ajay Paul. "3-D seismic velocity structure of the Kumaun–Garhwal (Central) Himalaya: insight into the Main Himalayan Thrust and earthquake occurrence." Geophysical Journal International 229, no. 1 (October 30, 2021): 138–49. http://dx.doi.org/10.1093/gji/ggab449.
Full textAbstract:
SUMMARY Objective assessment of seismic hazard and understanding of the Himalayan arc's tectonics requires detailed information on the crustal structure and geometry of the underthrusting Indian Plate beneath the Himalaya. Here, we present high-resolution 3-D P-wave velocity (Vp) and P-to-S-wave velocity ratio (Vp/Vs) images of the Kumaun–Garhwal Himalaya, a proposed potential region for the future great earthquake. We generate these images by inverting arrival times of 515 local earthquakes recorded by 41 broad-band stations during November 2006–June 2008. The tomographic images show a heterogeneous structure in the upper-mid crust. These images, along with available geophysical and geological information, indicate the presence of quartz-rich felsic rocks in the uppermost crust and the occurrence of saline-rich aqueous fluid/partial melt in the upper-mid crust. We propose that the Main Himalayan Thrust (MHT), having a flat-ramp-flat geometry, lies at the base of these fluid zones. The small- and moderate-to-strong-magnitude earthquakes are mainly confined to the fluid-rich zones along the MHT and quartz-rich rocks in the upper crust. Such an interpretation implies that the earthquake occurrence in the Kumaun–Garhwal Himalaya is largely controlled by the geometry of the MHT and crustal lithology.
49
Feng, Mei, Meijian An, James Mechie, Wenjin Zhao, Guangqi Xue, and Heping Su. "Lithospheric structures of and tectonic implications for the central–east Tibetan plateau inferred from joint tomography of receiver functions and surface waves." Geophysical Journal International 223, no. 3 (August 26, 2020): 1688–707. http://dx.doi.org/10.1093/gji/ggaa403.
Full textAbstract:
SUMMARY We present an updated joint tomographic method to simultaneously invert receiver function waveforms and surface wave dispersions for a 3-D S-wave velocity (Vs) model. By applying this method to observations from ∼900 seismic stations and with a priori Moho constraints from previous studies, we construct a 3-D lithospheric S-wave velocity model and crustal-thickness map for the central–east Tibetan plateau. Data misfit/fitting shows that the inverted model can fit the receiver functions and surface wave dispersions reasonably well, and checkerboard tests show the model can retrieve major structural information. The results highlight several features. Within the plateau crustal thickness is >60 km and outwith the plateau it is ∼40 km. Obvious Moho offsets and lateral variations of crustal velocities exist beneath the eastern (Longmen Shan Fault), northern (central–east Kunlun Fault) and northeastern (east Kunlun Fault) boundaries of the plateau, but with decreasing intensity. Segmented high upper-mantle velocities have varied occurrences and depth extents from south/southwest to north/northeast in the plateau. A Z-shaped upper-mantle low-velocity channel, which was taken as Tibetan lithospheric mantle, reflecting deformable material lies along the northern and eastern periphery of the Tibetan plateau, seemingly separating two large high-velocity mantle areas that, respectively, correspond to the Indian and Asian lithospheres. Other small high-velocity mantle segments overlain by the Z-shaped channel are possibly remnants of cold microplates/slabs associated with subductions/collisions prior to the Indian–Eurasian collision during the accretion of the Tibetan region. By integrating the Vs structures with known tectonic information, we derive that the Indian slab generally underlies the plateau south of the Bangong–Nujiang suture in central Tibet and the Jinsha River suture in eastern Tibet and west of the Lanchangjiang suture in southeastern Tibet. The eastern, northern, northeastern and southeastern boundaries of the Tibetan plateau have undergone deformation with decreasing intensity. The weakly resisting northeast and southeast margins, bounded by a wider softer channel of uppermost mantle material, are two potential regions for plateau expansion in the future.
50
Shirzad, Taghi, Marcelo Assumpcao, and Marcelo Bianchi. "Ambient seismic noise tomography in west-central and Southern Brazil, characterizing the crustal structure of the Chaco-Paraná, Pantanal and Paraná basins." Geophysical Journal International 220, no. 3 (December 3, 2019): 2074–85. http://dx.doi.org/10.1093/gji/ggz548.
Full textAbstract:
SUMMARY Surface wave analysis provides important information on crustal structure, but it is challenging to obtain accurate/robust models in aseismic regions because of the lack of local earthquake records. In this paper, interstation empirical Green's functions retrieved by ambient seismic noise in 75 broad-band stations from 2016 January to 2018 September were used to study crustal structure in west-central Brazil. Fast marching method was applied to calculate the 2-D surface wave tomographic maps, and local dispersion curves were estimated in the period range of 4–80 s for each geographic cell. 1-D damped least squares inversion method was then conducted to obtained shear wave velocity model. Finally, the average ($\tilde{\rm V}$S) of the calculated VSV and VSH quasi 3-D models were used to characterize the crustal structure. Besides the checkerboard test resolution, a stochastic test with the effect of errors in the dispersion curves and choice of inversion parameters were carried out to better evaluate model uncertainties. Our results show a clear relation between the sedimentary thickness and geological units with the shorter period tomographic maps. Agreement has also been observed in longer periods such as the clear N–S anomaly along the Asuncion and Rio Grande Arches representing the boundary between the Chaco-Paraná and the Paraná basins. A 3-D composite velocity model shows a crustal structure consisting of three main layers. Some differences in lower crustal properties were found between the Paraná and Chaco-Paraná basins, consistent with a recently postulated, gravity-derived Western Paraná suture zone. However, no high velocities along the SW–NE axis of the Paraná basin were found to confirm proposed underplating. At the eastern edge of the Pantanal basin, the thin crust seems to be associated with a very thin (or lack of) lower crustal layer, consistent with a recently proposed crustal delamination hypothesis for the formation of the Pantanal basin.