3D Gravitational Inversion Modelling of the South Dunedin Sub Basin, Otago, New Zealand
Total Page:16
File Type:pdf, Size:1020Kb
3D Gravitational Inversion Modelling of the South Dunedin Sub Basin, Otago, New Zealand Timothy John Lutter A thesis submitted for the degree of Masters in Geophysics at the University of Otago Dunedin New Zealand January 2018 Abstract 3D gravitational modelling of an infilled coastal sedimentary basin is rare, with almost all 3D gravity publications being based on mineral deposits related to hard rock mineral density features. The sediment-bedrock boundary under the South Dunedin tombolo, New Zealand provides a target to test the application of 3D inversion modelling on a small basin system. The infilled South Dunedin basin has never been fully characterised and is the perfect target for 3D inversion modelling due to a high contrast in density disparity between Holocene sedimentary materials and bedrock. Surveying was conducted using the University of Otago’s Worden gravimeter. Data were collected over 20 days covering about 6.4 km2 with the collection of 320 survey points. This survey was the first gravitational field work aimed at collecting data for 3D inversions by anyone from the Otago University’s Geophysical Team. The inversion software GRAV3D developed by the University of British Columbia-Geophysical Inversion Facility (UBC-GIF) was used for the 3D modelling. 2D modelling confirmed the bedrock sedimentary features seen in the 3D models. With the 3D model producing an insight into the bedrock interface that highlighted the potential of these features include paleo-valley and a potential offset of geological units usually seen as a result of faulting. The model shows that the paleo-channel bedrock interface was preferentially eroding down, rather than meandering across and paleo-surface. Resulting 3D models show the first view of the paleo-valley's sedimentary to bedrock interface clearly, with three main features unearthed for the first time: the ridge running down the centre of the valley, the symmetry of the valley sides and the identification of density variation in the Holocene sediments across the South Dunedin infill section. This small-scale variance, not seen in other models of the paleo-valley highlights the fact that more still can be learnt on what is beneath one of New Zealand’s main cities. Due to early errors within the early processing of data, that were caught only at the marking stage of the thesis, some features discussed above would have been affected and re-modelling from that stage onward is needed. In saying that all models that have been produced are still believed to have the same broad shape as discussed and will not have been affected much by the modelling process. Due to working commitments I have foregone the reprocessing of the data and will be satisfied that the first ever 3D gravitational modelling at Otago University was produced to the highest of standards. Further study is needed to determine the potential applications of 3D gravity modelling in sedimentary basins. With reprocessing of the models plus the confirmation of features seen in the models produced within this thesis could be completed by high-resolution seismic survey or drill holes into the South Dunedin Basin. Contents 1. Introduction ...................................................................................................................... 13 1.1 Gravity Surveying ...................................................................................................... 14 1.2 Gravitational Inversions ............................................................................................. 15 1.3 Global Importance ..................................................................................................... 15 1.4 South Dunedin Shallow Basin ................................................................................... 15 1.5 Regional Geology ...................................................................................................... 16 1.5.1 Otago Schist ....................................................................................................... 18 1.5.2 Sediments ........................................................................................................... 19 1.5.3 Dunedin Volcanic Group Synthesis................................................................... 23 1.5.3.1 Initial eruptive phase ................................................................................ 23 1.5.3.2 First eruptive phase .................................................................................. 24 1.5.3.3 Second eruptive phase ............................................................................. 24 1.5.3.4 Third eruptive phase ................................................................................ 24 1.5.3.5 Port Chalmers Breccia ............................................................................. 24 1.5.4 Development of Otago Harbour ............................................................................. 25 1.5.5 South Dunedin tombolo infill ................................................................................. 26 1.6 Previous Geophysical Work ...................................................................................... 26 1.7 Aims ........................................................................................................................... 29 1.8 Thesis Outline ............................................................................................................ 30 2. Data Collection ................................................................................................................ 31 2.1 The Worden Gravimeter ............................................................................................ 32 2.2 Survey Design ............................................................................................................ 34 2.3 Field Acquisition ........................................................................................................ 34 2.4 Data Reduction Process ............................................................................................. 36 2.4.2 Correction for Temporal Drift ........................................................................... 36 2.4.3 Meter Calibration ............................................................................................... 37 2.4.4 Absolute Gravity ................................................................................................ 37 2.4.5 Latitude Correction ............................................................................................ 37 2.4.6 Free Air Correction ............................................................................................ 38 2.4.7 Terrain Correction ............................................................................................. 38 2.4.8 Bouguer Correction ........................................................................................... 39 2.4.9 Resolving Regional Residual Anomalies .......................................................... 40 3. Data Reduction ................................................................................................................ 42 3.1 Collected Data ............................................................................................................ 43 3.1.1 Collected Gravity Data ...................................................................................... 43 3.2 Latitude Correction .................................................................................................... 44 3.3 Free Air Reduction ..................................................................................................... 45 3.4 Terrain/Hammer Cell Correction ............................................................................... 46 3.5 Bouguer Correction .................................................................................................... 47 3.6 Regional Background Gravitational Reduction ......................................................... 48 3.7 Anomalous Gravity Data ........................................................................................... 49 3.8 South Dunedin Gravitational Reduction Results ....................................................... 50 3.9 South Dunedin Gravity Anomalies ............................................................................ 51 4. 2D Forward Modelling ........................................................................................................ 56 4.1 GRAVCADW ............................................................................................................ 58 4.1.1 GravCadW Input Files ....................................................................................... 58 4.1.2 GravCadW Operating Specifications ................................................................ 58 1 4.2 Hillside Road Line A to A ........................................................................................ 59 4.3 Bayview Road Line B ................................................................................................ 60 4.4 Victoria Road Line C ................................................................................................. 62 1 4.5 Forbury Road Line D to D ........................................................................................ 63 4.6 Prince Albert Road Line E ......................................................................................... 65 4.7 Portsmouth Road Line F ...........................................................................................