Three-dimensional (3D) three-component (3C)shallow seismic refraction surveys across a shearzone associated with dryland salinity at the Spicers Creek Catchment, New South Wales, Australia Author: Nikrouz, Ramin Publication Date: 2005 DOI: https://doi.org/10.26190/unsworks/21772 License: https://creativecommons.org/licenses/by-nc-nd/3.0/au/ Link to license to see what you are allowed to do with this resource. Downloaded from http://hdl.handle.net/1959.4/20607 in https:// unsworks.unsw.edu.au on 2022-07-21
This document is posted to help you gain knowledge. Please leave a comment to let me know what you think about it! Share it to your friends and learn new things together.
Transcript
Three-dimensional (3D) three-component (3C)shallow seismicrefraction surveys across a shearzone associated with drylandsalinity at the Spicers Creek Catchment, New South Wales,Australia
Author:Nikrouz, Ramin
Publication Date:2005
DOI:https://doi.org/10.26190/unsworks/21772
License:https://creativecommons.org/licenses/by-nc-nd/3.0/au/Link to license to see what you are allowed to do with this resource.
Downloaded from http://hdl.handle.net/1959.4/20607 in https://unsworks.unsw.edu.au on 2022-07-21
Three-Dimensional (3D) Three-Component (3C) Shallow Seismic Refraction Surveys Across a Shear Zone Associated with Dryland Salinity at the Spicers
Creek Catchment, New South Wales, Australia
by
Ramin Nikrouz
B.Sc., Geology, IRAN M.Sc., Petrology, IRAN
A Thesis Submitted in Fulfillment of the Requirements for the Degree of
Doctor of Philosophy
School of Biological, Earth and Environmental Sciences The University of New South Wales
Sydney NSW 2052 Australia
April 2005
PLEASE TYPE THE UNIVERSITY OF NEW SOUTH WALES
Thesis/Project Report Sheet Surname or Family name:
First name:
Other name/s:
Abbreviation for degree as given in the University calendar:
School:
Faculty:
Title:
Abstract 350 words maximum: (PLEASE TYPE)
Declaration relating to disposition of project report/thesis I am fully aware of the policy of the University relating to the retention and use of higher degree project reports and theses, namely that the University retains the copies submitted for examination and is free to allow them to be consulted or borrowed. Subject to the provisions of the Copyright Act 1968, the University may issue a project report or thesis in whole or in part, in photostat or microfilm or other copying medium. I also authorise the publication by University Microfilms of a 350 word abstract in Dissertation Abstracts International (applicable to doctorates only). …………………………………………………………… Signature
……………………………………..……………… Witness
……….……………………...…….… Date
The University recognises that there may be exceptional circumstances requiring restrictions on copying or conditions on use. Requests for restriction for a period of up to 2 years must be made in writing to the Registrar. Requests for a longer period of restriction may be considered in exceptional circumstances if accompanied by a letter of support from the Supervisor or Head of School. Such requests must be submitted with the thesis/project report. FOR OFFICE USE ONLY
Date of completion of requirements for Award:
Registrar and Deputy Principal
THIS SHEET IS TO BE GLUED TO THE INSIDE FRONT COVER OF THE THESIS
N:\FLORENCE\ABSTRACT
To my wife Fariba and my son Amir
CERTIFICATE OF ORIGINALITY
I hereby declare that this submission is my own work and to the best of my
knowledge it contains no materials previously published or written by another
person, nor material which to a substantial extent has been accepted for the
award of any other degree or diploma at UNSW or any other educational
institution, except where due acknowledgement is made in the thesis. Any
contribution made to the research by others, with whom I have worked at
UNSW or elsewhere, is explicitly acknowledged in the thesis.
I also declare that the intellectual content of this thesis is the product of my own
work, except to the extent that assistance from others in the project’s design
and conception or in style, presentation and linguistic expression is
surveys were recorded over a shear zone at two sites associated with dryland
salinity in the Spicers Creek Catchment, near Dubbo in southeastern Australia.
The seismic data were recorded with the Australian National Seismic Imaging
Resources (ANSIR) 360-trace ARAM-24 seismic system and IVI MiniVibrator.
Dryland salinity occurs extensively throughout the Spicers Creek Catchment.
The high concentration of salt in the groundwater has led to a significant decline
in agricultural productivity and a reduction in native vegetation. Furthermore,
the saline groundwater in the surface soil has caused the destruction of the clay
and soil structure and as a result, large areas of the catchment have
experienced soil erosion and extensive alteration of the landscape.
The broad objective of this study was to use seismic refraction methods to map
in detail a shear zone, which is associated with the salination. Seismic
refraction methods were selected because of their potential ability to provide
greater lateral resolution of the narrow vertical shear zone, than is currently the
norm with electrical or electromagnetic methods. This situation was confirmed
with a number of resistivity depth images generated as part of this study.
The results of the seismic refraction surveys show that the shear zone occurs
as a narrow region with low seismic velocities and increased depths of
weathering. However, there were numerous ambiguities in generating
i
consistent parameters in the refractor between adjacent receiver lines.
Furthermore, these ambiguities were not resolved with the use refraction
tomography, because even major differences with the various starting models
still generated acceptable agreement with the traveltime data. Therefore, in
order to generate consistency between the results for the various receiver lines,
it was necessary to generate detailed starting models.
A major achievement of this study is the development of a smoothing method
which significantly removes the effects of the near surface irregularities.
Termed the GRM SSM (for statics smoothing method), it essentially generates
a time-depth model of the refracting interface for which the effects of the near
surface irregularities have been minimized, by taking an average of the time-
depth profiles for a range of XY values. The GRM SSM, which takes advantage
of the unique redundancy properties of the GRM computations, was a major
factor in deriving consistent detailed starting models for refraction inversion.
Furthermore, this consistency was achieved with both S-wave components, as
well as the P-wave results.
The major geological achievement, which was made possible with the GRM
SSM, was the demonstration that there were cross cutting features associated
with the major shear zone. Therefore, it appears that saline groundwater can
discharge at the surface where increased volumes of groundwater occur at the
intersection of different sets of shears. This model provides a useful
explanation for the discontinuous occurrence of salination along the major shear
zone.
ii
ACKNOWLEDGEMENTS
First and foremost, I would like to thank God for helping me to accomplish this
project which is going to affect my academic life for some time to come. This
accomplishment would not have been possible without the help and support of
my supervisor, Derecke Palmer, who acted as a source of inspiration and
practical scientific support. I am most grateful for his keen supervision,
continuous support and encouragement during my research study. The
geophysical knowledge which I have acquired as the result of the
accomplishment of this research project is indeed the result of his insight and
expertise. I have been inspired by his endurance and perseverance.
I would also like to acknowledge my indebtedness to Paul Lennox, my co-
supervisor, who saved no effort to keep me on track in my endeavors. A
special word of thanks also goes to one of my best friends and colleagues in
Australia, Andrew Spyrow, whose friendship I value most highly. He was
especially helpful and supportive in the many computing tasks associated with
this study. I also owe a word of gratitude to David Johnston and the ANSIR
team who made it possible for me to gather what was indeed the corner stone
of my project, the data of the study.
iii
On a different note, I would like to thank the officials in the Iranian Ministry of
Science, Research and Technology and the University of Urmia of which I am a
faculty member for the financial support they provided for my family and me.
Also, I need to acknowledge my sense of gratitude and indebtedness to my
wife, Fariba and my son, Amir. They served as a source of encouragement and
support for me and were commendably patient with me although due to the
requirements of my studies, I had less time to spend with them. I am also
thankful to my parents of whose encouragement and support I am always
appreciative.
I should also take this opportunity to express my most sincere thanks to my
fellow research students at UNSW and friends, Morteza Jami and Ahmad Reza
Mokhtari, for their support and encouragement throughout the study of my
research.
Last but not the least, I sincerely thank the staff of the School of Biological,
Earth, and Environmental Sciences at UNSW, for their continued support and
encouragement. I am also grateful to many others, whose names do not
appear here. Their help is greatly appreciated.
iv
CONTENTS
Abstract-------------------------------------------------------------------------------------- i Acknowledgements---------------------------------------------------------------------- iii Contents-------------------------------------------------------------------------------------- V Illustrations---------------------------------------------------------------------------------- xi
Chapter 1 Introduction--------------------------------------------------------------------------------- 1 1.1 - Recent Advances in the Resolution of Geophysical Data-------------- -- 1 1.2 - Three-Dimentional Shallow Seismic Refraction Methods----------------- 5 1.3 - Dryland Salinity---------------------------------------------------------------------- 8 1.4 - Detailed Mapping of Shear Zones---------------------------------------------- 10 1.5 - Objectives----------------------------------------------------------------------------- 12 1.6 - Problems Solved in this Study--------------------------------------------------- 15
Chapter 2 Dryland Salinity----------------------------------------------------------------------------- 17 2.1 - Summary------------------------------------------------------------------------------- 17 2.2 - Introduction---------------------------------------------------------------------------- 18 2.3 - Extent of Dryland Salinity in New South Wales------------------------------- 19 2.4 - Causes of Dryland Salinity--------------------------------------------------------- 23 2.5 - Dryland Salinity in the Spicers Creek Catchment---------------------------- 28 2.6 - Management of Dryland Salinity-------------------------------------------------- 33
Chapter 3 Shear Zones Characteristic------------------------------------------------------------ 43 3.1 - Summary------------------------------------------------------------------------------- 43 3.2 - Introduction---------------------------------------------------------------------------- 44 3.3 - Shear Zone Types------------------------------------------------------------------- 46 3.4 - Fluid Movement in Shear Zones------------------------------------------------- 51
v
3.5 - Shear Zone Evidences---------------------------------------------------------- 54 3.5.1 - Field Evidence------------------------------------------------------------------- 54 3.5.2 - Imaging Shear Zones with Geophysical Methods---------------------- 59
Chapter 4 Three Dimension Three Component Seismic Surveys-------------------- 61 4.1 - Summary---------------------------------------------------------------------------- 61 4.2 - The Seismic Refraction Method----------------------------------------------- 62 4.3 - Seismic Velocities in the Earth------------------------------------------------ 63 4.4 - Three Dimension Seismic Refraction Methods---------------------------- 65 4.5 - Three Component Seismic Refraction Surveys---------------------------- 68 4.6 - Shear Waves Versus P Waves------------------------------------------------ 70
7.4 - Field Procedure-------------------------------------------------------------------- 183 7.4.1 - Setting up the Survey Spread-------------------------------------------- 184 7.4.2 - Conducting the Survey----------------------------------------------------- 186
7.4.2.1 - Site 1---------------------------------------------------------------------- 188 7.4.2.2 - Site 2---------------------------------------------------------------------- 192
7.5 - Data quality in Site 1 and Site 2------------------------------------------------ 196 7.6 – Magnetic Data Acquisition------------------------------------------------------ 198 7.7 – Resistivity Data Acquisition----------------------------------------------------- 198
Chapter 8 Processing--------------------------------------------------------------------------------- 199 8.1 - Summary----------------------------------------------------------------------------- 199 8.2 - Introduction--------------------------------------------------------------------------- 200 8.3 - Processing of Magnetic and Resistivity Data-------------------------------- 200
8.6 - Procedure of the Data Processing---------------------------------------------- 213 8.6.1 - Pre-processing----------------------------------------------------------------- 213 8.6.2 - Refraction Convolution Sections (RCS)--------------------------------- 214 8.6.3 - First Arrival Picking----------------------------------------------------------- 220 8.6.4 - Traveltime Graphs------------------------------------------------------------ 223 8.6.5 - The GRM Refractor Velocity Analysis Function----------------------- 228 8.6.6 - Corrections for Surface Effects-------------------------------------------- 234 8.6.7 - Seismic Velocities in the Refractor--------------------------------------- 237 8.6.8 - The GRM Time-Depth and Depth Function------------------------------240 8.6.9 - Processing of the First Arrival Amplitudes------------------------------ 244 8.6.10 - Velocity Ratios--------------------------------------------------------------- 248 8.6.11 - Three-Dimensional Images---------------------------------------------- 249 8.6.12 - Traveltime Tomography--------------------------------------------------- 249
Chapter 9 Interpretation------------------------------------------------------------------------------ 253 9.1 - Summary------------------------------------------------------------------------------ 253 9.2 - Introduction--------------------------------------------------------------------------- 255 9.3 - Magnetic Data Interpretation----------------------------------------------------- 256 9.4 - Resistivity Data Interpretation---------------------------------------------------- 257 9.5 - 3D-3C Seismic Data interpretation---------------------------------------------- 265
9.5.1 - Refractor Seismic Images at the First Site----------------------------- 265
ix
9.5.2 - Refractor Seismic Images at the Second Site------------------------- 274 9.6 - Geology and Tectonic interpretation------------------------------------------- 282 Chapter 10 Conclusions------------------------------------------------------------------------------ 285 References-------------------------------------------------------------------------------- 290 Appendices Appendix A - Linear and Non-Linear Metod--------------------------------------- 307 Appendix B - Operation Reports------------------------------------------------------ 309 Appendix C - Specification for the Equipment------------------------------------ 312 Appendix D - Station and Shot Point Coordinates------------------------------ 315 Appendix E - Geophone Coordinates----------------------------------------------- 317 Appendix F - Comment Document From Observer------------------------------ 318 Appendix G - Assessment of Data Quality----------------------------------------- 322 Appendix H - Seismic Unix Shell Scripts and C Files--------------------------- 336 Appendix I - Trace Order--------------------------------------------------------------- 348 Appendix J - Refraction Convolution Images------------------------------------- 350 Appendix K - Refractor Seismic Velocities---------------------------------------- 362 Appendix L - Time-Depth and Depth Graphs------------------------------------- 374 Appendix M - Eavepath Eikonal Traveltime Inversion-------------------------- 390 Appendix N - Traveltime Tomography Errors------------------------------------- 395 Appendix O - Refractor Images------------------------------------------------------ 396 Appendix P - GRM Tomography Images------------------------------------------ 410 Appendix Q - GRM Statics for Shallow Refraction Data----------------------- 416 Q1 – Abstract---------------------------------------------------------------- 416 Q2 – Generating Detailed Refractor Model for Inversion--------- 417 Q3 – Non-Uniqueness in Determining Detailed Refractor Velocities-------------------------------------------------------------- 422 Q4 – “Statics” Corrections for Shallow Refraction Data – Model Study -------------------------------------------------------- 427 Q5 – “Statics” Corrections for Shallow Refraction Data – Case Study----------------------------------------------------------- 431 Q6 – Conclusions----------------------------------------------------------- 434 Q7 – References------------------------------------------------------------ 437 Q8 – Figure Captions------------------------------------------------------ 441 Q9 - Figures ----------------------------------------------------------------- 445
x
ILLUSTRATIONS FIGURES Figure 1 - 1: Improvements in the resolution of geological detail have been achieved by adopting a new generation of survey specifications and presenting the resulting grid as an image. This illustration is from an area of the Ebagoola 1:250 000 sheet in northern Queensland. Top: BMR survey flown (E-W) in 1973 at a mean terrain clearance of 150 m and a line spacing of 1500 m, presented as a contour map with contour interval 10 nT. Bottom: the same area flown (E-W) by Geoterrex in 1991 with 100 m terrain clearance and 400 m line spacing (Reeves, C.V., 1992)…………………………………………………………………2 Figure 1 - 2: Results from the 3-D Mt Bulga seismic refraction survey. Seismic velocities and interpreted faults are plotted over the contours of the time-depths in milliseconds. The bold arrows indicate the directions of the higher seismic velocity (Palmer 2001b)………………………………………………………………6 Figure 1 - 3: The standard 2-D refraction depth cross section obtained across a shear zone at Mt Bulga. The shear zone (between stations 54 and 60) exhibits a decrease in the seismic velocity and an increase in the depth of weathering (Palmer, 2001b)……………………………………………………………………….6 Figure 2 - 1: Dryland salinity hazards in New South Wales (taken from Bradd, et al. 1997). The area that is within each hazard category is approximately low 75,060,000 ha; moderate 4,250,000 ha; high 470,000 ha; and very high 310,000 ha…………………………………………………………………………….21 Figure 2 - 2: Dryland salinity risk in New South Wales 2000 (From National Land & Water Resources Audit, 2000)…………………………………………….22 Figure 2 - 3: Dryland salinity risk in New South Wales 2050 (From National Land & Water Resources Audit; 2000)…………………………………………….22 Figure 2 - 4: Map of the Spicers Creek catchment with the main town, road and drainage modified from Cobbora 1:50,000 topographic map……………………28 Figure 3 - 1: (A) Fault. (B) Fault zone. (C) Shear zone (From Hobbs et al, 1976)………………………………………………………………………………….46 Figure 3 - 2: Different types of shear zones from (Ramsay and Auber, 1983)….. …………………………………………………………………………………………47 Figure 3 - 3: Idealized active shear zones of Type I and II (right) and stress / shear strain curves for each material (left). Shaded areas are current, actively straining regions, at margins (m) or in interior (i) of zones of Type I and II, respectively. Stubby arrows indicate direction of migration over time of boundaries of straining regions. Ellipses represent finite-strain gradients (From Means, 1995)…………………………………………………………………………49 Figure 3 - 4: Finite shear strain vs. time curves (solid) for the three types of shear zones, for rock at the final zone margins (m) and in the interior (i). P (dashed line) represents some process-controlling parameter (like temperature or fluid chemistry) that is assumed to change over the shearing history. In Type II zones, the marginal and interior acquire their low-strain state at about the
xi
same time and p value (heavy dots). The marginal rocks, therefore, provide a good record of former low-strain state of the interior rocks. In zones of Types I and III, the corresponding dots are farther apart in time and p, so the marginal rocks provide a less reliable record of the former low-strain state of the interior rocks. The type III diagram has been drawn for the case where the marginal rocks begin straining at the same time as the interior rocks, but more slowly. In each diagram, the right-hand dotted line is drawn at the point along the m line where its slope (the marginal strain-rate) becomes zero. The horizontal position of this line, therefore, represents the time at which rock in the shear zone margins acquired its final state of (low) strain. The left-hand line represents the time at which rock in the shear zone interior passed through this same state of (low) strain (From Means, 1995)…………………………………………………..50 Figure 3 - 5: Different types of criteria used for the determination of shear sense in shear zones (From Ramsay and Auber, 1983)……………………….58 Figure 4 - 1: Snell’s Law…………………………………………………………..62 Figure 4 - 2: Seismic refraction wave raypaths…………………………………64 Figure 4 -3: Nine-component sections recorded by P-, SH-, and SV-geophones for P-, SH-, and SV sources. Without anisotropy the sections marked ? would be blank (From Tatham and McCormack, 1991)………………………………..69 Figure 5 - 1: Tectonic units within the Lachlan Fold Belt area (Modified from: Kingham, 1998)…………………………………………………………………….76 Figure 5 - 2: Structural map of the Lachlan Fold Belt in New South Wales (Suppel and Scheibner, 1990)……………………………………………………78 Figure 5 - 3: Tectonic model for the eastern Lachlan Fold Belt (Modified from Collins and Vernon (1994))………………………………………………………..83 Figure 5 - 4: The Sydney, Gunnedah and Bowen Basin in NSW (Tadros, 1993). ………………………………………………………………………………………..85 Figure 5 - 5: Lithostratigraphy of the Gunnedah Basin (Tadros, 1993)………87
Figure 5-6: Surat Basin Stratigraphy (From: www.minerals.unsw.gov.au/geosurvey/petrol)…………………………90 Figure 5 - 7: Locality map for the Dubbo 1:250,000 geological map sheet area, showing the six component 1:100,000 map sheets and geographical features (Meakin and Morgan, 1999)………………………………………………………91 Figure 5 - 8: Mean monthly maximum and minimum temperatures, mean monthly precipitation data and the average number of wet days for Dubbo and Dunedoo, pan evaporation data for Wellington is shown at the base of Figure (Schofield, 1998)……………………………………………………………………92 Figure 5 - 9: Topography, township and surface drainage map for the Dubbo and Cobbora 1:100,000 sheets……………………………………………………94 Figure 5 - 10: Schematic structure map of the Dubbo1:250 000 sheet area (Matson, 1975)………………………………………………………………………96 Figure 5 - 11: Simplified geological framework of the Dubbo 1:250,000 map sheet area (Meakin and Morgan, 1999)………………………………………….100 Figure 5 - 12: Dubbo 1:250,000 map sheet area showing the metamorphic zones (Meakin and Morgan, 1999)……………………………………………….104 Figure 5 - 13: The road map of the Dubbo area………………………………..106
Figure 5 - 14: Geological map of study area from Cobbora 1:100,000 geological map with the two sites marked. The various units are described in Table 1…108 Figure 5 - 15: Structural zones defined by Scott (1997) across the Dubbo geological map sheet………………………………………………………………118 Figure 5 - 16: Structural map of the thrust sheet on the Cowra Zone and Molong Zone of the Dubbo geological sheet (Modified from Meakin and Morgan, 1999)…………………………………………………………………………………122 Figure 5 - 17: Location of the geophysical data sets supplied by AGSO and the selected area of interest for the interpretation of this data (Schofield, 1998)..125 Figure 5 - 18: The HSV colour model – H, S, and V represent hue, saturation, and value of colour intensity, respectively. According to this model, for a colour of constant hue and saturation, if its value is decreased it darkens towards black. For a colour of constant hue and value, if the saturation is decreased it becomes whiter (Milligan and Gunn, 1977)……………………………………..131 Figure 5 – 19: Red-green-blue composite radiometric image for Dubbo area (Meakin and Morgan, 1999)……………………………………………………….134 Figure 5 - 20: Gamma-ray image of an area south of Dubbo, showing the geochemically different signatures of a number of minor intrusions in the area. Some intrusions have dominant U+Th signatures, others dominant K+U or K+Th, and yet others are high in all three elements (Jaques et al., 1997)…..135
Figure 5 - 21: TMI aeromagnetic image for the Dubbo region ( Modified from Meakin and Morgan, 1999)……………………………………………………….141 Figure 5 - 22: Point located data for the Dubbo and Ballimore regions (Schofield, 1998)……………………………………………………………………142 Figure 5 - 23: Imaged gravity data for the Dubbo region (Modified from Meakin and Morgan, 1999)…………………………………………………………………146 Figure 5 - 24: Regional features that can be identified from a comparison of the gravity and magnetic images (Schofield, 1998)…………………………………147
Figure 6 - 1: Model with a plane horizontal ground surface and a highly irregular refractor. The vertical and horizontal scales are equal (Palmer, 1980)………153 Figure 6 - 2: Traveltime curves derived from the model in Figure 6 - 1……..153 Figure 6 - 3: Coincident velocity analysis functions for XY-values from 0 to 30 m (with same reciprocal time), derived from the traveltime data in Figure 6 - 2…. ……………………………………………………………………………………….155 Figure 6 - 4: Stacked velocity analysis functions for XY-values from 0 to 30 m (with different reciprocal time), derived from the traveltime data in Figure 6-2… ……………………………………………………………………………………….155 Figure 6 - 5: Average velocity analysis graph, derived from the traveltime curves in Figure 6 - 2………………………………………………………………156 Figure 6 - 6: Linear versus non-linear refractor velocities from Figure 6 - 2. 157 Figure 6 - 7: Time-depth for XY-values from 0 to 30 m……………………….158
xiii
Figure 6 - 8: Depth section which shows particularly for the depression at 225 m and for the fault at 50 m………………………………………………………158 Figure 6 - 9: Model with irregular ground and refractor surfaces. The vertical and horizontal scales are equal (Palmer, 1980)………………………………159 Figure 6 - 10: Traveltime curves derived from the model in Figure 6 - 9….160 Figure 6 - 11: Velocity analysis functions for XY-values from 0 to 30 m derived from the traveltime data in Figure 7-10. The data for a 15 m XY-value are judged the best……………………………………………………………………161 Figure 6 - 12: Average velocity analysis graph, derived from the traveltime curves in Figure 6 - 10……………………………………………………………162 Figure 6 - 13: Refractor seismic velocity graph derived from Figure 6 – 10.163 Figure 6 - 14: Stacked time-depth for XY-values from 0 to 30 m derived from the traveltime data in Figure 6 - 10………………………………………………163 Figure 6 - 15: Depth sections calculated from time-depth using 15-m XY-values……………………………………………………………………………….163 Figure 6 - 16: Traveltime curves derived from the model in Figure 7-10 after one correction………………………………………………………………………164 Figure 6 - 17: First raw traveltime correction……………………………………165 Figure 6 - 18: Velocity analysis functions for XY-values from 0 to 30 m derived from the traveltime data in Figure 6 - 10 after one correction…………………166 Figure 6 - 19: Average velocity analysis graph, derived from the traveltime curves in Figure 15, after one correction…………………………………………167 Figure 6 - 20: Refractor seismic velocity graph derived from Figure 6 - 16, after one correction………………………………………………………………………168 Figure 6 - 21: Stacked time-depth for XY-values from 0 to 30 m derived from the traveltime data in Figure 6 - 10 after 1 iteration……………………………168 Figure 6 - 22: Depth sections calculated from time-depth using 15-m XY-values after 1 iteration……………………………………………………………………..169 Figure 6 - 23: Traveltime curves derived from the model in Figure 6 – 10 after four corrections…………………………………………………………………….170 Figure 6 - 24: Fourth time-depth correction graph…………………………….171 Figure 6 - 25: Velocity analysis functions for XY-values from 0 to 30 m derived from the traveltime data in Figure 6 - 10 after four corrections……………….171 Figure 6 - 26: Average velocity analysis graph, derived from the traveltime curves in Figure 6 - 23, after four corrections……………………………………172 Figure 6 - 27: Refractor seismic velocity graph derived from Figure 6 - 23, after four corrections…………………………………………………………………….172 Figure 6 - 28: Stacked time-depth for XY-values from 0 to 30 m derived from the traveltime data in Figure 6 – 10 after four corrections…………………….173
xiv
Figure 6 - 29: Depth sections calculated from time-depth using 15-m XY-values after 4 iteration……………………………………………………………………173 Figure 7 - 1: P-wave versus S-wave comparison…………………………….179 Figure 7 - 2: Geophones and shot points layout at site 1……………………191 Figure 7 - 3: Geophones and shot points layout at site 2……………………195 Figure 7 - 4: P-wave shot record………………………………………………..197 Figure 8 - 1: Arrangement of the blocks used in a model together with the datum points in the pseudosection that generate by RES2DINV software…202 Figure 8 - 2: An illustration of the ray paths of the velocity analysis function. Ray paths are shown as solid lines when they are added, broken lines when they are subtracted and as broken-dashed lines when they cancel (Palmer, 2003b)………………………………………………………………………………209 Figure 8 - 3: A schematic summary of the ray paths used in the computation of the time-depth function for GRM (Palmer, 2003b)…………………………….210 Figure 8 - 4: The process of refraction migration. The offset distance (XG and YG) is the horizontal distances between the point of emergence of the ray on the refracting interface and the point of detection at the surface (top). Inversion methods which do not explicitly recognise the offset distance suffer refractor smoothing. The GRM accommodates the offset distance by calculating different values of horizontal distance between forward and reverse receivers, then selecting the XY value for which the forward and reverse rays leave the refractor at the same point and arrive at different detector positions at the surface (bottom)……………………………………………………………………212 Figure 8 - 5: Vertical shot record. The receivers on lines 1, 2, 3 and 4 are numbered 100-128, 200-228, 300-328 and 400-428 respectively. The four receiver lines are wrapped continuously in a ‘W’ shape, where the receivers for lines 2 and 4 are in the reverse order……………………………………………215 Figure 8 - 6: The shot record in Figure 1 split into individual receiver lines. The receivers on lines 1, 2, 3 and 4 are numbered 100-128, 200-228, 300-328 and 400-428 respectively. The trace order of lines 2 and 4 have been reversed..217 Figure 8 - 7: Refraction convolution section for a vertical (left) and orthogonal horizontal (right) component showing the refractor structure and amplitudes proportional to the square of the head coefficients………………………… …218 Figure 8 - 8: The enhancement of S-wave in horizontal records using records of opposite polarity. Subtraction results in the sum of the S-wave components and cancellation of P-wave components (bottom left) whilst addition results in the opposite (bottom right)…………………………………………………………….219 Figure 8 - 9: Muting and killing traces can improve the quality of first-arrival picking. The original cross line horizontal component shot record with considerable noise and dead traces (left). The same shot record after killing and muting above the first-arrivals (right)…………………………………………….220
xv
Figure 8 - 10: The picked first-arrival traveltimes overlayed on the shot records. ………………………………………………………………………………………222 Figure 8 - 11: Survey geometry with the source line offset from the receiver line. For a source point s and receiver r, when X is less than Y or when they are comparable there is a significant distortion of the traveltime curve at that point. When X is much greater than Y the distortion is negligible…………………..224 Figure 8 - 12: Traveltime graph for shot points inside the receiver spread. Since the shot points are offset from the receiver line, direct arrival traveltimes are shifted up in time and have a concave upward shape………………………..225 Figure 8 - 13: Corrected traveltime graphs for shot points inside the receiver spread. By geometric proportioning the traveltimes are changed so that the first arrival at each shot point occurs at zero time………………………………….226 Figure 8 - 14: Traveltime graph for P-wave (site1, line1) before correction..226 Figure 8 - 15: Corrected traveltime for P-wave (site1, line1) after correction.227 Figure 8 - 16: Using phantoming to combine traveltime graphs. Each traveltime curve is shifted downwards by the average separation between itself and the base traveltime graph. The average of these graphs is then taken to obtain the combined traveltime graph……………………………………………………….228 Figure 8 - 17: Forward and reverse velocity analysis graph for line1, site1, P-wave…………………………………………………………………………………229 Figure 8 - 18: The velocity analysis function calculated for different XY values and stacked. The optimum XY value can be determined by observing for which XY value the others are symmetric about, in this case a XY value of about 5 m. ………………………………………………………………………………………230 Figure 8 - 19: Forward and reverse velocity analysis graph for XY=5, P-wave, line1, site1………………………………………………………………………….231 Figure 8 - 20: Velocity analysis residuals are obtained by subtracting the average of the velocity analysis functions from the individual velocity analysis functions. Variations between the velocity analysis functions for the different XY values are generally less than 1.5 ms…………………………………………..232 Figure 8 - 21: Velocity analysis residual for P- and S-waves, site1…………233 Figure 8 - 22: Velocity analysis graph that show near surface effect……….234 Figure 8 - 23: Velocity analysis graph before correction. Notice to the line trends for V2 and V3………………………………………………………………236 Figure 8 - 24: Velocity analysis graph after correction. Notice to the line trends for V2 and V3………………………………………………………………………237 Figure 8 - 25: Lateral variations in S-wave velocities (in m/s) within the refractor…………………………………………………………………………….239 Figure 8 - 26: Traveltime of a refracted arrival at the shot point (S) (Palmer, 2003b)………………………………………………………………………………241
xvi
Figure 8 - 27: The time-depth function calculated for different XY values and stacked. The optimum XY value can be determined by observing for which XY value small irregularities in the time-depth function are most pronounced, in this case a XY value of about 5 m…………………………………………………….242 Figure 8 - 28: time-depth residuals are obtained by subtracting the average of the time-depth functions from the individual time-depth functions. Variations between the time-depth functions for the different XY values are generally less than 2 ms……………………………………………………………………………242 Figure 8 - 29: Forward and reverse distance P-wave amplitudes for shots 5 & 14 with calculating 1/X, 1/X^2 and 1/X^3. The large geometric spreading component dominates. The effects of geometric spreading have been reduced. ……………………………………………………………………………………….247 Figure 8 - 30: The product of forward and reverse amplitudes for various offset shot pair. There are gross similarities in the shape of the amplitude products for all the shot pairs, suggesting that the variations are related to variations in the seismic velocities in the refractor…………………………………………………247 Figure 8 - 31: Final sections from WET tomography. Section generated using 1D starting model with multiple layers of assumed constant vertical velocity gradients (top). Section generated using 2D starting model derived from the GRM (bottom). Despite the large differences in the sections both fit the field data to acceptable accuracy………………………………………………………252 Figure 9 - 1: 2D contour map of magnetic data at site 1………………………258 Figure 9 - 2: 2D contour map of magnetic data at site 2………………………258 Figure 9 - 3: The resistivity result at Site1-Line1……………………………….263 Figure 9 - 4: The resistivity result at Site1-line2………………………………..263 Figure 9 - 5: The resistivity result at Site1-Line3……………………………….264 Figure 9 - 6: The resistivity result at Site2………………………………………264 Figure 9 - 7: The possible P-wave undetected layer. A change from partially saturated to completely saturated sediment above the refractor will result in an increase in P-wave velocity. If this layer of completely saturated sediment is relatively thin, then it will be undetectable in the P-wave traveltime graphs, resulting in the refractor being imaged shallower than the true depth. Since S-waves do not propagate through liquids, the degree of saturation in the sediment will have no effect on its propagation, thus the S-waves will image the refractor at its true depth…………………………………………………………..267 Figure 9 - 8: Summary of the refractor velocities at site 1. The refractor can be separated into three distinct regions of different seismic velocities…………..269 Figure 9 - 9: Cross-cutting features at site1…………………………………….270 Figure 9 - 10: Amplitude product for P-wave at site1………………………….272 Figure 9 - 11: Wave-type conversion and the GRM. The forward and reverse traveltimes employed by the GRM algorithms include one set derived from SV-P-SV-waves and the other set from pure SV-waves. When this occurs, the refractor is shifted laterally and is imaged at a greater depth…………………276
xvii
Figure 9 - 12: Summary of the refractor velocities at site 2. Refractor can be separated into four distinct regions of different seismic velocities……………277 Figure 9 - 13: Cross-cutting features at site2…………………………………..280 Figure 9 - 14: Distinction among the particle-displacement vectors associated with the three fundamental modes, P, SH, and SV, which compose vector-wavefield seismic data (Hardage et al., 2003)………………………………….281 Figure 9 - 15: Simplified regional geology map of the Hill End Trough (Modified from David et al., 2003)……………………………………………………………284 TABLES Table 1- 1: Areas (hectares) with a high potential to develop dryland salinity in Australia (From: National Land and Water Resources Audit, 2001)………….9 Table 2 - 1: Summary of the local, intermediate and regional systems (Coram, 1998)…………………………………………………………………………………26 Table 2 - 2: Evaluation attributes for different dryland salinity management options (from Coram et al., 2001)…………………………………………………38 Table 2- 3: Land management recommendations for the Yass Valley (modified from Nicoll and Scown, 1993)……………………………………………………..40 Table 4 - 1: Some applications of refraction methods………………………….66 Table 5 - 1: Description of different rocks type in the Spicers Creek area (From Cobbora 100,000 geological map)………………………………………………..109 Table 5 - 2: Proposed deformation events affecting rocks across the Dubbo 1:250,000 geological sheet Scott (Morgan, 1997)………………………………119 Table 5 - 3: Survey parameters for the Gilgandra and Dubbo airborne surveys………………………………………………………………………………126 Table 5 - 4: IAEA (International Atomic Energy Agency) recommended windows for conventional 3-channel airborne gamma-ray spectrometry (IAEA, 1991).129
Table 5 - 5: Interpretation of relative concentration of K. Th, and U in rocks based on the RGB image (Morgan, 1997)……………………………………….130
40
Table 7 - 1: Trace to station relationship at site 1………………………………190 Table 7 - 2: Trace to station relationship at site 2………………………………194 Table 7 - 3: The number of shot recorded in site 1 and site 2…………………196 Table 9 - 1: Minimum, maximum and average time-depth and depth for each wave-type……………………………………………………………………………266 Table 9 - 2: Minimum, maximum and average time-depth and depth for each wave-type…………………………………………………………………..………..274
xviii
PHOTOS Photo 2 - 1: An artesian bore discharging saline groundwater. The severity of the soil erosion is evident. Note the difference between the top of the concrete collar on the casing and the existing ground surface……………………………31 Photo 2 - 2: A photo taken to the left of photo 1 (the tree and post appearing in both). The degree of soil erosion caused by saline groundwater destroying the soil structure is evident, especially along the far edge of the accumulated groundwater……………………………………………………………………….…31 Photo 2 - 3: 50 m downstream from the artesian bore in photo 1. Salt completely covers the soil around the stream, indicating the streamwater’s high salinity levels…………………………………………………………………………32 Photo 2 - 4: Upstream from photo 3 the amount of salt covering the soil is even greater. Also notice the soil erosion along the banks of the stream………..….32 Photo 7 - 1: The ARAM24 CRU………………………………………………….177 Photo 7 - 2: IVI Mini-vibrator Model T-15000…………………………………..179 Photo 7 - 3: The vibrator operating in P-wave mode. The actuator (the blue cylinder above the base plate) is positioned vertically…………………………180 Photo 7 - 4: The vibrator operating in S-wave mode. The actuator is positioned horizontally and is perpendicular to the axis of the truck………………………180 Photo 7 - 5: GS-20DM 14 Hz three component geophone. A bubble spirit level (the white circle) is used to plant the geophone in the correct orientation. Each component is transmitted through a separate cable (marked red, white and yellow)……………………………………………………………………………….181 Photo 7 - 6: A cartridge of anzomex, a high frequency explosive composed of TNT and PENT……………………………………………………………………..182 Photo 7 - 7: The Pelton Shot Pro dynamite firing system was used to remotely control detonation……………………………………………………………….….182 Photo 7 - 8: A successful detonation of a small charge……………………….187 Photo 7 - 9: four receiver lines at site 1…………………………………………189 Photo 7 - 10: One of the receiver lines ready for recording. Pin flags are first surveyed. Three receiver cables are then laid out. Receivers are planted next to the pin flags and connected to the cables. The station units are then connected to the cables………………………………………………………………………..190 Photo 7 - 11: Minivib recording at site 2…………………………………………193 Photo 7 - 12: Geophones and cable deployment at site 2…………………….193 Photo 7 - 13: A successful detonation of a small charge at site2…………….194
xix
Chapter 1
Introduction
1.1 Recent Advances in the Resolution of Geophysical Data
In the last twenty-five years, there have been significant advances in the spatial
resolution of most geophysical sets of data. Whereas geophysical data were
once commonly acquired along widely spaced profiles, it is now more usual to
obtain geophysical data along numerous closely spaced traverses. As a result,
the spatial sampling is now comparable in both horizontal directions and the
greater resolution has greatly improved the interpretation of geological features.
Figure 1-1, taken from Reeves (1992), illustrates the effects of improvements in
the spatial sampling. The upper contour map was derived from data acquired in
1973 with a fluxgate magnetometer with accuracy of ±10 nanoTeslas and a line
spacing of nominally 1.5 km. The lower grey scale image was derived from
data recorded in 1991 with a cesium vapour magnetometer with a resolution of
±1 nanoTesla and a line spacing of 400 metres. Despite the resolution ±10
nanoTeslas of the fluxgate magnetometer, it is still not possible to recognize
features such as the folds on the left hand side or the dykes on the right hand
side, in the contour map, even though they can be readily observed in the grey
scale image. The contouring process has the properties of a low pass filter, and
does not permit the retaining of the higher frequency components.
1
Figure 1 - 1: Improvements in the resolution of geological detail have been achieved by
adopting a new generation of survey specifications and presenting the resulting grid as
an image. This illustration is from an area of the Ebagoola 1:250 000 sheet in northern
Queensland. Top: BMR survey flown (E-W) in 1973 at a mean terrain clearance of 150 m
and a line spacing of 1500 m, presented as a contour map with contour interval 10 nT.
Bottom: the same area flown (E-W) by Geoterrex in 1991 with 100 m terrain clearance and
400 m line spacing (Reeves, C.V., 1992).
2
However, it is clear that the line spacing of 1.5 kilometres is too large to permit
confident correlation of subtle anomalies on adjacent profiles, and as a result,
only a low resolution result is possible. Therefore, it can be argued that despite
its format as a two-dimensional map, the contour map has essentially the same
resolution as a series of one-dimensional profiles. It is essential to employ
significantly closer line spacings, in order to achieve a genuine 2-D image.
Perhaps one of the most spectacular demonstrations of the benefits of
improved spatial resolution in both horizontal directions, has been three-
dimensional (3-D) seismic reflection methods. In the past twenty-five years,
high resolution 3-D seismic reflection methods have revolutionized the
investigation of the sedimentary basins (Weimer and Davis, 1996). The
fundamental problem facing two-dimensional (2-D) seismic methods is the fact
that most geological targets are three-dimensional. In the past, two-
dimensional methods have attempted to counter this problem by locating lines
with the strike and dips of the major features in mind, minimizing but rarely
eliminating the effect of the third dimension (Brown 1996).
French (1974) demonstrated that 3-D migration eliminates many of the lateral
correlation ambiguities which are caused by “sideswipes” and “blind structures”.
Structure maps derived from 3-D migration processing give a true and precise
picture of 3-D models, whereas the same data when processed using 2-D
migration, results in the mapped structures being distorted in shape. However,
the distortions can be reduced by increased coverage and careful interpretive
techniques (French, 1974). As a result, to obtain a precise image from seismic
3
data it is necessary to employ sampling densities and processing methods
which recognize and accommodate the three dimensions.
Furthermore, the geological targets studied in the near surface surveys often
display considerable variations in cross-line and in-line directions in terms of
depth to and seismic velocities within (Palmer, 2001b). This makes the
adoption of 3-D methods inevitable.
Perhaps one of the most compelling demonstrations of the benefits of 3-D
geophysical methods has been the study by Nestvold (1992) which
demonstrates that 3-D seismic reflection results are often significantly different
from 2-D seismic reflection results. In many cases, 2-D methods can provide
an incorrect rather than an incomplete picture of the subsurface structure. The
differences can be explained by the fact that the geological targets can vary
significantly in both horizontal directions, thereby requiring the adequate spatial
sampling obtained with 3-D methods. In other words, the incorrect subsurface
structure can often be the result of spatial aliasing. Furthermore, most seismic
reflection processing routines, especially imaging (also known as migration),
perform much better in three dimensions because they recognize and
Per cent retention of native deep-rooted vegetation in high risk areasLong-term groundw ater level trendsLong-term stream salinity and salt load trendsExtent of land salinisedLand cover changeHydrogeological and soil characteristicsClimatic characteristics
Long-term groundw ater level trendsLong-term stream salinity and salt load trendsExtent of land salinisationLand cover changeHydrogeological and soil characteristicsClimatic characteristics
Long-term groundw ater level trendsLong-term stream salinity and salt load trendsExtent of land salinisationHydrogeological and soil characteristicsClimatic characteristics
Long-term stream salinity and salt load trendsExtent of land salinisationHydrogeological and soil characteristicsClimatic characteristics
Trea
tmen
t of c
ause
Prev
enta
tion
Land and surface w aters w ill remainfree from degradation from salinity
Land and surface w ater salinisationw ill diminish
Interception and disposal of salt, and reduction of groundw ater levels in transmissin zones
Land and surface w ater salinisationw ill diminish
Land and surface w ater salinisationw ill diminish
Reduce recharge to the groundw ater system by increasing vegetation w ater use through land management practices
Interception of w ater prior to infiltration or from groundw ater upgradient of discharge zone
Manage vegetation so deep drainagedoes not increase in high risk areas
38
Table 2- 2 (continued): Evaluation attributes for different dryland salinity management options (from Coram et al., 2001).
Managing saline discharge Long-term stream salinity and salt load trendsHydrogeological and soil characteristicsClimatic characteristics
Application of soil treatments (ameliorants) Extent of land salinisationSoil w ater characteristicsHydrogeological and soil characteristicsClimatic characteristics
Establishment of salt tolerant land cover Extent of land salinisationLand cover changeIncrease in production of biomass
Extent of land salinisationLand cover change
Optimistation of the use of non-saline resources Land cover change
Alternative use of saline land and w aterresources
Productivity from non-salt affectedland w ill increase
Productivity from salt affected landand w ater w ill increase
Prod
uctiv
e us
es o
f sal
ine
reso
urce
s
Land and surface w ater salinisationw ill diminish
Agriculture production from saline landis increased
Amel
iora
tion
of
sym
ptom
s
Productivity from salt affected landw ill increase
39
Table 2- 3: Land management recommendations for the Yass Valley (modified from Nicoll and Scown, 1993)
Land managementclass
Land use Description Recommended technique Attributes of recommendedland management technique
cropping country deep soils land capable of regular cultivationw ell structured maintain soil structure by crop rotation
protect area w ith w indbreaksavoid long, bare fallow s, avoid fallow s on f lood prone areasuse legume break crops to replace fallow smanage area to produce maximum vegetative production
rotational cropping country deep soils land capable of regular cultivationw eak structure maintain soil structure by crop rotation
prevent soil erosion by use of soil conservation earthw orksuse direct drill or reduced tillage practisesavoid long, bare fallow s, avoid fallow s on f lood prone areasestablish perennial pastures after cropping phaseprotect area w ith w indbreaks
good grazing country moderately deep soils agroforestry techniques utilise high WT and good soilsestablish w indbreaks and treelotsretain all existing timberfertilise according to soil tests, on a regular basismaintain a minimum 70% groundcover all yearspread stock w atering points to ensure even pasture grazingdo not graze pasture below 4cm in height
improved perennial pastures, potentially agroforestry as technology develops
LM 3
conservative farmingtechniques
LM 1
LM 2
minimum 4 years perennial pasture phase betw een crops
40
Table 2- 3 (continued): Land management recommendations for the Yass Valley (modified from Nicoll and Scown, 1993).
Land managementclass
Land use Description Recommended technique Attributes of recommendedland management technique
seasonal grazing country seasonally w aterlogged should not be cultivated due to the high risk of soil erosionretain w ater tolerant perennial speciesestablish perennial improved pasturesgraze w hen soils are firmmaintain even cover of vegetation, not less than 70%do not graze pasture below 4cm in heightremove stock w hen soils become w aterloggedfence out actively eroding drainage lines
poor grazing country soils highly erodible most suited to native pastures encourage re-seeding of desirable native specieslow available moisture use deep rooting perennial grasses
apply fertilisers on a regular basisavoid set-stocking natve pasturesmaintain minimum 70% ground cover all yearestablish w indbreaks and treelotsspread stock w atering points to ensure even pasture grazingdo not graze pasture below 4cm in heightaim tow ards 20% tree cover
fence out arearetain existing treesplant salt tolerant grasses, shrubs and treesimplement earthw orks as required, to stop soil erosionstraw mulch area for better germination of grasseslighlty graze for short periods onlydo not graze w hen w aterloggeddo not graze pasture below 5 cm in heightmaintain minimum 70% ground cover all year
most suited to w ater tolerant pastures. Not suitable for cultivation
suited to saltland agronomy, lands to be treated for saline discharge
41
Table 2- 3 (continued): Land management recommendations for the Yass Valley (modified from Nicoll and Scown, 1993).
Land managementclass
Land use Description Recommended technique Attributes of recommendedland management technique
steep, stony shallow soils retain mature forest cover and regrow thestablish block plantings of treesexclude livestock at critical periods of tree f low eringperiodically graze to reduce fire risksmaintain minimum 70% ground cover all yearselective thinning of trees is possible on slopes
non-agricultural land very steep slopes severe soil erosion hazardsrocky areas fence area to exclude livestock and feral animalsw ater bodies encourage native animals/habitat
encourage natural forest regeneration
best suited to natural vegetation, not suited to agricultural production
LM 8
suited to green timberproduction
best protected from soil erosion w ith trees, potentially high recharge areasLM 7
42
Chapter 3
Shear Zones Characteristic
3.1 Summary
Despite their significance, shear zones are difficult to map using conventional
geological methods. Shear zones can be very difficult to detect from the
surface, and therefore the use of surface mapping to define them is frequently
unreliable. Without the use of geophysical methods, shear zones are frequently
inferred from geological mapping, with little direct evidence of their existence.
Shear zones may act either as closed or open geochemical systems with
respect to movement of fluids and elements, irrespective of the size of the zone.
Shear zones are known to perform as preferred pathways for fluid circulation.
Little is understood about the 3D geometry (the relationships between length,
width and depth) of shear zones and in particular the depth of their root zones.
In the Spicers Creek Catchment magnetic, electrical resistivity and seismic
refraction methods were used to image the inferred shear zones. It was
anticipated that in this area, seismic refraction would provide the best results
since the propagation of seismic waves would not be affected by the high level
salinisation whilst also providing possibly greater resolution.
43
3.2 Introduction
Dennis (1967) showed that the term shear zone refers to a tabular region of
pervasively faulted rock, that is, a fault zone containing a very large number of
closely spaced and anastomosing fault surfaces. A shear zone is defined by
Bates and Jackson (1980) as “a tubular zone of rock that has been crushed and
brecciated by many parallel fractures due to shear strain”. Ramsay and Auber
(1983) defined a shear zone as a planar or curviplanar zone of high deformation
which is long relative to its width (length to width ratio greater than 5:1) and
which is surrounded by rocks showing a lower state of finite strain. Shear
zones may consist of a series of interleaving, anastomosing brittle faults and
crushed rock (cataclasite) formed near the surface, or of ductile shear zones
composed of mylonitic rocks produced by faulting at great depth (Hatcher,
1995). Shear zones represent regions of weak rock as well as sources of major
water flow.
There are some differences between fault, fault zone and shear zone. A fault is
a planar discontinuity between blocks of rock that have been displaced past one
another, in a direction parallel to the discontinuity (Figure 3-1A). A fault zone is
a tabular region containing many parallel or anastomosing faults (Figure 3-1B).
A shear zone, is a zone across which blocks of rock have been displaced in a
fault-like manner, but without prominent development of visible faults (Figure 3-
1c, Hobbs, 1976). Shear zones are thus regions of localized ductile
deformation, in contrast to fault zones that are regions of localized brittle
deformation. Unlike ordinary fault surfaces, shear zones commonly do not
44
display any discrete physical break (Davis and Reynolds, 1996). Another
distinction is that the normal component of displacement, which is negligible for
fault and fault zones, but may be appreciable for shear zone. The displacement
across a shear zone can be inclined at any angle, other than 90°, to the
boundaries of the zone.
Furthermore, little is known about the 3-D geometry of shear zones, especially
the depth of their root zones (Pili et al., 1997). The root zone is considered to
be deep when the outcrop length of vertical major shear zones is from a few
hundred to a few thousand kilometers because of their narrow parallel-sided
shape.
Foliation is a very common feature within the shear zones, and it is usually
inclined at an angle to the zone’s boundaries. This foliation is commonly curved
and rotates toward parallelism to the shear zone boundaries in the central
portion of the shear zone. In practice, however, although the obliquity of
foliation is relatively common, curvature of foliation within the shear zone is rare
and in many cases, foliation occurs parallel to the shear zone boundaries. The
foliation, which may be a cleavage or schistosity, represents a profound
modification of the rock’s fabric and as such imparts a mechanical anisotropy to
the rocks within the shear zones (Ramsay and Lisle, 2000).
45
Figure 3 - 1: (A) Fault. (B) Fault zone. (C) Shear zone (From Hobbs et al, 1976).
3.3 Shear Zone Types
Ramsay (1980) has classified the shear zones as brittle, brittle-ductile and
ductile (Figure 3-2). In ductile shear zones the deformation state varies
continuously form wall to wall through the zone, brittle shear zones where the
wall are separated are separated by a discontinuity on the fracture surface, and
various intermediate types known as brittle-ductile shear zones combining these
geometric features in different proportions. Brittle, brittle-ductile or ductile shear
zones are known in various lithologies and structural levels (Ramsay, 1980)
down to upper mantle (Vissers et al., 1991).
Ramsay (1980) attributes two conditions as geometric features of ideal shear
zones. The first is shear zones generally have parallel sides and the second is
the displacement profiles along any cross section through a shear zones should
46
be identical. However, some shear zones have complex geometries and have
subparallel margins and retain a fairly consistent thickness over much of their
length.
Figure 3 - 2: Different types of shear zones from (Ramsay and Auber, 1983).
Davis and Reynolds (1996) have classified shear zones in terms of general
characteristics (continuous or discontinuous); shear zones margins (parallel,
diverging and converging); shear zone sets (parallel, anastomosing and
conjugate); deflection and offset across shear zones (dextral, sinistral, normal,
reverse); plate tectonic setting of shear zones (ocean-continental convergence,
Figure 3 - 4: Finite shear strain vs. time curves (solid) for the three types of shear zones,
for rock at the final zone margins (m) and in the interior (i). P (dashed line) represents
some process-controlling parameter (like temperature or fluid chemistry) that is assumed
to change over the shearing history. In Type II zones, the marginal and interior acquire
their low-strain state at about the same time and p value (heavy dots). The marginal
rocks, therefore, provide a good record of former low-strain state of the interior rocks. In
zones of Types I and III, the corresponding dots are farther apart in time and p, so the
marginal rocks provide a less reliable record of the former low-strain state of the interior
rocks. The type III diagram has been drawn for the case where the marginal rocks begin
straining at the same time as the interior rocks, but more slowly. In each diagram, the
right-hand dotted line is drawn at the point along the m line where its slope (the marginal
strain-rate) becomes zero. The horizontal position of this line, therefore, represents the
time at which rock in the shear zone margins acquired its final state of (low) strain. The
left-hand line represents the time at which rock in the shear zone interior passed through
this same state of (low) strain (From Means, 1995).
50
3.4 Fluid Movement in Shear Zones
Fault zones can act both as barriers to fluid flow and as conduits along or
across which fluid flows ( Anderson et al., 1991, Anderson et al., 1994, Hooper,
1991, Knipe, 1993, Berg and Avery, 1995, Lopez and Smith, 1995, Heneberg,
1995, Sibson, 1996, Roberts et al., 1996, Yielding et al., 1997, Knipe, 1997,).
Faults and their movements play a key-controlling element in fluid flow systems
in geological settings. The direction of fluid flow through a rock mass is
governed by the maximum hydraulic gradient (not necessarily vertical), existing
permeability anisotropy (e.g. bedding, foliation), and superimposed structural
permeability (Sibson, 1996). Furthermore, the permeability of a fault zone
depends on additional factors such as juxtaposition of permeable units across
the fault, rock properties, stress field, the amount of sand, clay, or fault gouge in
the fault zone, the presence of material with high capillary entry pressure, and
the presence or absence of an open fracture network (Smith, 1980, Allan, 1989,
Miller, 1995, Sibson, 1996, Knipe, 1997). Philips (1991) demonstrated that
fractured rocks can have much larger permeabilities, and hence hydraulic
conductivity, than that of completely homogenous rocks. As a result,
permeability may vary by many orders of magnitudes, depending upon rock
type and the pressure-temperature conditions (Braun et al., 2003).
Shear zones have long been shown as preferential fluid pathways (Beach and
Fyfe, 1972). Shear zones may act both as closed or open geochemical
systems with respect to movement of fluids and elements, irrespective of the
size of the zone (Hatcher, 1995). Shear zones are known to perform as
51
preferred pathway for fluid circulation (Beach and Fyfe, 1972). Little is
understood about the 3D geometry (the relationships between length, width and
depth) of shear zones and in particular the depth of their root zones (Pili et al.,
1997).
For an understanding of the dynamics and rates of fluid motion around and
within faults, it is necessary to understand the behavior of hydrothermal
systems. Braun et al. (2003) presented the result of simple numerical
experiments of fluid flow around and within a permeable fault embedded in a
less permeable porous medium. They showed that fluid flow is controlled by
two timescales: and , where S is the specific storage
of the porous material, l the length of the fault, and and are the
hydraulic conductivities of the porous material and the fault, respectively.
Physically, the three conditions correspond to:
Ff KSl /2=τ MF KSl /2=τ
MK FK
If ( fττ < ), fluid rapidly travels through the fault driven by the initial vertical
pressure gradient (or hydraulic head);
If ( Ff τττ << ), the fluid flow in the fault leads to pressure changes in the rock at
both extremities of the fault that result in a decrease of the pressure gradient
between the extremities of the fault;
If ( Fττ < ), fluid flow in the matrix has adapted itself to the new permeability
structure and a steady-state pressure field has developed around the fault.
The first time constant fτ , is the time required by the fluid to travel inside the
fault. fτ is therefore a function of l and . The second time constant, FK Fτ , is
52
the time required by the fluid to reach a steady-state flow and associated
pressure gradient in the matrix in the vicinity of the fault. Fτ is therefore a
function of l and . For a 1 km-long fault, MK fτ varies between 30 years and
300 kyr, depending on the value of the hydraulic conductivity, whereas Fτ varies
between 3Myr and 30 Gyr. Braun et al. (2003) showed that the length of a fault
is the most important factor which determines the value of the two time
constants, fτ and Fτ .
It has been useful to find ‘steady-state flow’ in rocks, flow at constant stress,
constant temperature, constant strain rate and constant microstructure, so that
the flow law can be expressed in a strain-independent form that allows
separation of variables, thus (Rutter, 1999):
( ) ( ) ( )ST 321 ∫∫∫Α= σε
In which ε is strain rate, σ is flow stress, T is absolute temperature and is
some measure of rate controlling microstructure, which might be grain size. It is
widely recognized that fine-grained rocks can be weaker than their coarser-
grained counterparts at high temperature (Walker, et al., 1990), if grain size
sensitive flow processes were to be activated. As Rutter (1999) discussed,
rocks in shear zones are commonly tectonically reduced in grain size, an
observation which tends to reinforce the inference that shear zones are zones
of relative mechanical weakness.
S
53
3.5 Evidence of Shear Zones
Shear zones have certain characteristic to recognize them in the field, in thin
section, on geological maps and cross sections (Davis and Reynolds, 1996),
and geophysical methods. Geophysical methods are very useful when there is
no clear evidence of shear zones in the field. Shear zones exist at all scales
and their size depends on the conditions, such as brittle, ductile, or intermediate
conditions, under which they were formed.
3.5.1 Field Evidence
The following evidence for detecting shear zones (Ramsay and Auber, 1983;
Simpson and Schmid, 1983; Passchier and Simpson, 1986; Davis and
Reynolds, 1996) might be established in outcrops, in hand specimens or in thin
sections.
1 S-C band structure (Figure 3-5A): S-C fabrics consist of foliation (S-
surfaces) and shear bands (C-surfaces). The S and C surfaces are at an angle
to one another, indicating the sense of shear. In most S-C fabrics, the S-
surfaces are a typical foliation, being defined by flattened lenticular grains and
aggregates of grains. The oblique relationship of the two structures shows a
clear indication of the displacement sense. S-C fabrics can be divided into two
general types (Lister and Snoke, 1984). Type I are the dominant S-C fabric in
most mylonitic quartzofeldspathic rocks. Type II S-C fabrics are most common
in micaceous quartzite and some mylonites. In both types of S-C fabric, the C-
54
surfaces are commonly aligned parallel to the shear zone and are associated
with the same strain gradients as any continuous shear zone.
2 Rotation of porphyroblasts or porphyroclasts (Figure 3-5B): If the
large grains instead are new metamorphic grains that grew during or after
deformation, they are called porphyroblasts. Grains interpreted to be relics from
the protolith are called porphyroclasts because they commonly represent
fragments or clasts of original phenocrysts or detrital grains. In many
metamorphic rocks it is well known that large porphyroblasts can develop;
albite, garnet, staurolite are particularly useful for the determination of shear
sense. These crystals can roll like ball bearings and by matching the inclusion
trails they contain inside the crystal ( ) with those found outside ( ), the
rotation sense can be determined. Such crystals can grow as porphyroblasts
during the shear motion when they produce a “snowball structure”.
is es
3 σ -Structure (Figure 3-5C): This structure is found around
porphyroclasts, and relates to the form of the dynamically recrystallized
pressure shadow tail in relation to the porphyroclast. In this structure the
median line of the tail does not cross the tend of the average shistosity, and the
structure appears to form because the rotation rate of the tail is higher than that
of the porphyroclast. This structure is characteristic of shear zones with low
shear strains where the recrystallization rates are higher than the rotation rates.
4 δ - Structure (Figure 3-5D): This structure, like σ -structure, is
developed around porphyroclasts, but where the recrystallized pressure shadow
55
tail is strongly rotated by the porphyroclasts. This structure shows the
extremely curved and often embayed nature of the generally narrow tail, and
where the median line crosses the general trend of the schistosity. This type of
structure is particularly characteristic of origins of high shear strains where the
recrystallization rates are lower than the rotation rates.
5 Bookshelf sliding (Figure 3-5E): When broken crystals are subjected to
shear, the individual parts are rotated in the direction of shear, a feature which
sets up a contrary shear motion between the fragments much like a collapsing
set of books. Further shearing motion can lead to separation of the individual
fragments and the development of recrystallized pressure shadow zones
between the fragments.
6 Displaced crystals (Figure 3-5F): Where crystals possessing a well
developed crystal cleavage lie with this cleavage close to the shear plane, they
often become internally detached by gliding on this surface. In highly deformed
mylonites, for example, it is not uncommon to find individual mica fragments
connected by zonal films of phyllosilicates to form free “floating mica fish”.
7 Dynamic recrystallization (Figure 3-5G): Rocks that have been very
highly deformed by crystal plastic processes frequently show an intense
banding and schistosity. Within the more deformed crystal (especially with
quartz and calcite) the bands are crossed obliquely by small grains with a
marked individual alignment.
56
8 Preferred optical orientation (Figure 3-5H): Plastic flow of crystals
leads to a variety of different types of preferred alignment patterns of the optical
axes of the crystals in the rock aggregate. Many of these complications of
interpretation of preferred orientation patterns result from the nature of the
environmental conditions during deformation.
9 Protomylonite, mylonite, and ultramylonite: Mylonites are strongly
foliated metamorphic rocks that exhibit high ductile strain and incomplete
recrystallization or recovery. Mylonite contains 50 to 90 per cent matrix.
Mylonite is typical shear zones rock. The initial stages of mylonitization produce
a weakly to moderately mylonitic rock, called a protomylonite. Protomylonite
contains less than 50 per cent fine-grained matrix. Ultramylonite contains more
than 90 per cent matrix and less than 10 per cent relict grains. Some foliation in
mylonitic rocks is the expression of thin shear zones.
57
Figure 3 - 5: Different types of criteria used for the determination of shear sense in shear
zones (From Ramsay and Auber, 1983).
58
3.5.2 Imaging Shear Zones with Geophysical Methods
Despite their significance, shear zones can be difficult to map in detail using
conventional geological methods. However, their location can often be inferred
from topographic features or geophysical methods, such as electromagnetic
techniques. Shear zones can be difficult to detect from the surface, and
therefore the use of surface mapping to define them is often unreliable. Without
the use of geophysical methods, shear zones are frequently inferred from
geological mapping, often with little evidence of their existence either directly or
tectonically. It should be emphasized that detecting and studying shear zones
are not always a straightforward procedure.
Major faults and shear zones are routinely discovered by explorationists through
geophysical methods. Nowadays, integrated use of least two geophysical
methods with geological interpretation is important to obtain accurate geological
interpretation. For example Street and Engel (1990) showed that integrated use
of geophysical methods (magnetic and seismic refraction), drilling, and sample
analysis has hound geological controls on the location of dryland salinity in
Western Australia. The success of the electromagnetic and resistivity
geophysical methods in imaging shear zones lies in their ability to image these
narrow zones as regions of high conductivity. This increase in conductivity is
due to the development of secondary porosity and permeability from fracturing
which leads to an increase in groundwater saturation. The success of the
seismic refraction method lies in its ability to image the shear zone as a narrow
region of low seismic velocity and increased depths of weathering.
59
In the Spicers Creek Catchment magnetic, electrical resistivity and seismic
refraction methods were used to image the inferred shear zones. It was
anticipated that in this area, seismic refraction would provide the best results
since the propagation of seismic waves would not be affected by the high level
salinisation whilst also providing possibly greater resolution.
60
Chapter 4
Three-Dimension Three Component
Seismic Surveys
4.1 Summary
The objective of the seismic refraction method is to determine the velocity
distribution in the sub-surface. The seismic refraction method does have a
number of advantages over the reflection method, particularly in shallow
investigations such as in geotechnical, environmental and groundwater
applications. The seismic refraction methods seek to measure the spatial
variation of petrophysical parameters by using seismic velocity.
Many geological features in the subsurface are three-dimension (3-D) in nature,
and two-dimension (2-D) seismic section is a cross section of the 3-D seismic
wave field. Having a detailed image of the subsurface usually makes the
interpretation more reliable. The aim of using 3D-3C seismic data is to find
more information about rock properties with both P- and S-waves. The most
important parameters which can be extracted from P-wave and S-wave data are
the P-wave/S-wave velocity (Vp/Vs) ratio and Poisson’s ratio.
61
4.2 The Seismic Refraction Method
The foundation of the seismic methods is the theory of elasticity. The velocities
of the propagation of the compressional and shear seismic energy are functions
of the elastic moduli or constants, and the rock densities.
The seismic refraction method is based on the principle that when a seismic
wave (P or S wave) impinges upon a boundary across which there is a contrast
in velocity, then the direction of travel of that wave changes on entry into the
new medium. The amount of change of direction is governed by the contrast in
seismic velocity across the boundary according to Snell’s Law as shown in
Figure 4-1.
Figure 4 - 1: Snell’s Law.
When a seismic signal encounters a surface separating two media having
different elastic properties, it gives rise to reflected or refracted waves. While
most refracted events have not been reflected, most reflected events have been
refracted, because refraction occurs across any velocity interface in accordance
with the Snell’s Law.
62
4.3 Seismic Velocities in the Earth
The seismic velocity is generally used as a measure of earth elasticity, since P
and S seismic velocities can be derived from the elastic constants and the
density. Due to the variability of the near surface conditions, the seismic
velocities in shallow earth materials are also highly variable. The relevant
factors include major rock forming minerals, rock formation, the nature of the
cementing material, fracturing and weathering, structural deformation, etc. In
general, seismic velocities are less in unconsolidated materials.
Furthermore, the seismic velocities of earth materials usually increase with
depth, because of decreasing weathering, and increasing water saturation,
compaction and diagenesis with depth. The velocities of seismic waves in most
igneous and metamorphic rocks depend mainly on the elastic properties of the
minerals making up the rock material itself, since these rocks have little or no
porosity. In general, igneous rocks have seismic velocities which show a much
narrower range of variation than sedimentary or metamorphic rocks. The
average velocity for igneous rocks is higher than that for other types (Dobrin
and Savit, 1988).
Because the seismic velocities generally increase with depth, seismic energy is
systematically refracted away from the vertical, until it travels in a predominantly
horizontal direction. Furthermore, the horizontally traveling energy also radiates
energy back into the overlying layers. As a result, seismic energy which is
initiated at the surfaces propagates down into the subsurface and returns to the
63
surface at some greater distance. Figure 4-2 schematically shows the basics of
the seismic refraction method.
Figure 4 - 2: Seismic refraction wave raypaths.
It might be noted that the fundamental requirement for the seismic refraction
method to be a reality is the systematic increase of seismic velocities with
depth. By contrast, there is a systematic decrease in electromagnetic velocities
with depth and as a result, refraction methods are not possible with ground
penetrating radar methods.
The seismic refraction method has a number of advantages over the reflection
method, particularly in shallow investigations such as geotechnical,
environmental and groundwater applications. The seismic refraction method is
most effective in the near-surface for depths of less than several tens of metres.
64
On land, it can be difficult to obtain good seismic reflection results from the
near-surface. By contrast, the refraction method has a unique ability to provide
detailed velocity information on the deepest refractor enabling it to detect the
location of weathered and fractured zones in bedrock.
As a result refraction methods can obtain good results with shallow targets,
whereas reflection methods are more suited to deep targets. Refraction is now
rarely employed in oil exploration because of the large scale field operation
required (Dobrin and Savit, 1988). It should be noted that the processing
methods for shallow targets with refraction methods, such as the GRM, are
suitable and probably more reliable than reflection methods when there is are
lateral changes in seismic velocities of the refractor. Table 4-1 shows some
application of seismic refraction methods.
4.4 Three Dimension Seismic Refraction Methods
Perhaps the first three-dimensional 3D seismic methods were the fan shooting
refraction methods of the 1920s and 1930s. These 3D methods, which clearly
pre-date the 3D reflection methods of the last quarter of a century, demonstrate
the fundamental maxim that 3D targets require 3D seismic methods for
efficacious exploration. The fan shooting techniques later evolved into the
broadside shooting methods of the 1950s (Richards, 1959).
65
Table 4 - 1: Some applications of refraction methods.
Seismic Refraction Applications
Authors
Exploration of alluvial deposits and sediments
(Edge and Laby, 1931, Urquhart, 1956, pakiser and Black, 1957, Hobson and Grant, 1964, Paterson, 1965, Duguid, 1968, Young and Lucas, 1988, Biehler et al., 1991, Campos-Enriquez et al., 1997)
Groundwater studies, hydrology and hydrogeological investigations
(Bonini and Hickok, 1958, Warrick and Winslow, 1960, Levshin, 1961, Wiebenga and Jasson, 1962, Lennox and Carlson, 1967a, Lennox and Carlson, 1967b, Hasselstrom, 1969, Overmeeren, 1981, Haeni, 1986, Dobrin and Savit, 1988)
Investigation of engineering sites and geotechnical factors
(Shepard, 1939, Wood, 1940, Moore, 1952, Brown and Robertshaw, 1957, Welin, 1958, Hawkins, 1961, Bartlett, 1962, Drake, 1962, Stam, 1962, Svenson and Bowering, 1963, Bigelow, 1965, Irving, 1965, Lawson et al., 1965, Mallot, 1965, Atkinson, 1970, Cummings, 1979, Hatherly and Neville, 1986, Hatherly, 1986, Kilthy et al., 1986, Laa et al., 1991)
Rock fabric and anisotropy (Bamford and Nunn, 1979, Palmer, 2000c, Palmer, 2001a, 2001b, 2001d,2003b)
66
The development of 3D refraction methods in recent times has been restricted
to deeper crustal structures, and they are as much a result of the realities of
locating receivers in accessible locations, as with the acceptance of the 3D
geometry of the subsurface targets (See page 244 Palmer 1986).
The development of 3D refraction to shallow targets has not been widespread.
Bamford and Nunn (1979) used radial profiles to determine azimuthal
anisotropy. However, no attempt was made to resolve the dual challenge of
irregular refractor parameters and azimuthal anisotropy. Dean et al (2000)
applied refraction tomography to a set of data recorded across an alluvial
channel.
Palmer (2001a and 2003a) carried out a 3D refraction survey across a large
shear zone near the massive sulphide deposit at Mt Bulga, near Orange in
southeastern Australia. The survey showed that the 3D results can display a
significantly more complicated geological picture than the standard 2D depth
cross section. Not only did the 3D results reveal significant cross cutting
features (faults?), but the azimuthal anisotropy in one section of the survey was
orthogonal to that of the major structural elements. The significantly more
detailed structural model revealed by the 3D results would have considerable
value in modelling groundwater or contaminant fluid flow, designing grouting
programs, or even designing simple engineering constructions.
67
4.5 Three Component Seismic Refraction Surveys
The aim of using three-component coordinate system is to find more information
about rock properties with both P- and S-waves. Principally, propagation
velocity of P-waves is affected by both rock incompressibility and rigidity,
whereas, S-waves is affected by rock rigidity only. P-wave velocity is a function
of three separate rock properties, the bulk modulus, shear modulus and bulk
density, whilst S-wave velocity depends upon the shear modulus and bulk
density (Domenico and Donbom, 1986). Used separately these velocities can
be ambiguous indicators of rock lithology.
The motion associated with a wave arrival is a vector quantity, and measuring
its three orthogonal components requires the use of three orthogonal
geophones. Most three-component phones use one vertical phone (P-wave)
and two horizontal phones (SV- and SH-wave). A field-recording technique
sometimes adopted in fracture studies, which successively record from
compressional wave (P source) and shear waves (SH-, and SV- sources) with
three-component geophones; this is called nine-component recording because
separate sections can be generated for each of the three-geophone outputs for
each of the three sources (Sheriff and Geldart, 1995). If the sources are pure
and the earth isotropic, three of the nine sections would be conventional P-, SH-
, and SV-sections, two would show converted waves (P to SV and SV to P), and
the other four would be blank (Figure 4-3).
68
The coordinate system for measuring the three components of a P- or S-waves
can be rotated around any axis; in particular, they can be rotated into a natural
coordinate system. Such a coordinate rotation is sometimes called Alford
rotation or polarization filtering (Sheriff and Geldart, 1995).
Three-component recording can be employed to resolve fracture problems in a
3-D survey (Lewis et al., 1991). On the other hand, if the SH- and SV- sources
and geophones are not oriented parallel and perpendicular to the natural
orientation or azimuthal anisotropy is present, then shear wave splitting will
result in energy showing on all panels.
Conventional P-wave data
Conventional SH-wave data
Conventional SV-wave data
Mode conversion P to SV
Mode conversion P to SV
?
?
?
?
SV
SH
P
P SH SV
SOURCE MOTION
Figure 4 - 3: Nine-component sections recorded by P-, SH-, and SV-geophones for P-, SH-
, and SV sources. Without anisotropy the sections marked ? would be blank (From
Tatham and McCormack, 1991).
69
4.6 Shear Waves Versus P waves
P- and S-waves are considered as body waves because these two types can
propagate through the body of an elastic solid. P-waves are also known as,
longitudinal, primary, push, dilatational, irrotational or compressional waves.
Material particles oscillate about fixed points in the direction of wave
propagation by compressional and dilatational strain, similar to a sound wave.
S-waves are also known as the transverse, secondary, rotational, or shear
waves. The seismic velocity in rocks depends on many factors, including
The study area is located on the Dubbo 1:250, 000 sheet area (Figure 5-7).
This sheet includes six 1:100,000 maps (Dubbo (8663), Wellington (8632),
Cobbora (8733), Euchareena (8732), Gulgong (8833) and Mudgee (8832)). It is
also part of central western New South Wales and includes sedimentary,
volcanic and plutonic rocks of Early Paleozoic to Cainozoic age. The Dubbo
1:250,000 map sheet area covers the exposed north-eastern margin of the
Palaeozoic Lachlan Orogen, which is fringed by the western margin of the
Sydney Basin, the southern margin of the contiguous Gunnedah Basin and the
southern margin of the Surat Basin.
Figure 5 - 7: Locality map for the Dubbo 1:250,000 geological map sheet area, showing the six component 1:100,000 map sheets and geographical features (Meakin and Morgan, 1999).
91
5.4.2 Climate
Figure 5-8 shows the mean monthly maximum and minimum temperatures,
mean monthly precipitation data and the average number of wet days for Dubbo
and Dunedoo, and pan evaporation data for Wellington (Schofield, 1998). In
the Dubbo region rainfall is greatest in January, with an annual average of 584
mm reported on the Dubbo 1:100,000 sheet and 642 mm on the Cobbora
1:100,000 sheet. The annual distribution of rainfall is tilted toward the summer
months, but June is also a wet month. The climate is temperate with annual
mean temperatures ranging from 19 to 32 ˚C in summer, and 3 to 19˚C in winter
(Schofield, 1998).
Figure 5 - 8: Mean monthly maximum and minimum temperatures, mean monthly precipitation data and the average number of wet days for Dubbo and Dunedoo, pan evaporation data for Wellington is shown at the base of Figure (Schofield, 1998).
92
5.4.3 Physiography
The Dubbo 1:250,000 map area covers 15,600 square kilometers over parts of
the Shires of Cabonne, Wellington, Dubbo, Coolah, Mudgee, Merriwa and
Rylstone. The Western railway and several major roads, including the Mitchell
Highway, serve the area (Meakin and Morgan, 1999).
Topographically, the area falls from 740 m ASL in the east, to around 240 m
ASL near Dubbo (Figure 5-9). The topography in the Dubbo region has been
substantially reduced to a peneplain by erosion. This process mixes low-lying
lava flow relics with isolated topographic highs comprising plugs and earlier but
topographically inverted lava flows.
In the area there are three main drainage systems: Macquarie, Cudgegong and
Talbragar Rivers. The largest is the Macquarie River. The Macquarie River
flows from the western boundary of the study area to the NW through Dubbo.
The Talbragar flows in the north-west direction, through Dunedoo via Ballimore
to the Macquarie River. The two rivers link together 5 km north of Dubbo.
Talbragar River from east to west has 11 tributaries: Sandy Creek, Boomley
Creek, Spring Creek, Baragonumble Creek, Goan Creek, Spicers Creek, Back
Creek, Rocky Creek, Mitchells Creek, Jones Creek and Troy Creek. The
Gudgegong River flows in a southerly direction from the east side of Gulgong
town.
93
The Great Dividing Range extends along the eastern edge of the map area.
Surface elevation gradually decreases from as much as 1068 m ASL in the
south-east (near Harbon) to about 300 m ASL in the north-west of the map area
(Meakin and Morgan, 1999).
Figure 5 - 9: Topography, township and surface drainage map for the Dubbo and Cobbora 1:100,000 sheets.
94
5.4.4 Dubbo Area
The Dubbo 1:250,000 map sheet area covers the exposed north-eastern margin
of Lachlan Orogen. Meakin and Morgan (1999) divided the Lachlan Orogen
rocks into seven lithotectonic associations based on age, lithology and
depositional environment. These include the:
Early Ordovician to Late Ordovician generally deepwater, mafic to intermediate
volcanic and volcaniclastic rocks and intrusions;
?Early to Middle Ordovician quartz-rich turbidities;
late Early Silurian to Early Devonian shallow water to deepwater volcanic and
Sandstone) of the Surat Basin were deposited along the northern margin of the
area. They represent deposition due to erosion which commenced in the Late
Jurassic and continued into the Tertiary.
In the Cainozoic, deposits on the Dubbo 1:250,000 map are represented by
alluvial, colluvial and lacustrine sediments, basalt and duricrust. No formal
names have been applied to unconsolidated Quarternary units on the Dubbo
1:250,000 map sheet. Gravel deposition represents climatic change by high-
energy, braided river systems, which gradually diminished in energy with the
development of aridity in inland Australia. Tertiary basalt overlies and intrudes
all pre-Cretaceous rocks across the Dubbo 1:250,000 map sheet area, with
relict basalt preserved predominantly over the Mumbil Shelf and Capertee High
(Dubbo area and Liverpool Range), and also over the southern end of the Hill
End Trough (Meakin and Morgan, 1999).
102
5.4.6 Metamorphism
Figure 5-12 shows the metamorphic zones on the Dubbo 1:250,000 map sheet
area. There are four main metamorphic zones in the area (Pogson and
Watkins, 1998). The following is the four summary metamorphic zones:
M1 (blue): clay/carbonate, epidote, sericite, albite, with or without foliation.
M2 (green): prehnite/pumpellyite, normally without foliation.
M3 (purple): actinolite/biotite, with or without foliation.
MH (red): akin to the hornfels envelope surrounding an intrusion. This zone is
represented by the presence of metamorphic diopside or garnet, or other
phases such as cordierite, andalusite, or even other higher grade hornfels
assemblages.
The metamorphic map has been prepared from an assessment of the
metamorphic mineral assemblages of 1919 rocks corrected by the Geological
Survey of New South Wales from this sheet. Because of small data set, this
map is regarded as a preliminary attempt to reassess the metamorphic zones
published by (Smith, 1969). His work identified and showed the distribution of
six zones of differing metamorphic grade over the Buthurst-Dubbo region,
based on the presence or absence of key metamorphic minerals such as
carbonate, albite, chlorite, prehnite, pumpellyite, actinolite-tremolite and biotite.
Metamorphism can be seen adjacent to the intrusion on the Dubbo map. All the
Carboniferous plugs and plutons are altered to an M1 assemblage such as that
seen in the Yeoval Complex.
103
Figure 5 - 12: Dubbo 1:250,000 map sheet area showing the metamorphic zones (Meakin and Morgan, 1999).
104
5.5 Local Geology
5.5.1 Introduction
The study area is located in the Spicers Creek catchment within the Cobbora
1:100,000 map sheet area on the Dubbo 1:250,000 map (Figure 2-4). Spicers
Creek catchment is situated about 45 km from the city of Dubbo to the west, 20
km from the town of Cobbora to the east, 15 km from the town of Boomley to
the north and south 45 km to Wellington (Figure 5-13).
Cobbora region is situated in the north-eastern portion of the Lachlan Fold Belt
(LFB) and overlying marginal areas of the Gunnedah Basin and Surat Basin.
The basement rocks of the Lachlan Fold Belt have been regionally
metamorphosed and contain north-south trending faults, folds, and fractures.
The Gunnedah Basin Sequence includes a number of basic intrusions of
Mesozoic and Tertiary rocks. These are associated with massive extrusions of
Garrawilla Nolanic complex and the Liverpool, Warrumbungle and Nandewar
Ranges (Tadros, 1993). The Surat Basin contains 2500 metres of mainly
Jurrasic continental sedimentary rocks and lower Cretaceous marine beds
largely obscured by Cainozoic alluvium (Exon, 1976).
The area of study, is geologically complex due to the tectonic events it has
endured and the landscape evolution is closely related to the underlying
geology (Morgan and Jankowski, 2002). The Palaeozoic rocks range in age
from Ordovician through to Devonian. The depositional environments ranges
from marine to sub-marine (Schofield and Jankowski, 2000). Lithologies 105
include submarine felsic, rhyolitic tuffs, tuffaceous sedimentary rocks and latitic
lavas and intrusives (Meakin and Morgan, 1999). Permian to Triassic
sediments contain marine and fluvial sediments with some volcanics. Jurassic-
Cretaceous sediments consist of terrestrial-sourced shales, coals, and well-
sorted sands. Depositional environments include floodplains, meandering
streams and high energy braided streams. These sedimentary cover rocks dip
in a northwesterly direction and Tertiary basalts intrude the sediments as sills,
plugs, dykes, and relic basalt flows. Quaternary alluvial and collovial sediments
blanket drainage lines and low relief areas (Morgan and Jankowski, 2002).
Figure 5 - 13: The road map of the Dubbo area.
106
5.5.2 Climate
In the Spicers Creek catchment the highest temperatures are January with
maximum mean of 33˚C and the lowest temperatures are in July with a
minimum mean of 1.8˚C. The mean average rainfall for the catchment was
measured as 563 mm. Summer is the wettest season receiving 31% of the
annual rainfall (Jewel, 2000, Morgan and Jankowski, 2002).
5.5.3 Stratigraphy
The rock units that are represented in the Cobbora area are categorized into
those relating to basement and those that constitute the cover rocks. The
Ordovician to Late Devonian, mainly marine basement rocks were deformed
during compressive Middle Devonian and Middle Carboniferous orogenic
episodes, whilst the mainly terrestrial Late Carboniferous to Late Jurassic cover
rocks were deposited in an extensional tectonic regime.
Figure 5-14 shows the geological map around the area are scale of 1:100,000.
Table 5-1 is a summary of the different rock types found in the area.
Throughout the following sections the characteristics of the geological units
around the study area are documented (Meakin and Morgan, 1999).
107
Figure 5 - 14: Geological map of study area from Cobbora 1:100,000 geological map with the two sites marked. The various units are described in Table 1.
108
Table 5 - 1: Description of different rocks type in the Spicers Creek area (From Cobbora 100,000 geological map).
Symbol from
fig 2-14
Rock
Formation
/Group
Age
Composition
Qa
--- Quaternary Alluvial silt, clay and sand, variable humic content,
sporadic pebble to cobble-sized unconsolidated
conglomeratic lenses.
Rp Napperby
Formation
Triassic Siltstone thinly interbedded with fine-medium grained,
Silurian Rhyolitic to latitic lava, intrusive and tuff and
volcaniclastic sandstone.
θm Ungrouped
Ordovician
Intrusions
Late
Ordovician
Monzonite and diorite to more felsic syenite and
granite.
θco Oakdale
Formation
Ordovician Basalt, basaltic andesite, latite lava and intrusions,
volcaniclastic breccia, conglomerate, sandstone and
siltstone, minor allochthonous limestone.
109
5.5.3.1 Palaeozoic units
The Palaeozoic rocks range in age from Ordovician through to Devonian. The
early Ordovician sequence represents the oldest strata on the Dubbo 1:250,000
map sheet area and the study area.
5.5.3.1.1 Oakdale Formation (θco)
The Oakdale Formation is one of the most widespread units of the Cabonne
Group. The Formation takes it name from the Oakdale property, at the
southern end of the Oakdale Anticline (Meakin and Morgan, 1999). It forms a
broad, structurally complex, largely fault–bounded belt about 120 km in strike
length which extends from 35 km north of Wellington to 20 km south-east of
Orange.
Constituent units contain a number of allochthonous, and possibly rarely
autochthonous, limestone bodies, as well as limestone breccia and calcareous
sandstone. Primary volcanic rocks of the Oakdale Formation include
predominantly shoshonite, basalt to basaltic andesite and less commonly latite
lavas and high level intrusions. Volcanoclastic breccia and conglomerate are
common in the Oakdale Formation, particularly in the area to the north of
Wellington. The units are interbedded with sandstone turbidite packages and
primary volcanic layers.
110
The units are massive to coarsely bedded, clast- to matrix-supported, poorly
sorted and range from 3 m to over 20 m in thickness. This formation is
unconformably overlain in the main belt by late Early to Late Silurian Mumbil
Shelf strata, represented south of Wellington by the Dripstone Formation, and
north of Wellington by the Narragal limestone and Gleneski Formation of the
Mumbil Group.
In the Oakdale Formation, graptolite and brachiopod fauna have been
described by a number of researchers from the siltstone units (Meakin and
Morgan, 1999).
5.5.3.1.2 Ungrouped Ordovician Intrusions (θm)
Many small intrusive bodies of presumed Ordovician age occur within the
Oakdale Formation. None of these units have been isotopically dated. The
intrusions range from monzonite and diorite to more felsic syenite and granite in
composition and the outcrops range in size from less than 100 metres up to 900
metres in diameter.
Metamorphism by these intrusions is generally limited to weak to moderate
development of chlorite-sericite-epidote alteration of plagioclase and
ferromagnesian minerals, indicating a lower greenshist grade (Meakin and
Morgan, 1999).
111
5.5.3.1.3 Gleneski Formation (Sms)
The formation takes it name from the Gleneski property north of the Macquarie
River (Meakin and Morgan, 1999). The Gleneski Formation is comprised
predominantly of felsic to rhyolitic tuff and tuffaceous sedimentary rocks, with
lesser rhyolitic to felsic to latitic lava and intrusions, shale, and minor limestone.
Massive felsite, quartz felsite, quartz latite and rhyolite lava and intrusions
locally are intercalated with the volcanoclastic rocks of the Gleneski Formation.
These lithologies are interbedded on a variety of scales. Rare limestone and
calcareous sandstone are interbedded with volcanoclastic rocks of the Gleneski
Formation.
The massive and graded beds, erosional bases and local cross bedding in
clastic rocks of the Gleneski Formation suggest that they were deposited
subaqueously by turbidity currents. The magma source for the Gleneski
Formation was acid, but both primary volcanic and pyroclastic rocks, as well as
volcanic and lithic fragments in tuffaceous units, are deficient in quartz
phenocrysts (Meakin and Morgan, 1999).
5.5.3.1.4 Cuga Burga Volcanics (Dgc)
The name is derived from Cuga Burga Creek, which cuts through the formation
to the south-east of the Oakdale Anticline (Meakin and Morgan, 1999). The
Cuga Burga Volcanics contains no autochthonous fossils and age constraints
for the formation are provided mainly by conodont data for overlying and 112
underlying formations. It consists mainly of shoshonitic latitic volcanic and
volcanoclastic rocks. A thick sequence of pillow lavas is well-exposed in the
representative section of this units. A low radiometric character characterizes
the Cuga Burga Volcanics, although locally, the formation has a strong
potassium (pink) response.
The belt of Cuga Burga Volcanics separating the Yeoval Complex from the
Yennora Granite in the western area shows a complex intrusive history. Many
dykes are present, varying in composition from microdiorite and graphic
microgranite to dacite and rhyolite to basalt. The Cuga Burga Volcanics
contains no authochthonous fossils and age constraints for the formation are
provided mainly by conodont data for the overlying and underlying formations
(Meakin and Morgan, 1999).
The Cuga Burga Volcanics is characterized by a low radiometric character (dark
signature on radiometric images), although locally, such as in the fault-bounded
belt north of The Gap, the formation has a strong potassium (pink) response.
The formation has an overall anomalous magnetic high signature on total
magnetic intensity images. Igneous rocks of the formation have a high magnetic
susceptibility while volcaniclastic rocks are more variable (Meakin and Morgan,
1999).
113
5.5.3.1.5 Early Permian Undifferentiated (Pe)
This unit is largely exposed in a series of north-west trending lobes which
extend northwards under the Late Permian to Jurassic sedimentary rocks of the
Gunnedah Basin and the overlying Great Australian Basin. It contains
carbonaceous siltstone, quartz lithic sandstone, conglomerate coal lenses and
rare varves. This unit is somewhat better exposed in the Spicers Creek. In the
Spicers Creek the basal Early Permian sequence consists of poorly sorted,
sedimentary polymictic breccia and conglomerate which contains subangular to
rounded clasts of sandstone, phyllite and quartz which range up to 40 cm in
size.
The undifferentiated Early Permian, a possible distal equivalent to the Rylstone
Volcanics, correlates with the Late Carboniferous to earliest Permian, thick non-
marine fluvio-glacial conglomerates of the Talaterang Group in the southern
Sydney Basin, and the fluvio-lacustrine sediments of the Seaham Formation
(Meakin and Morgan, 1999).
5.5.3.1.6 Dunedoo Formation (Pd)
The Dunedoo Formation forms a relatively thin unit in the area of study. It
unconformably overlies the undifferentiated Early Permian (Pe) sedimentary
rocks and is disconformably overlain by Early Triassic sedimentary rocks. The
Dunedoo Formation is only moderately exposed at scattered localities, mainly in
114
creeks or occasionally on steeper slopes immediately below the more resistant
overlying Boulderwood Formation.
The formation consists of breccia, conglomerate and fine – to coarse- grained
sandstone and pebbly sandstone of quartzose to quartz-lithic and
quartzofeldspathic composition with a dominantly kaolinitic matrix. The lower
section of this formation is generally only poorly exposed. Surface exposures of
the Dunedoo Formation range up to 60 m on the Cobbora 1:100,000 map sheet
area. Common fining-up sequences in the Dunedoo Formation reflect fluvial
deposition (Meakin and Morgan, 1999).
5.5.3.2 Mesozoic units
Igneous rocks of Mesozoic age are scattered across Dubbo 1:250 000 map
sheet area, and are commonly associated with Permian to Triassic sedimentary
rocks of the Sydney and Gunnedah Basins. The Mesozoic rocks in the Spicers
Creek area are limited to the Boulerwood Formation and the Napperby
Formation. They cover of a large part of the north of area.
5.5.3.2.1 Boulerwood Formation (Rb)
The name is derived from Boulderwood homestead (Meakin and Morgan,
1999). It is typically exposed at the margins of slightly elevated plateaus
capped by the early Jurrasic Purlawaugh Formation. This formation ovelies the
115
Dunedoo Formation with probably disconformity or slight angular unconformity
and is conformably overlain by the Napperby Formation.
The Boulderwood Formation consists of very coarse-grained pebbly to
conglomeratic, lithic-quartz sandstone and grit, together with subordinate flaggy,
planar laminated, white, medium-grained quartz sandstone, claystone and
siltstone. Clasts are subangular to subrounded and predominantly quartz, with
minor fine silicic volcanic fragments and other siliceous metasedimentary
basement lithologies. The formation was deposited by a fluvial system.
Channel floor, channel bar and point bar sequences are common (Meakin and
Morgan, 1999).
5.5.3.2.2 Napperby Formation (Rp)
The name is derived from calcareous sandstone and laminated mudstone
exposed on Napperby station, near the Pine Vala well on the Coonabarabran-
Mullaley road (Meakin and Morgan, 1999). It conformably overlies the
Boulderwood Formation.
The formation consists of white, fine – to medium-grained, moderately to poorly
sorted lithic-quartz and quartzose, often flaggy, ferruginous sandstone with a
white clayey matrix thinly interbedded or interlaminated with grey siltstone and
minor conglomerate lenses. The Napperby Formation represents lacustrine-
delta deposition (Meakin and Morgan, 1999).
116
5.5.3.3 Cainozoic Unit
In the area of study the Cainozoic deposits are represented by alluvial deposits.
5.5.3.3.1 Alluvial Deposits (Qa)
No formal names have been applied to unconsolidated Quarternary units.
Alluvial deposits are geomorphological units. They can be regarded as falling
into three distinct morpho-stratigraphic types: actively depositing Holocene
zones; inactive Holocene/Pleistocene zones; and relict/elevated Pleistocene
zones. Generalized lithologies for the alluvium consist of silt and clay, with a
variable humic content.
5.5.4 Structure
The Spicers Creek region is situated in the east of the Dubbo 1:250,000 map
sheet and in the west of the Cobbora 1:100,000 sheet. It includes four of the
structural entities of the LFB; including the Cowra Trough, Molong High, Hill End
Trough and Capertee High (see Figure 5-10). Scott (1997) re-defined the
structural entities of the northeastern LFB into the four zones including the
Cowra-Molong Zone, Nindethana Zone, Hill End Zone and Capertee Zone
(Figure 5-15). He proposed there were three deformation events affecting rocks
across the map sheets. These are summarized in Table 5-2.
117
Site 2
Site 1Dubbo
Figure 5 - 15: Structural zones defined by Scott (1997) across the Dubbo geological map sheet.
D1 deformation in the late Middle Devonian was related to N-S compression
that produced east-trending folds and thrusts. The Ordovician and Silurian
strata were uplifted during this deformation because of thrusting them over early
Devonian sediments. This was followed by erosion as the elevated topography
was flattened. D2 deformation commenced in the early Carboniferous through
E-W shortening. It was a major deformation with an early compressional phase
which produced the dominant axial cleavage forming deformation. This
118
deformation varies in intensity style between the four structural zones as
described below.
Table 5 - 2: Proposed deformation events affecting rocks across the Dubbo 1:250,000 geological sheet Scott (Morgan, 1997).
Deformation Age Form Direction of Compression
D1 late Middle Devonian Folding and thrusting N-S compression
D2(major deformation) Early Carboniferous Cleavage, folding and E-W compression thrusting
D3 Early Carboniferous Kinking N-S compression
In the Cowra-Molong zone, the D2 deformation is represented by west-dipping
thrusts and east-dipping normal faults, with N-NE trending folds with near
vertical to west-dipping axial surfaces. The presence of small-scale synthetic
and antithetic strike-slip faults was probably associated with major strike-slip
faulting along the peripheries of the zone during the Early Carboniferous.
D2 deformation in Hill End Zone is not intense. Compression varied in
orientation from ENE-WSW in the north of the zone to ESE-WNW in the
southern part of the zone. This deformation is characterized by continuous,
shallowly plunging, open to tight folds with four-way dip closure. Folds have a
well-developed, axial plane cleavage.
119
In Nindethana Zone, D2 deformation is represented by E-NE dipping thrusts
and normal faults that cause repetition of the eastward-younging strata.
Deformation in the west of zone is higher than in the east part. Scott (Morgan,
1997) suggested that the zone was initially subjected to E-W compression
which produced north-trending faults. Subsequent ENE-WSW compression
resulted in dextral shear along these fault creating NW trending folds. In this
zone northwest trending folds have near vertical, east-dipping axial surfaces,
and are oblique to the north-trending faults.
D2 deformation in Capertee zone is characterised by west-dipping thrusts which
repeat the west-dipping stratigraphy. The direction of compression was ENE-
WSW, resulting in some dextral strike-slip faulting. Folds are open to tight with
an N-NE axial surface direction and have a well developed axial surface
cleavage.
D3 Early Carboniferous deformation involved N-S compression and
development of large scale easterly trending kink axes. There was no major
deformation after D3, but lineated faulting and cleavage development in some
parts the Wuuluman Granite and the Gulgong Granite. In the north of the map
sheet the sedimentary basin reflects the reactivation of the NNW trending faults,
with thick fault bound Permian outliers showing the significance post-Permian
faulting.
In the area of the study within the Molong Zone there are many faults and
thrusts striking north-south (Figure 5-16). These include the Nindethana Fault,
120
Narragal thrust, Neurea Fault and Macquarie fault. Throughout the following
sections the characteristic of the faults near the study area are documented.
5.5.4.1 Nindethana Fault
The Nindethana Fault is a major 50 km long thrust fault that extends from the
northern part of the Bathurst 1:250,000 map sheet area to the south-western
corner of outcrop of the Wuuluman Granite. The Nindethana Fault, for most of
its length, is localized within the Early Devonian Cunningham Formation.
However, west of Lake Burrendong, rocks east of the fault represented the Late
Silurian Mumbil Group and the overlying Early Devonian Cuga Burga Volcanics
in the cores and limbs of meridional anticlines (Meakin and Morgan, 1999).
Direct exposure of the Nindethana Fault occurs around Galwadgere (GR
692200 - 6383750) where it represents a 100 m wide shear zone, with strong
pyrite lineation. The dip of the fault is estimated to be 45-50˚ to the east
(Meakin and Morgan, 1999). The fault loses displacement northwards, so that
immediately south-west of the Wuuluman Granite, the footwall consists of Cuga
Burga Volcanics rather than the Cunningham Formation (Meakin and Morgan,
1999). The fault terminates to the north against a west-north-west trending
cross fault that forms part of a west-north-west corridor.
121
Figure 5 - 16: Structural map of the thrust sheet on the Cowra Zone and Molong Zone of the Dubbo geological sheet (Modified from Meakin and Morgan, 1999).
122
5.5.4.2 Narragal Thrust Sheet
The Narragal Thrust Sheet represents a combination of regional folds and
thrusts, and thus shows a structural transitional zone between the Molong Zone
to the west which is dominated by thrusts, and the Hill End Zone to the east that
is dominated by folds (Scott, 1997). This thrust sheet is bounded to the west
(below) by the east-dipping Narragal fault that marks the roof of Neurea Thrust
Sheet and to the east (end below) by the Nindethana Fault that also dips east.
The Narragal Thrust Sheet thus lies between two parallel thrusts (Meakin and
Morgan, 1999).
Cleavage in the Narragal Thrust Sheet dips steeply to the east and west and
swings in strike from NNW in the southern third of the map sheet area to NNE in
northern two third where it is approximately parallel in strike to fold axial traces.
The general obliquity with the meridonal bounding faults suggests a component
of left-lateral movement on them during deformation (Meakin and Morgan,
1999).
5.5.4.3 Macquarie Fault
The Macquarie Fault has a NNW trend and east dip. Its southern extent is
uncertain and it may disappear within the Garra Formation or link into a WNW
cross fault that cuts off the northern end of the Neurea Fault. This fault
juxtaposes Ordovician rocks in the hanging wall against the Late Devonian
Catombal Group in the Mount Arthur Syncline in the foot wall. The fault is
interpreted as one of a series of east-dipping thrusts with lesser displacement
123
and different geometry. The displacement history of these faults is uncertain.
They could represent either normal faults with varying displacement from north
to south or strike-slip fault (Meakin and Morgan, 1999).
5.5.4.4 Neurea Fault
The Neurea fault belongs to Neurea Thrust Sheet. The Neurea Thrust sheet is
bounded to the west by the Neurea Fault, which is interpreted to be an east-
dipping, oblique thrust. It is also bound to the east by the Narragal Fault,
interpreted to be an east-dipping thrust that juxtaposes Late Silurian Mumbil
Group and younger rocks in the Narragal Thrust sheet over Late Ordovician
Oakdale Formation (Meakin and Morgan, 1999).
In the Neurea Thrust Sheet some fold plunges are steep. These steep plunges
and the presence of strike-slip fibres/striae (plunging gently north) on steeply
east-dipping fault surface in the volcanics are interpreted in terms of strike-slip
movement on faults superimposed on vertical (east-side-up), dip-slip, thrust
movement.
124
5.6 Regional and Local Geophysics
5.6.1 Introduction
The Australian Geological Survey Organisation (AGSO) has acquired regional
geophysical data consisting of airborne radiometric and magnetic data and land
based gravity data. The geophysical data was provided as ERMapper grids,
one for both the Dubbo and Cobbora 1:100,000 sheets and a third grid for the
Gilgandra 1:250,000 sheet (Figure 5-17). The gravity data was provided as
Bouguer anomaly data corrected to a density of 2.67 tonnes/m.
Figure 5 - 17: Location of the geophysical data sets supplied by AGSO and the selected area of interest for the interpretation of this data (Schofield, 1998).
125
The radiometric and magnetic data for the Gilgandra Map Sheet was acquired
in 1979. Its resolution is significantly lower than that of the data from the Dubbo
Map Sheet which was flown as part of the National Geoscience Mapping
Accord (NGMA) in 1991. The survey parameters are listed in Table 5-3.
Table 5 - 3: Survey parameters for the Gilgandra and Dubbo airborne surveys.
Gilgandra Survey Dubbo Survey
Operator New South Wales Department New South Wales Department
Elements whose atomic nuclei contain the same number of protons but different
numbers of neutrons are called isotopes. They are forms of the same element
with different atomic weights. Certain nuclei disintegrate spontaneously
emitting α particles (helium nuclei He) and β particles (electrons and positrons).
This is the phenomenon of radioactivity. These emissions alter the nuclear
charge, α by –2, β (positron) by –1 and β (electron) by +1. As a result the
disintegrating nucleus is transformed into a nucleus of another element.
The natural decay of radioactive elements produces a variety of types of
radiation (alpha, beta, gamma) at specific energy levels. Only gamma-ray
radiation has sufficient energy to be useful for geological mapping or
exploration. There are a large number of radioactive minerals, but the more
common are potassium (K), uranium (U), and thorium (Th). The gamma
radiation is emitted when these elements decay gradually through time (see
Telford et al., 1976), for a description of the individual decay process and
series).
Radiometric surveying is employed in the search for radioactive deposits and
also for deposits associated with radioactive elements such as zirconium and
titanium (Kearey and Brooks, 2001). Radiometric surveys used in geological
mapping as different rock types can be recognized from their distinctive
signature (Pires and Harthhill, 1989). It is an important source of information for
127
soil, regolith and geomorphological studies, including their mineralogy, texture,
chemistry and style of weathering (Wilford et al., 1997).
As Jaques et al. (1997) demonstrated gamma-ray spectrometry is an important
tool for distinguishing igneous rock bodies because of composition and degree
of differentiation, as U, Th and K. These elements are critical elements involved
in magmatic differentiation processes.
5.6.2.2 Interpretation Methodology
Airborne radiometric methods measure gamma radiation emitted from the
Earth’s surface and it is recorded using a gamma-ray spectrometer. The most
diagnostically valuable of these for geological mapping correspond to potassium
(1.46 MeV), bismuth-214 (1.765MeV) from the uranium-238 decay series, and
thallium-232 (2.614 MeV) from the thorium-232 decay series (Horsfall, 1997).
When K, U and Th decay they create gamma radiation of differing energy, K
producing the lowest and Th the highest (Table 5-4).
The sensor detects the energy levels of the gamma spectrum, and the relative
contribution of each band is recorded as a separate channel along with a total
signal channel (Horsfall, 1997). Modern airborne gamma-ray data acquisition
for regional mapping consists of a multichannel spectrometer be able to
measure at least 256 channels of data in the energy range 0-3 MeV (Minty,
1997).
128
Table 5 - 4: IAEA (International Atomic Energy Agency) recommended windows for conventional 3-channel airborne gamma-ray spectrometry (IAEA, 1991).
Element analyzed Isotope used Gamma ray energy Energy window
Potassium K 1.46 MeV 1.370-1.570 MeV
Uranium Bi 1.76 MeV 1.660-1.860 MeV
Thorium Tl 2.61 MeV 2.410-2.810 MeV
The three bands of K, U and Th can be interpreted individually or combined as
a red-green-blue (RGB) composite image (Milligan and Gunn, 1977). The mix
of the three bands as the three primary colours produces a variety of colours,
which help in the interpretation of the relative concentrations of K (red), Th
(green) and U (blue) in the rock (Table 5-5). In creating images, the range of
each of the three colours is divided into 256 levels. For example, in RGB
images, the brightness of a particular colour equates to its relative abundance –
a dark red colour indicates a high level of potassium, with little or no thorium or
uranium. Furthermore, another factor for interpreting is using a cylindric
coordinate system in terms of colour variation perceived by the human eye as
the hue, saturation, and value model (HSV; Milligan and Gunn (1977)).
In the HSV model, hue indicates the combinations of primary colours (RGB),
value is the intensity (or energy of colour), and saturation is the relative lack of
white in the colour (Milligan and Gunn, 1977). Figure 5-18 shows the HSV
colour model that hue is measured around the vertical axis from 0˚ (red) to 360˚,
saturation varies from 0 on the vertical axis to 1 on the triangular surfaces of the
129
hexagonal, and value varies from 0 (black) through shades of grey along the
central axis to 1 (white) at the top. ( Milligan and Gunn, 1977).
Table 5 - 5: Interpretation of relative concentration of K. Th, and U in rocks based on the RGB image (Morgan, 1997).
K Th U Resultant colour
High Low Low Red
Low High Low Green
Low Low High Blue
High High Low Yellow
High Low High Purple
Low High High Light Blue
High High High White
An interpretation of gamma-ray data needs to incorporate many variables.
Besides the geometry and physical property contrasts of the radioactive source,
environmental effects can influence the interpretation of data. Environmental
factors include air temperature and pressure, soil moisture, non-radioactive
overburden, rainfall, vegetation and precipitation (Minty, 1997). Barren
overburden can reduce the radiation output from the Earth’s surface because of
its high density. In some areas, dense vegetation may shield the source of
radiation which would be equivalent to 50 m of air. The trunks of trees in dense
130
forests also reduce the effect of radiation from the ground. Changing
temperature and pressure can lead to a change in air density by up to 30 per
cent. Rainfall and soil moisture decrease the effect of radiation from the surface
as well (Minty, 1997)
Figure 5 - 18: The HSV colour model – H, S, and V represent hue, saturation, and value of colour intensity, respectively. According to this model, for a colour of constant hue and saturation, if its value is decreased it darkens towards black. For a colour of constant hue and value, if the saturation is decreased it becomes whiter (Milligan and Gunn, 1977).
131
5.6.2.3 Regional Interpretation
Figure 5-19 shows the RGB ternary radiometric image for the Dubbo area. It
shows that over the Dubbo 1:250,000 map sheet area the Ordovician volcanic
rocks are the generally low radioactivity, with slightly higher levels of potassium
parting a characteristic dark red to red appearance in RGB images. The
Silurian sedimentary rocks of the Cowra Trough show commonly of low to
moderate radioactivity, with thorium radiation dominant. This gives a
characteristic dark green anomaly in the RGB radiometric image. The
sedimentary rocks of the Hill End Trough represents the low to moderate levels
of radiation, with elevated thorium radiation. The Silurian sedimentary rocks of
the Capertee High appear in the RGB images as areas of dark green, indicating
low levels of radiation with slightly elevated levels of thorium radiation (Meakin
and Morgan, 1999).
The Devonian sedimentary rocks produce little radioactivity, resulting in dark
areas on the radiometric RGB image. For example, the Devonian Cuga Burga
Volcanics produces extremely low levels of radioactivity, giving a characteristic
black anomaly on RGB images. The Carboniferous intrusions over the Dubbo
area (such as the Gulgong Granite) exhibit moderate to high potassium
radioactivity with red colour in the RGB image. The northern outcrop of the
Gulgong Granite (East of Dubbo) shows high levels of potassium and thorium
radiation, giving a characteristic yellow colour in the RGB images. Permian to
Mesozoic sedimentary and Tertiary volcanic rocks generally emit low to
moderate levels of radioactivity.
132
In the north-eastern part of Dubbo 1:250,000 map sheet, the basalts may be
mapped from the radiometric images but in the north-western part, the generally
lower radioactivity levels emitted by the Great Australian Basin (Surat Basin)
sedimentary rocks makes it more difficult to distinguish between the
sedimentary rocks and the basalts (Meakin and Morgan, 1999).
As Jaques et al. (1997) showed there are some associated intrusive bodies in
the south of Dubbo ( probably Mesozoic age) with the development of alteration
and rare-metal mineralisation such as Zr, Hf, Nb, Ta, light REE and Y. The
gamma-ray image shows that each intrusion has a different geochemical
signature. Some have dominant U+Th signatures, others dominant K+U or
K+Th, and some are high in all three elements (Figure 5-20).
133
Figu
re 5
-19:
Red
-gre
en-b
lue
com
posi
te ra
diom
etric
imag
e fo
r Dub
bo a
rea
(Mea
kin
and
Mor
gan,
199
9).
134
Figure 5 - 20: Gamma-ray image of an area south of Dubbo, showing the geochemically different signatures of a number of minor intrusions in the area. Some intrusions have dominant U+Th signatures, others dominant K+U or K+Th, and yet others are high in all three elements (Jaques et al., 1997).
135
5.6.3 Aeromagnetic Methods
5.6.3.1 Theory, Application
Magnetic surveys can be carried out in the air from helicopters or fixed wing
aircraft, because large areas can be covered in great detail and relatively
inexpensively. The operation is known as airborne magnetic surveying or
aeromagnetic surveying. Magnetic survey can also be carried out on the ground
for more detail. The instrument used to measure the magnetic field is called a
magnetometer.
Airborne magnetic surveys are used extensively to map regional geology in
areas of poor outcrop. The usual procedure is to record data continuously
along equally spaced parallel flight lines covering the survey area. In Australia,
airborne surveys are flown routinely at line spacings of 400-500 m and in areas
where greater detail is required, flight-line spacing is reduced to 200 m. For
regional mapping the flight line direction is usually oriented orthogonal to the
predominant strike of the geology. For more specific applications, such as
mineral exploration targets, the flight-line separation and direction will be
selected to maximise the magnetic signature (Horsfall, 1997).
Ferromagnetic minerals in the earth’s crust are the source of anomalous
magnetic field. The major magnetic minerals are magnetite, titanmagnetite,
titanhematite, maghemite, pyrrhotite and native iron. These minerals can show
magnetic anomalies, because of their abnormally large magnetic susceptibilities
or because they have high remnant magnetizations. Of the magnetic minerals
136
that are found in the nature, magnetite is the most abundant (Grant, 1984).
Igneous rocks have a higher content of magnetic minerals than sediments and
can be mapped in sedimentary basins from magnetic data (Gunn, 1997a).
The contrasting proportions of these minerals in different crustal rocks produce
the magnetic anomalies that are the targets of magnetometer surveying.
Airborne magnetic surveys measure the total magnetic field that includes
induced magnetism and remnant magnetism. Any rock containing magnetic
minerals may possess both induced and remanent magnetizations (Gunn and
Dentith, 1997).
When a material is placed in a magnetic field it may acquire a magnetization in
the direction of the field which is lost when the material is removed from the
field. This phenomenon is known as induced magnetization or magnetic
polarization. The inherited magnetization remaining after removal of the applied
field is known as remanent, or permanent, magnetization (Kearey and Brooks,
2001).
The magnetic field decreases approximately as the inverse of the square of the
distance from the magnetic source. Therefore to record small variations in the
field, surveys must be flown close to the ground. The survey attitude is usually
in the range 80-60 m. Surveys with 200 m line spacing are typically flown 80 m
above ground while surveys with 100 m line spacing may be flown at 60 m
above ground (Horsfall, 1997).
137
5.6.3.2 Interpretation Methodology
Rock magnetism results from the magnetic minerals they contain. Magnetic
minerals also belong to iron-titanium-oxygen solid solution series from
magnetite (FeO) to ulvospinel (FeTiO) or the iron-sulphur group pyrrhotite
(FeS). Magnetite can occur as either primary magmatic minerals or as
metamorphic minerals. Basic igneous rocks are usually highly magnetic due to
their relatively high magnetite content. In general magnetite concentrations
decrease with increasing acidity in magmatic rocks (Kearey and Brooks, 2001).
Magnetic anomalies range in amplitude from a few tens of nanoteslas over
deep metamorphic basement to several hundred nanoteslas over basic
intrusions and may reach an amplitude of several thousand nanoteslas over
magnetite ores. The magnetic field for typical magnetic bodies in the
subsurface decays rapidly away from its source (Kearey and Brooks, 2001). As
a result shallow bodies are subjected by high wavenumber components,
creating narrower relatively low amplitude anomalies, whereas deep bodies
produce broader, smoother anomalies (Schofield, 1998).
Ferromagnetic rocks or sediments show positive magnetic susceptibility
whereas diamagnetic rocks or sediments show negative magnetic susceptibility.
For example pure salt comprising hatite, gypsum and anhydrite, is diamagnetic
and it has a negative magnetic susceptibility and as a result represent a
negative magnetic contrast relative to surrounding sediments (Gunn, 1997b).
138
The aim of magnetic interpretation is to obtain geological information from
magnetic survey data (Gunn, 1997b). The results of magnetic survey are
usually presented as a map, or an image, of the magnetic anomalies. The
magnetic anomaly is the product of the susceptibility of the rock and the size of
the body or rock. A colour total magnetic intensity (TMI) image provides a
useful method for interpreting geological units and their structure. Boundaries
between different geological units are determined on the basis of their magnetic
response. Common causes of magnetic anomalies include dykes, faulted,
folded, magnetite ore bodies and metamorphic basement rocks (Kearey and
Brooks, 2001).
There are some environmental factors that can influence magnetic anomalies
include temperature and weathering. It is unusual to find magnetite in
unmetamorphosed sedimentary rocks except in very minor amounts. Magnetite
is unstable in the low-temperature, highly oxidizing environment of chemical
weathering and sedimentation. Four factors effect the chemical weathering and
sedimentation of iron covsist of chemical weathering, leaching, transportation
and compaction and early diagenesis (Grant, 1984). Figure 5-21 presented the
total magnetic intensity (TMI) image for the Dubbo region.
5.6.3.3 Regional Interpretation
Sedimentary rocks over the Dubbo area show low to medium magnetic
susceptibility, but contain indications of internal north-south structure. In the
western part of Dubbo, the Silurian sedimentary rocks in the Cowra Trough 139
generally show low magnetic susceptibility. The magnetic susceptibility of the
Silurian volcanic rocks in this area is commonly low (range 10x10SI to
25x10SI). Silurian rocks on the Capertee High represent low magnetic
susceptibility (<20x10SI) that is reflected in low-amplitude magnetic anomalies
(Meakin and Morgan, 1999).
Devonian sedimentary rocks crop out over large portions of the Dubbo
1:250,000 map sheet area. They have low magnetic susceptibility (<20x10SI)
in relation to four structural zones including Cowra Zone (Cowra Trough),
Molong Zone, Hill End Zone and Capertee Zone (Meakin and Morgan, 1999).
The Carboniferous intrusions (e.g. Wuuluman Granite and Gulgong Granite) are
generally in the low to moderate range for magnetic susceptibility (between
20x10SI to 2000x10SI). The magnetic susceptibility values measured for the
Permian and Mesozoic sedimentary rocks of the Sydney Basin, Gunneda Basin
and Surat Basin are generally very low (<10x10SI). With low contrasts in
magnetic susceptibility, they do not produce significant magnetic anomalies
(Meakin and Morgan, 1999). However, the magnetic susceptibility values for
the Tertiary volcanic rocks are in the medium range (300x10SI to 1500x10SI).
Flows In the north-eastern part of Dubbo and in the west of Dubbo overlie
sedimentary rocks and older volcanic materials respectively, which have
moderate to high magnetic susceptibility values. Therefore, it is difficult to
distinguish the anomaly pattern of the Tertiary basalts from those of the
underlying rocks (Meakin and Morgan, 1999).
140
Survey area
Figure 5 - 21: TMI aeromagnetic image for the Dubbo region ( Modified from Meakin and
Morgan, 1999).
141
5.6.4 Gravity Methods
5.6.4.1 Theory, Application
Australian Geological Survey Organisation (AGSO) carried out numerous
traverses throughout the Dubbo and Ballimore region. The stations are on
average about 4 km apart. Figure 5-22 shows that the distribution of the
stations is not regular and accordingly parts of the area have dense coverage
while in others parts coverage is quite sparse. The gravity data are based on
the IGSN71 gravity datum and the ISOGAL84 scale.
Figure 5 - 22 Point located data for the Dubbo and Ballimore regions (Schofield, 1998).
142
Gravity surveying measures variations in the Earth’s gravitational field produced
by differences of density between subsurface rocks. The theory that explains
gravitational methods is based on Newton’s law. This law is expressing the
force of mutual attraction between two particles in terms of their masses and
separation. The force of attraction F between two masses m and m, that
dimensions are small with respect to the distance r between them, is given by
221
rmGmF =
where G is the Gravitational Constant.
The variation in gravity depends on lateral changes in the density of earth
materials in the vicinity of the measuring point (Dobrin and Savit, 1988). Gravity
method is a natural field method, which measure a bulk petrophysical property
such as density. Density depends on both composition and porosity. In
sedimentary rocks, porosity is the main cause of density variations and density
tends to increase with depth because of compaction, and with age due to
cementation. Igneous and metamorphic rocks have small porosity, and
composition is the main cause of density variation. In igneous rocks density
increase as acidity decreases.
Like radioactivity and magnetic methods, gravity is a natural source method in
which local variations in density of rocks near the surface cause small changes
in the main gravity field. Some applications of gravity surveying are in regional
Figure 6 - 1: Model with a plane horizontal ground surface and a highly irregular refractor. The vertical and horizontal scales are equal (Palmer, 1980).
First Arrival Times Graph
0
20
40
60
80
100
120
0 50 100 150 200 250 300 350
Distance (m)
Tim
e (m
s) Forward
Reverse
Figure 6 - 2: Traveltime curves derived from the model in Figure 6 - 1.
153
The velocity analysis functions for XY-values ranging from 0 to 30 are shown in
Figure 6-3. They have been stacked in order to clearly show the minor
differences between each function (Figure 6-4). Detailed study of the velocity
analysis functions indicates that the optimum XY-value is between 15 and 20 m,
thus an optimum XY-value of 20 m was chosen. For the XY-value of 20, the
velocity analysis data fall very close to two straight lines. Whilst for other XY-
values, the points are scattered about theoretically correct straight lines. For
the zero XY value, the real lateral change in the refractor velocity is obscured by
numerous fictitious velocity changes, which are related to the irregular refracting
interface.
The seismic velocities within the refractor were determined using a standard
“linear method” (see Appendix A). Figure 6-5 shows the average velocity
analysis function, which was computed by averaging the functions for a range of
XY values which were symmetrical about the optimum XY value of 20 m, i.e.
10m<XY<30m. A number of straight lines were fitted onto this function, with the
seismic velocities obtained by taking the reciprocal of the gradients of the lines.
In Figure 6-5, two straight lines were fitted onto the average velocity analysis
function indicating one lateral change of the seismic velocity within the refractor
at station 180. The gradients of these lines give seismic velocities of
approximately 4000 and 5000 m/s and they are within 1 percent of the correct
values.
154
Stacked Velocity Analysis Graph
0
20
40
60
80
100
120
0 100 200 300
Distance (m)
Stac
ked
Velo
city
Ana
lysi
s (m
s)XY=0
XY=5
XY=10
XY=15
XY=20
XY=25
XY=30
Figure 6 - 3: Coincident velocity analysis functions for XY-values from 0 to 30 m (with same reciprocal time), derived from the traveltime data in Figure 6 - 2.
Stacked Velocity Analysis Graph
0
20
40
60
80
100
120
0 100 200 300Distance (m)
Stac
ked
Velo
city
Ana
lysi
s (m
s) XY=0
XY=5
XY=10
XY=15
XY=20
XY=25
XY=30
Figure 6 - 4: Stacked velocity analysis functions for XY-values from 0 to 30 m (with different reciprocal time), derived from the traveltime data in Figure 6-2.
155
Average Velocity Analysis Graph
20
40
60
80
100
120
0 100 200 300
Distance (m)
Ave
rage
d V
eloc
ity A
naly
sis
(ms)
4000 m/s
5000 m/s
Figure 6 - 5: Average velocity analysis graph, derived from the traveltime curves in Figure 6 - 2.
In order to verify the seismic velocities obtained by the “linear method”, a “non-
linear method” was also used (see Appendix A). This method uses a series of
numerical differentiation formulas derived by differentiating Lagrangian
interpolation formulas (Beyer, 1975). The refractor velocities were calculated by
substituting the values of the average velocity analysis function into the
numerical differentiation formulas, to give the derivative at each point along the
function and then taking the reciprocal to give the refractor velocity at each of
these points. In Figure 6-6, the refractor seismic velocities obtained using both
the linear and non-linear methods are shown.
156
Refractor Seismic Velocities Graph
2000
3000
4000
5000
6000
7000
0 50 100 150 200 250 300 350
Distance (m)
Ref
ract
or S
eism
ic V
eloc
ity (m
/s)
3 point5 point7 pointlinear
Figure 6 - 6: Linear versus non-linear refractor velocities from Figure 6 - 2.
The time-depth functions for XY-values ranging from 0 to 30 are shown in
Figure 6-7. As with the velocity analysis function, they have been stacked in
order to more clearly show the minor differences between each XY value. All
the XY-values give essentially the same time-depths except where there are
changes in dip. At these points, use of the non-optimum XY values can
produce significant smoothing.
The resolution of the irregularities in the refractor is maximised with the use of
the optimum XY-value which in this example is 20 m. Depth section is shown
in Figure 6-8 which is particularly evident for the depression at 225 m and for
the fault at 50 m.
157
Stacked Time-depth Graph
0
5
10
15
20
25
30
35
0 50 100 150 200 250 300 350
Distance (m)St
acke
d Ti
me-
dept
h (m
s)
XY=0XY=5XY=10XY=15XY=20XY=25XY=30
Figure 6 - 7: Time-depth for XY-values from 0 to 30 m.
Figure 6 - 8: Depth section which shows particularly for the depression at 225 m and for the fault at 50 m.
158
6.5 Model 2: Irregular Refractor with Irregular Topography
The second model, shown in Figure 6-9, has the same very irregular refractor
interface with the lateral change in the seismic velocity as the first model, but
introduces an irregular topography. The first arrival traveltimes are shown in
Figure 6-10. A comparison with the traveltime graphs for the first model in
Figure 6-2 shows that the effects of near-surface variations are evident. Rapid
changes in topography, such as at 120 and 250 m, are shown as time-shifts in
both the forward and reverse traveltime graphs.
20 4060 80
100 120140
160180
200220 240
260 280300 320 340
20
0
-40
-20
20
0
-20
-40
1500 m/s
4000 m/s 5000 m/s
Rapid changes in topography
Figure 6 - 9: Model with irregular ground and refractor surfaces. The vertical and horizontal scales are equal (Palmer, 1980).
159
First Arrival Times Graph
0
20
40
60
80
100
120
0 50 100 150 200 250 300 350
Distance (m)
Tim
e (m
s)
ForwardReverse
Near surface effects
Figure 6 - 10: Traveltime curves derived from the model in Figure 6 - 9.
The stacked velocity analysis functions for XY-values ranging from zero to 30 m
are shown in Figure 6-11. Whilst none of the velocity analysis functions show
two distinct straight-line segments, an XY-value of 15 m was determined to be
the optimum. As with the first model, the refractor seismic velocities for the
second model were calculated using both the linear and non-linear methods. In
Figure 6-12, five straight lines were fitted on to the average velocity analysis
function, giving refractor seismic velocities of 4000, 3800, 8000, 3600, and 5000
m/s. It is evident that the velocity analysis has been adversely affected by the
variations in the near surface, particularly between 190 and 215 where the
velocity reaches an implausible 8000 m/s. Furthermore, the velocities are not
resolved with the non-linear method, shown in Figure 6-13, which returns an
erratic set of velocities.
160
Figure 6-14 shows the stacked time-depth functions for XY-values ranging from
0 to 30. The time-depths for all XY-values are highly irregular with the effects of
near-surface variations evident. This is particularly the case at 120 and 250 m
where the time-depths shift by 4 ms for the XY-value of zero and propagate
laterally across the time-depths for successive XY-values. Depth sections
calculated with XY-value of 15 m is shown in Figure 6-15. The computed
refractor surface, although following the general shape of the refractor in the
original model, is severely distorted. Without applying corrections for near-
surface effects, the inversion of the seismic refraction data has not generated
an accurate image of the subsurface.
Stacked Velocity Analysis Graph
0
20
40
60
80
100
120
0 50 100 150 200 250 300 350
Distance (m)
Stac
ked
Velo
city
Ana
lysi
s (m
s)
XY=0XY=5XY=10XY=15XY=20XY=25XY=30
Figure 6 - 11: Velocity analysis functions for XY-values from 0 to 30 m derived from the traveltime data in Figure 7-10. The data for a 15 m XY-value are judged the best.
161
Average Velocity Analysis Graph
V1=4000 m/s
V2=3800 m/s
V3=8000 m/s
V4=3600 m/s
V5=5000 m/s
20
40
60
80
100
120
0 50 100 150 200 250 300 350
Distance (m)
Ave
rage
Vel
ocity
Ana
lysi
s (m
s)
V1V2V3V4V5
Figure 6 - 12: Average velocity analysis graph, derived from the traveltime curves in Figure 6 - 10.
Figure 8 - 29: Forward and reverse distance P-wave amplitudes for shots 5 & 14 with
calculating 1/X, 1/X^2 and 1/X^3. The large geometric spreading component dominates.
The effects of geometric spreading have been reduced.
Amplitude Products
0
20
40
60
80
100
0 20 40 60 80 100 120 140distance
scale
d amp
litude
prod
uce (
%)
Shots 5 &14
Shots4&15
Shots1&18
Shots2&17
Shots 3&16
Figure 8 - 30: The product of forward and reverse amplitudes for various offset shot pair.
There are gross similarities in the shape of the amplitude products for all the shot pairs,
suggesting that the variations are related to variations in the seismic velocities in the
refractor.
247
8.6.10 Velocity Ratios
Shear wave (S-wave) data can be readily combined with the interpretation of
conventional seismic compressional-wave (P-wave) data. After determining the
seismic velocities for both P- and S-waves, it was possible to calculate the
velocity ratio Vp/Vs or Poisson’s ratio along each profile. Although Poisson’s
ratio (σ) is not a constant of proportionality between stress and strain, it can be
expressed in terms of elastic constants, and is often used in the same context
as the elastic constants (Tatham, 1985).
The ratio Vp/Vs is defined in terms of Poisson’s ratio (σ) and is given by
(Reynolds, 1997):
2/1
2/11
⎟⎠⎞
⎜⎝⎛
−−
=σ
σVsVp
where:
2/13/4
⎟⎟⎠
⎞⎜⎜⎝
⎛ +=
ρµkVp
( ) 2/1/ ρµ=Vs
k=Bulk modulus
µ = Shear (rigidity) modulus
ρ = density
Note that µ =0 for a fluid, as fluids cannot support shear-waves, and the
maximum value of Poisson’s ratio is 0.5; σ % 0.05 for very hard rocks, %0.45 for
loose, unconsolidated sediments, average % 0.25.
248
8.6.11 Three-Dimensional Images
Three-dimensional images were generated by using the results from the
processing of the first-arrival traveltimes and amplitudes across the four lines at
the first and second site. Image maps in Surfur (8) were used to generate
three-dimensional images.
8.6.12 Traveltime Tomography
Wavepath eikonal traveltime (WET) tomography was used to generate a
velocity-depth section of every line for each wave type. WET tomography
accommodates multiple signal paths contributing to each first break and is a
computationally efficient and geophysically robust method. Wavepath eikonal
traveltime inversion is a high frequency traveltime tomographic method. It is a
computationally efficient method, being an order of magnitude faster than wave-
equation traveltime inversion (yet comparable in effectiveness) since only
solutions to the eikonal equation are involved. It models multiple signal
propagation paths contributing to one first-arrival and as such is superior to
conventional ray-tracing tomography which is limited to the modelling of just one
ray per first-arrival (for a brief discussion on WET theory see Appendix M).
The first step was to generate starting model by using both one and two
dimensional inversion methods. The delta-t-V method was used to generate 1D
starting models (Gebrande and Miller, 1985). This method is a pseudo-
249
tomographic method that yields one-dimensional velocity profiles for each
common midpoint (CMP). This turning ray inversion method delivers
continuous depth vs. velocity profiles for all profile stations. The first step of the
WET inversion algorithm is to propose an initial slowness model, which may be
generated assuming constant vertical velocity gradients such as with 1D tau-p
inversion (see Appendix M), or assuming discrete velocity changes such as with
the generalised reciprocal method (Palmer, 2003a). The depth sections and
seismic velocities obtained with the GRM (as discussed previously) were used
to generate 2-D starting models.
The starting models were then refined with WET tomography until ideally the
modelled traveltimes matched the field traveltimes. A perfect match between
the traveltimes however was never realised, with the mean unsigned error
around 2 ms (Appendix N). This was because after a certain number of
iterations in the WET tomography algorithm, the amount of change between the
updated velocity fields of consecutive iterations decreased, that is the rate of
convergence between the modelled and field traveltimes decreased. At this
point the WET algorithm was terminated.
The differences between the final models generated using the 1-D and 2-D
starting models were significant. Furthermore the errors between them and the
field data were both comparable and within the bounds of acceptable accuracy.
This non-uniqueness is illustrated in Figure 6-31, the top and bottom sections
being derived from 1-D and 2-D starting models respectively. Both fit the field
data to acceptable accuracy. However the differences between the two are
250
large. A refracting interface is easily identifiable in the 2-D derived section at a
depth of 15 m. This is not the case with the 1-D derived section, with a
refracting interface tentatively being placed at a depth of 20 m along the 1000
m/s contour.
The second major difference between the sections is the depths to which they
are imaged. The 1-D derived section has been imaged to a depth of 80 m. The
distances between the offset source points and the receivers were of the same
order, and as such this depth seems unrealistically large. The 2-D derived
section has been imaged to a depth of 20 m which is more plausible for
refracted waves where the minimum source to receiver distances are greater
than four times the depth of the target.
251
0 20 40 60 80 100 120 140-80
-70
-60
-50
-40
-30
-20
-10
0
0 20 40 60 80 100 120 140-30
-20
-10
0
Figure 8 - 31: Final sections from WET tomography. Section generated using 1D starting
model with multiple layers of assumed constant vertical velocity gradients (top). Section
generated using 2D starting model derived from the GRM (bottom). Despite the large
differences in the sections both fit the field data to acceptable accuracy.
252
Chapter 9
Interpretation
9.1 Summary
Three-dimensional three-component shallow seismic refraction surveys at two
sites of dryland salinity in the Spicers Creek Catchment New South Wales have
verified that the shear zone exists as a relatively narrow region of decreased
seismic velocity and increased depths of weathering. Refraction tomography
indicated that the shear zone’s seismic expression continued below the
refracting interface and into the bedrock.
Three-dimensional surveying provided an improvement in the structural imaging
capabilities of the refraction method. Despite showing gross trends along the
orthogonal two-dimensional profiles, variations in the in-line direction suggest
that shear zone has considerable lateral variations parallel to their strike. The
validity of these in-line variations are emphasised by their consistency in the
results of all the wave types. Furthermore, these variations are not predictable
on the basis of orthogonal two-dimensional profiles.
A detailed analysis of the refractor seismic velocities and amplitudes at two
sites, show a number of linear features parallel to and cross-cutting the shear
253
zone (Nikrouz and Palmer, 2004). These lineaments, on basis of the Dubbo
geological map, have been interpreted as a series of recent faults which cut the
older shear zone. The discharge of saline groundwater appears to occur along
the intersection of these faults with the shear zone.
At the first site the results show the seismic geophysical anomaly of half a shear
zone. The velocity-depth sections show that the shear zone’s seismic
expression extends below the refracting interface into the bedrock. All the
wave-types have similar spatial distribution in both time-depth and depth
images and have been interpreted as the boundary between bedrock and
overlying sediment. Changes in refractor seismic velocities occur in both the
cross-line and in-line directions. In the in-line direction, the velocities for all
wave-types show a general decreases with the refractor being separated into
three main regions of different velocities. Furthermore, the seismic velocities in
the refractor in the cross-line direction suggest the presence of a number of
cross cutting and parallel lineaments which have been interpreted to be a series
of faults.
At the second site, the results show the seismic response of a shear zone as a
narrow region of low seismic velocity and increased depth of weathering. The
P-wave/S-wave velocity ratios, which is quite sensitive to pore fluid saturant and
lithology, and Poisson’s ratios (σ ) show distinctive high values over the region.
This has been attributed to an increase in secondary porosity and permeability
from fracturing, leading to an increase in groundwater saturation. At site2, as at
site1, an interpretation of the refractor images of the seismic velocities and
254
comparison with amplitudes show a number of lineaments parallel to and cross-
cutting the shear zone. These lineaments have been interpreted as a series of
recent faults that cut the older shear zone. Furthermore, these cross-cutting
features are consistent with the tectonic interpretation in the south part of the
study area. The discharge of saline groundwater occurs at the intersection of
these faults with the shear zone.
The shear zone’s existence at both sites confirms that there is a relationship
between the geological structure and the high levels of salinisation at the
discharge zones. Groundwater continuously discharges through shallow
alluvial and collovial aquifers along the faults. This forms a “saline plume”
which contaminates the soil down gradient from the seepage zone.
Furthermore, hydrogeology studies by Morgan and Jankowski (2004) have
supported this model of the salination process.
9.2 Introduction
Seismic exploration consists of three main stages: data acquisition, processing,
and interpretation (Yilmaz, 1988). The data acquisition and processing have
been discussed in the previous chapters. In this chapter the interpretation of
data is discussed. The interpretation is based on geophysical methods, which
include magnetic, resistivity and 3D-3C seismic refraction method, and
geological method which includes the interpretation of 1:100 000 Cobbora
geological map.
255
In the Spicers Creek Catchment magnetic, electrical resistivity and seismic
refraction methods were used to image the inferred shear zones. It was
anticipated that in this area, seismic refraction would provide the best results
since the propagation of seismic waves would not be affected by the high level
salinisation whilst also providing possibly greater resolution.
In general, the geophysical and geological interpretation show a consistent
result in terms of the being of shear zone and cross-cutting features at two sites
in the Spicers Creek Catchment.
9.3 Magnetic Data Interpretation
Figures 9-1 and 9-2 show the magnetic field in nanoteslas as a 2D contour plot
for the two sites. The magnetic results show reductions of approximately 2100
nanoteslas at site1 and 600 nanoteslas at site2 over the inferred zones of the
shear zones. The low magnetic values are inferred to be caused by the
enhanced oxidation of magnetite, the magnetite mineral, to hematite, goethite or
maghemite, resulting from fracturing of the rock associated with the shear zone.
Magnetite is unstable in the typically low temperature and highly oxidizing
environment of chemical weathering and sedimentation. Groundwater
containing CO causes the ferrous iron minerals to eventually break down with
Fe being slightly soluble in water containing CO . Meanwhile, the oxygen
fugacity sets the upper limit at which magnetite can exist as a stable compound.
2
+22
256
It is given by the fO2 vs T relationship that governs the oxidation of magnetite to
hematite (Grant, 1985):
)(32/1)(2 32243 hematiteOFeOmagnetiteOFe ⇔+
At the Spicers Creek Catchment, the magnetics show lows over the inferred
location of the shear zones. These magnetic lows are consistent with the
weathering of magnetite by the saline groundwater. Furthermore, the local
magnetic data is consistent with the airborne magnetic data obtained by AGSO
over the Dubbo 1:250,000 map (as discussed in chapter 5).
9.4 Resistivity Data Interpretation
Electrical resistivity survey methods have been employed for many decades to
image the near surface where they can provide useful information on geological
structure, lithologies and subsurface water resources. Resistivity methods are
widely employed in hydrogeological investigations (Kearey and Brooks, 2001),
because earth resistivity can be readily related to porosity and water quality.
2D resistivity imaging is a geophysical method which aims to generate a cross
section of the electrical properties of the subsurface by passing an electric
current along many different paths and measuring the associated potential
differences. This can be accomplished by connecting a linear array of
257
Figure 9 - 1: 2D contour map of magnetic data at site 1.
Inferred shear zone
Inferred shear zone
Figure 9 - 2: 2D contour map of magnetic data at site 2.
258
grounded electrodes through a multicore cable to a computer controlled
switching module and a resistivity meter. However, with the electrical resistivity
surveys at the Spicers Creek, the data were acquired with a traditional field
system using two current and two potential electrodes. A large number of
apparent resistivity measurements were made using differing combinations of
four electrodes along the array in a Wenner configuration. The apparent
resistivity measurements were used to plot pseudosections, which are
qualitatively related to the distribution of true resistivity. The pseudosections
are inverted using either a manual interactive approach or an automatic iterative
method to improve the detail and correctly scale the apparent resistivity image
giving a true resistivity-depth section (Griffiths and Barker, 1993).
The conductivity* of geological materials exhibits one of the largest ranges of all
physical properties (Reynolds, 1997). In a normal crust, from the surface to up
to 15 km depth, bulk conductivities are largely controlled by aqueous electrolytic
conduction (Ward, 1990). The major factors affecting the bulk electric
conductivity of soil or rock are (McNeill, 1990):
1. Porosity and permeability
2. Degree of saturation (fraction of pore space/fractures filled with moisture)
3. Presence of clays with moderate to high cation exchange capacity (CEC)
4. Conductivity of interstitial fluid
* Resistivity is the inverse of conductivity and has units of ohm-metres (Ωm). The following will refer to conductivity, which has units of siemens/metre (S/m). 259
Geological processes such as dissolution, faulting, shearing, and weathering
usually increase porosity and fluid permeability. Thus the result of such process
on a rock formation will be an increase in conductivity.
The fluid saturation of a rock formation is a direct influence on its conductivity.
A large proportion of pore spaces, fractures and other gaps filled with fluid in a
rock mass will enhance aqueous electrolytic conduction. In general, the
conductivity of a rock increases as the percent fluid saturation increases.
The presence of clay adds an additional component to conductivity, with clay
particles acting as separate conducting paths to the electrolyte path. Due to the
large cation-exchange capacity of clay, the cations of the diffuse layer are free
to move under the influence of an applied electric field. This leads to an
increase in the density of charge carriers and hence an increase in conductivity.
The salinity of interstitial fluid is perhaps the greatest factor determining the
conductivity of a rock. The ions that enable the flow of electric current in a fluid
result from the dissociation of salts, which occurs when salts are dissolved in
water. Since the quantity of charge that can be carried by an ion is finite, the
more ions available in a solution, the greater the total charge, and hence the
higher the conductivity. In general, a rock which contains water within its pores
will have a greater conductivity when the salinity of the water is high than when
it is low.
260
The resistivity of geological materials exhibits one of the largest ranges of all
physical properties, from m for native silver to m for pure
sulfur (Reynolds, 1997). In sedimentary rocks, the resistivity of the interstitial
fluid is probably more important than that the host rock. Archie (1942)
developed an empirical formula for the effective resistivity of a rock formation in
terms of the resistivity and volume of the pore water present as below:
Ω× −8106.1 Ω1610
wcb fa ρφρ −−=
Where φ is the porosity, the fraction of pores containing water of resistivity f
wρ and , and are empirical constants (a b c 5.0 5.2≤≤ a , , and
). Korvin (1982) has proposed a theoretical basis to account for Archie’s
Law. Salin groundwater may have a resistivity as low as 0.05 Ωm and some
groundwater can have resistivities in excess of 1000
5.23.1 ≤≤ m
2≈n
Ωm.
The 2D resistivity imaging results are shown in Figures 9-3 to 9-6. For each
Figure the top pseudosection is the measured apparent resistivity, whilst the
middle is the apparent resistivity starting model calculated by the inversion
software. The bottom is the resistivity-depth section resulting after three
iterations of the inversion program.
Evidence of the shear zone is not apparent in the resistivity-depth sections.
The sections do not show narrow vertical zones of low resistivity as would be
expected from the presence of the shear zone. Rather, they all show an
approximately layered subsurface with increasing resistivity with depth. The
very near surface does however show a degree of lateral variation in resistivity
261
such as at the second site (Figure 9-6) where there are very erratic resistivity
variations in the top 8 metres of the section. These near surface variations
represent geological noise which is prevalent where the overburden has
irregular resistivity or thickness.
A large resistivity contrast exists at line 2 of the first site (Figure 9-4) where
there is a change from 7 to 17 Ωm at an increasing depth along the line. This
contrast corresponds to the boundary between the bedrock and colluvium. The
depth to bedrock is not evident in the other sections despite all being imaged to
a depth of 20 m. This is possibly due to the inversion algorithm which, although
correcting the depth variations in resistivity, images sharp boundaries as
gradational. As a result resistivity contrasts appear less than true (Barker,
1997).
At the Spicers Creek Catchment, it was anticipated that resistivity methods
would not be the most suitable for imaging the shear zone. The widespread
occurrence and high level of salinisation in the catchment has resulted in a
highly conductive overburden. The geological noise produced by this
conductive overburden during resistivity surveys would obscure the conductivity
anomaly due to shear zone in the bedrock. In addition the resistivity image
sections show low resistivity zones in the area of the inferred shear zones.
However, the regions of low conductivity were quite extensive laterally and this
together with the limited penetration did not facilitate the precise delineation of
the shear zones.
262
Figure 9 - 3: The resistivity result at Site1-Line1.
Figure 9 - 4: The resistivity result at Site1-line2.
263
Figure 9 - 5: The resistivity result at Site1-Line3.
Figure 9 - 6: The resistivity result at Site2.
264
9.5 3D-3C Seismic Data Interpretation
Seismic data interpretation for the first and second sites is based on the plans
of the various refractor properties (shown in Appendix O) and velocity-depth
sections obtained from wavepath eikonal traveltime (WET) inversion of a two
dimensional starting model derived from the generalized reciprocal method
(GRM) (shown in Appendix P). Appendix O shows the different properties of
the refractor including time-depth, depth, seismic velocity, seismic velocity ratio,
Poisson’s ratio, amplitude product, and amplitude product ratio. Appendix P
includes GRM refractor seismic velocity images for both the P-wave and S-
waves for each line.
9.5.1 Refractor Seismic Images at the First Site
On the basis of the velocities, the refracting interface at the first site has been
interpreted as the boundary between the bedrock and the overlying sediment.
Refractor time-depth maps and depth maps for each wave-type are shown in
Appendix O. The time-depth and depth maps show a similar spatial distribution
and the depth to the refractor for all wave-types (P-, SV- and SH) and steadily
decrease in the increasing cross-line direction. The spatial distribution is most
prominent for SV- and SH- waves which have similar depth maps where the
refractor increases from a depth of 13 to 17 m. Table 9-1 shows the minimum,
maximum depth and time-depth for each wave-type.
265
In general, the refractor depths as imaged by the SH- and SV- waves across
the whole image are comparable. From cross-line 0 to 60, the refractor depths
as imaged by the P-waves and S-waves are similar. However, after cross-line
60, the rate of increase in depth to the S-waves refractor is greater than the P-
wave. By the end of profile, the S-wave images the refractor 3 m deeper than
the P-wave. This is especially evident on the time-depth and depth profiles for
the individual lines as show in Appendix L chapter 8.
Table 9 - 1: Minimum, maximum and average time-depth and depth for each wave-type.
Time-depth (ms) Depth (m)
Wave-
type
Minimum Maximum Minimum Maximum
P-wave -17 -22 12.5 15.5
SH- wave -37 -52 13.5 18.5
SV- wave -33 -55 11 19
The difference in depths that the P-waves and S-waves image the refractor can
be explained by the likely occurrence of an undetected layer consisting of
saturated sediments. The sediment overlying the bedrock changes from being
damp or partially saturated to completely saturated at a point along the bedrock
interface (shown in Figure 9-7). This increase in saturation represents a
genuine increase in P-wave velocity, from 1150 m/s to a likely 1750 m/s.
266
Bedrock
Partially saturated sediment
Saturated sediment
Vp1Vs1
Vp2 > Vp1Vs2 = Vs1
Vp3 > Vp2Vs3 > Vs2
Figure 9 - 7: The possible P-wave undetected layer. A change from partially saturated to
completely saturated sediment above the refractor will result in an increase in P-wave
velocity. If this layer of completely saturated sediment is relatively thin, then it will be
undetectable in the P-wave traveltime graphs, resulting in the refractor being imaged
shallower than the true depth. Since S-waves do not propagate through liquids, the
degree of saturation in the sediment will have no effect on its propagation, thus the S-
waves will image the refractor at its true depth.
However, the saturated layer is not thick enough to generate refracted P-wave
first arrivals, and therefore is not recognisable in the traveltime graphs. The
result of this undetected layer is to image the P-wave refractor shallower than
its true depth. Furthermore, in this case the problem of the undetected layer is
not overcome by the use of an average seismic velocity above the refractor
derived with the GRM. This is because the saturated sediment is not present
along the whole profile and also progressively changes thickness.
267
Since S-waves do not propagate through liquids, the degree of saturation in the
sediments has no influence on the velocity of these waves. Thus the presence
of the completely saturated layer along part of the profile will have no effect on
the S-wave image of the refractor. The layer of saturated sediment is present
where the refractor is imaged at differing depths by the P-waves and S-waves,
and is absent where the refractor is imaged at the same depth by both waves.
The refractor images for all wave-types are shown in Appendix O. Changes in
refractor seismic velocities (Figure O3) occur in both the cross-line and the in-
line directions. In the in-line direction, the velocities for all wave-types show a
general decrease with the refractor being separated into three main regions of
different velocities (Figure 9-8). The P-wave velocities are between 4700-5300
m/s between 0 and 30 m, around 4000 m/s between 30m and 120m and 3000
m/s between 120m and 140m. The constant in seismic velocities between the
regions is not as large for the S-waves, with both the SV- and SH-waves having
velocities greater than 2500 m/s in the first region, 1700 for the second, and
around 1400 m/s for the third region. The consistent zoning of refractor
velocities for all wave-types is compelling evidence for genuine lateral changes
within the bedrock. At the first site, the existence of half of a shear zone is
evident from the first 40 m in the in-line direction (line AB on Figure 9-8). Line
AB is coincident with a major fault in the Spicers Creek Catchment as is
discussed in the geology section.
268
B
AB
A B
A
Figure 9 - 8: Summary of the refractor velocities at site 1. The refractor can be separated
into three distinct regions of different seismic velocities.
An analysis of the seismic velocities in the refractor in the cross-line direction
suggests the presence of a number of cross cutting and parallel lineaments
which have been interpreted to be a series of faults (Figure 9-9). There is a
consistency in the patterns of the seismic velocities for all wave-types, which
suggests the existence of faults. These features can be seen such as the
velocity high at cross-line 35m and in-line 10m, which is also exist in all wave-
types; the velocity low at cross-line 40m and in-line 20m-30m. The similarities
269
in the spatial features and lineaments for all wave-types suggest that the
seismic velocities are true representations of the refractor, with minimal
artefacts from the inversion algorithm, or errors in the first arrival traveltimes.
V~4.9
V~5.3
V~3
V~3.5
V~4
V~2.2
V~2.5
V~2.7
Cross-line
In-li
ne
Cross-line
In-li
ne
Cross-line
In-li
ne
Figure 9 - 9: Cross-cutting features at site1.
270
Variations in P-wave/S-wave velocity (Vp/Vs) ratio and Poisson’s ratio are
associated with the distribution of crack or pore geometry and the subsequent
degree of fluid saturation, rather than with a change in the elastic constants of
the minerals comprising the matrix (Dohr, 1985). Increases in these ratios are
generally an indication of increasing porosity and permeability, and fluid
saturation. At the first site, the Vp/Vs ratio and Poisson’s ratio show a gross
increase in the cross-line direction. Vp/Vs ratio increases from around 1.4 to
over 2.2 (Appendix O, Figure O4). Poisson’s ratio ranges from -0.5 to 0.4
(Appendix O, Figure O5), which is within the theoretical limits of -1 for a solid of
very high rigidity and 0.5 for a fluid. The increase in Vp/Vs and Poisson’s ratio
can be attributed to an increase in fracturing and fluid saturation in the bedrock.
Overlaid on the gross variations in Vp/Vs ratio and Poisson’s ratio are features
related to changes in the SV- and SH-wave velocity. The S-wave velocity highs
at cross-line 5m, in-line 10m is represented by lows in Vp/Vs and Poisson’s
ratios, whilst the S-wave velocity low at cross-line 120 m to 140 m from in-line 0
to 20 m is represented by highs. The differences between the SV- and SH-
wave velocities are small as indicated by the SV-wave/SH-wave velocity ratio
which has an average of 1 and ranges from 0.75 to 1.25.
The amplitudes products are shown in Appendix O Figure O6. The amplitude
products are variable for P-waves and S-waves. However, there are similarities
with the lineaments on the amplitude products and the lineaments interpreted
from the seismic velocities, both showing a number of parallel and cross-cutting
features. Perhaps one features of interest is the wide zone of low amplitude
271
between cross-line 60m and 120 m for P-waves and SV and meanwhile cross-
line 60 m to 100 m for SH.
Figure 9-10 is the wireframe map of amplitude product with the geometric effect
correction being applied. The map can be divided into three regions. A region
of high amplitude product exists between cross-line between 0 to 60 m, a region
of low amplitude product exists between cross-line between 60m to 100 m, and
a region of high amplitude product exists between cross-line between 100 m to
140 m. The presence of these regions suggests that variations in amplitudes
are associated with changes in the refractor. A considerable degree of noise is
also evident on the amplitude products and is related to lateral variations in the
near surface layers. The low amplitude between cross-line 0 to 20 m and high
amplitude between cross-line 100 m to 120 m are possibly due to near surface
effects in the study area.
0.004
0.006
0.008
0.01
0.012
0.014
0.016
0.018
0.02
0.022
0.024
Amplitude Product (S1*S18) , P-wave, Site 1
Amplitude Product
Cross-line
In-line
Figure 9 - 10: Amplitude product for P-wave at site1.
272
As with the amplitude products, the amplitude product ratios show little
consistency between wave-types. However, a zone between cross-line 60 m to
100 m is similar to the amplitude low mentioned above. The cross-cutting
features are not obvious in amplitude product maps possibly due to near
surface effects.
Figures P1 to P3 (Appendix P) are velocity-depth sections obtained from
wavepath eikonal traveltime (WET) inversion of a two dimensional starting
model derived from the generalized reciprocal method (GRM) for all wave-
types.
The refractor depths for all wave-types are between 11 m to 13 m and are
similar to the depth sections as shown in Appendix L chapter 8. The velocities
for all wave-types show a general decrease in the increasing cross-line. The
refractor can be separated into three main regions of different seismic
velocities. The first region has P-wave velocities of between 5200 – 3600 m/s,
and around 2800 m/s and 2000 m/s for the second and third regions
respectively. For SH waves the first region has velocities of between 3600-
2200 m/s, and around 1600 and 800 m/s for the second and third regions
respectively. The distinction between regions is less evident for the SV waves
but lateral changes are obvious.
In general, the images show the lateral changes in refractor in term of velocity
for all type-waves. The presence of half a shear zone can be deduced from
these GRM tomography images. The existence of lateral changes is due to the
273
different geological strata and the presence of the fault at study area. Cross
cutting features are not obvious on these sections.
9.5.2 Refractor Seismic Images at the Second Site
Figures O8 and O9 (Appendix O) are the refractor time-depth and depth images
at the second site. All the wave-types have similar spatial distribution in both
images and have been interpreted as the boundary between bedrock and
overlying sediment. There is an increase in refractor depth in the middle of the
profile between cross-line 40 m to 100 m. The P-waves and SH-waves image
the refractor at a maximum depth of around 10.5 m between cross-line 40 m
and 90 m. The shallowest depths are at the ends of profile around cross-line 20
m and cross-line 120 m. The SV-waves image the refractor at a maximum
depth around 14 m between 30 m and 60 m, and are the shallowest at the ends
of profile around cross-line 15 m and 130 m. Table 9-2 shows the minimum,
maximum depth and time-depth for each wave-type.
Table 9 - 2: Minimum, maximum and average time-depth and depth for each wave-type.
Time-depth (ms) Depth (m)
Wave-
type
Minimum Maximum Minimum Maximum
P-wave -5 -12 4.5 10.5
SH- wave -10 -22 5.5 10.5
SV- wave -12 -26 6 14
274
The SH-wave depth image shows that the refractor has been shifted laterally by
20 m in the decreasing cross-line direction and is 3 m deeper than the other
wave-types. The lateral shift in the refractor depth is probably due to SV-P-SV
wave-type conversion (Figure 9-11) and the inconsistent use of the first-arrival
traveltimes in the GRM algorithms. SV-P-SV wave-type conversion occurs
when an SV-wave generated at the surface, propagates into the subsurface at a
velocity Vs. At the refracting interface, SV-wave is critically refracted and is
converted to a P-wave and travels along the refractor at a greater velocity Vp.
The head waves generated by this wave are then converted back to SV-waves
and travel at a velocity Vs to the surface. Thus waves which are thought to be
pure SV-waves are in fact SV-P-SV-waves. A lateral shift of the refractor may
result if the forward and reverse traveltimes used by the GRM algorithms has
one set derived from SV-P-SV-waves and the other set from pure SV-waves
(shown in Figure 9-11). Furthermore, the use of the faster SV-P-SV-wave
velocity in the depth conversion factor will result in the refractor being imaged at
a greater depth.
Unlike the first site, there is no considerable difference in the refractor’s depth
as imaged by the P-waves and SH-waves. This suggests that the layer of
partially saturated or completely saturated sediment is absent. The surface
velocities at the second site are greater than at the first, suggesting the absence
of the partially saturated layer.
275
A C D Bsurface
SVSVSV
SV refracting
PP
Figure 9 - 11: Wave-type conversion and the GRM. The forward and reverse traveltimes
employed by the GRM algorithms include one set derived from SV-P-SV-waves and the
other set from pure SV-waves. When this occurs, the refractor is shifted laterally and is
imaged at a greater depth.
The velocities for all wave-types are shown in Figure O10 Appendix O. As with
the first site, the change of velocities in the refractor occurs in both directions.
In the cross-line direction the refractor can be separated into four regions of
different seismic velocities, with a region of low seismic velocity occurring
between regions of high seismic velocities (Figure 9-12).
276
5
2.8
Figure 9 - 12: Summary of the refractor velocities at site 2. Refractor can be separated
into four distinct regions of different seismic velocities.
For all wave-types the boundary between the region of low and high seismic
velocity is quite distinct. The first region is between cross-line 0 to 20 m, the
second between 20 m to 50 m, the third between 50 m to 110 m, and the fourth
277
between 110 m to 140 m. The first region has P-wave velocities around 3000
m/s, and around 5000 m/s and 3800 m/s and 5400 m/s for the second and third
and fourth regions respectively. The seismic velocities for the SH-wave are
around 2000 m/s, 2300 m/s and 1500 m/s and 2500 m/s for the first, second,
third and fourth regions respectively. The seismic velocities for the SV-wave
are around 3500 m/s, 2800 m/s, 2000 m/s and 3000 m/s for the first, second,
third, and fourth regions respectively.
The velocity contrast at third region and the other three regions are smaller for
the P-waves than for S-waves. This is possibly due to the degree of fluid
saturation in the bedrock. The P-wave velocity contrast between unfractured
rock and fluid-saturated fractured rock is smaller than between unfractured rock
and dry fractured rock. Since they do not propagate through liquids, the velocity
of S-waves in fractured rock is the same no matter what the degree of fluid
saturation. The high degree of fluid saturation in the bedrock has led to a
decrease in the velocity contrast of the P-waves whilst not affecting the velocity
contrast of the S-waves.
The consistent zoning of refractor velocities for all wave-types is compelling
evidence for genuine lateral lithological changes within the bedrock especially
for the third mentioned region. These seismic velocity models together with the
depth images suggest a narrow shear zone at the third region (cross-line 50 m
to 110 m) characterized by low seismic velocities and increased depths of
weathering.
278
As with the first site, in the inline direction a detailed analysis of the seismic
velocities in the refractor suggests the presence of a number of cross cutting
and parallel lineaments which have been interpreted to be a series of faults
(Figure 9-13). There are a number of spatial similarities between all wave-types
velocities between cross-line 0 to 50 m, in-line 0 to 30 m, and between cross-
line 110 m to 140 m, in-line 0 to 30 m which seem to be related to the faults.
While The cross cutting features can be seen on both sides on the images for
all wave types, the SH-wave shows these features the best. This maybe due to
the direction of propagation SH-waves throughout the refractor. Figure 9-14
shows the direction of propagation for all wave-types.
There is a region of increased P-wave/SH-wave velocity ratio and associated
Poisson’s ratio between cross-lines 40 m and 110 m (Figure O11, O12,
Appendix O). The increase is considerable with a Vs/Vp ratio of 2.25 and
Poisson’s ratio of 0.35, contrasted with Vs/Vp ratio of 1.75 and Poisson’s ratio
of 0.3 in the surrounding regions. This increase correlates with the region of
low refractor velocity and is possibly related to an increase in porosity and
permeability, and a subsequent increase in fluid saturation.
The lateral variation in the SV-wave velocity can be resolved by examining the
velocity ratios. A P-wave/SV-wave velocity ratio high of greater than 2.25
occurs between cross-lines 20 m and 40 m, whilst the SV-wave/SH-wave
velocity ratio high of greater than 1.75 occurs between cross-lines 90 m and
110 m.
279
V~3.5
V~2.4
V~4.8
V~5.8
V~4.8
V~5.6
2.8
V~3.8 V~3.3
V~2.6 V~3.2
V~3.3
V~2.6
V~2.5
V~2
V~2.2
V~2.4
V~3
V~2
Cross-line
In-line
Cross-line
In-line
In-line
Cross-line
Figure 9 - 13: Cross-cutting features at site2.
280
Figure 9 - 14: Distinction among the particle-displacement vectors associated with the
three fundamental modes, P, SH, and SV, which compose vector-wavefield seismic data
(Hardage et al., 2003).
Figure O13 (Appendix O) shows the amplitude products for all wave-types.
There is a general increase in the increasing cross-line direction. Like the first
site there is no correlation between the amplitudes and the seismic velocities in
the refractor. However, there are similarities with a number of cross-cutting and
parallel faults in the seismic velocity profiles. The amplitude images show the
lateral changes in the cross-line related to the faults in the bedrock.
Variations in amplitudes related to changes in the bedrock may be resolved
using the P-wave/S-wave ratios of the amplitude product (O14). Features such
281
as the high between cross-line 40 m and 100 m in the P-wave/SV-wave ratio,
and between 50 m and 110 m in the P-wave/SH-wave ratio, are clearly
detectable. This region of high P-wave amplitude and low S-wave amplitude
may be due a high degree of fracturing causing the preferential attenuation of
S-waves with respect to P-waves as S-waves are more directly affected by
changes in rigidity modulus and density.
Figures P4 to P6 (Appendix P) show the GRM velocity-depth sections for all
wave-types. The low velocity zone can be seen for most of images including P-
waves (line1 and line4), SH-waves (line1, line3, line4) and SV-waves (line4).
There is no obvious low velocity zone for some images because of the limited
penetration of some of the wave-types into the refractor. The low velocity zones
support a region of increased porosity and permeability, and the subsequent
increases in the degree of groundwater saturation, that correlates with the
presence of the shear zone at the second site.
9.6 Geology and Tectonic Interpretation
As discussed in the geology chapter, the study area is situated in the mid-
Silurian to Early Devonian Hill End Trough in the eastern Lachlan Fold Belt
(LFB) and overlies marginal areas of the Gunnedah Basin and Surat Basin.
The Hill End Trough is bounded to the south by the Carboniferous Buthurst
Batholith and the Lachlan Transverse Zone (Glen and Walshe, 1999) and to the
east and west by Ordovician volcanics that are overlain by thin Silurian-Early
282
Devonian rocks and intruded in part by Carboniferous granites (David et al.,
2003). Figure 9-15 shows the simplified regional geology map of the Hill End
Trough (David et al., 2003).
The geology and tectonic interpretation shows a detailed model of the cross-
cutting features at the first and second site. They confirm that there is a good
relationship between the interpretation of seismic data and that of the geology
at the Spicers Creek Catchment. The detection of lineaments and cross-cutting
features through the interpretation of the seismic velocities and amplitude
products was made possible via the use of three dimensional seismic refraction
methods.
283
Figure 9 - 15: Simplified regional geology map of the Hill End Trough (Modified from
surveys were recorded over a shear zone at two sites associated with dryland
salinity in the Spicers Creek Catchment, near Dubbo in southeastern Australia.
The seismic data were recorded with the Australian National Seismic Imaging
Resources (ANSIR) 360-trace ARAM-24 seismic system and IVI MiniVibrator.
Dryland salinity occurs extensively throughout the Spicers Creek Catchment.
The high concentration of salt in the groundwater has led to a significant decline
in agricultural productivity and a reduction in native vegetation. Furthermore,
the saline groundwater in the surface soil has caused the destruction of the clay
and soil structure and as a result, large areas of the catchment have
experienced soil erosion and extensive alteration of the landscape.
The broad objective of this study was to use seismic refraction methods to map
in detail a shear zone, which is associated with the salination. Although the
shear zone had been recognized in airborne geophysical surveys and
subsequent regional geological mapping as the Nindethana Fault, the exact
geological factors which produced the discontinuous occurrences of saline
groundwater at the surface, were not known. Detailed ground magnetic surveys
carried out as part of this study, confirmed the location of the shear zone, but
285
they did not provide further explanation of the causes of salination. Therefore,
the detailed objectives were to map the extent of the shear zone as well as
other properties, such as the rock fabric.
Seismic refraction methods were selected because of their potential ability to
provide greater lateral resolution of the narrow vertical shear zone, than is
currently the norm with electrical or electromagnetic methods. This situation
was confirmed with a number of resistivity depth images generated as part of
this study. Furthermore, it was anticipated that the high electrical conductivities
of the near surface areas would have no significant impact on the depth of
investigation of the seismic energy. Finally it was also anticipated that the head
wave amplitude might provide a more detailed measure of the seismic velocities
as a function of azimuth and therefore of fracture porosity.
The 3D 3C survey was set out as four parallel lines 10 m apart with twenty-nine
3C geophones, 5 m apart, in each line. The source lines were located parallel
and orthogonal to the lines of geophones and passing through the geometric
centre of the receiver spread. Source points were nominally 10 m apart and
both P- and S-waves were generated.
The results of the seismic refraction surveys show that the shear zone occurs
as a narrow region with low seismic velocities and increased depths of
weathering. However, there were numerous ambiguities in generating
consistent parameters in the refractor between adjacent receiver lines.
Furthermore, these ambiguities were not resolved with the use refraction
286
tomography, because even moderate differences with the various starting
models still generated acceptable agreement with the traveltime data.
Therefore, in order to generate consistency between the results for the various
receiver lines, it was necessary to generate detailed starting models.
It was found that near surface irregularities of a very limited lateral extent were
the major sources of uncertainty in deriving detailed starting models. In the
Spicers Creek Catchment, the variations in the traveltimes caused by the near
surface irregularities, were frequently of a comparable magnitude to those
caused by the changes in depths and seismic velocities in the refractor. As a
result, the derivation of detailed parameters of the refractor using the
generalized reciprocal method (GRM) was subject to significant uncertainty
unless they are accommodated with the GRM SSM.
A major achievement of this study is the development of a smoothing method
which significantly removes the effects of the near surface irregularities.
Termed the GRM SSM (for statics smoothing method), it essentially generates
a time-depth model of the refracting interface for which the effects of the near
surface irregularities have been minimized, by taking an average of the time-
depth profiles for a range of XY values. Although the time-depths to the target
refractor show very little variation with changing XY distances, the traveltime
anomalies for the near surface irregularities migrate laterally through the time-
depths. Therefore, the average time-depth profile is largely free of the near
surface anomalies. When this averaged time-depth profile is subtracted from
that for zero XY, the result is essentially a set of corrections for the near surface
287
irregularities. These corrections are then subtracted from the field traveltimes,
and the GRM computations are then repeated with the corrected traveltimes.
The GRM SSM, which takes advantage of the unique redundancy properties of
the GRM computations, was a major factor in deriving consistent detailed
starting models for refraction inversion. Furthermore, this consistency was
achieved with both S-wave components, as well as the P-wave results.
The major geological achievement, which was made possible with the GRM
SSM, was the demonstration that there were cross cutting features associated
with the major shear zone. Therefore, it appears that saline groundwater can
discharge at the surface where increased volumes of groundwater occur at the
intersection of different sets of shears. This model provides a useful
explanation for the discontinuous occurrence of salination along the major shear
zone.
As is usually the case, increased numbers of recording channels would have
been beneficial. In particular, the four recording lines separated by 10 m did not
provide as high resolution in the cross-line direction as in the in-line direction.
Therefore, where cross cutting features are considered likely to occur and to be
significant, than re-deploying the recording spread in the orthogonal direction
provides an acceptable alternative to employing field systems with significantly
larger recording capacity.
The measurement of S-wave velocities provided a measure of Poisson’s ratio,
which is a useful elastic parameter. Although a detailed analysis of the head
288
wave amplitudes did not generate useful results, nevertheless, the fact that the
head wave amplitude is a function of the densities, as well as the seismic
velocities, suggests that the joint inversion of seismic refraction traveltimes and
head wave amplitudes, with detailed gravity profiles, should facilitate the
determination of both seismic velocity and density models, and in turn, the
derivation of elastic constants. It is likely that the detailed geotechnical
characterization of sites with 3D 3C seismic refraction methods, together with
detailed gravity profiles, could provide useful quantitative models for the
analysis of groundwater flow in fractured rock masses, as well as for the
traditional engineering construction site applications.
289
REFERENCES ACIL, Australia Pty Ltd, 1983, Causes, Extent and Effect of Salinity in Victoria:
A Report Prepared for the Salinity Committee of the Victorian Parliament, Melbourne, pp. 175.
Acworth, R. I. and Jankowski, J., 2001, Salt source for dryland salinity-evidence
from an upland catchment on the Southern Tablelands Of New South Wales: Aust. J.Soil.Res., 39, pp. 39-59.
Allan, U. S., 1989, Model for hydrocarbon migration and entrapment within
faulted structures: AAPG Bulletin, 73, pp. 803-811. Anderson, R. N., Flemings, P., Losh, S., Austin, J. and Woodhams, R., 1994,
Gulf of Mexico growth fault drilled, seen as oil, gas migration pathway: Oil and Gas J., 12, pp. 97-102.
Anderson, R. N., He, W., Hobart, M. A., Wilkinson, C. R. and Nelson, H. R.,
1991, Active fluid flow in the Eugene Island area, offshore Louisiana: The Leading Edge, 10, pp. 12-17.
Archie, G. E., 1942, The electrical resistivity log as an aid to determining some
reservoir characteristics: A.I.M.E., 146, pp. 389-406. Atkinson, T., 1970, Ground preparation by ripping in open pit mining: Mining
Magazine, 122, pp. 458-469. Bamford, D. and Nunn, K. R., 1979, In situ seismic measurements of crack
anisotropy in the carboniferous limestone of northwest England: Geophys. Prosp., 27, pp. 322-338.
Barker, R. D., 1997, Electrical imaging and its applications in engineering
investigations, In: McCann, D. M., Eddleston, M., Fenning, P.J., Reeves, G.M.(Ed): Modern geophysics in engineering geology, The Geological Society, London, pp. 37-43.
Barron, L. M., 1997, Igneous petrology overview of the Dubbo 1:250,000
geological sheet: Geological Survey of New South Wales, Report GS1997/080.
Bartlett, A. H., 1962, Refraction seismic survey in civil engineering: Australian
Road Research Board 1st Biennial Conference. Bates, R. L. and Jackson, J. A., 1980, Glossary of Geology: American
Geological Institute, U.S.A. Beach, A. and Fyfe, W. S., 1972, Fluid transport and shear zones at Scourie,
Sutherland: evidence for overthrusting?: Contrib.Mineral.Petrol., 36, pp. 175-180.
290
Bennett, G., 1999, 3-D seismic refraction for deep exploration targets: The Leading Edge, 18, pp. 186-191.
Berg, R. R. and Avery, A. H., 1995, Sealing properties of tertiary growth faults,
Texas Gulf Coast: AAPG Bulletin, 79, pp. 375-393. Beyer, W. H., 1975, Handbook of mathematical science, 5th edition: CRC
Press, Florida. Biehler, S., Ferguson, J., Baldridge, W. S., Jiracek, G. R., Aldren, J. L.,
Martinez, M., Fernandez, R., Romo, J., Gilpin, B., Braile, L. W., Hersey, D. R., Luyendyk, B. P. and Aiken, C. L., 1991, A geophysical model of the Espanol Basin, Rio Grande rift, New Mexico: Geophysics, 56, pp. 340-353.
Bigelow, N., 1965, Seismograph operations by marine state highway
commission: Highway Research Record No.81, Highway Research Board, Washington, pp. 9-15.
Bonini, E. W. and Hickok, E. A., 1958, Seismic refraction method in ground-
water exploration: AIME Trans., 211, pp. 485-488. Bourke, D. J., Hawke, J. M. and Scheibner, E., 1974, Structural subdivision of
the Great Australian Basin in New South Wales: Geological Survey of New South Wales, Quarterly Notes, 16, pp. 10-16.
Bradd, J. M. and Gates, G., 1996, Dryland salinity in New South Wales. A state
perspective: Department of Land and Water Conservation, Technical Report No. TS95.113.
Bradd, J. M., Milne-Home, W. A. and Gates, G., 1997, Overview of factors
leading to dryland salinity and its potential hazard in New South Wales, Australia: Hydrogeology Journal, 5, pp. 51-67.
Branagan, D. F. and Packham, G. H., 2000, Field geology of New South Wales:
Department of Mineral Resources New South Wales, Sydney, Australia. Braun, J., Munroe, S. and Cox, F., 2003, Transient fluid flow in and around a
fault: Geofluids, 3, pp. 81-87. Brown, A. R., 1996, Interpretation of three-dimensional seismic data: AAPG
Memoir, 42. Brown, P. D. and Robertshaw, J., 1957, Determining the thickness of
unconsolidated deposits overlying shallow mine working by seismic methods in methods and case histories in mining geophysics: Sixth Commonwealth Mining and Metallurgical Congress, Montreal.
291
Campos-Enriquez, J., Delgado-Redriguez, O., Chavez-Segura, R., Gomez-Contreras, P., Flores-Marquez, E. L. and Birch, F. S., 1997, The subsurface structure of the Chalco sub-basin (Mexico City) inferred from geophysical data: Geophysics, 62, pp. 23-35.
Carcione, J. M. and Cavallini, F., 2002, Poisson's ratio at high pore pressure.:
Geophysical Prospecting, 50, pp. 97-106. Cerveny, V. and Ravindra, R., 1971, Theory of seismic head waves: University
of Toronto Press. Cerveny, V., Molotkov, I.A., Psencik, I., 1977, Ray Method in Seismology:
Praha, Universita, Karlova. Chappell, B. W., White, A. J. R. and Hine, R., 1988, Granite provinces and
basement terranes in the Lachlan Fold Belt, Southeaster Australia: Australian Journal of Earth Sciences, 35, pp. 505-521.
Clarke, C. J., George, R. L., Bell, R. W. and Hatton, T. J., 2002, Dryland salinity
in south-western Australia: its origins, remedies and future research directions: Australian Journal of Soil Research, 40, pp. 93-113.
Collins, W. J. and Vernon, R. H., 1994, A rift-drift-delamination model of
continental evolution: Paleaozoic tectonic development of eastern Australia: Tectonophysics, 235, pp. 249-275.
Coney, P. J., 1992, The Lachlan belt of eastern Australia and Circum-Pacific
tectonic evolution: Tectonophysics, 214, pp. 1-25. Coram, J., 1998, National classification of catchments for land and river salinity
control: RIRDC Publication, No: 98/78. Coram, J., Dyson, P. and Evans, R., 2001, An evaluation framework for dryland
salinity: Department of Agriculture, Fisheries and Forestry Australia, Australian Government Publishing Service, Canberra.
Cramsie, J. N. and Hawke, J. M., 1984, Introduction In:Hawke, J. M. and
Cramsie, J. N.(Eds) Contribution to the geology of the Great Australian Basin in New South Wales: Geological survey of New South Wales, Department of Mineral Resources, Bulletin No.31, pp. 1-22.
Creelman, R., 2002, Reviewing the Science of Salinity: First Annual Symposium
on Salinity Assessment and Management, University of Sydney. Crook, K. A. W., 1980, Fore-arc evolution in the Tasman Geosyncline: the origin
of the southern Australia continental crust: Journal of the Geological Society of Australia, 27, pp. 215-232.
292
Cummings, D., 1979, Determination of depths an irregular interface in shallow refraction surveys using a pocket calculator: Geophysics, 44, pp. 1987-1998.
Dahlhaus, P. G., MacEwan, R. L., Nathan, L. and Morand, J., 2000, Salinity on
the southeastern Dundas Tableland, Victoria: Australian Journal of Earth Sciences, 47, pp. 3-11.
Dasios, A., McCann, C., Astin, T. R., McCann, D. M. and Fenning, P., 1999,
Seismic imaging of the shallow subsurface: shear-wave case histories: Geophysical Prospecting, 47, pp. 565-591.
David, V., 2003, Cross section four, Combo to Guntawang in Hill End Trough,
CD Geoscience Database: Geological Survey of New South Wales. David, V., Glen, R. A. and Spencer, R., 2003, Three-dimensional extension of
the Hill End Trough based on the modeling of the regional gravity data. In: Vassallo, J. J. G., R.A.(Ed), Evolution of the Hill End Trough: Hill End Trough CD Geoscience Data Package. Open File GS Report 2003/291, Geological Survey of New South Wales, Sydney, Australia, pp. 43-48.
Davis, H. D. and Reynolds, S. J., 1996, Structural geology of rocks and regions:
John Wiley & sons Inc., U.S.A. de Groot-Hedlin, C. and Constable, S., 1990, Occam's inversion to generate
smooth, two-dimensional models form magnetotelluric data: Geophysics, 55, pp. 1613-1624.
Deen, T. and Gohl, K., 2002, 3-D tomographic seismic inversion of a
paleochannel system in central New South Wales, Australia: Geophysics, 67, pp. 1364-1371.
Deen, T. J., Gohl, K., Leslie, C., Papp, E. and Wake-Dyster, K. D., 2000,
Seismic refraction inversion over a palaeochannel in the Lachlan fold belt, central NSW: Exploration Geophysics, 31, pp. 389-393.
Dennis, J. G., 1967, International tectonic dictionary: Am. Assoc. Pet. Geol.
Tusla, pp. 196. Dobrin, M. B., 1983, Introduction to geophysical prospecting: McGraw-Hill.Inc,
Japan. Dobrin, M. B. and Savit, C. H., 1988, Introduction to geophysical prospecting:
Domenico, S. N. and Danbon, S. H., 1986, Shear-wave technology in petroleum exploration - past, current, and future In: Dunbom, S. H. and Domenico, S. N.(Eds) Shear wave exploration, Society of Exploration Geophysicists, pp. 3-18.
Drake, C. L., 1962, Geophysics, "Geophysics" and engineering: Geophysics,
27, pp. 193-197. Duguid, J. O., 1968, Refracted determination of water table depth and alluvium
thickness: Geophysics, 33, pp. 481-488. Ebrom, D., McDonald, J. and Tatham, B., 1992, Sideswipe in a 3-D, 3-C model
survey: The Leading Edge, 11, pp. 45-51. Edge, A. G. and Laby, T. H., 1931, The principles and practice of geophysical
prospecting: Cambridge University Press, London. Edwards, L. S., 1977, A modified pseudosection for resistivity and induced
polarization: Geophysics, 42, pp. 1020-1036. English, P., 2002, Overview of salinity research at CSIRO land and water: First
Annual Symposium on Salinity Assessment and Management, University of Sydney.
Exon, N. F., 1976, Geology of the Surat Basin: Australian Government
Publishing Service, Canberra, French, S. W., 1974, Two-dimensional and three-dimensional migration of
model-experiment reflection profiles: Geophysics, 39, pp. 265-277. Gebrande, H. and Miller, H., 1985, Refraction seismology In: Bender, F.(Ed),
Applied geoscience, methods of applied geophysics and mathematical procedures in geosinces, F.Enke Publishing House, Stuttgart, pp. 226-260.
Glen, R. A. and Walshe, J. L., 1999, Cross Structure in the Lachlan Orogen
Australia: the Lachlan Transverse Zone example: Australian Journal of Earth Sciences, 46, pp. 641-658.
Glen, R. A., Walshe, J. L., Barron, L. M. and Watkins, J. J., 1998, Ordovician
convergent-margin volcanism and tectonism in the Lachlan sector of east Gondwana: Geology, 26, pp. 751-754.
Glen, R. A. and Watkins, J. J., 1999, Implications of Middle Devonian
deformation of the eastern part of the Hill End Trough, Lachlan Orogen, New South Wales: Australian Journal of Earth Sciences, 46, pp. 35-52.
Grant, F. S., 1985, Aeromagnetics, geology and ore environments, I.Magnetite
in igneous, sedimentary and metamorphic rocks, An overview: Geoexploration, 23, pp. 303-333.
294
Griffiths, D. H. and Barker, R. D., 1993, Two-dimensional resistivity imaging and modelling in areas of complex geology: Journal of Applied Geophysics, 29, pp. 211-226.
Gunn, P. J., 1997a, Regional magnetic and gravity responses of extensional
sedimentary basins: AGSO Journal of Australian Geology and Geophysics, 17, pp. 115-131.
Gunn, P. J., 1997b, Application of aeromagnetic surveys to sedimentary basin
studies: AGSO Journal of Australian Geology and Geophysics, 17, pp. 133-144.
Gunn, P. J. and Dentith, M. C., 1997, Magnetic responses associated with
mineral deposits: AGSO Journal of Australian Geology and Geophysics, 17, pp. 145-158.
Haeni, F. P., 1986, Application of seismic refraction methods in groundwater
modeling studies in New England: Geophysics, 51, pp. 236-249. Hasselstrom, B., 1969, Water prospecting and rock-investigation by the seismic
refraction method: Geoexploration, 7, pp. 113-132. Hatcher, R. D., 1995, Structure geology: principles, concepts, and problems:
Prentice-Hall, Inc., U.S.A. Hatherly, P. J., 1986, Attenuation measurements on shallow seismic refraction
data: Geophysics, 51, pp. 250-254. Hatherly, P. J. and Neville, M. J., 1986, Experience with the generalized
reciprocal method of seismic refraction interpretation for shallow engineering site investigation: geophysics, 51, pp. 255-265.
Hawke, J. M. and Cramsie, J. N., 1984, Contribution to the Geology of the Great
Australian Basin in New South Wales: Geological Survey of New South Wales, Bulletine 31, pp. 259.
Hawkins, L. V., 1961, The reciprocal method of routine shallow seismic
refraction investigations: Geophysics, 26, pp. 806-819. Heelan, P. A., 1953, On the theory of head waves: Geophysics, 18, pp. 871-
893. Helbig, K., 1986, Shear-waves-what they are and how they can be used
In:Danbom, S. H., Domenico, S.N.(Ed), Shear-wave exploration: Society of Exploration Geophysicists, pp. 19-36.
Heneberg, W. C., 1995, Steady state groundwater flow across idealized faults:
Water Resources Research, 31, pp. 1815-1820. Hobbs, B. E., 1976, An outline of structural geology: Wiley & Sons, Inc, pp.571. 295
Hobson, G. D. and Grant, A. C., 1964, Tracing buried river valleys in the Kirkland Lake area of Ontario with a hammer seismograph: Canadian Mining Journal, pp. 79-83.
Hooper, E. C. D., 1991, Fluid migration along growth faults in compacting
sediments: Journal of Petroleum Geology, 14, pp. 161-180. Horsfall, K. R., 1997, Airborne magnetic and gamma-ray data acquisition:
AGSO Journal of Australian Geology and Geophysics, 17, pp. 23-30. Hull, J., 1988, Thickness-displacement relationships for deformation zones:
J.Struct.Geol., 10, pp. 431-435. IAEA, 1991, Airborne gamma ray spectrometer surveying: International Atomic
Energy Agency, Technical Reports Series, No. 323. Irving, F. R., 1965, Seismic surveying methods, equipment and costs in New
York State: Highway Research Record No.81, Highway Research Board, Washington, pp. 2-8.
Jacobsen, T. and Adams, R. M., 1958, Salt and silt in ancient Mesopotamian
agriculture: Science, 128, pp. 1251-1258. Jaques, A. L., Wellman, P., Whitaker, A. and Wyborn, D., 1997, High-resolution
geophysics in modern geological mapping: AGSO Journal of Australian Geology and Geophysics, 17, pp. 159-173.
Jewel, C. M. and Associates, 2000, Hydrogeological study for Spicers Ck /
Talbragar catchment: Jewel CM and Associates Pty Ltd Rep No A59967.4, April 2000.
Kearey, P. and Brooks, M., 2001, An introduction to geophysical exploration:
Blackwell Scientific Publications. Kilthy, K. T., Norris, R. A., McLamore, W. R., Hennon and K.P.Euge, K., 1986,
Seismic refraction at Horse Mesa Dam: an application of the generalized reciprocal method: Geophysics, 51, pp. 266-275.
Kingham, R., 1998, Geology of the Murray-Darling Basin - Simplified
Lithostratigraphic Groupings: Australian Geological Survey Organisation, Record 1998/21.
Knight, M. J., Saunders, B. J., Williams, R. M. and Hillier, J., 1989, Geologically
induced salinity at Yelarbon, Border River area, New South Wales, Queensland: BMR Journal of Australian Geology and Geophysics, 11, pp. 355-361.
Knipe, R. J., 1993, The influence of fault zone processes diagenesis on fluid
flow In:Horbury, A. D. and Robonson, A. G.(Eds), Diagenesis and basin development studies in geology: AAPG Bulletin, pp. 135-151.
296
Knipe, R. J., 1997, Juxtaposition and seal diagrams to help analyze fault seals in hydrocarbon reservoirs: AAPG Bulletin, 81, pp. 187-195.
Korsch, R. J., Wake-Dyster, K. D. and Johnstone, D. W., 1993, Deep seismic
reflection profiling in the Gunnedah basin: regional tectonics and basin responses In: Swarbick, G. J. and Morton, D. J. (Eds), Proceedings NSW Petroleum Symposium: Petroleum Exploration Society of Australia, NSW Branch.
Korvin, G., 1982, Axiomatic characterization of the general mixture rule:
Geoexploration, 19, 267-276. Laa, A., Winit, Y., Chalong, B., Pansak, W., Somyot, H. and Ladda, W., 1991,
Integrated geothecnical and geophysical investigations of the groundwater dam at Ban Wa, case study: Khon Kaen, Thailand,
Lawson, C. E., Foster, W. R. and Mitchell, R. E., 1965, Geophysical equipment
usage in the Wiscons in highway commission organization: Highway Research Record No.81, Highway Research Board, Washington, pp. 42-48.
Lecomte, I., Gjoystdal, H., Dahle, A., Pedersen, O.C., 2000, Improving modeling
and inversion in refraction seismics with a first-order Eikonal solver: Geophysical Prospecting, 48, 437-454.
Lennox, D. H. and Carlson, V., 1967a, Geophysical exploration for buried
valleys in an area north of Two Hills, Alberta: Geophysics, 32, pp. 331-362.
Lennox, D. H. and Carlson, V., 1967b, Integration of geophysical methods for
groundwater exploration in the Prarie provinces, Canada: Geol. Surv. Canada, Economic Geology Report No. 26, pp. 517-535.
Levshin, A. L., 1961, Determination of groundwater level by the seismic
methods: Akad. Nauk. SSSR Izv. Ser. Geofiz, 9, pp. 857-870 (Eng. Trans.).
Lewis, C., Davis, T. L. and Vuillermoz, C., 1991, Three-dimensional
multicomponent imaging of reservoir heterogeneity, Silo Field, Wyoming: Geophysics, 56, pp. 2048-56.
Lewis, M. F., 1991, Lineaments and salinity in western Australia - carriers or
barriers?: International Hydrology and Water Resource Symposium, Perth-Australia.
Lister, G. S. and Snoke, A. W., 1984, S-C mylonites: Journal of Structural
Geology, 6, pp. 617-638. 297
Lopez, D. L. and Smith, L., 1995, Fluid flow in fault zones: analysis of the interplay of convective circulation and topographically driven groundwater flow: Water Resources Research, 31, pp. 1489-1503.
56, 645-653. Mallot, D. F., 1965, Application of geophysics to highway engineering in
Michigan: Highway Research Record No.81, Highway Research Board, Washington, pp. 20-41.
Martelat, J. E., Vidal, G., Lardeaux, J. M., Nicollet, C. and Rakotondrazafy, R.,
1995, Images spatiales et tectoniques profoundes des continents, l`exemple du Sud-Ouest de Madagascar: C.R.Acad.Paris, 321, pp. 325-332, (Eng. Trans.).
Matson, C. R., 1975, Mine data sheets to accompany metallogenic map, Dubbo
1:250 000 sheet: Geological Survey of New South Wales. Mawson, R. and Talent, J. A., 1999, The Devonian of eastern Australia:
stratigraphic alignments, stage and series boundaries, and the transgression - regression pattern reconsidered: Courier Forschungsinstitute Senckenberg, 225, pp. 243-270.
McNeill, J. D., 1990, The use of electromagnetic methods for groundwater
studies In:Ward, S. H.(Ed), Geotechnical and environmental geophysics, Volume I review and tutorial: Society of Exploration Geophysicists, Tulsa, pp. 191-218.
Meakin, N. S. and Morgan, E. J. c., 1999, Dubbo Geological sheet 1:250,000,
Means, W. D., 1995, Shear zones and rock history: Tectonophysics, 274, pp.
157-160. Mickus, K. L., Aiken, C. L. and Kennedy, W. D., 1991, regional-residual gravity
anomaly separation using the minimum-curvature technique: Geophysics, 56, pp. 279-283.
Miles, D. R., 1988, Three-component seismic data: Researcher's toy or
interpreter's tool?: The Leading Edge, 7, pp. 28-31. Miller, T. W., 1995, New insight on natural hydraulic fractures induced by
abnormally high pore pressures: AAPG Bulletin, 79, pp. 1005-1018. Milligan, P. R. and Gunn, P. J., 1977, Enhancement and presentation of
airborne geophysical data: AGSO Journal of Australian Geology and Geophysics, 17, pp. 63-75.
298
Minty, B. R. S., 1997, Fundamentals of airborne gamma-ray spectrometry: AGSO Journal of Australian Geology and Geophysics, 17, pp. 39-50.
Mitra, G., 1992, Deformation of granitic basement rocks along fault zones at
shallow or intermediate crustal levels In:Mitra, S. and Fisher, G. W.(Eds), Structural geology of fold and thrust belts: Johns Hopkins University Press, Baltimore, MD, pp. 123-144.
Moore, R. W., 1952, Geophysical methods adapted to highway engineering
problems: Geophysics, 17, pp. 505-530. Morgan, E. J. E., 1997, Geology and Metallogenesis of the Fifield-Peak Hill-
Wellington-Gulgong-Mudgee Region of New South Wales: Narromine and Dubbo 1:250,000 Geological Sheets (New South Wales Geological Survey, Sydney): New South Wales Department of Mineral Resources, Unpublished Report No.1997/080.
Morgan, K. and Jankowski, J., 2002, Determination of the source of salinity
within the dryland salinity affected Spicers Creek catchment, central west, NSW, Australia: Proceedings of the 27th Hydrology and Water Resources Symposium, Melbourne, Australia.
Morgan, K. and Jankowski, J., 2004, Saline groundwater seepage zones and
their impact on soil and water resources in the Spicers Creek Catchment, central west, New South Wales, Australia: Environmental Geology, 46, pp. 273-285.
Musgrave, A. W., 1970, Seismic refraction prospecting: Society of Exploration
Geophysics, pp.604. National Land and Water Resources Audit, 2001, Dryland salinity in Australia: a
summary of the National Land and Water Resources Audit’s Australian Dryland Salinity Assessment 2000: extent, impacts, processes, monitoring and management options: National Land & Water Resources Audit, ACT, Australia.
Nestvold, E. O., 1992, 3-D Seismic: Is the promise fulfilled?: The Leading Edge,
11, pp. 12-19. Nicoll, C. and Scown, J., 1993, Dryland salinity in the Yass Valley, processes
and management: Department of Conservation and Land Management, Australian Government Publishing Service, Canberra.
Nikrouz, R. and Palmer, D., 2004, 3-D 3-C Seismic Refraction Imaging of Shear
Zone Sources of Dryland Salinity: 17th geophysical conference and exhibition, 15-19 August 2004, Sydney – Australia.
O'Brien, P. N. S., 1967, The use of amplitudes in seismic refraction survey In:
Musgrave, A. W.(Ed), Seismic refraction prospecting, Society of Exploration Geophysicists, Tulsa.
299
Overmeeren, R. A. V., 1981, A combination of electrical resistivity, seismic refraction, and gravity measurements for groundwater exploration in Sudan: Geophysics, 46, pp. 1304-1313.
Packham, G. H., 1968, The lower and middle Palaeozoic stratigraphy and
sedimentary tectonics of the Sofala- Hill End- Euchareena region, NSW.: Proceedings of the Linnean Society of New South Wales., pp. 111-163.
Packham, G. H., 1999, Radiometric evidence for Middle Devonian inversion of
the Hill End Trough, Northeast Lachlan fold belt.: Australian Journal of Earth Sciences, 46, pp. 23-33.
Packham, G. H., 2003, Hill End Trough : Genesis and fill In:Vassallo, J. J. and
Glen, R. A.(Eds), Evolution of the Hill End Trough: Hill End Trough CD Geoscience Data package. Open File GS Report 2003/291: Geological Survey of New South Wales, Sydney, Australia, pp. 1-11.
pakiser, L. C. and Black, R. A., 1957, Exploration for ancient channels with the
refraction seismograph: Geophysics, 22, pp. 32-47. Palmer, D., 1974, An application of the time section in shallow seismic
refraction studies: Master's thesis, The University of Sydney, pp. 157. Palmer, D., 1979, What is the future for seismic refraction methods?: Bull.
Austral. SEG, 10, pp. 215-217. Palmer, D., 1980, The generalized reciprocal method of seismic refraction
interpretation: Society of Exploration Geophysicists. Palmer, D., 1981, An introduction to generalized reciprocal method of seismic
refraction interpretation: Geophysics, 46, pp. 1508-1518. Palmer, D., 1986, Refraction seismic, The lateral resolution of structure and
seismic velocity: Geophysical press. Palmer, D., 1990, The generalized reciprocal method - an integrated approach
to shallow refraction seismology: Exploration Geophysics, 21, pp. 33-44. Palmer, D., 1991, The resolution of narrow low-velocity zones with the
generalized reciprocal method: Geophysical Prospecting, 39, pp. 1031-1060.
Palmer, D., 1992, Is forward modelling as efficacious as minimum variance for
refraction inversion?: Exploration Geophysics, 23, pp. 261-266. Palmer, D., 2000a, Can amplitude resolve ambiguities in refraction inversion?:
Exploration Geophysics, 31, pp. 304-309. 300
Palmer, D., 2000b, Can new acquisition methods improve signal-to-noise ratios with seismic refraction techniques?: Exploration Geophysics, 31, pp. 275-300.
Palmer, D., 2000c, The measurement of weak anisotropy with the generalized
reciprocal method: Geophysics, 65, pp. 1583-1591. Palmer, D., 2001a, Amplitude "statics" in shallow refraction seismology: ASEG
15th Geophysical Conference and Exhibition, Brisbane. Palmer, D., 2001b, Digital processing of shallow seismic refraction data with the
refraction convolution section: PhD’s thesis, School of Geology, New South Wales, Sydney, Australia, pp. 256.
Palmer, D., 2001c, Imaging refractors with the convolution section: Geophysics,
66, pp. 1582-1589. Palmer, D., 2001d, Measurement of rock fabric in shallow refraction seismology:
Exploration Geophysics, 32, 307-314. Palmer, D., 2001e, Resolving refractor ambiguities with amplitudes:
Geophysics, 66, pp. 1590-1593. Palmer, D., 2001f, Shallow refraction seismology for the new millennium: ASEG
15th Geophysical Conference and Exhibition, Brisbane. Palmer, D., 2001g, A simple approach to 3-D shallow refraction seismology:
ASEG 15th Geophysical Conference and Exhibition, Brisbane. Palmer, D., 2003a, Application of amplitudes in shallow seismic refraction
inversion: ASEG 16th Geophysical Conferences and Exhibition, Adelaide.
Palmer, D., 2003b, Processing and interpreting shallow seismic refraction data
with the GRM & the CRM, Unpublished: University of New South Wales, pp. 154.
Palmer, D., Goleby, B. and Drummond, B., 2000, The effect of spatial sampling
on refraction statics: Exploration Geophysics, 31, pp. 270-274. Palmer, D., Nikrouz, R. and Spyrou, A., 2005, Statics Corrections for Shallow
Seismic Refraction Data: In prep.. Passchier, C. W. and Simpson, C., 1986, Porphyroclast systems as kinematic
indicators: Journal of Structural Geology, 8, pp. 831-843. Paterson, N. R., 1965, Seismic wave propagation in porous granular media:
Geophysics, 21, pp. 691-714. 301
Peck, A. J., Thomas, J. F. and Williamson, D. R., 1983, Salinity Issues: Australian Government Publishing Service, Canberra.
Philips, O. M., 1991, Flow and Reaction in Permeable Rocks: Cambridge
University Press, Cambridge, UK. Pili, E., Ricard, Y., Lardeaux, J. M. and Sheppars, S. M. F., 1997, Lithospheric
shear zones and mantle-crust connections: Tectonophysics, 280, pp. 15-29.
Pires, A. C. B., Harthhill, N., 1989, Statistical analysis of airborne gamma-ray
data for geological mapping purposes: Crixas-Itapaci area, Goias, Brazil: Geophysics, 54, pp. 1326-32.
Pogson, D. J. and Watkins, J. J., 1998, Buthurst Geological sheet SI/55-8:
Explanatory Notes: Geological Survey of New South Wales, Sydney. Powell, C. M., 1984, Ordovician to earliest Silurian: marginal sea and island arc,
In: Veevers, J. J.(Ed), Phanerozoic of earth history of Australia: Clarendon Press, Oxford, pp. 290-309.
Powell, C. M., Edgecombe, D. R., Henry, N. M. and Jones, J. G., 1976, Timing
of regional deformation of the Hill End Trough: a reassessment: Journal of the Geological Society of Australia, 23, pp. 407-421.
Power, P. E. and Devine, S. B., 1970, Surat Basin, Australia-subsurface
stratigraphy, history, and petroleum: AAPG Bulletin, 54, pp. 2410-2437. Qin, F., Luo, Y., Olsen, K., Cai, W., Schuster, G., 1992, Finite-difference
solution of the eikonal equation along expanding wavefront: Geophysics, 57, pp.478-487.
Ramsay, J. G., 1980, Shear zone geometry, a review: Journal of Structural
Geology, 2, pp. 83-99. Ramsay, J. G. and Auber, M., 1983, Modern structural geology: Academic
Press, London. Ramsay, J. G. and Lisle, R. J., 2000, Modern structural geology (volume 3):
Academic Press, UK. Reeve, R. C. and Fireman, M., 1967, Salt problem in relation to irregation.
In:Hagan, R. M., Haise, H. R. and Edminister, T. W.(Eds), Irrigation of agriculture lands: American Society of Agronomy, USA.
Reeves, C. V., 1992, New horizons for airborne geophysics mapping:
Exploration Geophysics, 23, pp. 273-280. Reston, T. J., 1990, Mantle shear zones and the evolution of the northern North
Sea basin: Geology, 18, pp. 272-275. 302
Reynolds, M. R., 1997, An introduction to applied and environmental geophysics: John Wiley & Sons, Inc., New York, USA.
Richards, T.C., 1959, Broadside refraction shooting: Geophysics, 24, 725-748. Roberts, S. J., Nunn, J. A., Cathles, L. and Cipriani, F. D., 1996, Expulsion of
abnormally pressured fluids along faults: Journal of Geophysical Research, 101, pp. 28231-28252.
Rutland, R. W. R., 1976, Orogenic evolution of Australia: Earth Science
Reviews, 12, pp. 161-196. Rutter, E. H., 1999, On the relationship between the formation of shear zones
and the form of the flow law for rocks undergoing dynamic recrystalization: Tectonophysics, 303, pp. 147-158.
Sasaki, Y., 1992, Resolution of resistivity tomography inferred from numerical
simulation: Geophysical Prospecting, 40, pp. 453-464. Scheibner, E., 1973, A plate tectonic model of the palaeozoic teconic history of
New Souht Wales: Journal of the Geological Society of Australia, 20, pp. 405-426.
Scheibner, E. and Basden, H. e., 1996, Geology of New South Wales-
Synthesis. Volume 1 Structural Framework: Geological Survey of New South Wales, Memoir Geology, 13(1), pp. 295.
Schofield, S. and Jankowski, J., 2000, Hydrogeology of Billimore Region,
Central New South Wales, Australia. In: O.Sililo et al. (Eds): Proc. Cong. IAH, Cape Town, Africa.
Schofield, S. K., 1998, The geology, hydrogeology and hydrogeochemistry of
the Ballimore Region central New South Wales, Phd thesis: University of New South Wales, pp. 453.
Schuster, G. T. and Quintus-Bosz, A., 1993, Wavepath eikonal traveltime
inversion, Theory: Geophysics, 58, pp. 1314-1323. Scott, J. H., 1973, Seismic refraction modeling by computer: Geophysics, 38,
pp. 271-284. Scott, M. M., 1997, Structural development of the Dubbo sheet In:Morgan, E.
J.(Ed), Geology and Metallogenesis of the Fifield-Peak Hill-Wellington-Gulgong-Mudgee Region of New South Wales: Narromine and Dubbo 1:250,000 Geological Sheets (New South Wales Geological Survey, Sydney): New South Wales Department of Mineral Resources, Report 1997/080.
Shepard, E. R., 1939, The seismic method of exploration applied to construcyion projects: The Military Engineer, 31.
Sheriff, R. E. and Geldart, L. P., 1995, Exploration Seismology: Cambridge
University Press, New York, USA. Sibson, R. H., 1996, Structural permeability of fluid-driven fault-fracture meshes:
Journal of Structural Geology, 18, pp. 1031-1042. Simpson, C. and Schmid, S. M., 1983, An evaluation of criteria to deduce the
sense of movement in sheared rocks: Geological Society of America Bulletin, 94, pp. 1281-1288.
Sjogren, B., 1980, The law of parallelism in refraction shooting: Geophysical
Prospecting, 28, pp. 716-743. Sjogren, B., 1984, Shallow refraction seismic: Chapman and Hall Ltd, New
York, USA. Smith, R.E, 1969, Zones of progressive regional burial metamorphism in part of
the Tasman Geosyncline, eastern Australia: Journal of Petrology, 10, 144-163.
Smith, D. A., 1980, Sealing and nonsealing faults in Louisiana Gulf Coast salt
basin: AAPG Bulletin, 64, pp. 145-172. Smith, D. N. and Ferguson, J. F., 2000, Constrained inversion of seismic
refraction data using the controlled random search: Geophysics, 65, pp. 1622-1630.
Stam, J. C., 1962, Modern developments in shallow seismic refraction
techniques: Geophysics, 27, pp. 198-212. Stockwell, J. W., Cohen, J. K. and Cohen, J., 2000, The New SU User's
Manual: The Society of Exploration Geophysicists, pp. 140. Street, G. J. and Engel, R., 1990, Geophysical Surveys of Dryland Salinity
In:Ward, S. H.(Ed), Geotechnical and environmental geophysics: Society of Exploration geophysicists.
Suppel, D. W. and Scheibner, E., 1990, Lachlan Fold Belt in New South Wales -
Regional Geology and Mineral Deposits In: Hughes, F. E.(Ed), Geology of the mineral Deposits of Australia and Papua New Guinea: The Australian Institute of Mining and Metallurgy, Melbourne, pp. 1321-1327.
Svenson, D. and Bowering, R. H., 1963, The engineering seismograph and its
use by the Snowy Mountain Authority: 4th Aust.-N.Z. Conf. Soil Mech. and Foundation Engineering.
304
Swain, C. J., 1992, The Kenya rift axial gravity high, a reinterpretation: Tectonophysics, 204, pp. 59-70.
Tadros, N. Z. E., 1993, The Gunnedah Basin, New South Wales: Geological
Survey of New South Wales, Memoir Geology 12, 649 pp. Taner, M. T., Wagner, D. E., Baysal, E. and Lu, L., 1998, A unified method for
the 2-D and 3-D refraction statics: Geophysics, 63, pp. 260-274. Tatham, R. H., 1982, Vp/Vs and lithology: Geophysics, 47, pp. 336-344. Tatham, R. H., 1985, Shear waves and lithology In:Helbig, K., Treitel, S.(Ed),
Seismic shear waves, Geophysical Press Limited, pp. 273. Tatham, R. H. and McCormack, M. D., 1991, Multicomponent Seismology in
Petroleum Exploration: Tusla: Society of Exploration Geophysicists. Telford, W. M., Geldart, L. P., Sheriff, R. E. and Keys, D. A., 1976, Applied
geophysics: Cambridge University Press, London. Toksoz, M. N., Cheng, C. H. and Timur, A., 1976, Velocities of seismic waves in
porous rocks: Geophysics, 47, pp. 621-645. Urquhart, D. F., 1956, The investigation of deep leads by the seismic refraction
method: Australia, Bureau of Mineral Resourses, 35. Vassalo, J. J., Glen, R. A. and David, V., 2003, Structural pathways and
locations of gold in the Hill End Trough, eastern Lachlan Orogen, NSW. In:Vassallo, J. J. and Glen, R. A.(Eds), Evolution of the Hill End Trough: Hill End Trough CD Geoscience Data package. Open File GS Report 2003/291: Geological Survey of New South Wales, Sydney, Australia, pp. 26-37.
Vidale, J., 1988, Finite-difference calculation of traveltimes: Bull.Sci.Soc.Am.,
78, 2062-2076. Vissers, R. L. M., Drury, M. R., Strating, E. H. H. and Van der Wal, D., 1991,
Shear zones in the upper mantle: a case study in an alpine type lherzolite massif: Geology, 19, pp. 990-993.
Walker, A. N., Rutter, E. H. and Brodie, K. H., 1990, Experimental study of grain
size sensetive flow of synthetic, host-pressed calcite rocks. In:Knipe, R. J. and Rutter, E. H.(Eds), Deformation Mechanisms, Rheology and Tectonics: Geol. Soc. London Spec. Publ., pp. 259-284.
Waltham, D., 1997, Why does salt start to move?: Tectonophysics, 282, pp.
117-128. Walton, G. G., 1972, Three-dimensional seismic method: Geophysics, 37, pp.
417-430. 305
Ward, S. H., 1990, Resistivity and induced polarisation methods In:Ward, S. H.(Ed), Geotechnical and environmental geophysics, Volume I review and tutorial: Society of Exploration Geophysicists, Tulsa, pp. 147-190.
Warrick, R. E. and Winslow, J. E., 1960, Application of seismic methods to a
groundwater problem in north-eastern Ohio: Geophysics, 25, pp. 505-519.
Weimer, P. and Davis, T. L., 1996, Applications of 3-D seismic data to
exploration and production: AAPG & SEG, Tulsa, Oklahoma, U.S.A. Welin, E., 1958, Seismic refraction survey of a hydro-electric plant in northern
Sweden, Geophysical surveys in mining, hydrological and engineering projects: E.A.E.G., pp. 262-270.
Werth, G. A., 1967, Method for calculating the amplitude of the refraction arrival
In:Musgrave, A. W.(Ed), Seismic refraction prospecting: Society of Exploration Geophysics, pp. 119-137.
Wiebenga, W. A. and Jasson, E. E., 1962, Geophysical exploration for
underground water: Burear of Mineral Resources, Australia, Record No.1962/172.
Wilford, J. R., Bierwirth, P. N. and Craig, M. A., 1997, Application of airborne
gamma-ray spectrometry in soil/regolith mapping and applied geomorphology: AGSO Journal of Australian Geology and Geophysics, 17, pp. 201-216.
Wood, A. E., 1940, Damsite surveying by seismograph: Engineering News-
Record, 124, pp. 438-441. Woodward, M., Rocca, F., 1988, Wave-equation tomography: 58th Ann.Internal.
Mtg., SEG, extended abstracts, 1232-1235. Xia, J., Miller, R. D., Park, C. B., Wightman, E. and Nigbor, R., 2002, A pitfall in
Yielding, G., Freeman, B. and Needham, D. T., 1997, Quantitative fault seal
prediction: AAPG Bulletin, 81, pp. 897-917. Yilmaz, O., 1988, Seismic data processing: Society of Exploration Geophysics. Young, R. A. and Lucas, J. E., 1988, Exploration beneath volcanics, Snake
River plain, Idaho: Geophysics, 53, pp. 444-452.
306
Appendix A: Linear and Non-Linear Method The linear method of refractor velocity determination can be a subjective
process since any number of straight lines could be justifiably fitted onto the
velocity analysis function. This was especially the case when the derivative of
the gradient changed rapidly, that is when function had a considerable concave
upwards or concave downward shape. Also, with the linear method, lateral
changes in the refractor velocity occur abruptly, over one receiver spacing.
In order to introduce some objectivity into the determination of the refractor
velocities, and to obtain smoother lateral changes, a non-linear method was
used as well. This method uses a series of numerical differentiation formulas
(1-3) to approximate the derivative at each point along the velocity analysis
function.
311
)1(0 6
1)(21 µδ−−≈ −TThT (1)
52112 30
1)88(121 µδ−+−+−≈ −− TTTT (2)
7321123 140
1)945459(601 µδ−−+−+−≈ −−− TTTTTT (3)
Formulas 1, 2 and 3 represent three-point, five-point and seven-point midpoint
formula respectively, where h is the node separation, Tn is the value of the
function at node n, Tn(1) is the first derivative of the function at node n, and µδn is
the order of the differences constituting the error term. The mid-point formulas,
together with a series of end-point formulas are derived by differentiating
Lagrangian interpolation formulas (Beyer, 1975).
307
308
With the non-linear method, refractor velocities were calculated by substituting
the values of the velocity analysis function into the different numerical
differentiation formulas, to give the derivative at each point along the function
and then taking the reciprocal to give the refractor velocity at each of these
points. For example using the three-point midpoint formula (1) the reciprocal of
the refractor velocity at a centre point T0(1) is approximately half the difference of
the velocity analysis function at the point to its right T1 and the point to its left T-1
divided by the point separation h. The velocities obtained using the different
formulas were similar, however the smoothness of the velocity trends tended to
increase as the number of points used in the formula increased.
Reference: Beyer, W. H., 1975, Handbook of mathematical science, 5th edition: CRC Press, Florida.
****Date: 2003/04/07 Time File # Comments 14:08:26 1 10 To 200 Hz 4 X 8 Second Sweeps 14:15:28 2 10-120 Hz 4 X 8 Seconds 14:21:58 3 10-250 Hz 4 X 8 Seconds 14:25:48 4 20-250 Hz 4 X 8 Seconds 14:29:18 5 30-250 Hz 4 X 8 Second Sweeps 14:33:21 6 20-250 Hz 9 X 8 Second Sweeps 14:40:01 7 20-250 HZ 9 X 8 Second Sweeps Correlate before stack 14:44:45 8 Start Production with 20-250 Hz 4 X 8 Second Sweep 16:05:53 26 Start of next line Sp 1 To 19 had Vertical swapped with horizontal 2 stations 199 and 201 on line 16:08:42 27 1 with 319 and 321 on line 3 16:19:30 28 Source Line 1 shot 20m north of SP 1 16:21:30 29 Source Line 1 shot 40M off SP 1 16:23:03 30 Source Line 1 shot 60M off SP 1 16:43:53 31 Back to Source line 2 Sp number reset to 25 this is the second shot on this line 17:48:41 49 Start of Source Line 3 18:20:14 64 Last Shot ****Date: 2003/04/08 ---------------------------------------------------------------------------------------------------- Time File # Comments 09:43:16 65 Start of Day 2 and Start of Source Line 4 10:14:47 73 Start of Line 5 10:52:59 87 Start of Line 6 14:44:49 90 Test Shot with One Det - Void Shot Record Length at 1500 Mil 14:46:16 91 Test Shot with One Det Record Length at 1500 Mil 15:19:27 92 Record Length at 1500 Mil - Explosives 15:45:34 93 Record Length at 1500 Mil - Explosives
318
Appendix F: Cont. 16:31:55 94 Record Length at 1000 Mil - Explosives ****Date: 2003/04/09 Time File # Comments 12:32:31 95 Start of S Wave Tests. Top of Mass on the left of the truck 12:15:33 114 Source line 1 the extra shots on the end 12:23:50 117 Mass Rotated 180 D. Top of Mass on the right of the truck 12:39:28 120 Back at the orginal start of line S1 12:41:14 121 Void file 120 13:12:38 140 Void file 138 14:12:33 141 Start of source line 1. Top of Mass at the front of the truck 14:57:46 163 Start of Line S1. Top of Mass at the back of the truck 15:33:36 181 Void file 180 15:58:39 186 Start of source line 2. Top of Mass at the back of the truck 16:22:25 201 Void file 199 16:33:08 203 Start of Source Line 2. Top of Mass at the front of the truck 17:06:31 219 Start of line S2. Top of Mass to the right of the truck 17:33:18 235 Start of Source S2. Top of Mass to the left of the truck. All references are from the back of the truck facing forward ****Date: 2003/04/10 Time File # Comments 09:21:10 251 Start of Diagonals. Top of Mass to the top of the truck. 09:24:34 252 A bit of wind noise 16:16:06 270 Top of Mass to the back of the truck. Shear Waves Start of Area 1. 3 Input Lines Used with 216 channels per line and 2 Aux Chan 16:17:38 271 High Cut filter increased to 328 16:18:57 272 Windy conditions 16:51:11 290 Start of next line. Top of Mass at the back of the truck 16:56:27 294 wind easing 17:16:02 305 Void File 304
319
Appendix F: Cont. ****Date: 2003/04/11 Time File # Comments 09:08:09 306 Start of next tests Top of Mass at the back of the truck 09:10:40 307 Windy Conditions 09:36:34 319 A few rain drops falling 10:10:14 321 Start tests again. Mass to the front of the truck 10:56:08 353 Receivers 162 on line 2 and 3 were clipped on the wrong line. Now swapped back 11:04:37 357 Light rain again 11:26:52 365 Very light rain 11:53:55 371 Start tests again. Top of Mass to the right of the truck 12:49:16 407 Still raining 14:56:46 409 Void file 408 - Rain Noise. Reshoot 38 later after some drying time ****Date: 2003/04/12 Time File # Comments 23:20:08 415 Start of tests. Top of mass on the left of the truck 10:49:19 465 Void files 403 To 414 Reshoot part of yesterday. Top of mass to the right of the truck 10:50:57 466 Reshoot part of yesterday. Top of mass to the right of the truck 11:21:20 481 Void file 480 missed sweep 11:26:29 483 Last shot for shear wave 11:49:58 484 Changed to 4 Aux Traces Explosive Shot 330g 11:55:52 485 Void File 484 Explosive Shot 330g 13:03:11 486 Void File 485 Explosive Shot 330g 13:24:03 487 Explosive Shot 220g 13:35:34 488 Explosive Shot 220g 13:45:18 489 Explosive Shot 110g 13:54:03 490 Explosive Shot 220g 14:01:24 491 Explosive Shot 330g
320
14:15:01 492 Explosive Shot 440g 14:24:51 493 Explosive Shot 220g 14:44:26 494 Last Explosive Explosive Shot 220g 15:03:10 495 Start of P wave tests 16:16:34 544 End of Spicer Program
321
Appendix G: Assessment of Data Quality
Site 1 Source: Vertical
Rating (1-5) File
number Source Point #
Source Line #
SourceOrientation(V, S1, S2)
SourceLocation(in/out
spread)
Distance to
NearestDetector V H1 H2 Overall
Use (Yes/No)
495 1 S2 V out 100 4 3 2 3 Yes 496 2 S2 V out 80 4 3 2 3 Yes 497 3 S2 V out 60 4 2 2 3 Yes 498 4 S2 V out 40 4 2 2 3 Yes 499 5 S2 V out 20 4 2 2 3 Yes 500 6 S2 V in 5 4 2 2 3 Yes 501 7 S2 V in 5 4 2 1 2 Yes 502 8 S2 V in 5 3 1 1 2 Yes 503 9 S2 V in 5 3 2 1 2 Yes 504 10 S2 V in 5 4 4 3 4 Yes 505 11 S2 V in 5 4 2 1 2 Yes 506 12 S2 V in 5 4 2 2 3 Yes 507 13 S2 V in 5 4 3 3 3 Yes 508 14 S2 V out 20 4 3 3 3 Yes 509 15 S2 V out 40 4 3 3 3 Yes 510 16 S2 V out 60 3 2 2 2 Yes 511 17 S2 V out 80 3 3 3 3 Yes 512 18 S2 V out 100 3 2 2 2 Yes 513 19 S2 V out 120 3 3 2 3 Yes 514 20 S2 V out 140 3 2 2 2 Yes 515 21 S2 V out 100 1 1 2 1 Yes 516 22 S2 V out 80 1 1 2 1 Yes 517 23 S2 V out 60 2 1 1 1 Yes 518 24 S2 V out 40 2 1 1 1 Yes 519 25 S2 V out 20 3 1 1 2 Yes 520 26 S2 V in 0 4 1 1 2 Yes 521 27 S2 V in 0 4 2 1 2 Yes 522 28 S2 V out 20 4 1 1 2 Yes 523 29 S2 V out 40 3 1 2 2 Yes 524 30 S2 V out 60 3 1 1 2 Yes 525 31 S2 V out 80 3 1 1 2 Yes 526 32 S2 V out 100 3 1 1 2 Yes 527 33 S2 V out 120 3 1 1 2 Yes 528 34 S2 V out 140 3 2 2 2 Yes 529 35 S2 V out 160 3 2 2 2 Yes 530 36 S2 V out 60 4 3 3 3 Yes 531 37 S2 V out 60 4 3 3 3 Yes 532 38 S2 V out 60 3 2 2 2 Yes 533 39 S2 V out 20 3 2 2 2 Yes
322
Cont. Site 1 Source: Vertical
Rating (1-5) File
number Source Point #
Source Line #
Source Orientation(V, S1, S2)
SourceLocation(in/out
spread)
Distance to
NearestDetector V H1 H2 Overall
Use (Yes/No)
534 40 S2 V in 0 3 3 3 3 Yes 535 41 S2 V in 0 3 3 3 3 Yes 536 42 S2 V out 20 3 2 2 2 Yes 537 43 S2 V out 20 3 2 1 2 Yes 538 44 S2 V in 0 4 2 3 3 Yes 539 45 S2 V in 0 3 3 2 3 Yes 540 46 S2 V out 20 2 2 2 2 Yes 541 47 S2 V out 30 4 2 2 3 Yes 542 48 S2 V out 20 4 3 3 3 Yes 543 49 S2 V out 20 4 4 3 4 Yes 544 50 S2 V out 30 4 2 2 3 Yes
Depth Section(P-wave), XY=5 - Site2, Line 4, S22-S19
0
5
10
15
20
0 20 40 60 80 100 120 140
Distance (m)
Dep
th (m
)
XY=5average
Time Depth (H2W) - Site2, Line 1, S22-S19
5
10
15
20
25
30
35
40
45
50
55
0 20 40 60 80 100 120 140
Distance (m)
Tim
e (m
s)
XY=-10
XY=-5
XY=0
XY=5
XY=10
XY=15
XY=20
384
Appendix L (cont.)
Time Depth (H2W) -Site2, Line 2, S22-S19
5
10
15
20
25
30
35
40
45
0 20 40 60 80 100 120 140
Distance (m)
Tim
e (m
s)
XY=-10
XY=-5
XY=0
XY=5
XY=10
XY=15
XY=20
Time Depth (H2W)- Site2, Line 3, S22-S19
15
20
25
30
35
40
45
0 20 40 60 80 100 120 140
Distance (m)
Tim
e (m
s)
XY=-10
XY=-5
XY=0
XY=5
XY=10
XY=15
XY=20
Time Depth(H2W) -Site2, Line 4, S22-S19
0
5
10
15
20
25
30
35
40
45
50
0 20 40 60 80 100 120 140
Distance (m)
Tim
e (m
s)
XY=-10XY=-5
XY=0XY=5
XY=10XY=15
XY=20
385
Appendix L (cont.)
Depth Section (H2W), XY=5 , Site2, Line 1, S22-S19
0
5
10
15
20
25
30
0 20 40 60 80 100 120 140
Distance (m)
Dep
th (m
)
XY=5average
Depth Section (H2W), XY=5 , Site2, Line 2, S22-S19
0
5
10
15
20
0 20 40 60 80 100 120 140
Distance (m)
Dep
th (m
)
XY=5
average
Depth Section (H2W), XY=5 ,Site2 , Line 3, S22-S19
0
5
10
15
20
0 20 40 60 80 100 120 140
Distance (m)
Dep
th (m
)
XY=5
average
386
Appendix L (cont.)
Depth Section (H2W), XY=5 ,Site2, Line 4, S22-S19
0
5
10
15
20
0 20 40 60 80 100 120 140
Distance (m)
Dep
th (m
)
XY=5
average
Time Depth (H1N)- Site2, Line 1, S23-S18
10
15
20
25
30
0 20 40 60 80 100 120 140
Distance (m)
Tim
e (m
s)
XY=-10
XY=-5
XY=0
XY=5
XY=10
XY=15
XY=20
Time Depth (H1N)- Site2, Line 2, S23-S18
8
13
18
23
28
0 20 40 60 80 100 120 140
Distance (m)
Tim
e (m
s)
XY=-10XY=-5
XY=0XY=5
XY=10XY=15
XY=20
387
Appendix L (cont.)
Time Depth (H1N) - Site2, Line 3, S23-S18
10
15
20
25
30
0 20 40 60 80 100 120 140
Distance (m)
Tim
e (m
s)
XY=-10XY=-5XY=0XY=5XY=10XY=15XY=20
Time Depth (H1N)- Site2, Line 4, S23-S18
10
15
20
25
30
0 20 40 60 80 100 120 140
Distance (m)
Tim
e (m
s)
XY=-10XY=0XY=5XY=10XY=15XY=20
Depth Section(H1N), XY=5 - Site2, Line 1, S23-S18
0
5
10
15
20
0 20 40 60 80 100 120 140
Distance (m)
Dep
th (m
) XY=5
average
388
Appendix L (cont.)
Depth Section (H1N) XY=5 - Site2, Line 2, S23-S18
0
5
10
15
20
0 20 40 60 80 100 120 140
Distance (m)
Dep
th (m
) XY=5
average
Depth Section (H1N), XY=5 - Site2, Line 3, S23-S18
0
5
10
15
0 20 40 60 80 100 120 140
Distance (m)
Dep
th (m
)
XY=5
average
Depth Section (H1N), XY=5 - Site2. Line 4, S23-S18
0
5
10
15
0 20 40 60 80 100 120 140
Distance (m)
Dep
th (m
)
XY=5
average
389
Appendix M: Wavepath Eikonal Traveltime Inversion
There are several methods of asymptotic tomography using the back projection
of traveltime residuals for the inversion of slowness fields in an acoustic
medium. These methods include ray-tracing (RT) tomography described by
Cerveny et al. (1977), the Woodward-Rocca (WR) method described by
Woodward and Rocca (1988), wave-equation traveltime (WT) inversion such as
described by Luo and Schuster (1991), and Schuster’s and Quintus-Bosz’s
(1993) wavepath eikonal traveltime (WET) inversion. It can be shown that
these methods are all derived from the general asymptotic back-projection
formula with different choices of the arbitrary weighting factor.
Wavepath eikonal traveltime inversion is so named because it computes
wavepaths using finite-difference solutions to the eikonal equation. Of the
asymptotic tomography methods, it is the most desirable in terms of
computational cost, robustness and effectiveness. It is generally better than the
WT and WR methods since it is an order of magnitude faster, yet is comparable
in effectiveness. It is also better than RT tomography since, whilst being
computationally comparable, it accounts for bandwidth and shadow effects in
the data by modelling multiple signal propagation paths contributing to one first-
arrival (which RT tomography does not).
The first step of the WET inversion algorithm is to propose an initial slowness
model, which may be generated assuming constant vertical velocity gradients
such as with 1D tau-p inversion, or assuming discrete velocity changes such as
390
with the generalised reciprocal method (Palmer, 2003). The eikonal equation
associated with this slowness field is then solved by a finite-difference method
such as that of Vidale, 1988 (discussed below). This results in a regular grid of
finite-difference first-arrival traveltimes within the medium which are subtracted
by the observed first-arrival traveltimes to give a grid of traveltime residuals.
These residuals are then used to evaluate the source weighting function of the
asymptotic back-projection formula at all points within the medium. The
slowness model is then updated. These steps are iteratively repeated until
convergence.
In order to use WET inversion, the eikonal equation associated with the
slowness field (1) must be solved for first-arrival traveltimes.
222
),( zxszT
xT
=⎟⎠⎞
⎜⎝⎛∂∂
+⎟⎠⎞
⎜⎝⎛∂∂
(1)
Vidale (1988) proposed an efficient method to do this by using finite-differences
to calculate the traveltimes on a regular grid. With this method a simple
equation, the traveltime estimator (2) is used to find the traveltime at a fourth
point on a grid cell knowing the traveltimes at the other three points (Figure 1,
left). By using an expanding square ring process, traveltimes are determined
along successive square rings centred at the source using those on the
previous ring until the traveltimes at all grid points inside the model are known
(Figure 1, right).
391
22 )()(2 ONMP tthstt −−±= (2)
s
O P
h
NM
Figure 1: Vidale’s eikonal equation solution (1988). The scheme for the
traveltime estimator, knowing the traveltimes at points O, M and N, the
traveltime at P can be estimated (left). The expanding square ring process,
traveltimes are determined along successive square rings centred at the
source point (right).
The major problem with Vidale’s method is that the expanding square ring
process is not entirely appropriate, being applied without considering the
surrounding slowness structure. Estimated traveltimes at a point may be
associated with supercritical incidence with respect to the slownesses. That is
causality is violated: “the time for the part of the raypath leading to a point must
be known before the time of the point can be found” (Vidale, 1988). As a
consequence the values inside the square root in the traveltime estimator
become negative resulting in an imaginary term in the estimator. This is readily
apparent when the slowness models have any considerable velocity contrasts.
392
The problem arising from solving the eikonal equation along an expanding
square ring can be largely eliminated by instead solving it along an expanding
wavefront (Qin et al., 1992). Such a strategy ensures that any point about to be
timed will have had its associated ray completely timed up to that point and
hence observing causality. This strategy gives accurate traveltimes in complex
velocity models and is stable in problem areas such as shadow zones and
waveguides. However as with the expanding square ring process this method
can still return negative square root arguments in extremely complex slowness
models.
Another way to overcome negative values inside the square root in the
traveltime estimator was proposed by Lecomte et al., 2000. Their method
achieves this by using five different traveltime estimators at each point and
taking the minimum estimated traveltime. Also, instead of using the expanding
square ring it uses an active point process where the central point is active and
may update the traveltimes at its eight surrounding points using different
estimators. The method can handle complex slowness models with strong
velocity contrasts. As well as this, its robustness allows for the recovery of
more than the first-arrival traveltime with later head waves able to be obtained
by masking fast deep layers with a high slowness mask.
Wavepath eikonal traveltime inversion is a high frequency traveltime
tomographic method. It is a computationally efficient method, being an order of
magnitude faster than wave-equation traveltime inversion (yet comparable in
effectiveness) since only solutions to the eikonal equation are involved. It
393
models multiple signal propagation paths contributing to one first-arrival and as
such is superior to conventional ray-tracing tomography which is limited to the
modelling of just one ray per first-arrival.
References:
Cerveny, V., Molotkov, I. A., Psencik, I., 1977, Ray Method in Seismology: Praha, Universita, Karlova.
Lecomte, I., Gjoystdal, H., Dahle, A., Pedersen, O.C., 2000, Improving modelling and inversion in refraction seismics with a first-order Eikonal solver: Geophysical Prospecting, 48, 437-454.
Palmer, D., 2003, Application of amplitudes in shallow seismic refraction inversion: 16th Conf., ASEG, extended abstracts.
Qin, F., Luo, Y., Olsen, K., Cai, W., Schuster, G., 1992, Finite-difference solution of the eikonal equation along expanding wavefronts: Geophysics, 57, 478-487.
Figure O12: Refractor Poisson’s ratio, site 2. P-wave/SV-wave (top), P-
wave/SH-wave (bottom).
407
0 20 40 60 80 100 120 140Cross-line
0
10
20
30
In-li
ne
20 40 60 80 100
0 20 40 60 80 100 120 140Cross-line
0
10
20
30
In-li
ne
20 40 60 80 100
0 20 40 60 80 100 120 140Cross-line
0
10
20
30
In-li
ne
20 40 60 80 100
Figure O13: Amplitude product, site 2. P-wave (top), SV-wave (middle),
SH-wave (bottom).
408
0 20 40 60 80 100 120 140Cross-line
0
10
20
30
In-li
ne
1 2 3 4
0 20 40 60 80 100 120Cross-line
1400
10
20
30
In-li
ne
1 2 3 4
0 20 40 60 80 100 120Cross-line
1400
10
20
30
In-li
ne
0.5
1.0
1.5
2.0
2.5
3.0
Figure O14: Amplitude product ratio, site 2. P-wave/SV-wave (top), P-
wave/SH-wave (middle), SV-wave/SH-wave (bottom).
409
Appendix P: GRM Tomography Images (Site1 & Site2)
GRM Velocity Images, Pwave, Site1
0 20 40 60 80 100 120 140
-15
-10
-5
0
400
1200
2000
2800
3600
4400
5200
0 20 40 60 80 100 120 140
-15
-10
-5
0
0 20 40 60 80 100 120 140
-15
-10
-5
0
0 20 40 60 80 100 120 140
-15
-10
-5
0
200
1100
2000
2900
3800
4700
200
110 0
2000
2900
3800
4700
300
1100
1900
2700
3500
4300
Line1
Line2
Line3
Line4
P1: GRM velocity-depth section obtained from WET inversion for P-wave at site 1.
410
GRM Velocity Images, h2N, Site1
0 20 40 60 80 100 120 140
-20
-15
-10
-5
0
0 20 40 60 80 100 120 140
-20
-15
-10
-5
0
0 20 40 60 80 100 120 140
-20
-15
-10
-5
0
0 20 40 60 80 100 120 140
-20
-15
-10
-5
0
100
800
1500
2200
2900
3600
0 800
1600
2400
3200
4000
0 800
1600
2400
0 800
1600
2400
Line1
Line2
Line3
Line4
P2: GRM velocity-depth section obtained from WET inversion for SH wave at site 1.
411
0 20 40 60 80 100 120 140
-20
-15
-10
-5
0
0 20 40 60 80 100 120 140-25
-20
-15
-10
-5
0
0 20 40 60 80 100 120 140-25
-20
-15
-10
-5
0
0 20 40 60 80 100 120 140-25
-20
-15
-10
-5
0
0 650
1300
1950
0 600
1200
1800
2400
3000
0 600
1200
1800
2400
0 800
1600
2400
3200
4000
4800
Line1
Line2
Line3
Line4
GRM Velocity Images, h1W, Site1
P3: GRM velocity-depth section obtained from WET inversion for SVwave at site 1.
412
GRM Velocity Images, Pwave, Site 2
0 20 40 60 80 100 120 140
-14
-12
-10
-8
-6
-4
-2
0
0 20 40 60 80 100 120 140
-15
-10
-5
0
0 20 40 60 80 100 120 140
-15
-10
-5
0
0 20 40 60 80 100 120 140
-15
-10
-5
0
400 900 1400
1900
2400
2900
3400
200
1100
2000
2900
3800
0 800
1600
2400
3200
4000
400
1000
1600
2200
2800
3400
4000
4600
Line1
Line2
Line3
Line4
P4: GRM velocity-depth section obtained from WET inversion for P-wave at site 2.
413
GRM Velocity Images, h2W, Site 2
0 20 40 60 80 100 120 140
-15
-10
-5
0
0 20 40 60 80 100 120 140
-15
-10
-5
0
0 20 40 60 80 100 120 140
-20
-15
-10
-5
0
0 20 40 60 80 100 120 140
-15
-10
-5
0
200
900
1600
2300
3000
3700
0 700
1400
2100
2800
3500
4200
4900
200 900 1600
2300
3000
200
800
1400
2000
Line1
Line2
Line3
Line4
P5: GRM velocity-depth section obtained from WET inversion for SH wave at site 2.
414
GRM Velocity Images, h1N, Site 2
0 20 40 60 80 100 120 140
-15
-10
-5
0
0 20 40 60 80 100 120 140
-14
-12
-10
-8
-6
-4
-2
0
0 20 40 60 80 100 120 140
-15
-10
-5
0
0 20 40 60 80 100 120 140
-14
-12
-10
-8
-6
-4
-2
0
0 700 1400
2100
0 700 1400
2100
0 700 1400
2100
2800
3500
0 700 1400
2100
2800
Line1
Line2
Line3
Line4
P6: GRM velocity-depth section obtained from WET inversion for SV wave at site 2.
415
APPENDIX Q: GRM Smoothing Statics for Shallow Refraction Data
Statics Corrections for Shallow Seismic Refraction Data.
by
Derecke Palmer1, Ramin Nikrouz2 and Andrew Spyrou3
Q1: Abstract
The determination of seismic velocities in refractors with near surface seismic
refraction investigations is an ill posed problem. Small variations in the
computed time parameters can result in quite large lateral variations in the
derived velocities, which are often artifacts of the inversion algorithms. Such
artifacts are usually not recognized or corrected with forward modeling.
Therefore if detailed refractor models are sought with model-based inversion,
then detailed starting models are required.
The usual source of artifacts in seismic velocities is irregular refractors. Under
most circumstances, the variable migration of the generalized reciprocal method
(GRM) is able to accommodate irregular interfaces and generate detailed starting
models of the refractor. However, where the very near surface of the earth is
also irregular, the efficacy of the GRM is reduced, and weathering corrections
can be necessary.
The standard methods for correcting for surface irregularities are usually not
practical where the very near surface irregularities are of limited lateral extent. In
such circumstances, the GRM smoothing statics method (SSM) is a simple and
1 School of BEES, UNSW, Sydney 2052, Australia. [email protected] 2 School of BEES, UNSW, Sydney 2052, Australia. [email protected] 3formerly School of BEES, UNSW, Sydney 2052, Australia, now Fugro Ground Geophysics Pty Ltd. [email protected]