Thermal E ﬀ ects of Magmatism on Surrounding Sediments and Petroleum Systems in the Northern O ﬀ shore Taranaki Basin, New Zealand

: In the past two decades, numerical forward modeling of petroleum systems has been extensively used in exploration geology. However, modeling of petroleum systems inﬂuenced by magmatic activity has not been a common practice, because it is often associated with additional uncertainties and thus is a high risk associated with exploration. Subsurface processes associated with volcanic activity extensively inﬂuence all the elements of petroleum systems and may have positive and negative e ﬀ ects on hydrocarbon formation and accumulation. This study integrates 3D seismic data, geochemical and well data to build detailed 1D and 3D models of the Kora Volcano—a buried Miocene arc volcano in the northern Taranaki Basin, New Zealand. It examines the impact of magmatism on the source rock maturation and burial history in the northern Taranaki Basin. The Kora ﬁeld contains a sub-commercial oil accumulation in volcanoclastic rocks that has been encountered by a well drilled on the ﬂank of the volcano. By comparing the results of distinct models, we concluded that magmatic activity had a local e ﬀ ect on the thermal regime in the study area and resulted in rapid thermal maturation of the surrounding organic matter-rich sediments. Scenarios of the magmatic activity age (18, 11 and 8 Ma) show that the re-equilibration of the temperature after intrusion takes longer (up to 5 Ma) in the scenarios with a younger emplacement age (8 Ma) due to an added insulation e ﬀ ect of the thicker overburden. Results of the modeling also suggest that most hydrocarbons expelled from the source rock during this magmatic event escaped to the surface due to the absence of a proper seal rock at that time.


Introduction
Basin modeling is a powerful tool for integrating geological and geophysical data over a range of scales and provides a means to test and refine concepts related to basin evolution. The term "basin modeling" is used here for dynamic forward modeling of geological and physical processes taking place in sedimentary basins over a certain geological time span [1]. Modeling incorporates deposition, subsidence and compaction, heat flow analysis, as well as hydrocarbon generation, expulsion and accumulation. Numerical forward modeling of petroleum systems is common practice in exploration, but areas with magmatic activity are often perceived as a high risk target for hydrocarbon exploration [2]. There are, however, multiple examples of petroleum systems associated with igneous intrusions all over the globe, some of which contain significant hydrocarbon accumulations (e.g., references [3,4]).
Volcanic activity includes multiple surface and subsurface processes that have an extensive impact on the thermal regime in the surrounding region. Subsurface volcanic processes such as sills/dikes emplacement and magma-chamber formation are of particular interest because the complex interaction between host rock and magmatic intrusion may have both positive and negative effects on all elements of a petroleum system. High heat flow resulting from magmatic activity, for example, may accelerate the thermal maturation of petroleum source rocks, whereas it may also lead to "overcooking" of organic matter and the loss of generated hydrocarbons [5,6]. It is, therefore, important to develop a good understanding of the architecture of a magmatic system, its age, and its geometric relationship to intruded strata, as well as associated heat flow variation through time when assessing a petroleum system. Interaction of magmatic intrusions (sills and dikes) with surrounding organic material may play a significant role in the timing of the hydrocarbon expulsion and accumulation, whereas the timing of emplacement and the size of intrusive body defines the volume of influence [7]. Hydrocarbon generation associated with strata rapidly heated by magmatic intrusion has been studied numerically for several regions [2,[8][9][10][11]. Modeling of magmatic intrusions is, however, rarely conducted in 3D petroleum systems models as their effects are generally considered to be local rather than regional in scale.
The offshore northern Taranaki Basin, New Zealand is an ideal location to study the influence of magmatism on a petroleum system due to a proven hydrocarbon presence within a well-preserved buried submarine volcanic center. To complement these data, exploration boreholes around the Kora volcano provide valuable information for petroleum system model calibration. The results of this study are relevant to other petroleum regions with magmatic activity.
This study focuses on detailed 1D and 3D numerical models of a buried Miocene arc volcano (Kora) to examine the impact of magmatism on source rock maturation and burial history in the northern Taranaki Basin. The geometries of the stratigraphic layers were interpreted from 3D seismic data. Recent studies have focused on the geomorphology, structure, and geometry of sills and dikes associated with the Kora volcano [12][13][14]. Previous numerical studies have focused on 1D thermal modeling and heat anomalies associated with shallow and deep seated intrusions of Kora and their impact on hydrocarbon generation [15]. However, a study that provides a full 3D perspective of magmatic thermal influence on maturation of organic matter is lacking.
The modeling results are used to answer the following scientific questions: • How does the magmatic intrusion affect the thermal maturation of the surrounding source rock? • What is the sensitivity of the results to the shape and distribution of volcanic intrusion? • What are the effects of a deep seated magmatic intrusion (e.g., magma chamber) in comparison to dikes and sills intruding the sedimentary sequence?
These questions, although widely discussed in the last decades based on 1D and 2D studies, remain not fully answered when the third dimension is added. Recent studies on the thermal effects of magmatic intrusion on hydrocarbon maturity show a strong relation between sill thickness and temperature, as well as maturation of the surrounding organic matter [16][17][18][19].
Industry 3D seismic reflection data tied to exploration wells were used to map the major stratigraphic units around Kora to define the geometry and shape of this sub-sea stratovolcano and underlying dikes and sill complexes. We use 3D numerical simulation to assess the influence of magmatic intrusion on the thermal maturation of organic matter and the petroleum system. The effects of deep seated intrusion, sill emplacement, and a system of dikes building a volcanic center are analyzed in this paper. We discuss local lateral changes in heat flow associated with different scenarios and its effect on hydrocarbon generation. To construct petroleum systems models we integrated previously published research data, core and well data, well completion reports, and geochemistry reports available online in the New Zealand Petroleum and Minerals Exploration Database, as well as further geological data from New Zealand Petroleum Exploration Data Pack 2018 provided by New Zealand Petroleum and Minerals.   [22] and Bischoff et al. [12]).

Figure 2.
Cretaceous-Cenozoic stratigraphic framework for the northern Taranaki Basin and interpreted horizons from this study. Principal source rock in the area (Waipawa Formation) is marked with a star (modified after King and Trasher, [20]; Stagpoole and Funnel, [15]). Not to scale. Taranaki Basin developed in the Late Cretaceous as Zealandia separated from the continental landmass of Gondwana [20,23]. The basin contains up to 11 km of Late Cretaceous to Recent strata that records continental breakup through to passive margin subsidence and the subsequent development of the modern plate boundary through New Zealand [20,24].
In its early phase, Late Cretaceous rifting in the Taranaki Basin was characterized by widely distributed NW-NNW to NE trending half-grabens developed in response to the Gondwana breakup and the opening of the Tasman Sea [23,25]. The predominantly Late Cretaceous extensional episodes developed large normal faults with thick terrestrial to marginal syn-rift sequences (2 to 3 km) [20]. In subsequent deformation phases, these Late Cretaceous faults have, in some cases, reactivated multiple times [24]. Extension waned in the Latest Cretaceous/Early Paleocene, resulting in a transition to a passive transgressive margin undergoing gradual thermal subsidence [20,23] [20,23]. The basin contains up to 11 km of Late Cretaceous to Recent strata that records continental breakup through to passive margin subsidence and the subsequent development of the modern plate boundary through New Zealand [20,24].
In its early phase, Late Cretaceous rifting in the Taranaki Basin was characterized by widely distributed NW-NNW to NE trending half-grabens developed in response to the Gondwana break-up and the opening of the Tasman Sea [23,25]. The predominantly Late Cretaceous extensional episodes developed large normal faults with thick terrestrial to marginal syn-rift sequences (2 to 3 km) [20]. In subsequent deformation phases, these Late Cretaceous faults have, in some cases, reactivated Geosciences 2019, 9, 288 5 of 28 multiple times [24]. Extension waned in the Latest Cretaceous/Early Paleocene, resulting in a transition to a passive transgressive margin undergoing gradual thermal subsidence [20,23] The development of the modern Australia-Pacific plate boundary through New Zealand began in the middle Eocene resulting in contraction along the eastern margin of the Taranaki Basin [26]. This contraction, characterized by reverse faulting and foreland basin development, resulted in uplift and erosion of Late Cretaceous to Paleocene strata along the eastern margins of the basin [20,27]. Maximum transgression in the basin occurred during the Oligocene and was followed by regression in the Early Miocene contemporaneous with initiation/reactivation of subduction along the Hikurangi margin to the east [20,24]. Arc volcanism along the Mohakatino Volcanic Centre developed in the northern Taranaki Basin from this time (ca. 18-20 Ma) and continued until ca. 5 Ma [22,28].
Contraction along the eastern basin margin continued until the Mid Miocene (ca. 12 Ma), when oblique subduction and steepening/roll back of the subducting Pacific plate resulted in extension across the northern part of the basin [20,21,29]. Arc volcanism and intra-arc rifting in the northern Taranaki Basin migrated southwards during the Late Miocene to Pliocene to its present location on the Taranaki peninsula (Mt. Taranaki) [21,22] (Figure 1C). From the Latest Miocene, uplift and exhumation of the central North Island and the southern Taranaki Basin resulted in high sedimentation rates, the rapid westward progradation of the shelf, and burial of stratovolcanoes of the northern Taranaki Basin [20,24,30]. The stratigraphic and tectonic framework for the northern Taranaki Basin is shown in Figure 2.
The Kora volcano belongs to a chain of buried arc-related stratovolcanoes and intrusions (Mohakatino Volcanic Centre) that developed in the northern Taranaki Basin from the Early Miocene [20,22,28]. This north-northeast-trending volcanic arc was subsequently dissected during rifting associated with the development of the Northern Graben [20,21,28]. Volcanic edifices of the Mohakatino Volcanic Centre cover about 20% of the Northern Graben area [20]. Radiometric dating and seismic reflection mapping tied to nearby exploration wells indicate that Kora developed between 19.5 and 5.5 Ma [22,28].
Petroleum exploration of the Mohakatino Volcanic Centre started in the 1960s with acquisition of 2D seismic data and drilling of the first well, targeting an off-shore volcanic center (Mangaa-1) in 1970 [12]. The first well into the Kora volcanic structure (Kora-1) was drilled in 1987 by Arco Petroleum NZ Ltd. The target of this drilling operation was a structural high within Eocene sandstones (Tangaroa Formation). However, an oil productive interval was discovered within Miocene volcanoclastics [28,31]. Therefore, based on the estimated hydrocarbon accumulations and reservoir rock properties, the Miocene volcanics became an attractive target for further exploration around the Kora volcano and a further three wells were drilled (Kora-2, Kora-3, and Kora-4).
Kora-1 well penetrated the Eocene to present day stratigraphic sequence represented by deep to shallow marine sediments and volcanoclastics, reaching a depth of 3421 m. Seismic reflection data and mapping show the Kora stratovolcano overlying deformed Early Miocene strata and onlapped by Pliocene-Recent strata (Figure 3; [28]). The oldest stratigraphic units encountered include several hundreds of meters of Eocene Turi Formation, which is predominated by shales with siltstones and sandstones. This unit is directly overlain by late Eocene-early Oligocene turbidite sandstones of the Tangaroa Formation submarine fan complex [32]. The Oligocene-early Miocene section is represented by deep-water limestones of Tikorangi and Taimana Formations. Series of pyroclastic deposits of andesite composition and volcanoclastics related to Miocene volcanic activity represent the equivalent of Mohakatino Formation in the region and possess good reservoir rock properties in several wells drilled around the Kora volcano. This volcanoclastic interval is underlain by up to 1000 m of claystones of Manganui Formation and overlain by a thick layer of progradational Giant Foresets Formation [33]. The main target of the drilling of the Kora-1/1A wells, the Eocene Tangaroa Formation was found to contain hydrocarbons [34], but in a too low permeability reservoir rock to be commercially productive. However, in the Kora-1A well, petroleum was produced from the Miocene volcanosedimentary section with an average flow rate of 668 barrels of oil per day for several days. The other wells either encountered oil shows in Miocene volcanics or were dry. During a drill stream test in well Kora-3, only formation water was produced and core saturation analysis of the well Kora-2 indicated the presence of only residual oil in the formation [34][35][36].
Geochemical analyses indicate that the oils from the Kora discovery were expelled from the Late Paleocene Waipawa Formation [37,38]. The Waipawa Formation is a widespread source rock in New Zealand basins, represented by about 2 to 80 m thick marine mudstone. Higher contributions of terrestrial (mainly woody) plant matter are indicated by higher land plant biomarkers and high phenol/naphthalene ratios [39]. The presence of a Late Paleocene source rock in the study region is confirmed by Ariki-1 well, where the Waipawa Formation is represented by a 12 m thick organic-rich marine shale located at depth of 4122-4134 m with total organic carbon (TOC) averaging 3.6 wt% [40]. A maximum TOC value of 11.5 wt% was derived for the samples from Ariki-1 well [41]. The Waipawa Formation is mainly characterized by terrigenous organic matter [39,41] but HI values are partly above 400 mg/g TOC falling into the field of type II and type II/III source rock. This is consistent  The main target of the drilling of the Kora-1/1A wells, the Eocene Tangaroa Formation was found to contain hydrocarbons [34], but in a too low permeability reservoir rock to be commercially productive. However, in the Kora-1A well, petroleum was produced from the Miocene volcano-sedimentary section with an average flow rate of 668 barrels of oil per day for several days. The other wells either encountered oil shows in Miocene volcanics or were dry. During a drill stream test in well Kora-3, only formation water was produced and core saturation analysis of the well Kora-2 indicated the presence of only residual oil in the formation [34][35][36].
Geochemical analyses indicate that the oils from the Kora discovery were expelled from the Late Paleocene Waipawa Formation [37,38]. The Waipawa Formation is a widespread source rock in New Zealand basins, represented by about 2 to 80 m thick marine mudstone. Higher contributions of terrestrial (mainly woody) plant matter are indicated by higher land plant biomarkers and high phenol/naphthalene ratios [39]. The presence of a Late Paleocene source rock in the study region is confirmed by Ariki-1 well, where the Waipawa Formation is represented by a 12 m thick organic-rich marine shale located at depth of 4122-4134 m with total organic carbon (TOC) averaging 3.6 wt% [40]. A maximum TOC value of 11.5 wt% was derived for the samples from Ariki-1 well [41]. The Waipawa Formation is mainly characterized by terrigenous organic matter [39,42] but HI values are partly above 400 mg/g TOC falling into the field of type II and type II/III source rock. This is consistent with general marine nature of the Waipawa Formation, but with a major contribution of terrestrial -derived woody material [39].

Input Data
The dataset for this study comprises a 3D seismic survey with an approximate area of 900 km 2 acquired by OMV and Argo in 2006, and well data for six exploration wells (Kora-1 and 1A, Kora-2, Kora-3, Kora-4, Ariki-1 and Tangaroa-1) ( Figure 1B). The 3D seismic data are in SEG normal polarity, 12.5 m × 25 m binning, and post-stack time migrated. Well data used for seismic interpretation and depth conversion include check-shot data, tops of the stratigraphic units, and well-log data. Well reports [33][34][35][36]43] and electric logs were used to provide additional input data for the numerical models.
Analysis of 3D seismic reflection data was carried out using reflectivity and variance/coherency (edge detect attribute) analysis. Seismic-to-well tie and an integration with stratigraphic data were carried out to create a framework for the initial seismic interpretation and a subsequent tectono-stratigraphic analysis. 3D seismic was tied to Ariki-1 well using an additional 2D seismic section.
Time-to-depth conversion was performed with Petrel 2016. Check-shot data were used to calculate interval velocities for the interpreted seismic units, and included in a final velocity model used for depth conversion (Table 1). Additionally, a high velocity interval of Miocene volcanics was introduced as a separate horizon. Uncertainties in the time-depth relationships (TDR) in the study area were analyzed by comparing the check-shots of all other wells with the TDR curve derived from synthetic seismograms at well Ariki-1. For the stratigraphic units below Kora volcano that are not present in Ariki-1 and Tangaroa-1 wells, interval velocities were assigned by extrapolating the TDR curve.

1D and 3D Modeling
The general workflow implemented in this study is shown on Figure 4.
To analyze the effect of volcanic activity in the Kora volcano region 1D, models were constructed for four wells (Ariki-1, Kora-1, Kora-4 and Tangaroa-1). In this study 1D models were used to simulate different scenarios for the timing of magmatic activity and to analyze the uncertainties of the models. To reconstruct the thermal and burial history of the area, well data were used that include the information about the depth of stratigraphic tops, formation age and its lithology. The cuttings description was used for defining the lithologies of the sediments and thickness of the stratigraphic units [33][34][35][36]43]. The input data used for the 1D model for Ariki-1 well are summarized in Table 2. Biostratigraphic ages, ages of intrusion emplacement, corrected borehole temperatures, and thermal maturation parameters were used to constrain the models. the information about the depth of stratigraphic tops, formation age and its lithology. The cuttings description was used for defining the lithologies of the sediments and thickness of the stratigraphic units [33][34][35][36]43]. The input data used for the 1D model for Ariki-1 well are summarized in Table 2. Thermal conductivity data for the lithologies in the model were introduced as reported in Funnel et al. [44] for Ariki-1 well. Physical rock properties such as permeability and compaction with depth were used as default settings for the specific lithologies from the software. Compaction calculation within the software is based on Athy's law.
Source rock kinetics for the Waipawa Formation as well as Rock-Eval data were available from the geochemical report for Ariki-1 well [41] (Table 3). Vitrinite reflectance measurements and borehole temperatures were used for calibration of the thermal history for each well [45].  Thermal conductivity data for the lithologies in the model were introduced as reported in Funnel et al. [44] for Ariki-1 well. Physical rock properties such as permeability and compaction with depth were used as default settings for the specific lithologies from the software. Compaction calculation within the software is based on Athy's law.
Source rock kinetics for the Waipawa Formation as well as Rock-Eval data were available from the geochemical report for Ariki-1 well [41] (Table 3). Vitrinite reflectance measurements and borehole temperatures were used for calibration of the thermal history for each well [45]. Table 2. Summary of lithostratigraphy intersected by Ariki-1 well in the northern Taranaki Basin. Thermal conductivity data adopted from Funnell et al. [44]. Lst = Limestone; Org = organic.

Thermal Calibration and Heat Flow
To reconstruct the thermal history of the study area, the knowledge of the present-day and paleo-heat flows is crucial. Present day heat flow was determined by calibration of the model to bottom-hole temperatures that were available from previous studies [45]. The paleo-heat flow from previous regional modeling [46] was used as a basis for the present study and modified to fit the modeled wells using the calibration data by implementing different scenarios of intrusion emplacement. Heat-flow input data and paleo-temperature used in the base case model and temperature history predicted for the Waipawa Formation are shown in Figure 4. The basal heat flow without taking into account magmatism decreases from 62 mW/m 2 at 80 Ma to a present day heat flow of 50 mW/m 2 ( Figure 5). The temperature history indicates the temperatures that a particular layer experienced at a certain point of time and primarily depends on the burial history and heat flow. The upper boundary condition is defined by the sediment/water interface temperature (SWIT), which depends on the paleowater depth (PWD), paleolatitude and climatic evolution [47]. The PWD in the northern Taranaki Basin was assigned based on Strogen et al. [48] and varies between 0 and 1500 m. SWIT was defined using an automatic function in PetroMod that is based on Wygrala [49].

Intrusion Modeling
Sub-volcanic intrusion was modeled for several stages during the Late Miocene. The timing for the Kora volcano activity and Mangaa volcano was adopted based on ages reported by Bergman et al. [28] and Giba et al. [22], who primarily used K-Ar dating and seismic stratigraphy as a method to determine the start and duration of the volcanic activity. Seismic stratigraphy analysis indicates that Kora volcano had been active between 12 and 8 Ma [22]. Results of K-Ar dating of andesitic tuffs from Kora wells give the age of 12.1 ± 2.6 Ma. The intrusion emplacement is a short lived event on the geological time scale and generally lasts no longer than several thousands to tens of thousands of years. However, volcanic activity on a regional scale might continue for several millions of years. In this study, the volcanic activity was modeled in two different ways. First the PetroMod "Intrusion" tool was used with magmatic intrusion being introduced to one or a number of layers. "Intrusion tool" for 1D allows to model a magmatic event by emplacing the "original lithology" by magmatic intrusion with corresponding parameters (type of intrusion, temperature etc.). In 3D, the "Intrusion" tool was used by replacing the host rock with a magmatic intrusion at the predefined area of a facies map. In this scenario (sub-) volcanic magmatic activity was introduced as andesitic type intrusion at 11 Ma to the model as described by Bergman et al. [28]. There is no direct evidence of the mechanism of intrusion emplacement at Kora, such as sill intrusions and/or dykes penetrated by a well. Therefore seismic reflection data showing high amplitude reflectors, interpreted as igneous intrusions, beneath Kora wells and widely distributed beneath the Kora volcano at the interpreted late Paleocene-Late Cretaceous layers were used as an input data source for the magmatic intrusion modeling ( Figures  6 and 7C). The base of layers P10-P00-K90 have been chosen for the emplacement of the intrusion. The intruded material is modeled to replace the original sedimentary rock with an andesitic melt, the thermal properties of which are shown in Table 4. We have modeled a second scenario of intrusion emplacement that represents a deep seated intrusion below the modeled section. In this scenario, the basal heat flow is dramatically increased at 11 Ma from 43 mW/m 2 to 200 mW/m 2 . The thermal effect of a deep seated intrusion gradually decreases from the basement to the upper layers.

Intrusion Modeling
Sub-volcanic intrusion was modeled for several stages during the Late Miocene. The timing for the Kora volcano activity and Mangaa volcano was adopted based on ages reported by Bergman et al. [28] and Giba et al. [22], who primarily used K-Ar dating and seismic stratigraphy as a method to determine the start and duration of the volcanic activity. Seismic stratigraphy analysis indicates that Kora volcano had been active between 12 and 8 Ma [22]. Results of K-Ar dating of andesitic tuffs from Kora wells give the age of 12.1 ± 2.6 Ma. The intrusion emplacement is a short lived event on the geological time scale and generally lasts no longer than several thousands to tens of thousands of years. However, volcanic activity on a regional scale might continue for several millions of years. In this study, the volcanic activity was modeled in two different ways. First the PetroMod "Intrusion" tool was used with magmatic intrusion being introduced to one or a number of layers. "Intrusion tool" for 1D allows to model a magmatic event by emplacing the "original lithology" by magmatic intrusion with corresponding parameters (type of intrusion, temperature etc.). In 3D, the "Intrusion" tool was used by replacing the host rock with a magmatic intrusion at the predefined area of a facies map. In this scenario (sub-) volcanic magmatic activity was introduced as andesitic type intrusion at 11 Ma to the model as described by Bergman et al. [28]. There is no direct evidence of the mechanism of intrusion emplacement at Kora, such as sill intrusions and/or dykes penetrated by a well. Therefore seismic reflection data showing high amplitude reflectors, interpreted as igneous intrusions, beneath Kora wells and widely distributed beneath the Kora volcano at the interpreted late Paleocene-Late Cretaceous layers were used as an input data source for the magmatic intrusion modeling (Figures 6  and 7C). The base of layers P10-P00-K90 have been chosen for the emplacement of the intrusion. The intruded material is modeled to replace the original sedimentary rock with an andesitic melt, the thermal properties of which are shown in Table 4. We have modeled a second scenario of intrusion emplacement that represents a deep seated intrusion below the modeled section. In this scenario, the basal heat flow is dramatically increased at 11 Ma from 43 mW/m 2 to 200 mW/m 2 . The thermal effect of a deep seated intrusion gradually decreases from the basement to the upper layers.    In the simulations for 1D models of Kora-1 and Kora-4 wells and the 3D models we assumed that that the magma temperature corresponds to the lower limit of andesite melts upon extrusion (950 °C). The heat pulse generated by the intrusion is highly dependent on the intrusion temperature and its thickness and significantly influences the overall output of the model. Convective heat transfer associated with volcanic activity such as hydrothermal fluid flows was not modeled in this study. In the simulations for 1D models of Kora-1 and Kora-4 wells and the 3D models we assumed that that the magma temperature corresponds to the lower limit of andesite melts upon extrusion (950 • C). The heat pulse generated by the intrusion is highly dependent on the intrusion temperature and its thickness and significantly influences the overall output of the model. Convective heat transfer associated with volcanic activity such as hydrothermal fluid flows was not modeled in this study.

Model Setup
A 3D model covering about 900 km 2 was built with a grid resolution of 250 × 250 m. This resolution allows for the construction of a detailed volcanic intrusion geometry and provides sufficient detail on the lithology distribution around the Kora volcano. Distribution of lithologies in the 3D model was based on interpolation between Kora wells. Lithologies for deeper stratigraphic units of Eocene to Late Cretaceous age were assigned based on the paleogeographic maps from Strogen [48]. The volcano edifice was introduced in the model as a uniform intrusion body in the shape of a cone with constant physical properties. The maps of stratigraphic units (see next chapter) were integrated with well data through lithological modeling and used as an input data for the 3D numerical model. A model simulating magmatic intrusion placed in basement rock was constructed using an increase in basal heat flow.

Seismic Interpretation in Kora 3D
Eighteen pre-, syn-and post-volcanic sequences were identified and mapped around the buried Kora volcano-sedimentary edifice in the offshore northern Taranaki Basin. These interpreted horizons were used as an input data for the 3D model. Out of eighteen interpreted horizons, thirteen represent key stratigraphic unit boundaries and six are interpreted within the Pliocene to recent progradational system of the Giant Foresets Formation (Figure 3 Figure 1C); (ii) secondary faults of primarily Eocene age ( Figure 7A); (iii) radial faults around the volcanic body ( Figure 1B). These radial faults surrounding the volcanic edifice are the result of cooling and collapse episodes that are typical for volcanic bodies [13,51]. Some of them are connected to the faults of the Northern Graben SE of the volcano ( Figure 1C) and have a higher offset.
Several seismic high amplitude zones were identified in the seismic data. High amplitudes located in the Giant Foresets Formation and in the shallow section directly above the volcanic edifice are interpreted to be related to gas-charged sediments ( Figure 7B). These are likely associated with biogenic gas generation or volcanic gases [13]. In the Paleocene and Eocene intervals, strong reflections occur in the vicinity of the Kora volcano-sedimentary edifice. These reflectors are interpreted as sill intrusions ("saucer sills") ( Figure 7C) and were also described by Bischof et al. [12] and Infante-Paez and Marfurt [13]. They differ in amplitude strength from shallower arcuate reflections that represent incised channels ( Figure 7D).

Seismic Horizons
The modern seafloor geometry is relatively smooth, with the slope gently deepening southeast and water depth increasing from 105 to 170 m ( Figure 8A).
The horizons representing the Giant Foresets Formation are significant for the understanding of the regional geology, variation in sedimentation rate and burial history. Pliocene horizons (N80 and N70) follow continuous positive reflectors. The reflection trend indicates that the sedimentary system prograded from southeast to northwest ( Figure 7B,C). The seismic units between the reflectors are characterized by large sigmoidal clinoforms with channel systems on the foresets.
The base Pliocene (N60) reflection is less continuous and in places displaced by normal faults and graben structures. These structures are typical for the region and can be attributed to the Late Miocene extension phase when the Northern Graben evolved.
Reflector N40 indicates the Middle to Late Miocene boundary ( Figure 8D). The unit between N40 and N60 representing the Miocene strata shows a partially continuous reflection pattern that becomes weaker toward the southwest of the seismic cube. As in the unit above, the reflectors are strongly affected by normal faulting, indicating the deepening of the Northern Graben southeast of the Kora volcano.
The Early Miocene (N10) reflector represents the top of the pre-magmatic activity sequence [12]. The unit below the reflector N10 is of continuous thickness of about 100 m. The reflector pattern is subparallel throughout the area. The unit shows less structural deformation in comparison to the units above, but several large faults are crossing the Kora volcano in the north and northwest. The horizons representing the Giant Foresets Formation are significant for the understanding of the regional geology, variation in sedimentation rate and burial history. Pliocene horizons (N80 and N70) follow continuous positive reflectors. The reflection trend indicates that the sedimentary system prograded from southeast to northwest ( Figure 7B,C). The seismic units between the reflectors are characterized by large sigmoidal clinoforms with channel systems on the foresets.
The base Pliocene (N60) reflection is less continuous and in places displaced by normal faults and graben structures. These structures are typical for the region and can be attributed to the Late Miocene extension phase when the Northern Graben evolved.
Reflector N40 indicates the Middle to Late Miocene boundary ( Figure 8D). The unit between N40 and N60 representing the Miocene strata shows a partially continuous reflection pattern that Horizons P50 (Eocene-Oligocene), P30 (Middle Eocene) and P10 (Paleocene-Eocene) were mapped based on the well tops from well Ariki-1. P50 reflector has a strong continuous positive response throughout the study area and shows a strong doming of strata underneath the Miocene volcano sedimentary edifice ( Figure 3) that is still present in the depth-converted surface ( Figure 8E). The units between these horizons are consistently subparallel and deformed by the major faults of the Northern Graben.
The Waipawa Formation was mapped between horizons P00 and P10. The approximate top Waipawa reflection (P10) corresponds to a negative reflection around well Ariki-1. Although there is no evidence from seismic data that the Waipawa Formation is present underneath the Kora volcano, the fact that it is related to a regional paleooceanographic event suggests a wide distribution in the northern Taranaki Basin [39,40]. The Top Cretaceous reflector P00 is also indicated in the Ariki-1 well ( Figure 3) as a SEG normal soft kick. This reflector is difficult to interpret throughout the 3D seismic, due to poor data quality below the Kora volcanics. In the central part of the dataset the Cretaceous horizon was interpolated, as it is crucial for the subsequent modeling to have a continuous surface.
Late Cretaceous unit tops K90 and K80 are represented by positive reflectors. Similar to the top Cretaceous (P00), the signal is partially lost in the deeper parts of the survey and underneath the volcanics. K90 was interpreted in Ariki-1 well; it has a SEG normal soft response. The continuous interpretation of the reflectors was difficult and, as for the Paleocene the top and bottom horizons (P10 and P00), had to be extrapolated underneath the volcanics.
The top acoustic basement is defined at well Ariki-1 as a low to medium amplitude positive reflection. The seismic to well tie from Ariki-1 to the Kora 3D seismic reflection survey correlates the basement along the same reflection loop, and the reflector becomes stronger in the faulted sequence and on basement highs ( Figure 8F). The acoustic basement is difficult to detect beneath the volcanic centers.

Thickness Variations
The Late Cretaceous syn-rift sequence discordantly overlays the faulted basement and fills the syn-rift structures. The thickness of this unit is dependent on the original basement topography and the characteristics and distribution the faults. The largest thicknesses are 3315 m east of Kora volcano and gradually decreases to the west where it reaches its minimum of about 740 m ( Figure 9A).

1D Modeling Results
Simulation results from 1D models provide information on the timing of oil and gas generation. Kora-1, Kora-4, Ariki-1 and Tangaroa-1 wells are good candidates for assessing the amount of heat produced by volcanos in the northern Taranaki Basin. 1D models provide detailed information about litho-stratigraphy and rock properties in the region. Tangaroa-1 is located northeast of the Kora volcano and was also drilled in the volcaniclastics; together with Ariki-1 well it was chosen for A prominent change in sediment thickness can be observed related to Late Miocene foreland basin evolution, where the sediment distribution was dramatically affected by volcano topography. The unit is thinning from about 400 to 0 m at the volcano center ( Figure 9C).
The progradational sequence evolved in the area during Pliocene time (N70-N80) and reaches its maximum thickness of 1560 m in the SSE, but it is only 880 m thick in the SSE ( Figure 9D).

1D Modeling Results
Simulation results from 1D models provide information on the timing of oil and gas generation. Kora-1, Kora-4, Ariki-1 and Tangaroa-1 wells are good candidates for assessing the amount of heat produced by volcanos in the northern Taranaki Basin. 1D models provide detailed information about litho-stratigraphy and rock properties in the region. Tangaroa-1 is located northeast of the Kora volcano and was also drilled in the volcaniclastics; together with Ariki-1 well it was chosen for comparison with the Kora wells. These wells intersect volcanoclastics and minor volcanic or intrusive rocks of Miocene age. Vitrinite data were used as reported in well reports [33][34][35][36]43] and were used for heat flow calibration.
The burial history reconstruction for Kora-4 well is shown on Figure 9A. About 4 km of subsidence occurred during Late Cretaceous-Paleocene rifting. Subsidence continued until the late Miocene, when the eastern margin of northern Taranaki Basin was inverted during development of the convergent plate boundary. Sedimentation rates markedly slowed from the late Eocene onwards due to sediment starvation of the basin and maximum transgression during the Oligocene and an increase in water depth to about 1400 m [48,52]. Sedimentation rates increased during the Miocene, and in particular after late Miocene inversion, when the shelf propagated across the Kora area and the Giant Foresets Formation was deposited. The highest sedimentation rates are documented for the last 5 Ma (Pleistocene-Pliocene) for all modeled wells, filling most of the previously created accommodation space. At the same time, the deepest burial was reached.
In Figure 10E, burial histories for the Tangaroa Formation are compared for the four wells. Evolution during the Paleogene is similar for all wells and low sedimentation rates prevail. Starting in the lower Miocene, the pattern becomes more complex and differences between the wells are more prominent. During the Plio-Pleistocene periods, Ariki-1 and Tangaroa-1 wells experienced even more burial than the two Kora wells.
The temperature overlay on the burial plot for Kora-4 well illustrates the thermal history of the stratigraphic column. Aside from the magmatic heat event at 11 Ma, there is a temperature increase during the last 5 Ma which is associated with the rapid subsidence of the basin due to high sedimentation rates during formation of the Giant Foresets propagation system during the Pliocene.
The modeled sub-volcanic intrusion emplacement as a consequence of the main volcanic activity leads to the temperature spike affecting both deeper and shallow layers. In the intrusion scenario-1, volcanic activity was modeled as subvolcanic intrusion into the Paleocene and Late Cretaceous, between layers P00 and K90 ( Figure 10B). The second scenario of intrusion emplacement represents a deep seated intrusion below the modeled section, where the basal heat flow was increased at 11 Ma from 53 mW/m 2 to 200 mW/m 2 . The thermal effect of a deep-seated intrusion is more gradual with depth than of a local intrusion in scenario 1 ( Figure 10C).
Both scenarios modeling intrusion emplacement provide a better fit to the calibration data than the scenario without considering magmatic activity. Without such local heat sources, vitrinite reflectance data in Kora-1 cannot be reproduced (Figures 11 and 12). In Kora-1 well a significant increase in vitrinite reflectance values occurs below 3000 m depth. For Kora-4 the deepest sample measured for vitrinite reflectance is at 2700 m depth. Below this depth, a pronounced increase of organic matter maturity is expected, as seen in Kora-1 ( Figure 11), but data are lacking.   The temperature overlay on the burial plot for Kora-4 well ( Figure 10A) illustrates the thermal history of the stratigraphic column using scenario 1. Sub-volcanic intrusion led to an increase of temperature in surrounding sediments including the Paleocene sediments of the Waipawa  The temperature overlay on the burial plot for Kora-4 well ( Figure 10A) illustrates the thermal history of the stratigraphic column using scenario 1. Sub-volcanic intrusion led to an increase of temperature in surrounding sediments including the Paleocene sediments of the Waipawa The more pronounced increases in vitrinite reflectance in Tangaroa-1 and Ariki-1 may also be due to local or hydrothermal heating. A recent example for such a system has been reported from ODP (Ocean Drilling Project) Site 859, where the mid-ocean ridge is subducted below forearc sediments leading to very high overall heat flows, strong local thermal effects [53] as well as maturation of organic matter [54].
The temperature overlay on the burial plot for Kora-4 well ( Figure 10A) illustrates the thermal history of the stratigraphic column using scenario 1. Sub-volcanic intrusion led to an increase of temperature in surrounding sediments including the Paleocene sediments of the Waipawa Formation. Compared to a scenario without any igneous heating event (blue line in Figure 11), this leads to a predicted heating of the Waipawa Formation to about 250 • C, compared to about 140 • C in the model without magmatic activity.
The model results demonstrate that magmatism has had profound effects on the thermal history of the Kora lithological column intersected by the well. Figure 13 shows a maturity overlay onto the burial history plot. The blue, green, red and yellow areas represent thermally immature, oil window, gas window, and over-mature zones, respectively. The oil and gas windows are influenced by the kinetic parameters and heat flow used. The striking feature in the maturity zone overlay is the sharp increase in maturity at the time of magmatic intrusion ( Figure 13). Figure 13A shows a change of the maturity in the scenario where several magmatic intrusions were placed in Paleocene and Cretaceous sediments. Prior to magmatic intrusion, the vitrinite reflectance values for the Waipawa Formation reached 0.73 %Ro. This parameter rapidly increases when organic matter gets heated by an intrusion located in the direct proximity to the organic rich sediment layer. By the time magmatic activity ceased, the organic matter of the Waipawa Formation became overmature (4.66 %Ro). The scenario with magmatic activity modeled in the basement indicates a smoother transition in maturity compared to the magmatic intrusion in the sediments scenario. Figure 12 B illustrates that the change of maturity for the Waipawa Formation starts at 13 Ma and gradually increases from 0.73 %Ro to 0.85 %Ro at 11.5 Ma, reaching its maximum of 2.2 %Ro after 11 Ma. Both scenarios indicate the location of the magmatic intrusion relative to the organic rich layer is the key parameter controlling the source rock maturation. The changes in temperature, and thus the geothermal gradient, allow rocks to pass through the liquid-and gas-generating maturity zones. Both scenarios of intrusion emplacement and deep seated intrusion show that the Paleocene Waipawa Formation becomes overmature in the vicinity of the intrusion.

3D Modeling Results
The 3D model integrating the geometries of the stratigraphic layers and facies distribution allows to study the lateral effects of sub-volcanic intrusions on the surrounding sediments. Facies modeling in Petrel software was carried out based on well log interpretation (gamma ray, spontaneous potential, and resistivity) and well cuttings descriptions. The lateral and vertical changes in lithology were considered. Facies distribution in Ariki-1 and its corresponding correlation with interpreted units is illustrated in Figure 3. The seabed to N80 interval is dominated by poorly consolidated silts. The N80 to N70 interval is defined by interbedded mudstones and siltstones, where the uppermost part is dominated by siltstones and the lower part of the interval by mudstones. A volcano-sedimentary interval, as interpreted for Kora-1 well from the gamma ray log and cuttings description, is identified for the well at the depth between 1800 m and 2480 m. The underlying interval penetrated by Kora-4 well includes Miocene mudstones with intercalated sandstone of varying thickness and Oligocene limestone. The sediments of the lowermost units are represented by a deep marine succession of marlstones, shales, and sandstones of the Tangaroa Fm. The Waipawa Formation is modeled as a uniform organic rich black shale layer of 10 m thickness. A transect through the study area showing facies distribution as implemented in the model is shown in Figure 14. Calibration followed the same principles as described for 1D. Three scenarios of intrusion timing were tested, at 18 varying thickness and Oligocene limestone. The sediments of the lowermost units are represented by a deep marine succession of marlstones, shales, and sandstones of the Tangaroa Fm. The Waipawa Formation is modeled as a uniform organic rich black shale layer of 10 m thickness. A transect through the study area showing facies distribution as implemented in the model is shown in Figure  14. Calibration followed the same principles as described for 1D. Three scenarios of intrusion timing were tested, at 18, 11, and 8 Ma, either with very high basal heat flows of 200 mW/m 2 during initial intrusion time or with replacement of sedimentary lithologies by hot andesite.  Figure 15A-C shows the differences in modeled thermal regimes within the study area at the location of Kora-1 well. All the models show a gradual increase of the temperature with time for the different layers and then an abrupt increase during magmatic activity. Temperatures drop rapidly shortly after igneous emplacement and then increase again significantly in the last 3-4 Ma before the present. The latter effect is related to the rapid subsidence associated with increased sedimentation rates during deposition of the Giant Foresets Formation. The temperature spike related to the intrusion (placed in the Paleocene and Late Cretaceous sediments) is the result of a short-lived event: a duration of magmatic activity of 1 Ma was assumed based on studies of magmatic activity by Baumgartner et al. [55]. It takes about 2-3 Ma until temperatures have re-equilibrated, which is strongly dependent on thermal conductivities of the sedimentary and magmatic rocks. Figure 15D shows the three different temperature scenarios for the Waipawa Formation with modeled ages of magmatic intrusion in the sediments. The re-equilibration of the temperature after intrusion takes longer (up to 5 Ma) in the scenarios with younger emplacement age due to an added insulation effect of the thicker overburden. The present-day temperatures ( Figure 16D) still show a variation across the buried Kora volcanic edifice. This is mainly due to the difference in lithologies and doming of the stratigraphic layer related to magmatic intrusion. The volcanic and volcanoclastic lithologies have a higher thermal conductivity, which results in a chimney effect (similar to salt domes) reducing the thermal insulation and hence the temperature at the Waipawa Formation. This observation is consistent with 1D modeling results.  In Figure 16A-C, modeled thermal maturity of the Waipawa Formation is shown for an intrusion timing of 11 Ma into Cretaceous to Paleogene rocks and three time steps: 12 Ma, 8 Ma, and present-day as well as present-day temperature distribution (Figure 16 D). Maturity close to the volcanic center is clearly determined by Miocene volcanic activity, whereas maturity at distances greater than 1-2 km is determined by late (Plio-Pleistocene) burial, reflecting present-day temperatures. magmatic intrusion in the sediments. The re-equilibration of the temperature after intrusion takes longer (up to 5 Ma) in the scenarios with younger emplacement age due to an added insulation effect of the thicker overburden. The present-day temperatures ( Figure 16D) still show a variation across the buried Kora volcanic edifice. This is mainly due to the difference in lithologies and doming of the stratigraphic layer related to magmatic intrusion. The volcanic and volcanoclastic lithologies have a higher thermal conductivity, which results in a chimney effect (similar to salt domes) reducing the thermal insulation and hence the temperature at the Waipawa Formation. This observation is consistent with 1D modeling results.
In Figure 16A-C, modeled thermal maturity of the Waipawa Formation is shown for an intrusion timing of 11 Ma into Cretaceous to Paleogene rocks and three time steps: 12 Ma, 8 Ma, and present- Figure 16. A-C Maturity change at the unit of Waipawa Formation before and after intrusion. D-Present day temperature at Waipawa Formation. The Waipawa Formation has been drilled in the Ariki-1 well, where it is represented by a 12 m thick layer of organic rich marine shale. The hydrocarbon potential of the source rock was evaluated based on the geochemical analysis of source rock samples in that well [37]. Rock-Eval measurements of Waipawa Formation samples in Ariki-1 have yielded TOC values of 11.5 %, a hydrogen index (HI) of 351 mg/gTOC and Tmax of 421 • C [41]. In order to model hydrocarbon generation in the Kora area, custom kinetics for the Waipawa Formation measured on Ariki-1 samples were used [41]. Geochemical analysis and the Kora discovery indicate that Waipawa Formation in the northern Taranaki Basin is oil-prone [41].
The model results demonstrate that both the high sedimentation rates during the last 5 Ma and Miocene magmatic activity have a profound effect on the thermal history of the lithological column around Kora volcano. Hydrocarbon generation was accelerated by the local temperature increase in direct proximity to any sills and dikes. The thermal impact of the igneous intrusion extends laterally more than 5 kilometers from the intrusion center, but maturation effects at distances greater than about by the volcanic intrusion became overmature, as indicated by the yellow color on the map (>4 %Ro). This effect rapidly decreases from the center of intrusion with vitrinite reflectance values changing from 3.8 to 2.0 %Ro within about 1000 m from intrusion ( Figure 16B) and not visible at more than 2000 m distance. Close to volcanic activity the temperatures at the level of the Waipawa Formation is predicted to have increased during the intrusion emplacement from about 120 • C to 260 • C and returned to the pre-emplacement level within approximately 2 Ma. The thermal effect of magmatic intrusion in the basement with a local increase in basal heat flow has significantly less impact on the thermal maturation of the source rock in comparison to magmatic intrusion in the Paleocene -Cretaceous sediments.
Modeling results suggest that the Waipawa Formation was primarily immature around the Kora volcano before the volcanic activity started. Expulsion/generation of hydrocarbons for three modeled scenarios of Kora-1 and Kora-4 wells is shown in Figure 17. For Kora-4 well, the scenarios with deep seated intrusion in the basement and magmatic intrusion in sediments show the same results, with a rapid increase in maturation at 11 Ma associated with high temperature. In the base case scenario for Kora-1 and Kora-4 wells (Figure 17), just 50% of kerogen conversion in Waipawa source rock takes place for Kora-4 and less than 5% in Kora-1 at 11 Ma. Modeling results of different age of volcanic activity suggest that hydrocarbons locally generated by the magmatic activity would have been expelled from the Waipawa Formation prior to the deposition of a proper seal rock. Absence of seal and/or insufficient overburden would likely have resulted in a loss of most of the hydrocarbons generated. However, the Waipawa Formation that was not affected by the volcanic activity still retained generative potential. Hydrocarbons have likely been expelled from the source rock during the later burial associated with deposition of the Giant Foresets Formation in the deeper area of the Taranaki Basin and migrated to the Kora structure.
around Kora volcano. Hydrocarbon generation was accelerated by the local temperature increase in direct proximity to any sills and dikes. The thermal impact of the igneous intrusion extends laterally more than 5 kilometers from the intrusion center, but maturation effects at distances greater than about 1 to 2 km are overlapped by later deep burial and related temperature increase ( Figure 16). The model with magmatic intrusion in sediments shows that in areas that were not affected by volcanic activity, vitrinite reflectance values are in the range between 0.5 and 1.2 %Ro. The source rock directly affected by the volcanic intrusion became overmature, as indicated by the yellow color on the map (>4 %Ro). This effect rapidly decreases from the center of intrusion with vitrinite reflectance values changing from 3.8 to 2.0 %Ro within about 1000 m from intrusion ( Figure 16B) and not visible at more than 2000 m distance. Close to volcanic activity the temperatures at the level of the Waipawa Formation is predicted to have increased during the intrusion emplacement from about 120 °C to 260 °C and returned to the pre-emplacement level within approximately 2 Ma. The thermal effect of magmatic intrusion in the basement with a local increase in basal heat flow has significantly less impact on the thermal maturation of the source rock in comparison to magmatic intrusion in the Paleocene -Cretaceous sediments.
Modeling results suggest that the Waipawa Formation was primarily immature around the Kora volcano before the volcanic activity started. Expulsion/generation of hydrocarbons for three modeled scenarios of Kora-1 and Kora-4 wells is shown in Figure 17. For Kora-4 well, the scenarios with deep seated intrusion in the basement and magmatic intrusion in sediments show the same results, with a rapid increase in maturation at 11 Ma associated with high temperature. In the base case scenario for Kora-1 and Kora-4 wells (Figure 17), just 50% of kerogen conversion in Waipawa source rock takes place for Kora-4 and less than 5% in Kora-1 at 11 Ma. Modeling results of different age of volcanic activity suggest that hydrocarbons locally generated by the magmatic activity would have been expelled from the Waipawa Formation prior to the deposition of a proper seal rock. Absence of seal and/or insufficient overburden would likely have resulted in a loss of most of the hydrocarbons generated. However, the Waipawa Formation that was not affected by the volcanic activity still retained generative potential. Hydrocarbons have likely been expelled from the source rock during the later burial associated with deposition of the Giant Foresets Formation in the deeper area of the Taranaki Basin and migrated to the Kora structure.

Sensitivity Analysis
The following paragraphs are meant to provide an idea on how the sedimentary system in the Taranaki Basin reacts to changes of basal heat flow due to a deep seated intrusion and on changes in

Sensitivity Analysis
The following paragraphs are meant to provide an idea on how the sedimentary system in the Taranaki Basin reacts to changes of basal heat flow due to a deep seated intrusion and on changes in magmatic intrusion thickness placed in the sedimentary column. Testing the sensitivity of the model to heat flow changes and proximity to the organic rich layer helps to quantify the effects of these parameters on temperature and maturity of the surrounding sediments.
Heat flow in volcanically active regions in the direct proximity to the volcano reach up to 500 mW/m 2 . The heat flow values increase near a volcano [56]. We modeled five scenarios with a basal heat flow of 100, 200, 300, 400 and 500 mW/m 2 . The other parameters were not changed. The resulting difference in temperature and maturity are shown in Figure 18A. The Kora 1 data fit best to a high basal heat flow of 300 mW/m 2 , while at greater distance from the volcanic center lower heat flows are expected.

Sensitivity Analysis
The following paragraphs are meant to provide an idea on how the sedimentary system in the Taranaki Basin reacts to changes of basal heat flow due to a deep seated intrusion and on changes in magmatic intrusion thickness placed in the sedimentary column. Testing the sensitivity of the model to heat flow changes and proximity to the organic rich layer helps to quantify the effects of these parameters on temperature and maturity of the surrounding sediments.
Heat flow in volcanically active regions in the direct proximity to the volcano reach up to 500 mW/m 2 . The heat flow values increase near a volcano [56]. We modeled five scenarios with a basal heat flow of 100, 200, 300, 400 and 500 mW/m 2 . The other parameters were not changed. The resulting difference in temperature and maturity are shown in Figure 18A. The Kora 1 data fit best to a high basal heat flow of 300 mW/m 2 , while at greater distance from the volcanic center lower heat flows are expected.
Recent studies indicate that the thickness and proximity of magmatic intrusion to organic matter rich layers are major factors responsible for the maturation and petroleum generation [16,57]. To test the sensitivity of the thickness and distribution of the magmatic intrusion we assumed two magmatic sills both with andesitic magma. We used three scenarios changing the thickness of these sills from 25 to 50 and 75 m each, that clearly show the intrusion thickness has an enormous influence on temperature and maturity of surrounding sediments, in particular between the two intrusions ( Figure 18B).  Recent studies indicate that the thickness and proximity of magmatic intrusion to organic matter rich layers are major factors responsible for the maturation and petroleum generation [16,57]. To test the sensitivity of the thickness and distribution of the magmatic intrusion we assumed two magmatic sills both with andesitic magma. We used three scenarios changing the thickness of these sills from 25 to 50 and 75 m each, that clearly show the intrusion thickness has an enormous influence on temperature and maturity of surrounding sediments, in particular between the two intrusions ( Figure 18B).

Conclusions
We present the results of a seismo-stratigraphic interpretation and 1D and 3D basin modeling, that is used for exploring the effect of magmatic events on the thermal history and petroleum generation around the Kora volcano. The modeling study represents the first investigation of potential interaction of the Waipawa Formation and magmatic intrusion in 3D in the northern Taranaki Basin, New Zealand. The thermal history of the study area was constrained using vitrinite reflectance data from four wells, with two of them intersecting the flank of the Kora structure, and two located within 15 km from it, as well as knowledge of the tectonic settings. Distinct scenarios of the intrusion emplacement and its timing provide information on the timing of hydrocarbon expulsion and illustrate the interaction of hydrocarbon generation and widespread magmatic activity during the Late Miocene time.
The modeling results show that the Waipawa Formation was thermally immature before the magmatic activity started in the late Miocene. Modeled magmatic intrusion in the Paleocene and Late Cretaceous sediments (P00 and K90) at 11 Ma caused rapid heating of the surrounding sediments and maturation of the regional source rock (Waipawa Formation). Depending on the proximity of the source rock to the intrusion, the Waipawa Formation became overmature or was in the gas window directly after the magmatic activity ceased. The scenario with magmatic intrusion placed in the basement (basal heat flow increase) shows similar results. A sensitivity study involving variations of the intrusion age (18 Ma,11 Ma and 8 Ma) indicates that Waipawa Formation would experience higher temperatures if magmatic activity occurs at 8 Ma due to the coupling of two events strongly influencing the thermal history: magmatic activity, and rapid subsidence related due to deposition of Giant Foresets Formation. The same magmatic event modeled at 11 and 18 Ma has less effect on the thermal regime at the level of source rock because the system cooled down more rapidly.
3D modeling results show that the lateral effect of magmatic activity has a local effect on the maturation of the source rock. The thermal aureole caused by volcanic activity at the level of the source rock gradually decreases from its center influencing the area in a radius of 5 km, with a strong impact being restricted to 1-2 km.
The modeling results suggest that most of the hydrocarbons expelled from the Waipawa Formation due to local magmatic heat escaping to the surface, because at the time of expulsion there was no well-developed seal rock that could prevent hydrocarbon leakage ( Figure 19). reflectance data from four wells, with two of them intersecting the flank of the Kora structure, and two located within 15 km from it, as well as knowledge of the tectonic settings. Distinct scenarios of the intrusion emplacement and its timing provide information on the timing of hydrocarbon expulsion and illustrate the interaction of hydrocarbon generation and widespread magmatic activity during the Late Miocene time.
The modeling results show that the Waipawa Formation was thermally immature before the magmatic activity started in the late Miocene. Modeled magmatic intrusion in the Paleocene and Late Cretaceous sediments (P00 and K90) at 11 Ma caused rapid heating of the surrounding sediments and maturation of the regional source rock (Waipawa Formation). Depending on the proximity of the source rock to the intrusion, the Waipawa Formation became overmature or was in the gas window directly after the magmatic activity ceased. The scenario with magmatic intrusion placed in the basement (basal heat flow increase) shows similar results. A sensitivity study involving variations of the intrusion age (18 Ma, 11 Ma and 8 Ma) indicates that Waipawa Formation would experience higher temperatures if magmatic activity occurs at 8 Ma due to the coupling of two events strongly influencing the thermal history: magmatic activity, and rapid subsidence related due to deposition of Giant Foresets Formation. The same magmatic event modeled at 11 and 18 Ma has less effect on the thermal regime at the level of source rock because the system cooled down more rapidly.
3D modeling results show that the lateral effect of magmatic activity has a local effect on the maturation of the source rock. The thermal aureole caused by volcanic activity at the level of the source rock gradually decreases from its center influencing the area in a radius of 5 km, with a strong Figure 19. Petroleum system chart for the study area.
Author Contributions: Conceptualization, A.K., R.L and K.K.; methodology, A.K., S.B., K.K. and H.S.; software, A.K., S.B. and K.K.; validation, A.K., K.K., and S.B.; formal analysis, A.K., R.L. and K.K.; investigation, A.K.; resources, A.K., K.K., H.S.; data curation, A.K., R.L., K.K., H.S and S.B.; writing-original draft preparation, A.K.; writing-review and editing, K.K., R.L., H.S. and R.L.; visualization, A.K.; supervision, R.L. and K.K.; project administration, A.K., R.L. and K.K. All authors confirmed they have contributed to the intellectual content of this paper and have met the following four requirements: (a) significant contributions to the conception and design, acquisition of data, or analysis and interpretation of data, (b) drafting or revising the article for Author Contributions: Conceptualization, A.K., R.L and K.K.; methodology, A.K., S.B., K.K. and H.S.; software, A.K., S.B. and K.K.; validation, A.K., K.K., and S.B.; formal analysis, A.K., R.L. and K.K.; investigation, A.K.; resources, A.K., K.K., H.S.; data curation, A.K., R.L., K.K., H.S and S.B.; writing-original draft preparation, A.K.; writing-review and editing, K.K., R.L., H.S. and R.L.; visualization, A.K.; supervision, R.L. and K.K.; project administration, A.K., R.L. and K.K. All authors confirmed they have contributed to the intellectual content of this paper and have met the following four requirements: (a) significant contributions to the conception and design, acquisition of data, or analysis and interpretation of data, (b) drafting or revising the article for intellectual content, (c) final approval of the published article, and (d) agreement to be accountable for all aspects of the article, which ensures that questions related to the accuracy or integrity of any part of the article are appropriately investigated and resolved.
Funding: RWTH Aachen University, Germany and the New Zealand Ministry of Business, Innovation and Employment through the GNS Science-led research program on New Zealand petroleum source rocks, fluids, and plumbing systems (contract C05X1507) supported this project.