A Seismic Investigation into Structural Controls on ...

5 downloads 0 Views 5MB Size Report
calculated from available well log checkshots provided by Arco Petroleum in 1992, International. 157. Department Dallas in 1984 and NZ Oil and Gas in 1984, ...
1

A Seismic Investigation into Structural Controls on Submarine

2

Channel Distribution in the Taranaki Basin, New Zealand

3

Aisling Scully1

4

1

5

ABSTRACT

School of Biological, Earth and Environmental Sciences, University College Cork, Cork, Ireland

6

Submarine or abyssal channels are formed primarily from turbidity currents in continental

7

shelf margins. Still poorly understood, these channel systems can span up to hundreds of

8

kilometers and are economically appealing, being one of the main sediment transport mechanisms

9

in deep water settings. Using high quality 3-D seismic reflection data, a series of these channels

10

were mapped in the Parihaka Basin, offshore New Zealand. The Parihaka Basin has an aerial extent

11

of 1520 km2, with a multiphase deformational history and is transected by the Parihaka Fault. This

12

complex fault is approx. 50 km in length and has undergone several phases of activation leading

13

to substantial throw displacements. Seismic coherence was used to record channel thalweg

14

positions, at depth intervals of 25 ms (two-way travel time). The main parameters recorded were

15

channel point location, sinuosity, vertical stacking and channel density. Supporting data collected

16

included sediment thickness maps and fault throw measurements. In typical marine half-graben

17

structures, the areas of maximum sediment supply and drainage are positioned mainly by the

18

footwall scarp or in transfer/relay zones. While the data collected supported this, the results also

19

showed concentration of channel points at sites of rupture associated with oblique deformation,

20

particularly in the opposing sides of the main breached transfer zone. This was repeated numerous

21

times in similar locations throughout the study area. Overall the results displayed a predictable

22

pattern of sedimentation, which was directly controlled by this type of en echelon fault system,

23

with concentration of sediment buildup and channel distribution at localized areas of subsidence

24

and uplift, particularly in the zones of linkage between segments.

25

INTRODUCTION

26

One of the most important zones being targeted for hydrocarbon exploration at present is

27

the deep-water continental slope margin. It is clear that a general understanding of the dynamics

28

in these offshore areas has been greatly enhanced by 3D seismic investigation (Posamentier,

29

2004). Increased detail in current studies shows far more complexity in these systems than

30

previously thought and new methodologies have been constructed in order to classify factors such

31

as morphology and depositional patterns in these submarine environments (McHargue et al., 2001;

32

Deptuck et al., 2007).

33

Sediment transport in continental shelf margins is largely provided by turbidity currents

34

that migrate suspended material downslope (Bouma, 2004). The individual features of each

35

deposit are constrained by physical parameters such as the steepness of the slope, the composition

36

of the sediment and confining physical barriers within the basin topography (Posamentier and

37

Walker, 2006). These transport flows lead to the development of a variety of characteristic features

38

typical of such a setting, with submarine or abyssal channels and canyons covering vast distances.

39

These channel systems can be either straight or highly sinuous depending on the gradient of the

40

area, with higher sinuosity occurring in flatter terrain (Deptuck et al., 2003). Insight into the

41

evolution of these abyssal channel systems, which may contain economically viable sand rich

42

reservoirs, is of critical importance and can be reliant on a variety of controls. Studies have shown

43

that tectonic activity can have a dominant influence on the patterns and rates of deposition and

44

supply of sediment, particularly in the instance of fault related subsidence (Bouma, 2004; Lin et

45

al., 2004; Shepard and Emery, 1973).

46

Finding a consistent and repeatable array of sediment distribution in relation to specific

47

areas on a particular fault system is an important topic for investigation. This is due to the

48

prediction of potential hydrocarbon reservoirs being the main goal of exploration companies, with

49

focus being placed on accurately identifying the optimal location to drill. Intrabasin transfer zones

50

such as the one seen in the study area are often the site of hydrocarbon fields and classification of

51

them has important implications for predicting reservoir location and size (Gawthorpe and Hurst,

52

1993). This study aims to use high quality seismic reflection data and associated attributes in order

53

to map the sediment transport surrounding the oblique normal system seen in the Parihaka Fault

54

and determine the level of influence the fault system has on drainage and topography locally. It 55

also aims to ascertain any predictable

56

localities of sediment build-up that

57

relate to the en echelon structure by

58

constraining the various segments

59

connected to the transfer relay zones.

60

GEOLOGICAL BACKGROUND

61

The study area, which is known

62

as the Taranaki Basin, is situated

63

approximately 25 km offshore from the

64

north-western

65

Plymouth, New Zealand (Figure 1),

66

encompassing ~1520 km2 of the North

67 Figure 1. Overview of the study area, highlighted in the black 68 polygon. Red and green polygons indicate previous

Taranaki Graben (Veritas, 2005). This

petroleum exploration zones. Modified from (Veritas, 2005).

boundary

of

New

entire region is divided into three

69

separate sections, the Western Stable Platform, the Southern Inversion Zone and the Eastern

70

Mobile Belt, where the North Taranaki Graben is situated. At present there is extension occurring

71

in the Northern and Eastern sections, with compression happening in the South (King and

72

Thrasher, 1996; Giba et al., 2010). The entire North Graben is divided down the centre by an

73

extensive volcanic complex called the Mohakatino Volcanic Centre, comprising a series of

74

submarine stratovolcanoes and intrusive bodies that permeate the region (King and Thrasher,

75

1996). The Basin has changed over its duration from a post-rift basin during the Paleogene to an

76

active foreland-basin in the Neogene. It is characterised largely by the main prograding sequence

77

present, which is called the Giant Foresets Formation, a sedimentary section up to several

78

kilometres in thickness (Stern and Davey, 1989).

79

The Taranaki Basin has undergone three main deformational phases, the first being a period

80

of extension during the Late Cretaceous (approx. 84- 50 ma), driven mainly by the dissolution of

81

Gondwana (Giba et al., 2010). The second phase was a period of compression between the Late

82

Eocene up recent times (approx. 40-0 ma) and second extensional period spanning the Late

83

Miocene to Recent (approx. 12-0 ma), largely initiated by subductive processes at the Hikurangi

84

Trough and Pacific plate interaction zone, which caused back arc rifting to occur (Stern and Davey,

85

1989). There was an overlap of compression and rifting between the Late Miocene to recent,

86

attributed to subduction driven differences in tectonic behaviour on a regional scale (Giba et al.,

87

2010).

88

THE PARIHAKA FAULT

89

The Parihaka Fault is the most prominent feature within the seismic survey, dividing the

90

area from the centre up to the northern section (see Figure 2). The fault itself is a complex structure

91

Figure 2. The sequence of images shown the left display the three main stages of activity in the Parihaka Basin, with Cretaceous extension show in image 1 and the Early Pliocene oblique reactivation shown in image 2. Current extension is shown in image 3. The sequence to the right displays the stages of fault linkage from hard at depth up to soft linkage in shallower horizons. Modified after (Giba et al., 2012). 92

93

comprised of four individual fault segments, ranging between 8 and 15 km in length, soft linked

94

in shallower sections by varying degrees of intensity and hard linked at depth. The fault has

95

undergone several stages of activation and displacement and its unique en-echelon pattern is

96

attributed to re activation of the fault during the Early Pliocene in a different direction to the

97

original movement (Giba et al., 2012). It is the primary structural system that borders the western

98

edge of the Taranaki Basin and has a total length of approx. 50 km. The fault trace displays a range

99

of strikes, dependent on which horizon is viewed and its activation history. It is near-continuous

100

at the Late Cretaceous level, with throws of approx. 1800 m and striking NNW to SSE (Giba et

101

al., 2012). Changing at the Early Pliocene level, the four individual fault segments are noticeably

102

separate and the fault strike in this location is approx. N-S to NE-SW, with vertical displacements

103

of up to 1450 m (Giba et al., 2012). After the Miocene, both breached and unbreached relays divide

104

the fault segments due to oblique reactivation of the original fault. These are hard linked to the

105

basement fault, which strikes in an N-S direction (Giba et al., 2012).

106

Giba et al. (2012) suggested that the relay zones of the Parihaka Fault support a coherent

107

fault model, where the rapid formation of both fault segments and relay ramps occur geologically

108

instantaneously, usually originating from a connection to a continuous basement fault at depth

109

(Walsh et al., 2003). The Parihaka Fault reflects this, where it displays soft linkage at shallower

110

horizons in the en echelon shaped sections, which combine at depth to the preceding unbroken

111

fault trace at (Figure 2). The creation of these distinct fault segments was caused by reactivation

112

around 3.7 Ma, at a 45-degree angle to the original extension during the Cretaceous. Although

113

fault growth was instantaneous, further growth was accommodated by later displacement and

114

rotation of the relay zones (Giba et al., 2012). The sedimentary record during the Late Cenozoic

115

was exceptionally preserved due to high sedimentation rates, leading to the blanketing and

116

conservation of the growth history of the main faults. When displacement analysis was undertaken

117

on this fault, it was seen that movement post- Miocene was concentrated in the segments located

118

at the north of the survey. The NE- SW fault segments underwent the most extension as they were

119

the most optimally aligned to the direction of movement (Giba et al., 2012). Similar orientation in

120

the fault strikes between the Cretaceous and Pliocene fault segments has produced an overall 3D

121

fault shape that varies between sigmoidal and straight, with straight fault sections occurring in

122

areas with particularly close alignment (Giba et al., 2012). Due to the sigmoidal shape that occurs

123

in certain areas of the fault, the need for accommodation of volumetric strain has led to the

124

formation of some smaller antithetic faults adjacent to a couple of the fault segments located in the

125

hanging wall (Giba et al., 2012).

126

METHODS

127

3D Seismics

128

3D seismic reflection consists of firing compressed air shots which generate sound waves

129

that penetrate the subsurface. This seismic energy is produced in the form of wavelets, which

130

reflect off different materials at a specific time interval. This reflected acoustic data is then

131

processed into a three dimensional volume, which can be used to create images in both horizontal

132

and vertical directions and then utilised for geological interpretation (Kearey et al., 2002).

133

The Parihaka seismic survey was shot in the Taranaki Basin using Veritas Viking II

134

between January 12th and February 24th, 2004 (Veritas, 2005). A total of eight streamers, each

135

4500 meters in length and separated by 100 meters, were used in combination with alternating dual

136

sources, laterally spaced at 50 meters. This combination resulted in a 60-fold subsurface

137

illumination area of 400 meters by 2250 meters for each shot. The source arrays were fired

138

alternately at 25.0 meter intervals, during which time 6.1440 seconds of seismic reflection data

139

Figure 3. Image A displays a representative cross section through the inline of the survey area, complete with dated horizons and fault interpretation lines. Image B shows the coherence cube created of the study area, displaying the survey boundary and the four main fault segments that cut across it.

140

were recorded at 2 milliseconds sampling rate. This data was then processed and filtered to create

141

a seismic volume.

142

The seismic attribute generated from this volume was a coherency cube. This was produced

143

using Petrel 2014 by Schlumberger. Coherency is the measurement of similar adjacent traces

144

which are located within a specific area by changing the volume from continuity to discontinuity.

145

This makes it useful for interpretation of any disruptions within the volume which may represent

146

features such as faults, channel points, sedimentary structures or unconformities. It can also be

147

used in order to observe stratigraphic changes such as parallel or chaotic bedding and mapping the

148

lateral extent of a certain reflector.

149

Data Analysis

150

All interpretation, calculations of area and modelling were carried out using Petrel 2014

151

Software package by Schlumberger. The fault interpretation tool was used to map each individual

152

fault segment and two main horizons were mapped using the picking tool, selecting the strongest

153

reflectors above and below the section where channels were located. These two main reflectors

154

will be referred to as the Top Chanel Horizon and the Bottom Channel Horizon. The seafloor and

155

basement horizons were also mapped in order to accurately date the various layers and features.

156

Data from previously completed depth conversion graphs (Giba et al., 2012), which were

157

calculated from available well log checkshots provided by Arco Petroleum in 1992, International

158

Department Dallas in 1984 and NZ Oil and Gas in 1984, provided a time to depth reference for

159

the seismic lines (Fig 3A). The biostratigraphic dating of each horizon, which was completed by

160

Giba in 2012, was utilised and compared to the corresponding TWTT lines in every seismic section

161

produced, in order to date the different layers (Figure 3A). A seismic coherence cube was created

162

(Figure 3B) and the polygon feature was used in order to create channel points. This was carried

163

out by analysing a plan view of the coherence image and working at regular intervals of 25 ms

164

between the depths of -864 ms and 1248 ms, the first and last appearance of channels in the images.

165

The channels were traced by placing points at a fixed distance along each channel thalweg. The

166

thalweg is the central, deepest section of the submarine channel. The point data, coordinates and

167

horizon grids were converted into a shapefile and exported into ARC GIS for analysis.

168

ARC Map was used for further analysis of channel morphology. The channel length was

169

measured by tracing each individual channel along its thalweg and measuring the entire distance

170

(Appendix II). Valley length was also recorded, a parameter defined as the straight line distance

171

between final points at each side of the channel (Flood and Damuth, 1987). These two datasets

172

(Appendix II) were recorded in order to determine the sinuosity of each channel, which was

173

calculated by dividing channel length by valley length. Meander wavelength was measured by

174

counting the number of meanders in each channel, then calculating the ratio of number of meanders

175

to valley length.

176

The supporting maps that were attached to the seismic survey were modified using ARC

177

Map. A channel point density map was also created using this software by exporting the shapefile

178

from Petrel and calculating the seismic point data using GIS.

179

RESULTS

180

Overview

181

The survey field is separated into four zones for analysis, displayed in Figure 4B. Zone A

182

is the topmost section of the fault, encompassing both hanging wall and footwall and highlights

183

the area around the main unbreached relay ramp. The seismic section S-A (Appendix I) is a

184

representative slice through this feature. Zone B1 and Zone B2 span the central section of the

185

survey, with B1 displaying the footwall features and B2 the hanging wall. The seismic slice S-B

186

(Appendix I) is a section through the fault which transects this part. Zone C is the bottom sector

187

of the survey, with the seismic slice S-C (Appendix I) displaying features throughout the fault in

188

this location.

189

Channel point distribution on the footwall of the Parihaka Fault is displayed in Figure 4A,

190

at depths of between -864 ms to -1248 ms, measured in Two Way Travel Time (TWTT). The

191

submarine channels initiate at number of locations along the main fault boundary and fan out in a

192

general N-W trend. Channel initiation occurs at both the inside and outside bend of the sigmoidal

193

S shape of the fault trace. There is also some channel point data seen within the fault link zones

194

and horsetail splay features. The channel systems on this side of the fault are widely spread and

195

display moderate sinuosity, spanning distances of up to 200 kilometres. Directly east of the

196

Parihaka Fault, channel point data on the hanging wall (Figure 4C) can be found between depths

197

of -1248 to -1600 ms in TWTT. Submarine channel initiation in this sector of the survey area

198

appears much more concentrated to the central section of the main fault, with three distinctive

199

localities heavily populated by points. There is no visible trend of channel direction in the hanging

200

wall as the data points are so tightly packed. However, the overall profile appears to be a

201

hemispherical pattern, spreading from the centre of the main fault.

202

Zone A

203

This section is located to the Northern half of the survey area and the main focus is the

204

relay ramp feature. As seen in the seismic section S-A (Appendix I), which is vertical profile

205

through the relay itself, the channel points can be observed on both sides of the central ramp. There

206

is an element of drag displayed in the beds that feed into the ramp and the channel points in these

207

horizons follow this trend accordingly. The channels exit on the hanging wall in the same horizon

208

as the footwall, however, the throw appears to be moderate due to rotation of the beds within the

209

ramp, also highlighted in the seismic profile.

210

In the coherence image displayed in Figure 4A, the channel points plotted on the footwall

211

show a general west- southwest trend, and are highly sinuous. The channels in this sector can be

212

seen feeding into the relay ramp, however there is also a channel initiating from the fault block

213

behind the main fault line which appears to bisect intrusive bodies. Point data on the hanging wall

214

propagates from one distinct area, where the channels cross the relay ramp and follow the same

215

general trend as the foot wall. Again the channels are sinuous and concentrated from the central

216

point of the C shape of the fault block.

217

Zone B1 and Zone B2

218

Zones B1 and B2 cover the footwall and hanging wall sectors of the centre of the survey

219

area. The seismic slice S-B in Appendix I is a vertical profile through both Zone B1 and Zone B2

220

across the main fault. In the seismic profile S-B (Appendix I) that intersects both Zones B1 and

221

B2, there is a clear difference between sediment thicknesses on both sides. A number of layers on

222

the footwall pinch out as they near the fault surface, whereas on the hanging wall there are many

223

extra layers of deposition between the top and bottom horizon. There is limited vertical stacking

224

of channels on the footwall side of approx. 200 ms thickness and they are situated directly adjacent

225

to the fault. On the hanging wall, vertical channel accretion is double that of the footwall with

226

vertical thickness of 700 ms. The throw in this central fault segment is 250 ms in difference

227

between footwall and hanging wall, particularly on the bottom horizon.

228

In Zone B1, which is the footwall portion of the central area, channels initiate in a number

229

of distinct locations. The fault in this section comprises a number of fault segments, connected by

230 Figure 4. Submarine channel point data is shown displaying extents and trends of channels in the survey area. A) Channel distribution located on the footwall. B) Overview of the area, showing Zones A, B1, B2 and C. Also displays the channel trends compared to individual fault strikes and the central hanging wall depocentre. C) Channel distribution on the hanging wall. 231

232

relay ramps, observed in plan view and in the seismic profile. In Figure 4A, there are three separate

233

fault segments each joined by a breached relay ramp with channel initiation occurring at each

234

linked zone. The channels in this area are moderately sinuous and spread out, with large visible

235

branching. The trend of channel direction varies depending on the fault trace strike. In Segment 1,

236

which has an N-S strike, the channels are running at a 45-degree angle to the south-west. However,

237

in Segments 2 and 3, which both have an N/NE-S/SW strike the channels also travel at a 45-degree

238

angle from the main fault strike, trending directly west. Zone B2, or the hanging wall section of

239

this area (Figure 4C), has a much more concentrated distribution of channels. Where the fault

240

segment traces display a sigmoidal profile, channel point initiation occurs mainly on the curvature

241

maximum of the “S” shape. This is repeated in three separate locations in this section. The channel

242

points in this location are highly sinuous and tightly packed near the area of initiation. There is

243

also a large quantity of channels in this area, as seen both in the plan view and in the seismic

244

profile.

245

Zone C

246

In the seismic profile for Zone C, S-C (Appendix I), throw between footwall and hanging

247

wall is minimal. The fault profile itself is much more steeply dipping in this section compared to

248

the other zones. Channel distribution on the footwall side is small, only visible directly adjacent to

249

the fault surface. Channel distribution on the hanging wall side is more prominent with repeated

250

vertical stacking occurring in this part of thicknesses up to 500 ms.

251

In Zone C, there is one main fault segment which displays the sigmoidal profile in plan

252

view. Channel point initiation on the footwall occurs mainly on the inside bend of this outline,

253

whereas on the hanging wall it occurs on the opposite bend of this sigmoidal fault trace. The

254

channels are concentrated largely on the hanging wall side and follow an E-W trend. There is also

255

some channel initiation occurring in the horsetail splay features at the very bottom of the main

256

fault. Channel initiation that is visible in the central area of the hanging wall is occurring adjacent

257

to the intrusive body which is situated to the south east of the coherence image in Figure 4A.

258

The channel point density map created for the entire area displays a number of features

259

(Figure 5B). Overall trends suggest three main lines of heavy population which transect the survey

260

area. Each is adjacent to the areas where the fault segments are connected, either by soft linked

261

relay ramps or where they have fully been breached and classified as transfer zones. The red and

262

yellow zones in the channel density map highlight the areas of greatest point concentration. This

263

is also directly related to vertical stacking as the channel density calculation takes into account

264

every point across all depth layers. There is a significant concentration of channels in the central

265

area of the hanging wall, this is also the location with the most vertical stacking, although some

266

can be seen on the foot wall as well.

267

An isochron sediment thickness map was also created (Figure 5A), with the top channel

268

horizon and bottom channel horizon being used as parameters. The constraining horizons selected

269

above and below the submarine channel appearance in the seismic imagery, are dated at 2.0 and

270

3.4 Ma respectively, during the Pliocene-Pleistocene syn-rift sequence and were used in order to

271

calculate sediment thickness. When compared with the channel point density map, the isochron

272

displayed a similar trend, with thickness maximums focused in the hanging wall across all three

273

zones (Figure 5B).

274

Figure 5. A) Ischron sediment thickness map produced of the survey area, calculated between the Top Channel Horizon and the Bottom Channel Horizon. B) Channel point density map created displaying concentration of channel points and vertical stacking. 275

276

Channel Sinuosity

277

When longitudinal profiling was carried out on a representative selection of channels, in

278

both the hanging wall and footwall, a number of trends were observed (Figure 6A). In the footwall,

279

sinuosity varied depending on proximity to the fault. In Zones B1 and C, sinuosity was classified

280

as straight (1.5) at the furthest distance

282

away from the fault. The average sinuosity in zone B1 was 1.18 next to the fault zone and 1.9 in

283

the section to the west of survey area. Averages in Zone C showed a similar gradient, moving from

284

1.2 near the fault and 1.8 away from it. This varied slightly in Zone A, which displayed some

285

moderate sinuosity in channels located immediately next to the fault and transfer zone, measured

286

at 1.46 and 1.8, fault adjacent and non-fault adjacent respectively. The sinuosity variation between

287

different depths in the footwall was negligible.

288

In the hanging wall, however, these trends were reversed (Figure 6B). Instead, sinuosity

289

was lower away from the scarp area, becoming less sinuous as proximity to the fault lessened. The

290

most notable characteristic of the parameter found in the hanging wall section was the variance

291

with depth. Channels changed morphology completely, moving from straight to sinuous to

292

meandering as depth decreased (Figure 7). There was little variation within zones in the hanging

293

wall, however, most channels were concentrated around the central B1 Zone on this half of the

294

study area.

295 296

Overall, the trends varied completely between foot wall and hanging wall sections, displaying opposing sinuosity measurements and behaviour.

297 Figure 6. Sinuosity profiles for the hanging wall and footwall channels. A) Overview showing different zones of sinuosity. B) Hanging wall distribution over depth changes. C) Footwall distribution based on zone and proximity to the fault system.

298 299 300

Figure 7. Changes in hanging wall sinuosity throughout a variety of depths. A) Sinuosity at

301

deepest level, 1584 ms, average of 1.10. B) Sinuosity at mid-level, 1464 ms, average of 4.50. C)

302

Sinuosity at shallowest level, -1320 ms, average of 10.85.

303

304

DISCUSSION

305

Central Hanging Wall Depocentre

306

Standard normal fault growth models suggest that displacement is highest in the centre of

307

a fault and decreases gradually towards the tips (Walsh and Watterson, 1987, 1988). It has also

308

been observed that values for throw or rates of localized depression and uplift may change along

309

the length of a fault (Gawthorpe et al., 1994). While that concept of fault growth is quite basic,

310

this study displays a strong correlation to that idealized type of model, particularly on the hanging

311

wall side. Offset between horizons varies from zone to zone, with the largest throw recorded in the

312

central Zone B1 at 250 ms, while measurements taken in the north section Zone A and the south

313

section Zone C are minimal in comparison, although in Zone A it can be attributed to intense

314

hanging wall deformation. An isochron sediment thickness map produced between the top and

315

bottom channel horizons (Figure 5A) displays deposition maximums or red zones in the central

316

section of the main fault, with thickness gradually decreasing along the length of the fault towards

317

the ends, indicating the accommodation space available is controlled by fault activity. This pattern

318

is consistent with displacement analysis undertaken on this fault system by Giba et al., in 2012,

319

which indicates throw maximums in the centre of the fault, decreasing at greater distance away

320

from the centre. This has implications correlating throw displacements to sediment thickness, as

321

both follow the same trend when moving from fault centre to tips. A 3D model created of the

322

bottom channel horizon, complete with structural contours (Figure 8A), clearly shows a large

323

depression in the centre of the hanging wall and provides a visual approximation of the topography

324

at this level. The extent of this basin and volume of sediment accumulation observed are

325

constrained by the fault type present and its related activity, as is typical of marine half-graben

326

systems (Leeder and Gawthorpe, 1987).

327

While significant subsidence occurred during the Cretaceous–Eocene extension period,

328

attributing to the total displacement figure of 1800 m (Giba et al., 2012), this study has focused on

329

the syn-rift deposition occurring during the Pliocene reactivation period. It is during this time

330

frame that the submarine channel system was formed, providing insight into the tectono-

331

sedimentary facies type present adjacent to this fault. This study shows a significant concentration

332

of channels in the central hanging wall area, with two localized zones of focus near the two main

333

transfer zone/relay ramps as seen in the channel density map produced. This has implications for

334

the type of drainage occurring in this system, as it supports a model put forward by Leeder and

335

Gawthorpe in 1987 of a continental half-graben with transverse through drainage. Several key

336

indicators support this type of facies model, firstly, the two main sites of channel concentration

337

surrounding the transfer zone and relay ramp, suggesting drainage initiating between the en

338

echelon segments and moving down into the sub-basin, while the presence of transverse migrating

339

channels, which produced both vertical and lateral channel stacking throughout the sequence, was

340

observed repeatedly throughout the seismic images. Sinuosity profiles completed on the hanging

341

wall channel systems at a range of depths support both the normal fault model and the facies model

342

proposed and provide evidence to help constrain slope and direction of flow in the hanging wall.

343

Channel sinuosity profiles were classified as straight at the base of the sequence, flowing down off

344

the footwall scarp and into the basin, while alternative drainage was recorded from channels which

345

emerged from the east of the survey area and ran into the central depression zone adjacent to the

346

fault. This provides two key details with regards the morphology of the sub-basin; that the slope

347

was relatively steep coming down off the fault scarp, as sinuosity decreases with increased slope.

348

349 350 351

Figure 8. Layer images of the Bottom Channel Horizon and Top Channel Horizon. A) Bottom

352

Channel Horizon illustrating the hanging wall depocentre. B) Top Channel Horizon showing the

353

en echelon morphology on footwall and structural dip for the area.

354

It also represents an element of initial steepness in the hanging wall dip slope. This would support

355

the coherent fault model suggested by Giba et al., in 2012, as it subscribes to the idea of

356

instantaneous fault growth, with maximum displacement occurring over a brief period of time, and

357

a gradual levelling out of the hanging wall slope, a concept which might justify such steep hanging

358

wall dip followed by rapid infilling. That the channels increase in sinuosity as the sequence moves

359

up depth, is consistent with the evolution of these types of submarine channel systems, changing

360

from immature, linear channels in high gradients to mature, meandering channel systems with a

361

complex geometry (Gee and Gawthorpe, 2007) and displaying vertical accretion as seen in this

362

study. This is due to slope levelling out as the sequence thickens, a concept which is familiar in

363

the infilling of half-graben hanging wall structures (Gawthorpe and Leeder, 2000).

364

Transfer Zone/ Relay Ramp Footwall Controls

365

Displacement accommodation between individual fault segments is often aided by the

366

formation of intrabasin transfer zones or relay ramps, particularly where the fault shape follows an

367

en echelon side stepping pattern (Leeder and Gawthorpe, 1987). These intrabasin transfer zones

368

have a strong influence on sedimentation trends on a localised scale and may act as a pathway for

369

turbidity currents or footwall deposits, especially in transverse drainage facies models (Gawthorpe

370

and Hurst, 1993). Individual transfer zone morphologies are poorly constrained at the subsurface

371

level (Gawthorpe and Hurst, 1993), as the majority of seismic surveys are imaged parallel to the

372

strike of these areas (Gibbs, 1989).

373

In this study area, there are two separate transfer zones within the main fault system. The

374

first transfer zone is located in Zone A, and can be classified as an unbreached relay ramp. The

375

submarine channel data which was recorded surrounding the ramp, shows a clear convergence of

376

footwall drainage moving down through the relay ramp and feeding into the hanging wall sub-

377

basin. The second transfer zone, located in Zone B, also acts as a conduit for sediment flow into

378

the hanging wall and is fully breached. Both of these factors are consistent with the transverse

379

drainage facies model (Leeder and Gawthorpe, 1987) suggested in this study already from hanging

380

wall observations. This theory is backed up by the results of the channel point density map which

381

indicates a build-up of sediment in two distinct lines of concentration focused around each transfer

382

zone from across the footwall.

383

There are a number of notable trends found within this footwall channel data which may

384

have implications for the en echelon segment morphology. The first element for consideration, is

385

how the flow direction of an individual channel follows a repeatable pattern depending on the

386

particular point of initiation along the sigmoidal “S” fault trace. This infers separate zones of

387

subsidence and uplift which are directly controlled by the fault profile. In the lower half of each

388

en echelon segment, channels initiate in the central internal point of curvature, then run off at

389

approx. 45° to the fault strike, this angle does not change throughout the entire fault system, and

390

when the fault strikes vary from south to north, the direction of channel movement maintains the

391

same angle and varies accordingly. As fluid flow of any type migrates downslope taking the path

392

of maximum gradient, this behaviour in has implications for topographical high zones within the

393

segments of this fault type, suggesting an area of uplift is situated in the lower half of the en echelon

394

unit. Previous research indicates that certain structures observed in these typical relay ramp

395

transfer zones hold similarities to strike slip fault systems, particularly the alternating zones of

396

extension and compression (Reading, 1980; Gawthorpe and Hurst, 1993), an idea that appears to

397

be supported by the behaviour of the fault system examined in this study. In the upper section of

398

the same fault, where channels initiate at the point of maximum external curvature, the angle is

399

observed to be between 85°-90° from the fault strike and channels flow in the opposite direction

400

to the lower half. This implies that the topography in this half of the en echelon segment slopes

401

down towards the relay transfer zone. This result corresponds with the characteristic elevation

402

decrease towards the ramp area associated with these type of intrabasin transfer zones (Gawthorpe

403

and Hurst, 1993).

404

Sinuosity profiles completed on the channels in the footwall section show some

405

corresponding trends. General patterns in sinuosity show low to zero sinuosity directly adjacent to

406

the fault segments and moderate to high sinuosity as distance from the fault gradually increases

407

towards the west. This type of progression provides evidence of footwall backtilting, which is

408

frequently observed in other half graben systems, however the direction of channel flow indicates

409

the oblique nature of this particular fault controls the orientation of this slope. When the channels

410

which are fault adjacent are examined in greater detail, there are slight variations in sinuosity

411

between the upper and lower halves of the en echelon segments. Channels which run off the lower

412

half of the unit, display zero to low sinuosity until they reach a moderate distance away from the

413

fault. Alternatively, the channels which propagate in the upper half of the segment show low to

414

moderate sinuosity as they feed down into the transfer relay ramps. This substantiates the

415

hypothesis that there are localised zones of interlinked uplift and depression associated with this

416

style of fault. The variation in sinuosity also helps constrain the different levels of elevation on

417

opposing sides of the en echelon segment, inferring a steeper slope which runs at 45° to fault strike,

418

dipping south from the inside bend of the “S” and a gentler slope which runs between 85°-90° to

419

fault strike dipping north from the outside bend of the “S” and leading into the transfer/relay zone.

420

There is a slightly higher grade of sinuosity in Zone A which does not fit in with the rest of the

421

model, however this is most likely attributed to the fact that the northern section of the fault system

422

underwent the least amount of displacement, due to the angle of extension being sub-optimal in

423

this area. This factor led to repeated rotation of beds in the relay ramp and significant hanging wall

424

deformation (Giba et al., 2012), which may account for the lower gradient slope and offset in this

425

sector.

426

By combining submarine channel flow direction with sinuosity gradient it is possible to

427

postulate a model of the morphology of this fault system and consequently how it might control

428

sediment supply to the surrounding areas.

429

CONCLUSIONS

430

Outlined are the main conclusions drawn from this study;

431



Submarine channels and other fault scarp and transfer zone drainage observed on the

432

hanging wall are directed towards the sub- basin located at the centre of the fault system,

433

where the accommodation space was created by fault activity and sediment thickness maps

434

correlate with fault throw, decreasing along the fault length from middle to tip.

435



Sinuosity profiles in the hanging wall sub-basin indicate an evolution of the channel system

436

present, from immature, high gradient channels to mature, meandering systems as depth

437

decreases. This type of high sinuosity profile and vertical channel accretion combined with

438

en echelon fault transfer drainage indicates a facies model of a continental half graben with

439

transverse through drainage.

440



Channel drainage patterns in the footwall collaborate with sinuosity profiles in order to

441

constrain the topography produced by the en echelon fault shape. The data supports the

442

concept of intra basin highs and lows associated with transfer zones and provides an

443

indication of the associated slopes of each.

444 445



Sinuosity gradient in the footwall provides evidence of footwall back tilting, the orientation of which is controlled by the oblique fault system.

446

While sedimentation patterns on a regional level are controlled by overall topography and

447

tectonics, the trends found in the Parihaka Basin are influenced locally by the main oblique normal

448

fault. Accommodation space created by the original faulting in the area has resulted in a significant

449

hanging wall depocentre, while the complex zones of uplift and subsidence associated with the en

450

echelon fault system has controlled distribution surrounding the transfer zones and relay ramps

451

between segments. Syn-rift abyssal channel development on the footwall, initiated during the

452

oblique reactivation of the fault system, follows a predictable and repeatable pattern, with channel

453

initiation and flow direction being strongly influenced by specific areas in each fault zone.

454

APPENDIX Table of Contents Appendix I: Seismic Sections A) Seismic cross section through the complete survey area, inline vertical section from N-S. B) Representative Seismic Slice S-A, vertical section through the fault at Zone A. C) Representative Seismic Slice S-B, vertical section through the fault at Zone B1 and B2. D) Representative Seismic Slice S-C, vertical section through the fault at Zone C. Appendix II: Channel Measurement Data A) Representative image of the type of channel that was mapped, complete with the set of measurements taken. B) Footwall Channel Measurement Data C) Hanging Wall Measurement Data

APPENDIX I Seismic Cross Sections A) Seismic Section Inline Area

Appendix I (A). Seismic cross section from N-S encompassing the whole area. Annotations include the main features and the relative ages of the layers. B) Seismic Section S-A

Appendix I (B). Representative cross section through the main fault in Zone A of the study area, complete with annotations of the main features.

C) Seismic Section S-B

Appendix I (C). Representative cross section through the main fault in Zone B1 and B2 of the study area, complete with annotations of the main features. D) Seismic Section S-C

Appendix I (D). Representative cross section through the main fault in Zone C of the study area, complete with annotations of the main features.

APPENDIX II A) Representative Channel Image

Appendix I (A). Example of the submarine channels mapped in survey area. Annotations include the relevant parameters that were measured and an example of a thalweg.

B) Footwall Channel Data

CHANNEL NUMBER

FAULT SECTION

SURVEY ZONE

DEPTH

CHANNEL LENGTH

VALLEY LENGTH

-936 -936 -936 -936 -936 -936 -936 -936 -936 -936 -936 -936 -936 -936 -936 -936 -936 -936 -936 -936

m 1449.84 3210.56 479.19 595.84 2119.77 2051.05 3757.40 1384.24 1367.36 2581.28 3174.34 2506.35 7229.85 4259.88 8834.38 1189.98 2155.08 1326.68 1847.44 1660.55

m 1377.29 3005.80 469.94 575.98 1443.76 1442.09 2651.66 1060.73 1238.57 2298.13 2919.40 1513.93 3962.57 3425.77 7507.27 1033.11 1953.30 1181.98 1468.51 1530.45

1 2 3 4 5 6 4 5 6 7 8 9 10 11 12 13 14 15 16 17

FW FW FW FW FW FW FW FW FW FW FW FW FW FW FW FW FW FW FW FW

ZONE C C C C C C C C C C C C B B B A A A A A

FAULT /NON ADJACE NT FA FA FA FA NA FA FA FA FA FA FA FA NA FA NA FA FA FA FA FA

18 19 20 21 22 23 24 25 26 27 28 29 30

FW FW FW FW FW FW FW FW FW FW FW FW FW

C C B B B B B B B A A A A

FA NA FA FA FA NA NA FA NA FA FA FA FA

-984 -984 -984 -984 -984 -984 -984 -984 -984 -984 -984 -984 -984

4413.23 10901.95 4863.14 6084.16 3012.41 9643.29 10821.82 1726.64 11313.20 5455.10 3977.03 4704.75 5664.31

3929.30 7346.25 4141.95 4708.56 2242.29 6303.12 4979.55 1661.07 2250.96 4806.92 3720.17 3165.12 917.71

31 32 33

FW FW FW

C C B

FA FA FA

-1008 -1008 -1008

7854.60 6821.28 4497.21

5040.61 4921.51 3336.28

CHANNEL NUMBER

FAULT SECTION

SURVEY ZONE

DEPTH

CHANNEL LENGTH

VALLEY LENGTH

FAULT /NON ADJACE NT NA

-1008

m 11397.54

m 8574.52

34

FW

ZONE B

35 36 37 38 39 40

FW FW FW FW FW FW

B B A A A A

NA NA FA FA FA NA

-1008 -1008 -1008 -1008 -1008 -1008

17497.46 12787.83 2497.01 5133.71 5375.37 5692.16

9199.64 9908.92 2264.57 4231.23 4543.58 3042.36

41 42 43 44 45 46 47 48 49 50 51 52

FW FW FW FW FW FW FW FW FW FW FW FW

C C C B B B B B A A A A

FA NA NA FA NA FA NA FA FA FA FA FA

-1032 -1032 -1032 -1032 -1032 -1032 -1032 -1032 -1032 -1032 -1032 -1032

5823.83 15376.24 4100.16 8920.92 16851.57 6066.54 11788.78 5806.75 4148.00 4612.17 2001.57 3026.41

5647.75 9803.60 1289.86 7876.85 10438.99 5696.20 7029.97 5000.76 3644.08 3475.55 1449.82 2654.62

53 54 55 56 57 58 59 60

FW FW FW FW FW FW FW FW

C C B B B B A A

FA FA NA FA FA FA FA FA

-1056 -1056 -1056 -1056 -1056 -1056 -1056 -1056

6112.42 5701.24 15280.74 4501.72 9666.93 6366.20 3653.53 1738.59

3385.60 4457.92 8847.03 3945.30 8345.79 5704.61 2960.86 1603.64

CHANNEL NUMBER

1 2 3 4 5

1.05 1.07 1.02 1.03 1.47

>1.15 STRAIGHT * * * *

1.15-1.5 SINUOUS

*

MEANDER

WAVE LENGTH

1 4 0 1 6

0.00073 0.00133 0.00000 0.00174 0.00416

1.15 STRAIGHT 6 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47

1.42 1.42 1.30 1.10 1.12 1.09 1.66 1.82 1.24 1.18 1.15 1.10 1.12 1.26 1.09 1.12 1.48 1.17 1.29 1.34 1.53 2.17 1.04 5.03 1.13 1.07 1.49 6.17 1.56 1.39 1.35 1.33 1.90 1.29 1.10 1.21 1.18 1.87 1.03 1.57 3.18 1.13 1.61 1.07 1.68

1.15-1.5 SINUOUS * * *

MEANDER

WAVE LENGTH

4 3 2 3 8 3 5 12 3 7 2 2 1 1 0 3 8 5 4 1 14 11 1 15 6 1 4 7 5 6 4 12 10 13 1 7 4 9 2 10 5 5 15 2 7

0.00277 0.00113 0.00189 0.00242 0.00348 0.00103 0.00330 0.00303 0.00088 0.00093 0.00194 0.00102 0.00085 0.00068 0.00000 0.00076 0.00109 0.00121 0.00085 0.00045 0.00222 0.00221 0.00060 0.00666 0.00125 0.00027 0.00126 0.00763 0.00099 0.00122 0.00120 0.00140 0.00109 0.00131 0.00044 0.00165 0.00088 0.00296 0.00035 0.00102 0.00388 0.00063 0.00144 0.00035 0.00100

1.15 STRAIGHT 48 49 50 51 52 53 54 55 56 57 58 59 60

1.16 1.14 1.33 1.38 1.14 1.81 1.28 1.73 1.14 1.16 1.12 1.23 1.08

1.15-1.5 SINUOUS *

MEANDER

WAVE LENGTH

4 3 6 1 3 4 3 17 4 6 3 4 1

0.00080 0.00082 0.00173 0.00069 0.00113 0.00118 0.00067 0.00192 0.00101 0.00072 0.00053 0.00135 0.00062