Kevin C. Hill

The Porgera Transfer Zone (PTZ) is a major crustal and probably lithospheric structure across Papua New Guinea recording >50 km offset of ophiolites and very different patterns of geology and topography on either side. In the Late Jurassic, the PTZ probably separated oceanic crust and thick Jurassic Om shales to the west from a continental promontory to the east. During the Late Miocene to Recent orogenesis, the differential compression of these features is interpreted to have created a dextral strike slip fault across the fold belt with pull-apart basins at sites of fault relays. This facilitated the ascent of intrusions and mineralization at Porgera. The acquisition of high-resolution LIDAR data semi-regionally around the Porgera Gold Mine greatly improved interpretation of the regional geology and particularly the recognition of normal faults. By correlating with sparse dip data and paly-dated samples, it was possible to create stratigraphic sections and interpret structural cross-sections using the LIDAR data. As the area involved strike–slip offsets, it was important to construct sections in multiple orientations in order to interpret the 3D geology. Both dips and fault orientation could be directly inferred from the LIDAR data such that sections could be constructed orthogonally to them. A balanced, restored and forward-modelled cross-section illustrates the interaction between thrust faults and normal faults during compression and that it was synchronous with the development of a pull-apart basin. A semi-regional 3D geological model, which was developed mainly from the LIDAR data, supports the hypothesis of inversion of the thick Om beds to the west before or during compression of the continental promontory to the east resulting in dextral strike–slip offsets across the PTZ. A jog or relay in the faults occurred and caused a pull-apart collapse basin to develop in the area of the Porgera mine. Similar pull-apart graben, or negative flower structures, were detected nearby and may be areas for future exploration.

The acquisition of regional, 1 metre resolution LIDAR in the PNG Highlands combined with 3D modelling in MoveTM has revolutionized our understanding of the evolution of the Porgera gold mine and areas of new potential. The new 3D model demonstrates active pull-apart basins along a regional transfer during ongoing fold and thrust deformation. When overlain on regional aeromagnetic data, new potential is revealed. The Porgera gold mine is one of the richest in the world and lies in a wide valley at an elevation of 2800m surrounded by mountains up to 4000m high in the middle of the PNG Highlands. It lies along a major lineament, the Porgera Transfer Zone (PTZ), cutting across the orogen that is associated with a 50 km offset of the ophiolite belt. To the NW of the mine there is an extensive belt of low- to high-grade metasediments that formed in deep water during Jurassic rifting, but were metamorphosed during Eocene subduction to the north. Subduction ceased due to accretion of the Sepik Terrane which caused inversion and mild erosion of the Porgera area. The Late Oligocene onset of wrenching in the Mobile Belt to the north placed that area into extension, emplacing metamorphic core complexes, and led to regional subsidence in the Early Miocene. Collision of the margin with the Melanesian Arc in the Middle Miocene caused Late Miocene to Pliocene orogenesis creating the broad mountain belt that we see today, that is still active, as shown by the 2018 MW 7.5 earthquake. Regional mapping of the area using high resolution LIDAR in association with limited field mapping, analysis of gravity and magnetics data and drilling of 300-500m deep core-holes has allowed development of a detailed 3D structural model. The Porgera valley is constrained laterally by the Eastern Boundary and Western Boundary dextral strike-slip faults that lie along the PTZ cutting across the orogeny. However, whilst the position and nature of the faults is clear, due to bending, fracturing and offset of major synclines and anticlines, the faults do not crop-out as significant through-going features. Rather, the dextral offset of basement is soft-linked to the Tertiary limestones at surface due to intervening thick, ductile Cretaceous shales. Fold and thrust structures are ubiquitous in the Tertiary limestones, but the youngest structural features are large extensional faults, particularly around the Porgera valley. The intrusive underlying the Porgera ore-body was emplaced at 6.0+-0.3 Ma, the time of maximum compression during orogenesis. At this time the Eastern Boundary fault and Western Boundary fault were both active allowing a pull-apart basin to form locally along the PTZ at the south-eastern boundary of the Jurassic metasediments. This enabled emplacement of the intrusive. The gold-bearing fluids from the intrusive and metasediments were brought up through the Mesozoic muds to the contact with the Paleogene carbonates where the Porgera ore body was emplaced.

Definition of the basement architecture within the Ceduna Sub-basin is poorly understood due to the depth of the basin and limited deep crustal geophysical datasets. This study uses seismic interpretation of BightSPAN™ deep crustal seismic lines combined with 2D forward models of the gravity to define the basement architecture in the Ceduna Sub-basin. Gravity forward models identify two depocentres with maximum depths of 25 km overlying thinned continental crust. Syn-rift sedimentary packages are several km thick and underlie up to 10 km of postkinematic Cretaceous deltaic sequences. Regional Moho models contain only sparse data points offshore and isostatic residual gravity data suggests substantial local variation in the depth to Moho. This local variation in the Moho surface is interpreted as boudinage of a weak lower crust along crustal shear zones. A gradual increase in basement density towards the continent-ocean transition suggests a diffuse rather than distinct boundary between continental and oceanic crust. The results of this work are used to constrain the tectonic evolution, particularly in terms of the subsidence history and uplift and erosional pulses recorded on the seismic data in the Ceduna Delta (Hill et al., this volume) and then used as input for finite element numerical models of rifting (Farrington et al., this volume).

A regional, high-quality seismic traverse from the coast to oceanic crust across the Bight Basin has been assembled and interpreted in detail, then balanced, restored, decompacted and replaced at paleo-water depths. The Late Cretaceous Ceduna Delta developed above a Late Jurassic-Early Cretaceous rift basin in three stages punctuated by significant pulses of uplift and erosion across areas >100 km wide and with up to 1 km of erosion. The Cenomanian White Pointer delta prograded into deepening water and hence underwent gravitational collapse. This was terminated in the Santonian when the Antarctic margin was pulled out from below, thus supplying heat to a remnant thicker outer margin crust causing doming and erosion. Importantly, this established the saucer-shaped geometry of the Ceduna Delta that persisted throughout its development, so that any hydrocarbons generated in the southern half of the basin would have migrated towards this outer margin high. The Tiger Formation was deposited in shallow water in a full rift basin prior to breakup which was followed by regional thermal subsidence. The Hammerhead delta developed on the newly formed passive margin, but was terminated by another pulse of uplift and erosion, perhaps associated with a Paleogene change in plate motion at the end of the Cretaceous. Finite element modelling of this proposed tectonic evolution will test its validity and predict hydrocarbon generation and migration through time.

This paper outlines a Finite Element modelling workflow using the software package ‘Underworld’ that is applied to sedimentary basin undergoing multiple rifting events. The evolution of the resulting strain rate, viscosity structure, temperature field, and sedimentary structures can be tracked. Application to the Ceduna Sub-basin is discussed.

A crustal-scale, fully restored section across the PNG orogenic belt reveals the Oligocene to Recent compressional deformation of the margin. The northern end of the section comprises the Landslip Metamorphics, an accreted continental terrane, separated from the main part of the Fold Belt by the Jurassic April Ultramafics and Om Metamorphics, interleaved with Eocene volcanics which together constitute an accretionary prism. Existing maps show that the suture is overlain by distal Miocene sediments indicating Oligocene docking and probable compression prior to Early Miocene subsidence. The latter is consistent with Early Miocene extension in PNG and the emplacement of metamorphic core complexes in the Sepik area, but is also related to dynamic topography causing subsidence of the whole northern margin of Australia. Neogene compression commenced around 12 Ma with ~70km shortening in the Om terrane and ~38km shortening in the Fold Belt. Existing thermochronology data indicate shortening of ~12mm/year from 12-4 Ma, but only 2.5mm/year from 4-0 Ma, consistent with a change in structural style in the Fold Belt from thrust to more ductile, fold-dominated deformation. The model also requires substantial thickening of the continental crust beneath the Muller Ranges, here represented by ‘basement’ underthrusting. Gravity modelling indicates the presence of sedimentary graben up to 10 km deep beneath the Fold Belt, which were strongly inverted, such as beneath the Lavani Valley. A key issue is when this inversion occurred, in the Oligocene or Pliocene, as this has a significant influence on the timing of hydrocarbon generation and migration.

Southeast Asia experienced a complex tectonic and geodynamic history related to the subduction of the eastern Tethyan ocean basins, resulting from the long-term convergence between the Indo-Australian, Eurasian, and Pacific plates since Pangea breakup. The complex collage of continental and island arc terranes can be reconstructed into an estimated ancient arrangement using plate tectonic reconstruction approaches based on a synthesis of continental and marine geological and geophysical data. We use the open-source and cross-platform software GPlates (www.gplates.org) to refine the evolution of the eastern Neo-Tethys since the latest Jurassic rifting episodes along northern Gondwana. We apply the resulting plate motions to drive numerical models of mantle flow in order to predict the evolving mantle structure. New Guinea’s northward motion over subducted slabs, related to the Sepik back-arc basin and the Maramuni subduction system, resulted in long-term flooding of the margin since ~20 Ma, despite falling long-term global sea levels. The Sundaland continental promontory experienced dynamic uplift in the latest Cretaceous to Eocene times due to the accretion of the Woyla Arc at ~80 Ma, leading to slab breakoff and a temporary interruption of subduction. However, renewed subduction along the Sunda margin resulted in renewed dynamic subsidence from ~30 Ma, which was amplified by regional basin rifting events. In addition, the sinking Sunda slab likely triggered a mantle slab avalanche, resulting in a counterintuitive combination of contemporaneous basin inversion and strong dynamic subsidence from ~15 Ma. The evolution of the eastern Tethyan oceanic gateway provides an important framework for understanding the role of plate tectonics in controlling long-term oceanic circulation and climate, as well as shedding light on the complex interplay between deep Earth and surface processes in driving basin formation and evolution. These results provide new avenues for reconciling stratigraphic and tectonic processes, as well as contributing new approaches for basin analysis and hydrocarbon exploration.

Apatite fission track (AFT) and vitrinite reflectance (VR) data for early Miocene outcrops from the Waitemata Basin reveal that the basin sequence was subjected to shallow burial before denudation. AFT results suggest that the total sediment thickness within the basin was ≤1 km and maximum paleotemperatures during burial never exceeded c. 60°C. Statistical analyses of the detrital AFT ages distinguish four dominant sources of sediment supply: contemporaneous volcanism; metagreywacke rocks of the Waipapa Group; the Northland Allochthon; and an unidentified source south of the basin. The apatite and zircon fission track results from the Waipapa Group rocks (Gondwana Terrane) adjacent to the basin suggest two discrete phases of accelerated cooling: the first during the early Cretaceous (c. 117 Ma) and the second during the mid Cretaceous (c. 84 Ma). These events probably reflect key stages in the tectonic development of the New Zealand microcontinent during the Cretaceous period, the earlier event being related to the climax of compressional deformation (Rangitata Orogeny) and the latter to extensional tectonism associated with the opening of the Tasman Sea. Waipapa Group rocks now exposed at the surface cooled from maximum paleotemperatures of c. 250°C at an estimated rate of c. 180–36°C/m.y., involving substantial denudation.