Issue
Bull. Soc. géol. Fr.
Volume 188, Number 5, 2017
Petroleum source rocks
Article Number 30
Number of page(s) 18
DOI https://doi.org/10.1051/bsgf/2017188
Published online 10 November 2017

© V. Crombez et al., Published by EDP Sciences 2017

Licence Creative CommonsThis is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

1 Introduction

In the past decade, with the development of shale plays, there was a renewed interest in describing and characterizing in greater details fine-grained and organic-rich deposits (Bruns et al., 2015; Chalmers and Bustin, 2012; Jarvie 2012a, b; Løseth et al., 2011; Passey et al., 2010). Recent studies show that the location of organic accumulations results from the interplay of three controlling factors: primary productivity, preservation and dilution (Bohacs et al., 2005; Slatt and Rodriguez, 2012) and that these factors are strongly impacted by the regional geodynamic settings (Trabucho-Alexandre et al., 2012). Understanding the complex and dynamic interactions between these controlling factors requires a full 3D quantification of the primary productivity, the dilution and the preservation of organic matter at basin scale. Stratigraphic forward modeling provides a powerful tool to achieve this goal.

The aim of the present work is to use the recent development in stratigraphic forward modeling of organic-rich rocks (DionisosFlow and DORS) (Chauveau et al., 2013; Granjeon, 2014) to better understand the impact on organic accumulation of the spatial and temporal evolution of the primary productivity, sedimentation rates and oxygen levels at basin-scale in the Montney and Doig Formations.

Our work focuses on the Lower and Middle Triassic Montney and Doig Formations from the Western Canada Sedimentary Basin. Recent studies (Crombez et al., 2016a) proposed a 3D stratigraphic architecture and conceptual models for the formation of organic rich heterogeneities in the Montney and Doig Formations. Based on these previous works, the present study aims at testing different scenarios of primary productivity, basin restrictions and organic matter transportation through stratigraphic modeling, in order to predict the location of organic-rich heterogeneities and discuss allocyclic controls on primary productivity, dilution and preservation.

2 The Montney and Doig Formations

The Montney and Doig Formations are present in the foreland basin of the Canadian Cordillera in the subsurface of Alberta and British Columbia (Fig. 1a). These Formations were deposited during the Early and Middle Triassic (Golding et al., 2015b; Orchard and Zonneveld, 2009) on the western margin of Pangaea and recent studies suggest a fore-arc basin geodynamic setting (Beranek and Mortensen, 2011; Rohais et al., 2016). At the time of the deposition, the basin was located in latitudes inferior to 30°N (Fig. 1b). These low latitudes imply a dry and arid climate over the basin during the deposition of the studied interval (Davies et al., 1997). The Lower and Middle Triassic strata mainly consist of siltstone and fine-grained sandstone deposited in wave dominated environments (Armitage, 1962; Davies et al., 1997; Zonneveld et al., 2010). Detailed analysis of the sedimentology highlighted an interval of turbiditic deposits (Moslow and Davies, 1997) and an interval of phosphatic shale deposits (Golding et al., 2014). Recent work showed that the studied interval was deposited in a marine basin with a maximum estimated bathymetry of 300 m (Crombez, 2016; Crombez et al., 2016a). Study on the organic content shows that these formations can be considered as a type II/III marine source rock with actual hydrogen index up to 600 mgHC/gTOC and TOC up to 12 wt% (Riediger, 1997; Riediger et al., 1990). In the Montney and Doig Formations, two intervals present significant organic matter accumulation, the upper Montney (initial TOC of non-thermicaly mature rocks (TOCini) up to 8 wt% and average TOCini of 2.2 wt%) and the basal Doig phosphates zone (TOCini up to 14 wt% and average TOCini of 1.1 wt%) (Crombez et al., 2016b).

Recent studies on the stratigraphic architecture of these formations (Crombez, 2016; Crombez et al., 2016a) show that the Montney and Doig/Halfway Formations respectively represent the expression of two second order stratigraphic cycles (sequence A and B). Moreover, the Montney Formation (sequence A) is subdivided in three third order sequences (sequence 1, 2 and 3) that respectively represent the second order TST (transgressive systems tract), HST (highstand systems tract) and FSST (falling stage systems tract); whereas the Doig and Halfway Formations (sequence B) were only subdivided into 2 system tracts: the transgressive system tract 4 (TST4) and the highstand system tract 4 (HST4, Fig. 2) that are interpreted as second order TST and the HST-FSST. The sequence boundary between the two second order cycles is marked by a major erosion that is the result of a second order lowstand sea-level (Fig. 2). In the Montney Formation, preserved strata only present sediment inputs from the East and the South-East whereas in the Doig and Halfway Formations sediment inputs are also coming from the West (Fig. 2).

Nutrient supply in this interval is interpreted to be linked to both continental (sequence A) and marine sources (sequence B). Connection to upwelling cells, and marine nutrient supply in sequence B are interpreted by Crombez et al. (2016b) as resulting from the regional geodynamic evolution. Over this period, anoxic layer development is interpreted to be linked to basin restriction (sequence A3) or to high primary productivities (TST4).

In the Lower and Middle Triassic strata, two significant organic accumulations are present respectively in the sequence 3 and in the TST4. A recent study (Crombez et al., 2016b) proposed three conceptual models (Fig. 2) of basin types to explain the organic matter accumulations in the Montney and Doig Formations:

  • type 1: an oxic basin (sequences 1, 2 and HST-FSST 4), with no restriction and a moderate primary productivity that lead to limited accumulation of organic matter (TOCini < 3 wt%);

  • type 2: a restricted basin (sequence 3), with high basin restriction and moderate productivity that lead to moderate concentration of organic matter in the central part of the basin (1 wt% < TOCini < 5 wt%);

  • type 3: a high primary productivity basin (TST4), with no restriction and an increased primary productivity that lead to accumulation of very rich interval on the margins of the basin (TOCini up to 14 wt%).

The sequence 3 and the TST4 contain organic-rich deposits and they represent a perfect case study to investigate the regional controls on primary productivity, dilution and anoxia.

thumbnail Fig. 1

A. Location of the study area (faults from Reed et al., 2005, Triassic subcrop edges from Edwards et al., 1994; basin limits from Wright et al., 1994). B. Paleogeographic map during the lower Triassic over the Western Canada sedimentary basin. Qs: sedimentary fluxes; Ab.: Alberta; B.C.: British Columbia; Y.: Yukon; N.T.: Northern Territories.

thumbnail Fig. 2

Stratigraphic framework of the Lower and Middle Triassic of the Western Canada Sedimentary Basin (modified from Crombez, 2016, eustatic variations from Hardenbol et al., 1998).

3 Data and methods

3.1 Data from stratigraphy

The isopach, bathymetric and paleogeographic maps used as input data for the model were derived from Crombez (2016). The depositional architecture developed in this study resulted from the correlation of over 400 wells, using the sequence stratigraphic method and concepts defined in Hunt and Tucker (1992) and Catunaenu et al. (2009). It provided maps for the area where the Lower and Middle Triassic strata are still present (Fig. 1). A reconstruction of the missing parts of these formations was carried out, based on our understanding of the regional stratigraphic architecture and of the paleogeographic and geodynamic (Beranek and Mortensen, 2011; Crombez, 2016; Rohais et al., 2016) evolution of the basin.

3.2 Numerical modeling

The aim of this study is (1) to simulate the stratigraphic architecture of the Montney and Doig Formation and quantify the controlling factors (accommodation, sediment supply and transport parameters), (2) to test the effects of different primary productivity values and basin restriction conditions on organic matter distribution and (3) to test different scenarios of the organic matter transport.

3.2.1 Stratigraphic modeling principles

DionisosFlow is a process-based multi-lithology diffusion driven forward modeling tool developed by IFPEN (Granjeon, 2014). Based on a diffusion algorithm that simulates the transport of different class of sediments at basin scale, the model allows for testing the relative impact of numerous parameters on the basin fill architecture (e.g. sediment flux, location of the sources, subsidence history, eustatic variations) (Csato et al., 2014; Gvirtzman et al., 2014; Hawie et al., 2015; Leroux et al., 2015). Diffusion of the sediments in DionisosFlow mainly depends on two factors: (1) the gravity: the conversion of potential energy to kinetic energy, (2) the water: the transport capacity of the water, that are calibrated using the diffusion coefficients Kigravi Kiwater respectively linked to the slow creeping transport and water driven displacement. As this diffusion only account for transportation induced by gravity or water fluxes, it can be completed by a transport equation that incorporates wave induced diffusion (associated with  Kiwave). It is important to notice that theses equations are only relevant over long time period and at basin scale, DionisosFlow cannot simulate single events (e.g. floods, crevasse splay, storm), sedimentary architectural elements (e.g. levees, dunes…) or hydrodynamic processes.

Recently, a new model, DORS, has been plugged into DionisosFlow. The DORS model aim at simulating the production, exportation, transportation and burying of organic matter (Chauveau et al., 2013, 2016). As a matter of simplification, the numerical modeling does not aim at representing the complex nutrient transportation and delivery to the basin, but only the primary productivity variations resulting from these processes. In modern oceans, at the scale of the present study chlorophyll concentration (Sigman and Hain, 2012) tends to show that regardless of the nutrient sources (continental runoff or upwelling cells), the primary productivity is correlated to the distance to shore. In DORS, the primary productivity is therefore defined as a function of the distance to the shore, decreasing toward the basin center. Once produced, the organic matter is degraded while settling through the water column, following the law proposed by Martin et al. (1987). Then the organic matter is transported along the sea floor as silt particles. Lastly, the organic matter is buried and degraded using burial efficiency laws (Burdige, 2007). Here the preservation of the OM is a function of the sea-bottom oxygen concentration, and sedimentation rates. This equation globally represents the different processes (chemical and biological) that occurs on the sea bed (Chauveau et al., 2013, 2016).

3.2.2 Model generalities and input data

The simulated area is 1000 km by 500 km with a 10 × 10 km2 grid cells. It is centered on the present-day Montney and Doig subcrop area (Fig. 1). The simulation starts at the base of the Lower Triassic (252.2 Ma) and ends at the top the Middle Triassic [237 Ma, in the present work chronostratigraphic boundaries are taken from Cohen et al. (2013)]. It covers a time interval of 15.2 My with 0.2 My time step (76 time steps).

The model takes into account five sediment classes: sand, silt, shale, bioclastic carbonate and organic matter. In order to constrain the simulation, six paleo water depth/topography maps and five thickness maps were used. The first five paleo water depth maps and four thickness maps were used to define subsidence. The last elevation and thickness maps were used to compute the compaction at present day using compaction laws as a function of burial depth for each sediment class.

The eustatic variations used in the model is presented on the Figure 3. In order to fit the stratigraphic evolution, the eustatic curve from Hardenbol et al. (1998) was modified (Crombez, 2016, long term and medium term on Fig. 3) in order to simulate the regional subsidence/uplift due to the tectonic activity associated with the geodynamic evolution of the basin.

Lastly, the location and rates of the sediment inputs were defined based on paleogeographic reconstructions from Crombez (2016) and Rohais et al. (2016). During the deposition of Montney Formation (sequence A), the sediments mainly come from the eastern, southern and northern part of the model, whereas during the deposition of the Doig-Halfway Formations (sequence B), the sediments mainly come from the East, South and West (Fig. 4). The diffusion coefficients (K) used in the present study are reported in Table 1. These coefficients are calibrated in order to reproduce realistic depositional slopes and sediment partitioning at basin scale (Alzaga-Ruiz et al., 2009; Burgess et al., 2006; Granjeon, 2014).

thumbnail Fig. 3

Eustatic variations used in the stratigraphic modeling. The composite curve is based on Hardenbol et al. (1998) and on two user defined relative sea-level curves that represent the impact of regional uplifts and subsidence.

thumbnail Fig. 4

Location of the sedimentary inputs of the stratigraphic model. Note the turnover between sequence A and sequence B.

Table 1

Diffusion coefficients used in DionisosFlow.

3.2.3 Tested scenarios

Three scenarios corresponding to the three aforementioned basin types were developed to investigate the impact of diffusion coefficient variation on organic matter distribution (Tab. 2). In these scenarios, the diffusion coefficient of the organic matter (KOM) was successively set equal to diffusion coefficients of sand (Ksand), silt (Ksilt) and shale (Kshale). This screening test aimed at selecting the optimal diffusion coefficient to be used in the rest of the study.

In a second step, scenarios were designed to test the impact of primary productivity and basin restriction on organic matter accumulation (Tab. 3). In these scenarios, the primary productivity is a function of the distance to the shore: end-values are defined at the shoreline and at 150 km from the shore. The primary productivity is interpolated in between these two extreme values and is not affected by the bathymetry or the emerged area topography. This gradient aims at representing the trend of the primary productivity along continental margin. Primary productivity values defined in the models range from 225–112.5 gC/m2/a to 112.5–50 gC/m2/a. In DORS, the water stratification due to the basin restriction is simulated through a basin restriction coefficients (Rwater). This coefficient affects in association with the bathymetry the oxygen renewal in marine environments. It ranges from Rwater = 1 (open marine) to Rwater = 0.05 (highly restricted) (Chauveau et al., 2013, 2016). The objective of this work is to estimate and discuss the controls on primary productivity and basin restriction conditions required to reproduce the organic accumulation and oxygen conditions described in the three basin types.

Table 2

Scenarios tested in order to study the impact of organic matter diffusion on organic accumulations.

Table 3

Scenarios tested in order to study the organic accumulation in the different basin type.

4 Modeling results

4.1 Best-fit scenario

4.1.1 3D stratigraphic forward model

The best-fit scenario is presented on Figure 5. On this figure, the color scale represents the variation of a synthetic gamma ray, which is computed based on the relative sediment proportion GR = 4 × %shale + 2 × %silt + 0.1 × %sand). As this synthetic gamma ray is based on sediment proportion, it will not show high gamma ray values that are induced by authigenic uranium enrichment. It will therefore not represent the high gamma ray values linked to high uranium concentration in sequence 3 and TST4 (Crombez et al., 2016b).

The model of Figure 5 reproduces the trends of the stratigraphic architecture of the Montney and Doig Formations highlighted in previous studies: (1) the maximum backstep of the shoreline during the early sequence 2; (2) a long regression from the highstand 2 up to the end of sequence 3; (3) the unidirectional filling of the sedimentary basin during sequence 1, 2 and 3 and the bi-directional filling during sequence 4. Moreover, the Figure 5 highlights typical sedimentary geometries present in the Montney and Doig Formations: (1) prograding clinoforms; (2) sandy turbiditic interval (Davies et al., 1997; Moslow and Davies, 1997).

Errors between the best-fit scenarios and the case study are presented on Figures 6 and 7. The Figure 6 shows a relative error that is generally inferior the 30% in the Montney and inferior to 50% in the Doig and Halfway. The increased error in this last interval is linked to the non-calibration of the Halfway thickness that is not considered as a potential organic rich rock. However, thickness in the Doig Formation was better calibrated and therefore resulted in a lower relative error (lower part of sequence B, Fig. 7). In addition to the maps, Figure 7 illustrates the comparison between real data and pseudo-wells extracted from the simulation at specific well locations. On the wells, the depths of the stratigraphic surfaces are from Crombez (2016), on the synthetic wells the surface are placed over the time step that corresponds to the age of the sequence boundary.

Sedimentation rates have a direct impact on the preservation and the dilution of the organic content (Tyson, 2001), it is therefore important that for each sequence, at basin scale, thickness errors are inferior to 25 m. On this figure, it is apparent that some of the sedimentological details present on the well logs are not represented in the model. Despite these differences, the best fit model successfully reproduces overall stratigraphic trends and sediment repartition; with reasonable thickness errors. It can therefore be used to further investigate the controls on organic matter accumulations.

thumbnail Fig. 5

Cross sections in the best fit model. Here the color is related to an estimated gamma ray based on the lithology proportion. This figure highlights the typical geometries of the Montney Formation (A) and the Doig Formation (B).

thumbnail Fig. 6

Thickness and errors maps for each modeled sequence. It is apparent that relative errors are low (< 30%) and therefore acceptable.

thumbnail Fig. 7

Comparison between synthetic wells extracted from the best fit model and real well from the same location in the Western Canada Sedimentary Basin. Note the thickness similarities between numerical simulated data and real data.

4.1.2 Organic matter transport

The results of the experiments summarized in Table 2 are presented on Figure 8. This figure shows similar impact of the diffusion coefficient changes on the three basin types: the organic matter tends to concentrate in the center of the basin. On Figure 8, the concentration of the organic matter in the central part of the basin is more pronounced in the type 3 configuration and weaker in type 1 configuration. The Figure 9 presents a dip-oriented profile of organic matter concentration (location on Fig. 8) together with associated oxygenation and bathymetry for each configuration.

The scenarios with high diffusion show a relative depletion (compared to low diffusion scenarios) of organic matter in areas located at less than 140 km of the shoreline and a concentration of the organic deposits beyond 140 km. The amplitude of organic matter variations is higher in type 3 (TOCini up to 4.5 wt%) than in type 2 (TOCini up to 2.5 wt%) and type 1 (TOCini up to 0.5 wt%) basins. However, in the deepest part of the basin, the relative increase of organic matter concentration is close to 70% in both type 2 and type 3 basins, which shows that the relative impact of an increase of organic matter diffusion is the same regardless of the primary productivity. Here, the increase of the transport coefficient induced a displacement of the maximum accumulation area (Fig. 9): between Ksand and Kshale, the maximum accumulation area, shifted 50 km seaward. Lastly, on the three configurations, the Figure 9 shows no change in the bathymetry profile and only minor change in the oxygen concentration.

Various studies focused on the sedimentary processes that drive the formation of mudstone and others fine-grained sediments (Ghadeer and Macquaker, 2011; Macquaker and Bohacs, 2007; Macquaker et al., 2010; Schieber and Southard, 2009; Schieber et al., 2007). All these studies highlight that mudstone and others fine-grained sedimentary rocks get mainly deposited in environments that still present unidirectional or oscillatory currents. In marine source rocks, the organic fraction that reaches the water sediment interface is mainly composed of faecal pellets, flocculates or other aggregates (Ghadeer and Macquaker, 2011) called the “marine snow” (Alldredge and Silver, 1988). Therefore, the organic fraction presents on the sea bed is composed of various particles with various sizes and densities (Alldredge and Silver, 1988; Wakeham and Lee, 1993). Recent study (Schieber et al., 2007) suggested that the organic sediments may have similar transport characteristics than fine-grained sandstone, but it does not take into account the wide variety of particles that compose the sedimentary organic matter.

In DORS, the organic matter is considered as a sediment class with its own transport coefficient. In this study, we tested various diffusivity coefficients in order to properly reproduce the organic matter accumulations (from Ksand to Kshale, Figs. 8 and 9). In the center of the basin, the organic matter concentration was approximately 70% higher when using Ksand than when using Kshale, regardless of the primary productivity or oxygenation of the water (Fig. 9). Although DORS does not reproduce explicitly the flocculation or aggregation processes, we found that using a transport coefficient of silt to fine-grained sand is a good approximation for modeling basin-scale distribution of organic matter.

thumbnail Fig. 8

Organic content maps extracted from the model in order to test the influence of the diffusion coefficient on the organic matter distribution in the different basin types. On this figure a weak concentration in the deepest part of the basin in apparent. HP: high productivity; MP: medium productivity; LP: low productivity; HD: high diffusion; MD: medium diffusion; LD: low diffusion.

thumbnail Fig. 9

TOCini, benthic oxygen and bathymetry profile in the three basin types. This figure shows a significant enrichment in organic matter in the deepest part of the basin with the increase of the diffusion coefficient. Note that this concentration of the organic matter has no impact on the benthic oxygen profile.

4.2 Sediment inputs variation

The values of the sediment and water fluxes used in the best-fit model are presented in the Figure 9A. Average sediment inputs are higher in the Montney Formation (sequence A: 28.103 km3/Ma) than in the Doig-Halfway Formations (sequence B: 5.103 km3/Ma). Within the Montney Formation, the total sediment input increases from sequence 1 up to sequence 3 (from 16.103 km3/Ma to 43.103 km3/Ma) and the limit between the Montney and the Doig Formation is marked by a sharp drop of total sediment inputs. In sequence A, the relative proportion of sediment delivered from the three sources (N, E, S) is not constant (Fig. 9A). During sequence 1, the maximum sediment input is coming from the North, whereas in sequence 2 the maximum sediment input is coming from the East. In the Montney Formation, the sediments coming from the South only contribute 15 to 20% of the total sediment supply, whereas the source from North or from East provides up to 60% of the sediments inputs. During the deposition of the Doig Formation, sources from West represent more than 40% of the sediment input whereas the sources from East and South respectively carry 35% and 25% of the sediments. Similar trends are also visible in water discharge (Fig. 9B), the Montney Formation (sequence A) shows water discharge from 1.103 m3/s to 1.105 m3/s. The Montney/Doig boundary presents a major turnover with a sharp drop in the total water discharge (6.102 m3/s).

4.3 Primary productivity and restriction variations

In order to investigate the impact of the primary productivity and basin restriction on anoxia and organic matter distribution, nine scenarios were tested on the best fit model. Bottom water oxygen concentration and organic matter maps were extracted from the model at intervals representative of type 1, type 2 and type 3 basins. The bathymetries corresponding to those intervals are presented on Figure 11. It shows that the bathymetry in type 1 and 2 are deeper than in type 3. It also shows a major change in the physiography of the basin: type 1 and 2 intervals only present sediment supply from the East whereas the interval of type 3 presents sediment supply from the East and the West. The Figures 1214 present the benthic oxygen and organic matter concentration map resulting from the sensibility analysis.

thumbnail Fig. 11

Bathymetric maps of the interval representative of type 1 (A), type 2 (B) and type 3 (C) basins. On this figure, it is apparent that the bathymetry in type 3 basin is shallower than in type 1 and 2.

4.3.1 Type 1 basin

On Figure 12, it is apparent that all scenarios with high restriction and the model with medium restriction and medium to high primary productivity present anoxic conditions in the central part of the basin. On Figure 12, only LP-0.15 (low organic productivity – basin restriction coefficient = 0.15 ∼ moderately restricted) and the three low restriction scenarios present oxic conditions. Here, only scenarios with high or medium primary productivity, associated with medium to high restriction present organic accumulation higher than 1.5 wt% (with a maximum of 4 wt% in HP-0.05 scenario). The Figure 12 shows that the model with MP-0.15 and HP-1 present the same distribution of organic matter than the simulation with MP-1 and LP-0.05. Among the nine scenarios tested, only the three models with low restriction resulted in a fully oxic environment with low organic accumulations (Fig. 12).

4.3.2 Type 2 basin

The Figure 13 shows similar results than Figure 12 due to the similar shape of the basin (see Fig. 10): the three models with high restriction present anoxia as well as the scenarios with medium restriction associated with medium to high primary productivity. Oxic conditions are prevalent in models with low restriction and in the model with moderate restriction and low productivity (LP-0.15 scenario). The organic content of the nine scenarios showed in Figure 13 are significantly lower than the one displayed in Figure 12 because of the higher sedimentation rates and the resulting dilution of organic matter (Fig. 9). On Figure 13 only two scenarios with high restriction and high to moderate productivity (HP-0.05 and MP-0.05) result in anoxic deposits and TOCini up to 1 wt% in scenario HP-0.05.

thumbnail Fig. 10

Total sediment inputs (A) and fluvial discharge (B) in the model. It shows that an important decrease of both sediment inputs and fluvial discharge occurs at the end of sequence A.

4.3.3 Type 3 basin

On Figure 14 only two scenarios show evidences of anoxia (HP and MP-0.05) and two others show local anoxia and widespread dysoxia (HP-015 and LP-0.05). In type 3 basin (Fig. 14), the range of organic richness is larger than in type 1 and 2 basins (Figs. 12 and 13). Here, the maximum TOCini reaches 10 wt% in the HP-0.05 scenario. Four scenarios (MP-0.05, HP-0.15, MP-0.15 and HP-1) result in significant organic deposits (> 2 wt% up to 10 wt%) and three only present limited organic accumulations (MP and LP-1 and LP-0.15). In all the scenarios, organic rich areas are located along the Western margin of the basin. Among the 9 scenarios only two present the organic distribution and oxygen condition consistent with our definition of a type 3 basin: HP-0.05 and MP-0.05.

thumbnail Fig. 13

TOCini and benthic oxygen maps from the stratigraphic modeling of type 2 basin. On this figure, only the scenarios with restriction present anoxic conditions. Here only highly restricted basin with high productivity present significant organic accumulation. HP: high productivity; MP: medium productivity; LP: low productivity; MD: medium diffusion.

5 Discussions

5.1 Regional controls on primary productivity, dilution and preservation

5.1.1 Evolution of sediment inputs

Dilution of organic matter by non-organic particles is considered to be one of the main controlling factors of source rock accumulation (Bohacs et al., 2005; Katz, 2005; Tyson, 2001). In the present study, the use of DionisosFlow provides a way to investigate the sediment supplies variations through time at basin scale (Hawie et al., 2015; Leroux et al., 2015). The relations between water discharge, sedimentation load, drainage area and topography of modern rivers has been extensively studied (Dai and Trenberth, 2002; Milliman and Meade, 1983; Milliman and Syvitski, 1992; Syvitski et al., 2003; Syvitski et al., 2005).

The Figure 15 illustrates the relation between the drainage area and the water discharge from more than 160 rivers (Dai and Trenberth, 2002). On this figure is plotted the average water discharge of the sediment sources used in the best-fit model. As no major climatic changes from dry to humid is known to occur during the Lower and Middle Triassic on Western Canada (Davies, 1997a; Davies, 1997b; Golonka et al., 1994; Hallam, 1985; Sellwood and Valdes, 2006), the variations of water discharge are most likely linked to change of the drainage area. During sequence A, the increase of water discharge and sediment supply from the North, East and South (Fig. 10) suggests an expansion of the drainage areas from the continent (Fig. 15) (Golding et al., 2015a). At the transition between sequence A and B, a sharp drop of water discharge and sediment supply suggest a significant reduction of the continental drainage area to the east (Milliman and Syvitski, 1992). However, the drop of water discharge (divided by 370) is of larger magnitude than the drop of sediment supply (divided by 15), resulting in a higher sediment concentration in rivers during sequence B. This may be explained by a topographic rise of this area, resulting in a steeper topography and higher erosion rates of the drainage area during sequence B. Concomitantly, the exhumation of sediment sources from the proto Canadian Cordillera to the west, further indicates significant paleogeographic and geodynamic transformations at that time.

Edwards et al. (1994) show the occurrence of Triassic strata in two basins separated by the Sweetgrass arch (Kent and Christopher, 1994): (1) the Alberta basin and (2) the Williston basin. In the Alberta basin, recent publications suggest that the Triassic deposition took place in a fore-arc basin (Beranek and Mortensen, 2011; Rohais et al., 2016). The DionisosFlow model presented in this study suggests that during the sequence A (Lower Triassic) a significant part of the sediment supply derived from the continent, whereas during the sequence B (Middle Triassic) sediments came from the proto Canadian Cordillera and from the Canadian continent, but with a narrower and steeper drainage area. In the Williston basin, detailed stratigraphy (Cumming, 1956) shows that the continental Lower Watrous Formation was deposited during the Middle and Upper Triassic and rest upon a major unconformity (Carlson, 1968; Fowler and Nisbet, 1985; Kent, 1994). Here, the deposition of continental strata on top of an unconformity highlights a major change of the regional subsidence rates.

In our case study, the geodynamic evolution seems to have a significant impact on the sources and sinks of Lower and Middle Triassic strata. During sequence A, the Triassic strata are sourced by a drainage area that may extend eastward the Canadian Shield (Carlson, 1968; Kent, 1994). During sequence B due to the evolution of the proto Canadian Cordillera (Rohais et al., 2016), the basin is most-likely sourced by the cordillera and to a lesser extent by the continent. Here, the reduction of the drainage is interpreted to be linked to the rise of a structural high that will split a widespread drainage area into two catchment areas: one flowing to the Alberta basin and a second flowing to the Williston basin. The rise of this topographic high will increase the erosion and induce a drop in the sediment supply and water discharge due to the reduced drainage area. In the same time, the geodynamic evolution of the proto Canadian Cordillera affected the subsidence rates in Western Canada by: (1) reducing the subsidence rate in the Alberta Basin, (2) provoking the rise of a relative structural high between the Alberta basin and the Williston basin, that reduced the sediment and water delivery to the Alberta basin and (3) increasing the subsidence in the Williston basin, allowing the deposition of continental red beds.

Huc et al. (2005) highlighted the first order control of the geodynamic settings on source rock accumulation because of their link to volcanic activities increases that affect primary productivity. The present study shows another important control of the geodynamic and physiographic evolution on the sedimentary delivery in the basin at the scale of second order stratigraphic cycles:

  • in sequence A, the extensive drainage area resulted in high sedimentation rates (> 150 m/Ma) that caused a high dilution of the organic matter (Fig. 13). In anoxic environment, this dilution results in a decrease of the organic concentration whereas in oxic environment it significantly increases the preservation (Tyson, 2001);

  • in sequence B, the reduced drainage area resulted in low sedimentation rates (< 50 m/Ma). In this scenario, anoxia is required to generate significant organic matter accumulation (Fig. 14).

thumbnail Fig. 15

Relationship showing an increase of the water discharge with drainage area (data from Dai and Trenberth, 2002). This figure shows a major decrease of the drainage area of the continental sources at the end of Lower Triassic.

5.1.2 Evolution of primary productivity and preservation of the organic matter in the Montney and Doig Formations

In sedimentary basins, anoxia can be due to high productivity (e.g. the Peruvian margin or the Namibian margin), to a restriction of the basin (e.g. the Black Sea) or to a combination of both (Arthur and Sageman, 1994; Brumsack, 2006; Sageman and Lyons, 2003). In the Montney and Doig Formations, three types of basin are present (Crombez et al., 2016b):

  • an oxic basin;

  • a restricted basin;

  • a high productivity basin.

In these three types of basin, the organic matter distribution is not the same:

  • in oxic basin, little organic matter is preserved (TOCini < 3 wt%);

  • in restricted basin, moderate accumulation of organic matter is preserved in the central part of the basin (TOCini < 5 wt%);

  • in high productivity basin, large accumulation of organic matter occurs on the edge of the basin (TOCini < 14 wt%).

The sensitivity tests performed with the new DORS model (Figs. 1214) allows for a quantification of organic productivity in each basin types:

  • type 1: the results of the modeling show that even a moderate restriction induces dysoxia (Fig. 12). The results also show that a productivity of 150–75 gC/m2/y is too low to generate 1.5 wt.% of TOCini, and that a productivity of at least 225–125 gC/m2/y associated with basin restriction is needed to decrease the oxygen content in the bottom water (Fig. 12). In order to obtain the correct organic content described by previous studies, the primary productivity must be increased to 250–125 gC/m2/y which is a moderate rate for coastal areas (Saliot, 1994) and is not high enough to induce dysoxia;

  • type 2: the modeling results show that for primary productivity higher than 150–75 gC/m2/y associated with strong and moderate restriction, the basin becomes anoxic (Fig. 13). As for type 1, the primary productivity in our simulations was underestimated. An increase to 250–125 gC/m2/y is therefore necessary to fit with the organic content in type 2 basin;

  • type 3: the results of the stratigraphic modeling (Fig. 14) show that even the highest primary productivity (225–112.5 gC/m2/y) does not induced dysoxia if not associated with strong restriction. Crombez et al. (2016b) suggested the possibility for a restriction in type 3 basin, which appears to be necessary. The results of the modeling show that if a moderate restriction with a production of 225–112.5gC/m2/y generates anoxia locally, it mostly results in dysoxia.

As a perspective, additional scenarios could be tested: in type 1 and 2 basins primary productivity values were underestimated in our simulations and in type 3 basin higher primary production scenarios could to be used (400–200 and 350–175 gC/m2/y). If these new scenarios do not induce anoxia in the basin, then scenarios with strong restriction will be required for type 3 basin. Crombez et al. (2016b) did not discuss a possible increase of the primary productivity between sequence 1, 2 and 3. The present study emphasizes the increase of sediment supply from sequence 1 to sequence 3. If the sediment sources remain the same throughout sequence A, the increase of the sediment delivery will therefore induce an increase of the total volume of nutrient delivered. This increase of nutrient delivery may induce increase of primary productivity from the sequence 1 to sequence 3. Moreover, Figure 14 shows that with a primary productivity of 225–112.5 gC/m2/y, TOCini values are up to 10 wt.% which is close to the value published in previous works. Therefore, the increase of the primary productivity to 400–200 gC/m2/y might be too high inducing an over-accumulation of organic matter.

thumbnail Fig. 12

TOCini and benthic oxygen maps from the stratigraphic modeling of type 1 basin. On this figure, only the scenarios with no restriction present oxic conditions, they also show low organic accumulation. HP: high productivity; MP: medium productivity; LP: low productivity; MD: medium diffusion.

thumbnail Fig. 14

TOCini and benthic oxygen maps from the stratigraphic modeling of type 3 basin. Note that anoxia cannot be induced by high primary productivity alone. HP: high productivity; MP: medium productivity; LP: low productivity; MD: medium diffusion.

5.2 Impact of the geodynamic evolution on the organic distribution in sedimentary basin

A previous study on the stratigraphic architecture of the Montney and Doig Formations show that this interval is composed of two second order sequences (Crombez, 2016). In these two sequences, organic rich intervals are not located in the same systems tracts: in sequence A organic rich layers are located in the FSST, whereas in the sequence B they are located in the TST.

5.2.1 Dilution variation

In the present work, the stratigraphic modeling emphasized a major drop of the fluvial discharge and sediment inputs at the boundary between sequence A and sequence B. This reduction is interpreted to be linked to a reduction of the drainage area caused by the relative rise of an arch on the continent. In the Lower and Middle Triassic strata of the Alberta basin, it is apparent that the regional geodynamic evolution of Western Canada and particularly the evolution of the Proto Canadian Cordillera (Beranek and Mortensen, 2011; Rohais et al., 2016) controls the dilution rates in sequence A and sequence B. Based on our simulations, we can estimate that during sequence A the average sediment fluxes was 23.103 km3/Ma whereas during sequence B it was only 5.103 km3/Ma. This important drop of detrital input resulted in drastically different dilutions between the two second order sequences: organic particles in sequence A were 4.5 times more diluted than in sequence B.

5.2.2 Primary productivity variation

Crombez et al. (2016b) concluded that in sequence A primary productivity is controlled by continental nutrients delivery whereas in sequence B, the high productivity of the TST is linked to marine nutrients brought through upwelling. In this study the primary productivity profile are assumed to have similar trends but different values. The present work highlighted an increase of the sediment fluxes from sequence 1 to sequence 3. In sequence A, as the primary productivity is controlled by continental nutrient inputs (Crombez et al., 2016b) the sediment fluxes increase probably induces a small primary productivity increase. In this interval, the potential moderate increase of the primary productivity is therefore connected to the regional geodynamical evolution of the area that controlled the size of the drainage area. In sequence B, the increase of nutrient inputs is interpreted to be related to upwelling activities (Crombez et al., 2016b). The sequence boundary between sequence A and sequence B is interpreted to be associated with a major rearrangement of the basin paleogeography (Rohais et al., 2016). Consequently, the high primary productivity in TST4 is linked to the regional geodynamic evolution that connected the basin to upwelling cells. In sequence B, the fact that high productivity only occurs in the TST suggests that during the geodynamically induced transgression, adequate water exchange made it possible to deliver the deep marine nutrient in the basin. Wignall and Newton (2001) proposed a model of “transgressive nearshore black shale” where high primary productivity is controlled by transgressive trends. The controls on the organic accumulation in the TST4 of our model are probably more complex as it encompasses longer time scale and thicker deposits.

5.2.3 Dynamics of anoxia

In the present case study we propose two different models for generating anoxia. In the FSST of sequence A, the anoxia is interpreted to be triggered by a major basin restriction (Crombez et al., 2016b) whereas this study shows that the anoxia in the TST of sequence B is linked to both primary productivity and a moderate basin restriction. In both models, the development of anoxia is associated to the regional geodynamic evolution of western Canada:

  • the major regression at the end of sequence A is induced by a regional tectonic uplift;

  • the physiography of the basin during sequence A and sequence B is linked to the geodynamical contex;

  • the primary productivity that generate anoxia in the TST of sequence B is linked the geodynamic evolution between sequence A and B.

This interpretation contrasts with the model of Algeo et al. (2008) that shows a major control of the climate and water delivery on anoxia, but is more similar to the model proposed by Tyson et al. (1979) and Wignall (1991) where anoxia is controlled by the basin physiography.

6 Conclusions

Based on previous works on the Montney and Doig Formations and on a 3D stratigraphic forward modeling this study shows that:

  • a major reduction of the continental sources drainage area occurs between sequence A and B. It induces a major drop in the sediment inputs during the deposition of the sequence B. This shift is interpreted to result from a major geodynamic change that occurs in western Canada at the Lower/Middle Triassic boundary (Rohais et al., 2016). During sequence A, the increasing sediment supply results in higher nutrient delivery and organic matter dilution;

  • an increase of the basin restriction is needed in type 3 basin (with high primary productivity) to trigger anoxia. Here both restrictions are induced by the basin physiography that is linked to the geodynamic evolution of the area;

  • the estimated productivity profile is 250–125 gC/m2/y in type 1 (oxic basin) and 2 (restricted basin) and 400–200 gC/m2/y in the type 3 basin (Doig phosphates zone). These values are respectively characteristic of passive margin and upwelling areas. Here additional scenarios are required to test the accuracy of these values;

  • in this case study, the distribution of the organic rich layer is controlled by the regional geodynamical evolution: it controls the sediments fluxes variations, the basin restriction and the nutrient delivery to the basin. In the present study, we show that geodynamic have an impact on marine (upwelling) and continental (runoff) nutrient delivery. These conclusions are concordant with Trabucho-Alexandre et al. (2012) but drastically differ from the transgressive model of Slatt and Rodriguez (2012). In the present study, as no climatic changes occur during the deposition of the Lower and Middle Triassic strata, geodynamic framework is the first control on source rock development, the transgressive and regressive trends only play a second role.

Acknowledgements

Acknowledgments go to IFP Energies nouvelles and the ISTeP for funding this project and to all those who have contributed to it: IFP Technologies (Canada) Inc; and Alberta Energy regulator and the Geological Survey of Canada, BC Oil and Gas Commission for the access to core data and samples. Authors also would like to thank the anonymous reviewer for his helpful comments.

References

  • Algeo TJ, Heckel PH, Maynard JB, Blakey R, Rowe H. 2008. Modern and ancient epeiric seas and the super-estuarine circulation model of marine anoxia. In: Pratt BR, Holmden C, eds. Dynamics of Epeiric seas: sedimentological, paleontological and geochemical perspectives. St. John's, Canada: Geological Association of Canada, Special Publication, vol. 48, pp. 7–38. [Google Scholar]
  • Alldredge AL, Silver MW. 1988. Characteristics, dynamics and significance of marine snow. Progress in Oceanography 20: 41–82. DOI: 10.1016/0079-6611(88)90053-5. [CrossRef] [Google Scholar]
  • Alzaga-Ruiz H, Granjeon D, Lopez M, Seranne M, Roure F. 2009. Gravitational collapse and Neogene sediment transfer across the western margin of the Gulf of Mexico: Insights from numerical models. Tectonophysics 470: 21–41. DOI: 10.1016/j.tecto.2008.06.017. [CrossRef] [Google Scholar]
  • Armitage JH. 1962. Triassic Oil and Gas Occurrences in Northeastern British Columbia, Canada. Bulletin of Canadian Petroleum Geology 10: 35–56. [Google Scholar]
  • Arthur MA, Sageman BB. 1994. Marine Black Shales: Depositional Mechanisms and Environments of Ancient Deposits. Annual Review of Earth and Planetary Sciences 22: 499–551. DOI: 10.1146/annurev.ea.22.050194.002435. [CrossRef] [Google Scholar]
  • Beranek LP, Mortensen JK. 2011. The timing and provenance record of the Late Permian Klondike orogeny in northwestern Canada and arc-continent collision along western North America. Tectonics 30: 1–23. DOI: 10.1029/2010TC002849. [CrossRef] [Google Scholar]
  • Bohacs KM, Carroll AR, Mankiewicz PJ, Miskell-Gerhardt KJ, Schwalbach JONR, Wegner MB, et al. 2005. Production, destruction, and dilution-the many paths to source-rock development. In: Harris NB, ed. The Deposition Of Organic-Carbon-Rich Sediments: Models, Mechanisms and Consequences. Tulsa, USA: Special Publications of SEPM, vol. 82, pp. 61–101. [CrossRef] [Google Scholar]
  • Brumsack H-J. 2006. The trace metal content of recent organic carbon-rich sediments: Implications for Cretaceous black shale formation. Palaeogeography, Palaeoclimatology, Palaeoecology, 232: 344–361. DOI: 10.1016/j.palaeo.2005.05.011. [CrossRef] [Google Scholar]
  • Bruns B, Littke R, Gasparik M, Van Wees J-D, Nelskamp S. 2015. Thermal evolution and shale gas potential estimation of the Wealden and Posidonia Shale in NW-Germany and the Netherlands: a 3D basin modelling study. Basin Research 28: 2–33. DOI: 10.1111/bre.12096. [CrossRef] [Google Scholar]
  • Burdige DJ. 2007. Preservation of organic matter in marine sediments: controls, mechanisms, and an imbalance in sediment organic carbon budgets ? Chemical reviews 107: 467–85. DOI: 10.1021/cr050347q. [CrossRef] [Google Scholar]
  • Burgess PM, Lammers H, Van Oosterhout C, Granjeon D. 2006. Multivariate sequence stratigraphy: Tackling complexity and uncertainty with stratigraphic forward modeling, multiple scenarios, and conditional frequency maps. AAPG bulletin 90: 1883–1901. [CrossRef] [Google Scholar]
  • Carlson CE. 1968. Triassic-Jurassic of Alberta, Saskatchewan, Manitoba, Montana, and North Dakota. AAPG Bulletin 52: 1969–1983. [Google Scholar]
  • Catuneanu O, Abreu V, Bhattacharya JP, Blum MD, Dalrymple RW, Eriksson PG, et al. (2009). Towards the standardization of sequence stratigraphy. Earth-Science Reviews 92: 1–33. DOI: 10.1016/j.earscirev.2008.10.003. [CrossRef] [Google Scholar]
  • Chalmers GRL, Bustin RM. 2012. Geological evaluation of Halfway–Doig–Montney hybrid gas shale–tight gas reservoir, northeastern British Columbia. Marine and Petroleum Geology 38: 53–72. DOI: 10.1016/j.marpetgeo.2012.08.004. [CrossRef] [Google Scholar]
  • Chauveau B, Granjeon D, Huc AY. 2013. Depositional Model of Marine Organic Matter Coupled With a Stratigraphic Forward Numerical Model (DIONISOS): Application to the Devonian Marcellus Formation (Abstr.). In: AAPG Hedberg Conference Petroleum Systems: Modeling the Past, Planning the Future abstracts. American Association of Petroleum Geologists. [Google Scholar]
  • Chauveau B, Granjeon D, Ducros M. 2016. 3D numerical stratigraphic model for basin scale modeling of the organic matter deposition in a marine environment: Application to the Natih formation (Late Cretaceous, Oman) (Abstr.). In: AAPG Geosciences Technology Workshop abstracts. American Association of Petroleum Geologists. [Google Scholar]
  • Cohen KM, Finney SC, Gibbard PL, Fan J-X. 2013. The ICS International Chronostratigraphic Chart. Episodes 36: 199–204. [Google Scholar]
  • Crombez V. 2016. Petrofacies, sedimentologie et architecture stratigraphique des roches riches en matiere organique : Étude multi-approches des Formations Montney et Doig (Trias inferieur et moyen, Alberta − Colombie Britannique, Canada). France : Thèse de doctorat de l'Université de Paris VI, 225 p. [Google Scholar]
  • Crombez V, Rohais S, Baudin F, Euzen T. 2016a. Facies well-log patterns, geometries and sequence stratigraphy of a wave dominated margin: insight from the Montney Formation (Alberta, British Columbia, Canada). Bulletin of Canadian Petroleum Geology 64 (4): 516–537. [CrossRef] [Google Scholar]
  • Crombez V, Baudin F, Rohais S, Riquier L, Euzen T, Pauthier S, et al. 2016b. Basin scale distribution of organic matter in marine fine-grained sedimentary rocks: Insight from sequence stratigraphy and multi-proxies analysis in the Montney and Doig formations. Marine and Petroleum Geology 83: 382–401. DOI: 10.1016/j.marpetgeo.2016.10.013. [CrossRef] [Google Scholar]
  • Csato I, Catuneanu O, Granjeon D. 2014. Millennial-scale sequence stratigraphy: numerical simulation with Dionisos. Journal of Sedimentary Research 84: 394–406. [CrossRef] [Google Scholar]
  • Cumming AD. 1956. The Watrous strata in Saskatchewan (Abstr.). First International Williston Basin Symposium abstracts. North Dakota Geological Society & Saskatchewan Geological Society 165–169. [EDP Sciences] [Google Scholar]
  • Dai A, Trenberth KE. 2002. Estimates of Freshwater Discharge from Continents: Latitudinal and Seasonal Variations. Journal of Hydrometeorology 3: 660–687. DOI: 10.1175/1525-7541(2002)003<0660:EOFDFC>2.0.CO;2. [CrossRef] [Google Scholar]
  • Davies GR. 1997a. The Triassic of the Western Canada Sedimentary Basin: Tectonic and Stratigraphic Framework, Paleogeography, Paleoclimate and Biota. Bulletin of Canadian Petroleum Geology 45: 434–460. [Google Scholar]
  • Davies GR. 1997b. Aeolian Sedimentation and Bypass, Triassic of Western Canada. Bulletin of Canadian Petroleum Geology 45: 624–642. [Google Scholar]
  • Davies GR, Moslow TF, Sherwin MD. 1997. The Lower Triassic Montney Formation, West-Central Alberta. Bulletin of Canadian Petroleum Geology 45: 474–505. [Google Scholar]
  • Edwards DE, Barclay JE, Gibson DW, Kvill GE, Halton E. 1994. Triassic Strata of the Western Canada Sedimentary Basin. In: Mossop GD, Shetsen I, eds. Geological Atlas of the Western Canada Sedimentary Basin. CSPG and Alberta Research Council 259–273. [Google Scholar]
  • Fowler CMR, Nisbet EG. 1985. The subsidence of the Williston Basin. Canadian Journal of Earth Sciences 22: 408–415. [CrossRef] [Google Scholar]
  • Ghadeer SG, Macquaker JHS. 2011. Sediment transport processes in an ancient mud-dominated succession: a comparison of processes operating in marine offshore settings and anoxic basinal environments. Journal of the Geological Society 168: 1121–1132. DOI: 10.1144/0016-76492010-016. [CrossRef] [Google Scholar]
  • Golding ML, Orchard MJ, Zonneveld J-P, Henderson CM, Dunn L. 2014. An exceptional record of the sedimentology and biostratigraphy of the Montney and Doig formations in British Columbia. Bulletin of Canadian Petroleum Geology 62: 157–176. [Google Scholar]
  • Golding ML, Mortensen JK, Ferri F, Zonneveld J-P, Orchard M. 2015a. Determining the provenance of Triassic sedimentary rocks in northeastern British Columbia and western Alberta using detrital zircon geochronology, with implications for regional tectonics. Canadian Journal of Earth Sciences 53: 140–155. DOI: 10.1139/cjes-2015-0082. [CrossRef] [Google Scholar]
  • Golding ML, Orchard MJ, Zonneveld J-P., Wilson NSF. 2015b. Determining the age and depositional model of the Doig Phosphate Zone in northeastern British Columbia using conodont biostratigraphy. Bulletin of Canadian Petroleum Geology 63: 143–170. [CrossRef] [Google Scholar]
  • Golonka J, Ross MI, Scotese CR. 1994. Phanerozoic paleogeographic and paleoclimatic modeling maps. In: Embry AF, Glass DJ. eds. Pangea: Global Environments and Resources . CSPG Special Publications 17: 1–47. [Google Scholar]
  • Granjeon D. 2014. 3D forward modelling of the impact of sediment transport and base level cycles on continental margins and incised valleys. In: Martinius AW, Ravnas R, Howell JA, STEEL RJ, Wonham JP. eds. From Depositional Systems to Sedimentary Successions on the Norwegian Continental Margin. International Association of Sedimentologists special publication 46: 453–472. [Google Scholar]
  • Gvirtzman Z, Csato I, Granjeon D. 2014. Constraining sediment transport to deep marine basins through submarine channels: The Levant margin in the Late Cenozoic. Marine Geology 347: 12–26. [CrossRef] [Google Scholar]
  • Hallam A. 1985. A review of Mesozoic climates. Journal of the Geological Society 142: 433–445. DOI: 10.1144/gsjgs.142.3.0433. [CrossRef] [Google Scholar]
  • Hardenbol J, Thierry J, Farley M, Jacquin T, De Graciansky PC, Vail PR. 1998. Mesozoic and Cenozoic sequence chronostratigraphic framework of European Basins. In: Graciansky PC, Hardenbol J, Thierry J, Vail PR, eds. Mesozoic and Cenozoic sequence stratigraphy of European basins. Special Publications of SEPM 60: 3–13. [Google Scholar]
  • Hawie N, Deschamps R, Granjeon D, Nader FH, Gorini C, Müller C, et al. 2015. Multi-scale constraints of sediment source to sink systems in frontier basins: a forward stratigraphic modelling case study of the Levant region. Basin Research 29 (S1): 418–445. DOI: 10.1111/bre.12156. [CrossRef] [Google Scholar]
  • Huc AY, Van Buchem FSP, Colletta B. 2005. Stratigraphic Control on Source-Rock Distribution: First and Second Order Scale. In: Harris NB, ed. The Deposition Of Organic-Carbon-Rich Sediments: Models, Mechanisms and Consequences. Special Publications of SEPM 82: 225–242. [Google Scholar]
  • Hunt D, Tucker ME. 1992. Stranded parasequences and the forced regressive wedge systems tract: deposition during base-level fall − reply. Sedimentary Geology 95: 147–160. DOI: 10.1016/0037-0738(94)00123-C. [CrossRef] [Google Scholar]
  • Jarvie DM. 2012a. Shale Resource Systems for Oil and Gas: Part 1–Shale-gas Resource Systems. In: Breyer JA, ed. Shale reservoirs–Giant resources for the 21st century. AAPG Memoir 97: 69–87. DOI: 10.1306/13321446M973489. [Google Scholar]
  • Jarvie DM. 2012b. Shale Resource Systems for Oil and Gas: Part 2–Shale-oil Resource Systems. In: Breyer JA, ed. Shale reservoirs—Giant resources for the 21st century. AAPG Memoir 97: 89–119. DOI: 10.1306/13321447M973489. [Google Scholar]
  • Katz BJ. 2005. Controlling Factors on Source Rock Development–A Review of Productivity, Preservation, and Sedimentation Rate. In: Harris NB, ed. The Deposition Of Organic-Carbon-Rich Sediments: Models, Mechanisms and Consequences. Special Publications of SEPM 82: 7–16. [Google Scholar]
  • Kent DM. 1994. Paleogeographic evolution of the cratonic platform-Cambrian to Triassic. In: Mossop GD, Shetsen I, eds. Geological Atlas of the Western Canada Sedimentary Basin. Calgary: CSPG and Alberta Research Council, pp. 69–86. [Google Scholar]
  • Kent DM, Christopher JE. 1994. Geological history of the Williston Basin and Sweetgrass arch. In: Mossop GD, Shetsen I, eds. Geological Atlas of the Western Canada Sedimentary Basin. Canada: CSPG and Alberta Research Council, pp. 421–429. [Google Scholar]
  • Leroux E, Rabineau M, Aslanian D, Gorini C, Bache F, Moulin M, et al. 2015. Post-rift evolution of the Gulf of Lion margin tested by stratigraphic modelling. Bulletin de la Société Géologique de France 186: 291–308. DOI: 10.2113/gssgfbull.186.4-5.291. [CrossRef] [Google Scholar]
  • Løseth H, Wensaas L, Gading M, Duffaut K, Springer M. 2011. Can hydrocarbon source rocks be identified on seismic data? Geology 39: 1167–1170. DOI: 10.1130/G32328.1. [CrossRef] [Google Scholar]
  • Macquaker JHS, Bohacs KM. 2007. On the Accumulation of Mud. Science 318: 1734–1735. [CrossRef] [Google Scholar]
  • Macquaker JHS, Bentley SJ, Bohacs KM. 2010. Wave-enhanced sediment-gravity flows and mud dispersal across continental shelves: Reappraising sediment transport processes operating in ancient mudstone successions. Geology 38: 947–950. [CrossRef] [Google Scholar]
  • Martin JH, Knauer GA, Karl DM, Broenkow WW. 1987. VERTEX: carbon cycling in the northeast Pacific. Deep Sea Research Part A. Oceanographic Research Papers 34: 267–285. DOI: 10.1016/0198-0149(87)90086-0. [CrossRef] [Google Scholar]
  • Milliman JD, Meade RH. 1983. World-Wide Delivery of River Sediment to the Oceans. The Journal of Geology 91: 1–21. [CrossRef] [Google Scholar]
  • Milliman JD, Syvitski JPM. 1992. Geomorphic/Tectonic Control of Sediment Discharge to the Ocean: The Importance of Small Mountainous Rivers. The Journal of Geology 100: 525–544. [CrossRef] [Google Scholar]
  • Moslow TF, Davies GR. 1997. Turbidite Reservoir Facies in the Lower Triassic Montney Formation, West-Central Alberta. Bulletin of Canadian Petroleum Geology 45: 507–536. [Google Scholar]
  • Orchard MJ, Zonneveld J-P. 2009. The Lower Triassic Sulphur Mountain Formation in the Wapiti Lake area: lithostratigraphy, conodont biostratigraphy, and a new biozonation for the lower Olenekian (Smithian). Candian Journal of Earth Sciences 46: 757–790. [CrossRef] [Google Scholar]
  • Passey QR, Bohacs K, Esch WL, Klimentidis R, Sinha S. 2010. From Oil-Prone Source Rock to Gas-Producing Shale Reservoir − Geologic and Petrophysical Characterization of Unconventional Shale Gas Reservoir (Abstr.). In: CPS/SPE International Oil & Gas Conference and in China abstracts. Society of Petroleum Engineers. DOI: 10.2118/131350-MS. [Google Scholar]
  • Reed JC, Wheeler JO, Tucholke BE, Stettner WR, Soller DR. 2005. Decade of North American Geology Geologic Map of North America–Perspectives and explanation. Boulder, USA: The Geological Society of America, 32 p. DOI: 10.1130/DNAG-CSMS-v1. [Google Scholar]
  • Riediger CL. 1997. Geochemistry of Potential Hydrocarbon Source Rocks of Triassic Age in the Rocky Mountain Foothills of Northeastern British Columbia and West-Central Alberta. Bulletin of Canadian Petroleum Geology 45: 719–741. [Google Scholar]
  • Riediger CL, Brooks PW, Fowler MG, Snowdon LR. 1990. Lower and Middle Triassic source rocks, thermal maturation, and oil-source rock correlations in the Peace River Embayment area, Alberta and British Columbia. Bulletin of Canadian Petroleum Geology 38: 218–235. [Google Scholar]
  • Rohais S, Crombez V, Euzen T, Baudin F. 2016. The Montney-Doig-Halfway Formations from Western Canadian Sedimentary Basin (WCSB): Passive margin, Back-Arc or Fore-Arc geodynamic setting? (Abstr.). In: Geoconvention: Optimizing Resources abstracts. Canadian Society of Petroleum Geologists, Canadian society of exploration geophysicists & Canadian Well Logging Society. [Google Scholar]
  • Sageman BB, Lyons TW. 2003. Geochemistry of fine-grained sediments and sedimentary rocks. Treatise on geochemistry 7: 115–158. [Google Scholar]
  • Saliot A. 1994. Biochimie organique marine. Oceanis 20: 197. [Google Scholar]
  • Schieber J, Southard JB. 2009. Bedload transport of mud by floccule ripples-Direct observation of ripple migration processes and their implications. Geology 37: 483–486. [CrossRef] [Google Scholar]
  • Schieber J, Southard JB, Thaisen K. 2007. Accretion of mudstone beds from migrating floccule ripples. Science 318: 1760–1763. [CrossRef] [Google Scholar]
  • Sellwood BW, Valdes PJ. 2006. Mesozoic climates: General circulation models and the rock record. Sedimentary Geology 190: 269–287. DOI: 10.1016/j.sedgeo.2006.05.013. [CrossRef] [Google Scholar]
  • Sigman DM, Hain MP. 2012. The biological productivity of the ocean. Nature Education Knowledge 3(6): 1–16. [Google Scholar]
  • Slatt RM, Rodriguez ND. 2012. Comparative sequence stratigraphy and organic geochemistry of gas shales: Commonality or coincidence ? Journal of Natural Gas Science and Engineering 8: 68–84. DOI: 10.1016/j.jngse.2012.01.008. [CrossRef] [Google Scholar]
  • Syvitski JP, Peckham SD, Hilberman R, Mulder T. 2003. Predicting the terrestrial flux of sediment to the global ocean: a planetary perspective. Sedimentary Geology 162: 5–24. DOI: 10.1016/S0037-0738(03)00232-X. [CrossRef] [Google Scholar]
  • Syvitski JPM, Vörösmarty CJ, Kettner AJ, Green P. 2005. Impact of Humans on the Flux of Terrestrial Sediment to the Global Coastal Ocean. Science 308: 376–380. [CrossRef] [Google Scholar]
  • Trabucho-Alexandre J, Hay WW, De Boer PL. 2012. Phanerozoic environments of black shale deposition and the Wilson Cycle. Solid Earth 3: 29–42. DOI: 10.5194/se-3-29-2012. [CrossRef] [EDP Sciences] [Google Scholar]
  • Tyson RV. 2001. Sedimentation rate, dilution, preservation and total organic carbon: some results of a modelling study. Organic Geochemistry 32: 333–339. DOI: 10.1016/S0146-6380(00)00161-3. [CrossRef] [EDP Sciences] [Google Scholar]
  • Tyson RV, Wilson RC, Downie C. 1979. A stratified water column environmental model for the type Kimmeridge Clay. Nature 277: 377–380. [CrossRef] [Google Scholar]
  • Wakeham S, Lee C. 1993. Production, Transport, and Alteration of Particulate Organic Matter in the Marine Water Column. In: Engel M, Macko S, eds. Organic Geochemistry. Springer US 145–169. DOI: 10.1007/978-1-4615-2890-6_6. [Google Scholar]
  • Wignall PB. 1991. Model for transgressive black shales ? Geology 19: 167–170. DOI: 10.1130/0091-7613(1991)019. [CrossRef] [Google Scholar]
  • Wignall PB, Newton R. 2001. Black shales on the basin margin: a model based on examples from the Upper Jurassic of the Boulonnais, northern France. Sedimentary Geology 144: 335–356. DOI: 10.1016/S0037-0738(01)00125-7. [CrossRef] [Google Scholar]
  • Wright GN, Mcmechan ME, Potter DEG. 1994. Structure and Architecture of the Western Canada Sedimentary Basin. In: Mossop GD, Shetsen I, eds. Geological Atlas of the Western Canada Sedimentary Basin. CSPG and Alberta Research Council 25–40. [Google Scholar]
  • Zonneveld J-P., Macnaughton RB, Utting J, Beatty TW, Pemberton SG, Henderson CM. 2010. Sedimentology and Ichnology of the Lower Triassic Montney Formation in the Pedigree-Ring/Border-Kahntah River Area, Northwestern Alberta and Northeastern British Columbia. Bulletin of Canadian Petroleum Geology 58: 115–140. [CrossRef] [Google Scholar]

Cite this article as: Crombez V, Rohais S, Baudin F, Chauveau B, Euzen T, Granjeon D. 2017. Controlling factors on source rock development: implications from 3D stratigraphic modeling of Triassic deposits in the Western Canada Sedimentary Basin, Bull. Soc. géol. Fr. 188: 30.

All Tables

Table 1

Diffusion coefficients used in DionisosFlow.

Table 2

Scenarios tested in order to study the impact of organic matter diffusion on organic accumulations.

Table 3

Scenarios tested in order to study the organic accumulation in the different basin type.

All Figures

thumbnail Fig. 1

A. Location of the study area (faults from Reed et al., 2005, Triassic subcrop edges from Edwards et al., 1994; basin limits from Wright et al., 1994). B. Paleogeographic map during the lower Triassic over the Western Canada sedimentary basin. Qs: sedimentary fluxes; Ab.: Alberta; B.C.: British Columbia; Y.: Yukon; N.T.: Northern Territories.

In the text
thumbnail Fig. 2

Stratigraphic framework of the Lower and Middle Triassic of the Western Canada Sedimentary Basin (modified from Crombez, 2016, eustatic variations from Hardenbol et al., 1998).

In the text
thumbnail Fig. 3

Eustatic variations used in the stratigraphic modeling. The composite curve is based on Hardenbol et al. (1998) and on two user defined relative sea-level curves that represent the impact of regional uplifts and subsidence.

In the text
thumbnail Fig. 4

Location of the sedimentary inputs of the stratigraphic model. Note the turnover between sequence A and sequence B.

In the text
thumbnail Fig. 5

Cross sections in the best fit model. Here the color is related to an estimated gamma ray based on the lithology proportion. This figure highlights the typical geometries of the Montney Formation (A) and the Doig Formation (B).

In the text
thumbnail Fig. 6

Thickness and errors maps for each modeled sequence. It is apparent that relative errors are low (< 30%) and therefore acceptable.

In the text
thumbnail Fig. 7

Comparison between synthetic wells extracted from the best fit model and real well from the same location in the Western Canada Sedimentary Basin. Note the thickness similarities between numerical simulated data and real data.

In the text
thumbnail Fig. 8

Organic content maps extracted from the model in order to test the influence of the diffusion coefficient on the organic matter distribution in the different basin types. On this figure a weak concentration in the deepest part of the basin in apparent. HP: high productivity; MP: medium productivity; LP: low productivity; HD: high diffusion; MD: medium diffusion; LD: low diffusion.

In the text
thumbnail Fig. 9

TOCini, benthic oxygen and bathymetry profile in the three basin types. This figure shows a significant enrichment in organic matter in the deepest part of the basin with the increase of the diffusion coefficient. Note that this concentration of the organic matter has no impact on the benthic oxygen profile.

In the text
thumbnail Fig. 11

Bathymetric maps of the interval representative of type 1 (A), type 2 (B) and type 3 (C) basins. On this figure, it is apparent that the bathymetry in type 3 basin is shallower than in type 1 and 2.

In the text
thumbnail Fig. 10

Total sediment inputs (A) and fluvial discharge (B) in the model. It shows that an important decrease of both sediment inputs and fluvial discharge occurs at the end of sequence A.

In the text
thumbnail Fig. 13

TOCini and benthic oxygen maps from the stratigraphic modeling of type 2 basin. On this figure, only the scenarios with restriction present anoxic conditions. Here only highly restricted basin with high productivity present significant organic accumulation. HP: high productivity; MP: medium productivity; LP: low productivity; MD: medium diffusion.

In the text
thumbnail Fig. 15

Relationship showing an increase of the water discharge with drainage area (data from Dai and Trenberth, 2002). This figure shows a major decrease of the drainage area of the continental sources at the end of Lower Triassic.

In the text
thumbnail Fig. 12

TOCini and benthic oxygen maps from the stratigraphic modeling of type 1 basin. On this figure, only the scenarios with no restriction present oxic conditions, they also show low organic accumulation. HP: high productivity; MP: medium productivity; LP: low productivity; MD: medium diffusion.

In the text
thumbnail Fig. 14

TOCini and benthic oxygen maps from the stratigraphic modeling of type 3 basin. Note that anoxia cannot be induced by high primary productivity alone. HP: high productivity; MP: medium productivity; LP: low productivity; MD: medium diffusion.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.