Two-dimensional Velocity Model Research Articles

Implications on crustal structure from the South Finland Coastal (SOFIC) deep seismic sounding profile

We present results from a deep seismic sounding (DSS) experiment carried out along the southern coast of Finland in summer 2015. Data used in the survey derived from industrial blasts recorded by temporary project stations and permanent network stations. The western 220 km part of the 450 km long Southern Finland Coastal (SOFIC) profile runs along the Uusimaa belt (UB) in the 1.7–1.9 Ga Southern Finland subprovince (SFS) of the Svecofennian domain, while the 170 km part in the east crosses the 1.62–1.65 Ga Wiborg rapakivi batholith (WRB). The farthest 60 km cross a geologically diverse area consisting of supracrustal rocks and granitoids of the Saimaa area (SA), an eastern extension of the SFS. Our results show that the Moho boundary depth varies significantly, from ca. 52–54 km below UB to 40–45 km below WRB. All three crustal layers (upper, middle, and lower) have their maximum depth in the contact zone between UB and the WRB. Below WRB, a lower crust with Vp ~6.7–6.9 km/s is observed. High velocity lower crust was observed below UB (Vp ~7.2 km/s) and possibly below SA (Vp ~7.35 km/s). The modelling was based on ray tracing, using the extrapolation of seismic wave arrival times with the help of travel times predicted from a one-dimensional velocity model. The resulting two-dimensional velocity model partly relies on data from the intersecting DSS profiles and supports previous observations of the lithospheric structure of southeastern Fennoscandia.

A Comprehensive Study of Local, Global, and Combined Optimization Methods on Synthetic Seismic Refraction and Direct Current Resistivity Data

Most geophysical inversions face the problem of non-uniqueness, which poses a challenge in the mapping and delineation of the subsurface anomalies. To tackle this challenge, a combined local and global optimization approach is considered for jointly inverting two-dimensional direct current resistivity (DCR) and seismic refraction (SR) data that aim to estimate the corresponding physical model parameters. In this combined approach, the output of the local optimization method is used to determine the search space and tuning parameters for the global optimization algorithm. The multi-objective genetic algorithm (non-dominated sorting genetic algorithm) was utilized to jointly optimize the objective functions of two different methods. Because the genetic algorithm is a population-based optimization method, it requires numerous forward calculations. To deal with the expected high computational cost associated with this approach, parallel computing was utilized for the forward function evaluations to reduce the run time of the entire process. The proposed approach was tested using synthetic two-dimensional resistivity and velocity models that had three different types of anomalies (dyke, positive, and combined positive and negative). The results showed an improvement in the anomaly delineation in the output of the combined local and global optimization method compared with the local optimization method. Additionally, similar synthetic models were tested using only the single objective global optimization algorithm (conventional global optimization), which showed promising anomaly delineation.

Study on Flow Velocity Distribution in Open Channel With Flexible Vegetation

Ecological management of river channels is a hot topic for current sustainable development and flow measurement of ecological river is an important part. In this article, a flow velocity distribution model of the channel containing flexible vegetation is constructed from the vegetation riverbed theory and the bursting phenomenon to reveal the microscopic mechanism of the flow velocity distribution in the upper layer of vegetation. In the vegetation riverbed law, the effect of flexible vegetation is evaluated by the mixed length formula. The bursting phenomenon law considers the influence of the channel sidewalls on the flow and a two-dimensional velocity model is established by introducing the concept of average turbulence structure. The mechanism of the downward shift of the maximum flow velocity point on the channel sidewall is explained. The verification of the calculated velocity profiles is carried out based on data obtained in laboratory experiments. The results show that the combination of the two models can well describe the velocity distribution of the whole channel. At the end, the phenomenon of flow velocity zoning in open channel is discussed, which provides a solution for flow measurement in ecological channel.

Inversion of Love waves in earthquake ground motion records for two-dimensional S-wave velocity model of deep sedimentary layers

We propose a new waveform inversion method to estimate the 2D S-wave velocity structure of deep sedimentary layers using broadband Love waves. As a preprocessing operation in our inversion scheme, we decompose earthquake observation records into velocity waveforms for periods of 1 s each. Then, we include in the inversion only those periods for which the assumption of 2D propagation holds, which we propose to determine through a principal component analysis. A linearized iterative inversion analysis for the selected Love wave segments filtered for periods of 1 s each allows a detailed estimation of the boundary shapes of interfaces over the seismic bedrock with an S-wave velocity of approximately 3 km/s. We demonstrate the effectiveness of the technique with applications to observed seismograms in the Kanto Plain, Japan. The differences between the estimated and existing velocity structure models are remarkable at the basin edges. Our results show remarkable differences from previous existing structural models, particularly near the basin edges while being in good agreement with the surface geology. Since a subsurface structure at a basin edge strongly affects the earthquake ground motions in a basin with the generation of surface waves, our method can provide a detailed model of a complex S-wave velocity structure at an edge part for strong ground motion prediction.

Estimation of shallow subsurface S-wave velocity structure in urban area based on inversion of apparent dispersion curve

AbstractDue to its efficiency, convenience, non-destructive nature and strong anti-interference capability, the microtremor survey method (MSM) has found wide applications in urban geological surveys. The spatial autocorrelation method is diffusely applied to extract the dispersion curves from microtremor signals, which needs to satisfy the assumption that the energy of the fundamental Rayleigh wave is dominant. However, for layered media containing a layer with a significant low- or high-velocity contrast, this assumption is distinctly incorrect for certain frequency ranges. We present a processing methodology comprising the extraction and inversion of the apparent dispersion curves based on extended spatial autocorrelation method and fast simulated-annealing algorithm. We analyse synthetic microtremor signals for three selected geological models, and then compare the S-wave velocity structures estimated from their inversions with the actual models. Subsequently, a filed data example is given to detect the shallow stratigraphic structures in Guangzhou city, China, in which the new MSM was used. The estimated two-dimensional S-wave velocity model provided an accurate description of the thickness and depth of the strata in the study area, based on a priori information. Moreover, the S-wave velocity structures estimated from the MSM and the results from the drilling match very well, indicating that MSM is a reliable geophysical technique in urban geological surveys. Combined with available borehole information, MSM can be a very robust and effective method for detecting the shallow three-dimensional velocity structures in an urban area.

Locating Mine Microseismic Events in a 3D Velocity Model through the Gaussian Beam Reverse-Time Migration Technique.

Microseismic (MS) source location is a fundamental and critical task in mine MS monitoring. The traditional ray tracing-based location method can be easily affected by many factors, such as multi-ray path effects, waveform focusing and defocusing of wavefield propagation, and low picking precision of seismic phase arrival. By contrast, the Gaussian beam reverse-time migration (GBRTM) location method can effectively and correctly model the influences of multi-path effects and wavefield focusing and defocusing in complex 3D media, and it takes advantages of the maximum energy focusing point as the source location with the autocorrelation imaging condition, which drastically reduces the requirements of signal-to-noise ratio (SNR) and picking accuracy of P-wave arrival. The Gaussian beam technique has been successfully applied in locating natural earthquake events and hydraulic fracturing-induced MS events in one-dimensional (1D) or simple two-dimensional (2D) velocity models. The novelty of this study is that we attempted to introduce the GBRTM technique into a mine MS event location application and considered utilizing a high-resolution tomographic 3D velocity model for wavefield back propagation. Firstly, in the synthetic test, the GBRTM location results using the correct 2D velocity model and different homogeneous velocity models are compared to show the importance of velocity model accuracy. Then, it was applied and verified by eight location premeasured blasting events. The synthetic results show that the spectrum characteristics of the recorded blasting waveforms are more complicated than those generated by the ideal Ricker wavelet, which provides a pragmatic way to evaluate the effectiveness and robustness of the MS event location method. The GBRTM location method does not need a highly accurate picking of phase arrival, just a simple detection criterion that the first arrival waveform can meet the windowing requirements of wavefield back propagation, which is beneficial for highly accurate and automatic MS event location. The GBRTM location accuracy using an appropriate 3D velocity model is much higher than that of using a homogeneous or 1D velocity model, emphasizing that a high-resolution velocity model is very critical to the GBRTM location method. The average location error of the GBRTM location method for the eight blasting events is just 17.0 m, which is better than that of the ray tracing method using the same 3D velocity model (26.2 m).

Traveltime tomography and prestack depth migration for vertical seismic profiling of an angle-domain walkaway on a complex surface

Walkaway VSP cannot obtain accurate velocity field, as it asymmetrically reflects ray path and provides uneven coverage to underground target, thereby presenting issues related to imaging quality. In this study, we propose combining traveltime tomography and prestack depth migration for VSP of an angle-domain walkaway, in a bid to establish accurate two-dimensional and three-dimensional (3D) velocity models. First, residual curvature was defined to update velocity, and an accurate velocity field was established. To establish a high-precision velocity model, we deduced the relationship between the residual depth and traveltime of common imaging gathers (CIGs) in walkaway VSP. Solving renewal velocity using the least squares method, a four-parameter tomographic inversion equation was derived comprising formation dip angle, incidence angle, residual depth, and sensitivity matrix. In the angle domain, the reflected wave was divided into up- and down-transmitted waves and their traveltimes were calculated. The systematic cumulative method was employed in prestack depth migration of a complex surface. Through prestack depth migration, the offset-domain CIGs were obtained, and dip angle was established by defining the stack section horizon. Runge—Kutta ray tracing was employed to calculate the ray path from the reflection point to the detection point, to determine the incident angle, and to subsequently calculate the ray path from the reflection point to the irregular surface. The offset-domain residual depths were mapped to the angle domain, and a new tomographic equation was established and solved. Application in the double complex area of the Tarim Basin showed the four-parameter tomographic inversion equation derived in this paper to be both correct and practical and that the migration algorithm was able to adapt to the complex surface.

An efficient waveform inversion using the common mid-point gather in the wavenumber-space-time domain

As full waveform inversion (FWI) requires large computation time, a variety of techniques have been suggested to reduce the computational burden. In this study, we use wavenumber-space-time domain modelling, which directly generates common mid-point (CMP) gathers, to implement the FWI algorithm. The modelling technique, which is suitable for laterally invariant velocity models, synthesises CMP gathers efficiently by using limited wavenumber components, and therefore allows reduced computation time for FWI. To consider the accuracy as well as the efficiency of FWI, the Gauss-Newton method using the approximate Hessian matrix is chosen in this study. Rather than using all of the wavenumber components, our FWI algorithm can be accelerated by using only a few components. The wavenumber components can be chosen through an analysis of the residual wavefields. To validate the usefulness of our method, we first use a one-dimensional (1D) velocity model. From the 1D model example, we note that our FWI algorithm can be successful if given a reliable initial velocity model and sufficient data with a long offset distance. Even though our algorithm is valid for only horizontally layered velocity models, we apply our algorithm to a two-dimensional (2D) velocity model with lateral velocity variations. Through the 2D velocity example, we confirm that our FWI can be used to estimate subsurface structures with dipping interfaces if the dips are moderate and the structures can thus be considered to be locally flat.

A 2‐D tomographic model of the Juan de Fuca plate from accretion at axial seamount to subduction at the Cascadia margin from an active source ocean bottom seismometer survey

AbstractWe report results from a wide‐angle controlled source seismic experiment across the Juan de Fuca plate designed to investigate the evolution of the plate from accretion at the Juan de Fuca ridge to subduction at the Cascadia margin. A two‐dimensional velocity model of the crust and upper mantle is derived from a joint reflection‐refraction traveltime inversion. To interpret our tomography results, we first generate a plausible baseline velocity model, assuming a plate cooling model and realistic oceanic lithologies. We then use an effective medium theory to infer from our tomography results the extent of porosity, alteration, and water content that would be required to explain the departure from the baseline model. In crust of ages >1 Ma and away from propagator wakes and regions of faulting due to plate bending, we obtain estimates of upper crustal hydration of 0.5–2.1 wt % and find mostly dry lower crust and upper mantle. In sections of the crust affected by propagator wakes we find upper estimates of upper crustal, lower crustal, and upper mantle hydration of 3.1, 0.8, and 1.8 wt %, respectively. At the Cascadia deformation front, we find that the amount of water stored at uppermost mantle levels in the downgoing JdF plate is very limited (<0.3 wt %), with most of the water carried into the subduction zone being stored in the oceanic crust.

Long-range cortical connections give rise to a robust velocity map of V1

This paper proposes a two-dimensional velocity model (2DVM) of the primary visual cortex (V1). The model’s novel aspect is that it specifies a particular pattern of long-range cortical temporal connections, via the Connection Algorithm, and shows how the addition of these connections to well-known spatial properties of V1 transforms V1 into a velocity map. The map implies a number of organizational properties of V1: (1) the singularity of each orientation pinwheel contributes to the detection of slow-moving spots across the visual field; (2) the speed component of neuronal velocity selectivity decreases monotonically across each joint orientation contour line for parallel motion and increases monotonically for orthogonal motion; (3) the cells that are direction selective to slow-moving objects are situated in the periphery of V1; and (4) neurons in distinct pinwheels tend to be connected to neurons with similar tuning preferences in other pinwheels. The model accounts for various types of known illusionary perceptions of human vision: perceptual filling-in, illusionary orientation and visual crowding. The three distinguishing features of 2DVM are: (1) it unifies the functional properties of the conventional energy model of V1; (2) it directly relates the functional properties to the known structure of the upper layers of V1; and (3) it implies that the spatial selectivity features of V1 are side effects of its more important role as a velocity map of the visual field.

Ductile Gap between the Wenchuan and Lushan Earthquakes Revealed from the Two-dimensional Pg Seismic Tomography

A high-resolution two-dimensional Pg-wave velocity model is obtained for the upper crust around the epicenters of the April 20, 2013 Ms7.0 Lushan earthquake and the May 12, 2008 Ms8.0 Wenchuan earthquake, China. The tomographic inversion uses 47235 Pg arrival times from 6812 aftershocks recorded by 61 stations around the Lushan and Wenchuan earthquakes. Across the front Longmenshan fault near the Lushan earthquake, there exists a strong velocity contrast with higher velocities to the west and lower velocities to the east. Along the Longmenshan fault system, there exist two high velocity patches showing an “X” shape with an obtuse angle along the near northwest-southeast (NW-SE) direction. They correspond to the Precambrian Pengguan and Baoxing complexes on the surface but with a ~20 km shift, respectively. The aftershock gap of the 2008 Wenchuan and the 2013 Lushan earthquakes is associated with lower velocities. Based on the theory of maximum effective moment criterion, this suggests that the aftershock gap is weak and the ductile deformation is more likely to occur in the upper crust within the gap under the near NW-SE compression. Therefore our results suggest that the large earthquake may be hard to happen within the gap.

Study on the limitations of travel-time inversion applied to sub-basalt imaging

Abstract. The difficulties of seismic imaging beneath high velocity structures are widely recognised. In this setting, theoretical analysis of synthetic wide-angle seismic reflection data indicates that velocity models are not well constrained. A two-dimensional velocity model was built to simulate a simplified structural geometry given by a basaltic wedge placed within a sedimentary sequence. This model reproduces the geological setting in areas of special interest for the oil industry as the Faroe-Shetland Basin. A wide-angle synthetic dataset was calculated on this model using an elastic finite difference scheme. This dataset provided travel times for tomographic inversions. Results show that the original model can not be completely resolved without considering additional information. The resolution of nonlinear inversions lacks a functional mathematical relationship, therefore, statistical approaches are required. Stochastic tests based on Metropolis techniques support the need of additional information to properly resolve sub-basalt structures.

The Fossil Accretionary Wedge of the Bay of Biscay: Critical Wedge Analysis on Depth-Migrated Seismic Sections and Geodynamical Implications

AbstractThe accretionary wedge of the Bay of Biscay is an east-west compressive belt buried under recent sediments of the abyssal plain at the north Iberian margin. This structure formed through the partial closure of the previously extended Biscay basin during the Cenozoic north-south collision between Europe and Iberia, the same collision that produced the Cantabrian-Pyrenean range on land. Three north-south seismic sections have been prestack depth migrated, showing a narrow-tapered wedge (7°–8°) whose internal structure corresponds to a set of south-dipping thrusts converging toward a basal decollement. There are differences along strike within the wedge: thrust spacing, the dip of the basal thrust, and the thickness of the sediments at the trough augment toward the east, increasing its overall size. The two-dimensional velocity models obtained through migration analysis reflect values between 2000 km/s at the sea floor (4500 m) and 5000 km/s at 12-km depth. The syntectonic package thickness varies fr...

Seismic velocities on the Nova Scotian margin to estimate gas hydrate and free gas concentrations

This article provides new constraints on gas hydrate and free gas concentrations in the sediments at the margin off Nova Scotia. Two-dimensional (2-D) velocity models were constructed through simultaneous travel-time inversion of ocean-bottom seismometer (OBS) data and 2-D single-channel seismic (SCS) data acquired in two surveys, in 2004 and 2006. The surveys, separated by ∼5 km, were carried out in regions where the bottom-simulating reflection (BSR) was identified in seismic reflection datasets from earlier studies and address the question of whether the BSR is a good indicator of significant gas hydrate on the Scotian margin. For both datasets, velocity increases by 200–300 m/s at a depth of approximately 220 m below seafloor (mbsf), but the results of the 2006 survey show a smaller velocity decrease (50–80 m/s) at the base of this high-velocity layer (310–330 mbsf) than the results of the 2004 survey (130 m/s). When converted to gas hydrate concentrations using effective medium theory, the 2-D velocity models for both datasets show a gas hydrate layer of ∼100 m thickness above the identified BSR. Gas hydrate concentrations are estimated at approximately 2–10% for the 2006 data and 8–18% for the 2004 survey. The reduction in gas hydrate concentration relative to the distance from the Mohican Channel structure is most likely related to the low porosity within the mud-dominant sediment at the depth of the BSR. Free gas concentrations were calculated to be 1–2% of the sediment pore space for both datasets.

Crustal seismic structure and depth distribution of earthquakes in the Archean Kuusamo region, Fennoscandian Shield

Two-dimensional crustal velocity models are derived from passive seismic observations for the Archean Karelian bedrock of north-eastern Finland. In addition, an updated Moho depth map is constructed by integrating the results of this study with previous data sets. The structural models image a typical three-layer Archean crust, with thickness varying between 40 and 52 km. P wave velocities within the 12–20 km thick upper crust range from 6.1 to 6.4 km/s. The relatively high velocities are related to layered mafic intrusive and volcanic rocks. The middle crust is a fairly homogeneous layer associated with velocities of 6.5–6.8 km/s. The boundary between middle and lower crust is located at depths between 28 and 38 km. The thickness of the lower crust increases from 5–15 km in the Archean part to 15–22 km in the Archean–Proterozoic transition zone. In the lower crust and uppermost mantle, P wave velocities vary between 6.9–7.3 km/s and 7.9–8.2 km/s. The average V p/ V s ratio increases from 1.71 in the upper crust to 1.76 in the lower crust. The crust attains its maximum thickness in the south-east, where the Archean crust is both over- and underthrust by the Proterozoic crust. A crustal depression bulging out from that zone to the N–NE towards Kuusamo is linked to a collision between major Archean blocks. Further north, crustal thickening under the Salla and Kittilä greenstone belts is tentatively associated with a NW–SE-oriented collision zone or major shear zone. Elevated Moho beneath the Pudasjärvi block is primarily explained with rift-related extension and crustal thinning at ∼2.4–2.1 Ga. The new crustal velocity models and synthetic waveform modelling are used to outline the thickness of the seismogenic layer beneath the temporary Kuusamo seismic network. Lack of seismic activity within the mafic high-velocity body in the uppermost 8 km of crust and relative abundance of mid-crustal, i.e., 14–30 km deep earthquakes are characteristic features of the Kuusamo seismicity. The upper limit of seismicity is attributed to the excess of strong mafic material in the uppermost crust. Comparison with the rheological profiles of the lithosphere, calculated at nearby locations, indicates that the base of the seismogenic layer correlates best with the onset of brittle to ductile transition at about 30 km depth. We found no evidence on microearthquake activity in the lower crust beneath the Archean Karelian craton. However, a data set of relatively well-constrained events extracted from the regional earthquake catalogue implies a deeper cut-off depth for earthquakes in the Norrbotten tectonic province of northern Sweden.

Two-dimensional velocity model of the crust beneath the South Cambay Basin, India from refraction and wide-angle reflection data

SUMMARY A deep seismic sounding study was carried out in western India, in the southern part of the Cambay Basin, along a 230-km-long profile between the cities of Mehmadabad (22 ◦ 50 � N, 72 ◦ 46 � E) and Billimora (20 ◦ 46 � N, 72 ◦ 58 � E) during the year 1976–1977. Seismic refraction, as well as near- and wide-angle reflection, data were obtained from various boundaries from the sedimentary section to the Moho. Kaila et al. published models based on analogue versions of these data. In this study, we digitized, reprocessed, modelled, and reinterpreted the data from the Mehmadabad–Billimora seismic survey using modelling techniques that were not available at the time that Kaila et al. published their results. Our extensive modelling and interpretation of a large number of seismic record sections in south Cambay Basin reveal three Tertiary sedimentary subbasins: the Mehmadabad–Jambusar basin, the Jambusar–Kosamba basin, and the Kosamba–Billimora basin. These basins are vertically offset relative to the each other along high-angle faults, and the maximum depth (∼9.5 km) to the granitic/Proterozoic basement (5.8–6.0 km s −1 ) in the region occurs in the Jambusar–Broach region. Near-surface sedimentary rocks within the subbasins are underlain by Deccan Trap volcanics (4.75–5.2 km s −1 ), which have also been vertically offset throughout the study area. Although a prominent traveltime delay in first-arrivals, which is often indicative of low-velocity layers, is not observed on many of the seismograms, our modelling of secondary reflected arrivals indicate the presence of lower-velocity Mesozoic sediments below the Deccan Trap volcanics within the subbasins. The upper crust (6.0–6.3 km s −1 ), which is about 15 km thick, also appears to be stratified in places, with a low-velocity zone (5.5 km s −1 ) in the depth range of 12.5–14.5 km between the cities of Nadiad and Kathana. However, this low-velocity zone tapers and merges with the top of lower crust (6.5–6.8 km s −1 ) south of the city of Jambusar. The lower crust consists of two layers, an upper layer (6.5–6.8 km s −1 ) and a lowermost layer (7.2–7.5 km s −1 ), with the discontinuity between the two lower crustal layers occurring at 23–25 km depth. We interpret the higher velocity (7.2–7.5 km s −1 ) lowermost crustal layer as rocks that were magmatically underplated to the crust during rifting associated with the large-scale extrusion of the Deccan volcanics. The Moho discontinuity (8.0 km s −1 ) is about 31–33 km deep between the cities of Mehmadabad and Jambusar and increases to 37 km from south of the city of Ankleshwar. The large thickness of Tertiary sediments in the Cambay Basin and a relatively thin crust, as compared to the surrounding regions, suggests that rifting in the region occurred from the Cretaceous until the Tertiary.

Probable low-angle thrust earthquakes on the Juan de Fuca–North America plate boundary

In 2004, two clusters of earthquakes occurred in the central part of the Cascadia forearc, which displays several characteristics indicative of along-strike and downdip variations in plate coupling. Moment tensor analysis for the main shock in each cluster indicates that both events, which had magnitudes of 4.9 and 4.8, were compatible with low-angle thrust motion on fault planes dipping 6°–15° to the east, consistent with the plate boundary dip of ∼12°. By tracing rays through a high-resolution two-dimensional crustal velocity model to match arrival times of secondary arrivals (pP and PmP), we estimate that the source depth was 9–11 km for the M4.9 event and 15–17 km for the M4.8 event. We conclude that these earthquakes probably represent seismogenic thrust motion on the nominally locked or transitional part of the Cascadia megathrust, updip of where episodic tremor and slip (ETS) have been documented. Continued high-resolution observations of seismicity and improved crustal models are needed to confirm whether apparent temporal correlations between ETS and continued seismic activity in these clusters indicate stress transfer along the megathrust.

Quality control of compacted grounds using seismic velocities

ABSTRACTThe speed‐up plaque load experiment and a seismic refraction survey were used in quality‐control tests of compaction in the runway grounds of the Süleyman Demirel Airport, Isparta, located to the north of Burdur Lake, Turkey. The straight‐line correlation between plaque load experiments and seismic velocities ( and ) has shown that compaction is sufficient when both the and values exceed 590m/s and 379m/s, respectively. In addition, two‐dimensional velocity models for both P‐ and S‐waves map out the level of compaction along the profiles.

Geophysical evidence for a deep crustal root beneath the Yilgarn Craton and the Albany-Fraser Orogen

Using a marine seismic reflection survey as a source for onshore and offshore seismic refraction recorders, high quality refraction data were recorded to a maximum offset of 400 km. A two-dimensional velocity model for Line 19 of Geoscience Australia survey 280 in the Bremer Sub-basin and its onshore extension was derived by forward modelling using iterative ray tracing. Our velocity model indicates the presence of a deep crustal root located under the southern margin of the Yilgarn Craton and the Albany-Fraser Orogen. Comparison of the gravity response derived from this model with the gravity profile extracted from regional grids shows reasonable correlation. However, fine tuning of the seismic model is required to better define the complex geometry of the crustal root. Deep crustal roots are common under Proterozoic mobile belts that surround Archaean cratons. However, these results also indicate the presence of a crustal root under the southern margin of the Yilgarn Craton. This is consistent with the view that Proterozoic taphrogenesis penetrated tens of kilometres into the Yilgarn Craton.

Tau-p velocity imaging of regolith structure

The tau-p velocity imaging method, first developed for obtaining the velocity field from marine multichannel seismic data, has been applied to refracted waves from the regolith in regional seismic reflection surveys on land. The technique converts travel time picks from the refracted wavefield into two-dimensional velocity models, by transforming from time-offset into the tau-p domain. Each arrival is mapped individually, and the ?true? velocity and position of the ray turning point is obtained by considering reversed raypaths. Thus the data are transformed directly into a depth or two-way time image of the subsurface displayed in seismic velocity. The method is extremely fast and involves no interpretive steps or iteration. Ideal datasets contain the refracted wavefield sampled densely and equally in the shot and receiver domains. It was therefore decided to test the application of the method for mapping the velocity structure of a portion of regolith and to compare the results with those obtained using more conventional methods. The area chosen for study was part of a regional seismic reflection line across the Lachlan River palaeo-valley in central NSW. The data set consisted of the first break picks for 240 channels with receivers spaced every 40 m and vibration points every 40 m. The velocity images were produced as both time and depth sections and compared with one and two layer refractor models from a refraction tomographic approach. A low velocity region on the image corresponds to the deepest part of the refractor model, interpreted as the thickest part of the palaeo-valley. Bedrock velocity variations are also mapped and agree with the changes along the (lowest) refractor velocity profile. While further tuning may be required for land work, the technique has the advantage that velocities can be directly imaged and potentially related to regolith physical properties.