Detection of Late Cretaceous eustatic signatures using quantitative biostratigraphy James S. Crampton† GNS Science, P.O. Box 30-368, Lower Hutt, New Zealand
Poul Schiøler The Geological Survey of Denmark and Greenland, Ø. Voldgade 10, DK-1350, Copenhagen K., Denmark
Lucia Roncaglia Danish Oil and Natural Gas A/S, Agern Alle 24-26, DK-2970 Hørsholm, Denmark ABSTRACT This study tested the potential for highresolution quantitative biostratigraphic methods to distinguish tectonic from supposedly eustatic sequences. The study employed automated multidimensional graphic correlation using constrained optimization (program CONOP). Our analysis was based on 15 onshore sections in eastern New Zealand spanning the Coniacian to Maastrichtian stages (Upper Cretaceous), and 398 wellconstrained lowest- or highest-occurrence records derived from 245 dinoflagellate species. The age of the resulting composite section was calibrated using six dated bioevents; two of these events were tied to geochemical or paleomagnetic datums. Following calibration, the composite section had an average relative temporal resolution of less than 130 k.y. through the study interval, although the succession of events was not dated to this level of precision in an absolute sense. Although the CONOP composite will not resolve unconformities or condensed intervals that occur across all studied sections, such “global” unconformities can be identified using clusters of events that are placed at a single composite level. Fifteen major event clusters were identified, 10 of which coincided closely with sequence boundaries identified in the Northern Hemisphere; the probability of this coincidence arising by chance alone is ~8%. Although these results must be regarded as provisional, they suggest that high-resolution quantitative stratigraphy can provide a potent tool for the identification and correlation of stratigraphic sequences. Furthermore, our findings add significant †
E-mail:
[email protected].
support for recent studies that have argued for the presence of Late Cretaceous eustasy. Keywords: biostratigraphy, correlation, eustasy, Late Cretaceous, New Zealand, quantitative analysis, sequence stratigraphy, unconformities. INTRODUCTION Distinguishing tectonic from eustatic signatures within stratigraphic successions remains problematic (e.g., Pedley and Frostick, 1999, and associated papers). This issue is of particular interest for the Upper Cretaceous system, for which the evidence of high-frequency, largemagnitude global sea-level changes, and the implied growth and decay of continental-scale ice sheets, remains equivocal (e.g., Abreu et al., 1998; Barrera and Savin, 1999; Miller et al., 1999, 2003, 2004, 2005a, 2005b; Price 1999; Stoll and Schrag, 2000; Huber et al., 2002). In part, the problem reduces to one of resolution in biostratigraphic and sequence stratigraphic correlation: arguments for eustasy require demonstration of global synchroneity of third- and higher-order sequences (for relevant discussions of this problem, see Miall, 1992; Miall and Miall, 2004). This can be extremely difficult using conventional correlation tools, especially in regions of stratigraphic or structural complexity. In this paper, we employ the newly developed quantitative stratigraphic approach of automated multidimensional graphic correlation using constrained optimization (CONOP; Kemple et al., 1995) to derive a very high-resolution biostratigraphic framework for lithologically monotonous, but structurally and stratigraphically complex, Coniacian to Maastrichtian successions in eastern New Zealand. Using this framework, we differentiate numerous local unconformities from those of regional extent and develop an age
model that allows correlation of our regional unconformities with sequence boundaries recognized in the Northern Hemisphere. Although the results presented here must be considered provisional, they seem to provide a global-scale test of hypothesized Late Cretaceous eustasy. During the Late Cretaceous, New Zealand was positioned at ~70°S (Sutherland et al., 2001). Coniacian to Maastrichtian strata of eastern New Zealand are inferred generally to have been deposited during thermal subsidence on the passive paleo-Pacific margin of New Zealand, following its separation from Gondwanaland (e.g., Carter, 1988; Ballance, 1993; King et al., 1999). The studied successions lie close to the base of the first-order transgressive megasequence that characterizes New Zealand Late Cretaceous and Paleogene stratigraphy. Cretaceous subduction in eastern New Zealand is likely to have ceased by about Cenomanian times, although this interpretation remains contentious (Field et al., 1997, p. 121, and references therein). Regardless of this, Coniacian to Maastrichtian formations are dominated by fine-grained siliciclastics or pelagic carbonates and show little conspicuous evidence for significant syndepositional tectonism. The present study was focused on strata immediately overlying a regional onlap surface that started forming at the same time as initiation of Tasman Sea spreading (MacPherson, 1948; King et al., 1999). In stratigraphic updip locations, this surface has been called the Waipounamu erosion surface, and it is strongly diachronous, recording perhaps 50 m.y. of fluvial erosion, marine transgression, and planation (LeMasurier and Landis, 1996). Here, we concentrate on fully marine strata deposited basinward; it is these that preserve the most complete sedimentary record of geological events during the Coniacian to Maastrichtian stages.
GSA Bulletin; July/August 2006; v. 118; no. 7/8; p. 975–990; doi: 10.1130/B25826.1; 9 figures; Data Repository item 2006121.
For permission to copy, contact
[email protected] © 2006 Geological Society of America
975
Crampton et al. Fifteen outcrop sections were studied in detail (Fig. 1). Three of these lie within the northern part of the Canterbury basin (Field et al., 1989); the rest are distributed throughout the onshore extent of the East Coast basin (Field et al., 1997). Sections in the East Coast basin were selected, in part, to provide a series of three transects perpendicular to the paleo–continental margin, although it is possible that the original relationships have been disrupted substantially by Neogene deformation (see following). Sections from the Canterbury basin were included in an effort to resolve possible basinwide tectonic events in the East Coast basin. Because calcareous macro- and microfaunas are typically depauperate, poorly preserved, and/or of low diversity in the studied formations (e.g., Crampton et al., 2000), the present biostratigraphic analysis is based instead on diverse, abundant and well-known dinoflagellate floras (Wilson, 1984; Schiøler and Wilson, 1998; Roncaglia et al., 1999; Schiøler et al., 2000, 2001, 2002). Throughout this paper, in order to facilitate comparisons with previous studies, we have calibrated our biostratigraphic events to the global geochronological scale of Gradstein et al. (1995). We note that this scale differs somewhat from the recently proposed time scale of Gradstein et al. (2004), although these differences are minor in the studied interval (6 and/or mean penalties >50 (Fig. DR3A). These events were given weights of 0 in the input file used to generate the final composite section. Zero-weighted events do not influence the coexistence constraint in CONOP 9, and they do not contribute to the penalty score. In other words, they were effectively eliminated from the analysis. These events are indicated in the composite section (Table DR2, see footnote 1), but are omitted from all other plots and outputs based on the composite. RESULTS The CONOP Composite Section The CONOP composite section is presented in Figures DR4–DR17 and Table DR2 (see footnote 1). Because the individual events in the composite are too numerous to label and for ease of comparison between logs, the vertical scale of the composite is expressed as percentage units above the base. This convention deliberately separates the task of integrating biostratigraphic and field data from the explicit inference of an age model for the composite section (undertaken later). In the composite, the 398 reliable or well-constrained events are distributed across 182 event levels. In Figures DR4–DR17 (see footnote 1), the composite is shown against each of the studied sections; one example section, Ben More, is shown in Figure 4. Note that the “event clusters” marked on these plots are explained below. Typically, 25%–50% of expected taxa are present at each sampled horizon in each section. In Figure 4 and the plots of Figures DR4–DR17, convergence of many composite event levels to a single horizon in a given section, and corresponding steps in the event-depth plot, can indicate either true relative condensation of the section at that point and/or apparent condensation due to wide sample spacing or structural disruption. Bundling of composite events due to wide sample spacing results because CONOP locates events only at actually sampled horizons (see Methods)—wide sample spacing in a particular section relative to the composite inevitably causes several composite events to be clustered at sample horizons on either side of a sample gap. With these caveats in mind, all significant composite event convergences and bundles in each section have been examined against all available evidence to determine whether they mark true intervals of sedimentary condensation or are artifacts of sampling, faulting, or the method. Disconformities, unconformities, and condensed intervals indicated in the central
columns of Figures 4 and DR4–DR17 are those inferred to be real. From these results, it can be seen that CONOP has correctly identified all unconformities or disconformities recorded in the field. In addition, it has detected a large number of other possible disconformities or highly condensed intervals, many of which are supported by field observations. Using the Ben More section as an example (Fig. 4), later in this paper we cite corroborating field and other evidence for the unconformities and condensed intervals that have been identified using the CONOP composite. Detailed interpretation of all sections and an explanation of inferred depositional patterns, however, will be undertaken elsewhere. Age Interpretation To a first approximation, the spacing of events in the CONOP composite section is likely to be proportional to their true spacing in time (Sadler and Cooper, 2004). In order to correlate precisely the composite with the chronometric time scale, it is necessary to calibrate the composite using events of known age. Sadly, there are few suitable events in the Coniacian to Maastrichtian interval in New Zealand, and calibration to the local and international time scales is rather problematic (Crampton et al., 2000; Cooper, 2004). Events used here to correlate the composite section to the international time scale of Gradstein et al. (1995) are listed next in ascending stratigraphic order. These correlations are based either on direct geochemical or paleomagnetic ties in New Zealand, or on other Southern Hemisphere mid- to high-latitude bioevents that have themselves been correlated using a range of paleontological and geochemical criteria (refer to the cited publications). Only taxa that are believed to be well known taxonomically are used; some events, such as the lowest occurrence of Isabelidinium cretaceum s.s., are excluded because of uncertainties regarding consistency of identification and taxonomic concept. 1. The lowest occurrence (LO) of the dinoflagellate Conosphaeridium striatoconum, which is inferred to lie close to the base of the Coniacian stage in Australia (Young and Laurie, 1996) and in the Northern Hemisphere (Williams et al., 2004). In New Zealand, this event apparently lies low in the local Teratan stage (G.J. Wilson, March 1995, personal commun.; Crampton et al., 2001) and just above the base of the Coniacian according to the correlation of Cooper (2004). This event has been identified in three of the studied sections. Herein, we adopt an age of 88.5 Ma for the LO of C. striatoconum.
2. The LO of the dinoflagellate Nelsoniella aceras, which is known to lie close to the C34/ C33 magnetochron boundary in New Zealand (Crampton et al., 2000). The age of this boundary is somewhat contentious (see Crampton et al., 2000, and references therein); here we assume that it is at the Santonian-Campanian boundary and adopt an age of 83.5 Ma, following Gradstein et al. (1995). The LO of Nelsoniella aceras was identified in nine of the studied sections. 3. The LO of the dinoflagellate Cerodinium diebelii, which has been inferred to lie at the base of the Upper Campanian in Antarctica and the Southern Ocean (Crame et al., 1991; Mohr and Mao, 1997; Pirrie et al., 1997) and, more recently, to lie within the middle Campanian in southern mid-latitudes (Williams et al., 2004). In Australia, there are conflicting reports of the age of this event (see references in Roncaglia et al., 1999; Young and Laurie, 1996). The LO of C. diebelii was identified in three of the studied sections. Herein, following Williams et al. (2004), we assume a mid-Campanian age of 78.5 Ma for this event. 4. The LO of the dinoflagellate Isabelidinium korojonense, which lies somewhere in the middle Campanian in Australia and Antarctica (Helby et al., 1987; Pirrie and Riding, 1988; Marshall, 1990b; Young and Laurie, 1996; Pirrie et al., 1997). This event was identified in nine of the studied sections. An age corresponding to the mid-point of the stage is adopted herein: 77.5 Ma. 5. The LO of the dinoflagellate Canninginopsis bretonica, which is inferred to lie at the base of the Upper Campanian in Antarctica, Australia, and southern mid-latitudes (Marshall, 1990a; Crame et al., 1991; Riding et al., 1992; Pirrie et al., 1997; Williams et al., 2004). This event was identified in 11 of the studied sections. An age of 76 Ma is adopted herein. 6. The highest occurrence (HO) of the dinoflagellate genus Odontochitina, which is inferred to lie at the base of the Upper Maastrichtian in Australia, Antarctica, and the Southern Ocean (Helby et al., 1987; Pirrie and Riding, 1988; Riding et al., 1992; Young and Laurie, 1996; Mohr and Mao, 1997). This event was identified in 14 of the studied sections. An age of 69.5 Ma is adopted herein. 7. The LO of the dinoflagellate Trithyrodinium evittii, which is known to lie close to the precisely located Cretaceous-Paleocene boundary in New Zealand (Wilson, 1987; Strong et al., 1995; Willumsen, 2000; Cooper, 2004). This event is known from just one of the studied sections. Its placement, however, is well constrained by numerous associated events that define a major event cluster at the CretaceousCenozoic boundary (see the following).
Geological Society of America Bulletin, July/August 2006
981
Crampton et al. Ben More
Composite section
P30/93091662 93381703 outcrop continues
Thickness (m) 260
40
250
Herring Formation
240
230
220
cm )m
(
?
210
200
190
Minor fault repeats section
170
? G P
HF
cm G cm G
160
30
?
0
G
Composite section (percentage units)
180
dm (G)
150
140
? 130 dm
0
120
110
90
80
Paton Formation
100
20
cmdm
70
0 (G)
60
50 cmdm
40 0
?
30
20
cmdm
0
10
(G)
10
0
Mst
Sst Cgl
0
20
40
60
Composite section (percentage units)
982
Geological Society of America Bulletin, July/August 2006
Figure 4. Log of the Ben More section. The figure comprises four parts: On the left is the measured section based on field data. Lithostratigraphic and graphic symbols used on the log are explained in Figure DR4 (see text footnote 1). On the right of the figure is the relevant part of the composite section derived from constrained optimization program (CONOP); “global” event clusters identified later in this paper are indicated here by black triangles against the composite scale. The two central parts of the figure show events of the CONOP composite projected into the measured section—shown both as an eventdepth plot and, to the right of that, as bundled composite event levels. Not all sampled horizons have associated events. For clarity, events that cluster at the upper- and lowermost samples—i.e., events that most probably lie outside the range of the section—have been omitted. Unconformities or condensed intervals marked in the central column are inferred from all available evidence—from field data, from significant local convergence and bundling of biostratigraphic events relative to the CONOP composite, and from “global” event clusters. Not all local event convergences are interpreted as unconformities. Instead they may mark intervals of low sedimentation rate, faults, or result simply from wide sample spacing relative to the spacing of events in the composite. Mst—mudstone; Sst— sandstone; Cgl—conglomerate.
Quantitative biostratigraphy and Cretaceous eustasy 60
7
65
6
70
Age (Ma)
The correlation between the position of these seven points in the composite section and their inferred ages is high (Pearson’s product-moment correlation coefficient, r = −0.967), but there is still significant scatter in the data (Fig. 5). In part, this scatter results from the arbitrary placement of biostratigraphic range ends at substage boundaries that, for most cited studies, represent the finest practical units of chronostratigraphic resolution. Also, of course, scatter may result from true diachroneities in the timing of particular events between regions. With these considerations in mind, we regard point number 5, the LO of Canninginopsis bretonica, to be an outlier and exclude it subsequently from the construction of the age model. Based on the remaining six points, r = −0.990. In Figure 5, a straight line is fitted to these six data points. Because the last event, the LO of T. evittii, is considered to be a precise marker of the Cretaceous-Paleocene boundary, the regression is constrained to pass through this point. The regression line defines our best estimate of an age model for the composite section. In using a straight line, the model assumes that the spacing of events in the composite is indeed proportional to their spacing in time. From the available data, there is no reason to suspect a nonlinear relationship for the data set as a whole (but see following). Based on this age model, the 182 composite event levels yield an average relative temporal resolution of ~128 k.y. If one considers just the latest Coniacian to earliest Maastrichtian interval, then the average relative temporal resolution is 102 k.y. These figures refer just to the relative resolving power of the composite section; as discussed already, precision in absolute dating is limited by comparatively large uncertainties in our age model. Using this age model, Figure 6 shows the composite against the global geochronological time scale of Gradstein et al. (1995). Also shown are the New Zealand dinoflagellate zones and subzones of Wilson (1987), Schiøler and Wilson (1998), Roncaglia et al. (1999), and Cooper (2004). The older zonal boundary index events occur in their expected order in the composite sequence, and the zones can be identified unambiguously. In contrast, some of the younger zones cannot be identified unambiguously because index events are apparently out of order in the composite. This problem may reflect any of the factors that contribute to the correlation problem, discussed already. In essence, during the conventional process of creating a zonation, complex biostratigraphic data are filtered in a subjective but intelligent way. The CONOP composite, on the other hand, is an objective ordering that follows simple mechanical rules. For a data set of this size, however, there is no
75 5 4 3 80
2 85
1 90
0
20
40
60
80
100
Composite section (percentage units)
Figure 5. Construction of the age model used to calibrate the composite section. Numbered points shown are the selected dinoflagellate event datums described in the text; all other event levels in the composite section are shown as vertical lines above the x axis. Error bars on inferred ages are conservative and based on possible age limits discussed in the text. Error bars on the height in the composite section, based on the best-fit intervals, are all negligible (i.e., 0 or smaller than the plotted symbols). Point 5, representing the lowest occurrence of Canninginopsis bretonica, is regarded as an outlier and is not used in the construction of the age model. The correlation coefficient, r, for the remaining six points is −0.990. The linear regression is constrained to go through the youngest point (point 7, LO of Trithyrodinium evittii). The ages of all other event levels are determined by interpolation through the regression line. The equation of this line is: y = −0.286x + 91.633.
single, unique best composite, but a range of equally best-fit solutions. There is, therefore, a stochastic element to the final CONOP composite sequence. For these reasons, it is unsurprising that there is some disagreement between event orders inferred during conventional zonation and those inferred using CONOP, as noted previously by Sadler and Cooper (2004). Because of this problem, in Figure 6 the bases of the Isabelidinium korojonense and I. pellucidum zones are placed at the lowest co-occurrences of three key or characteristic taxa identified in Roncaglia et al. (1999). Using the same convention, it is not possible to separate the overlying Alterbidinium acutulum and Manumiella druggii zones;
resolution of this problem is beyond the scope of the present paper. Identification of “Global” Event Clusters As noted previously, without independent age control, CONOP cannot resolve “global” changes in depositional rate that affect all sections in a sample. Thus, for example, unconformities that are present in all sections will not be resolved in the individual sections. We would expect, however, that given relatively uniform event-level spacing in time, “global” unconformities would be marked by event clusters in the composite section—the placement of
Geological Society of America Bulletin, July/August 2006
983
(part)
Chlamydophorella ambigua
Turonian
(part)
Cymososphaeridium benmorense
Conosphaeridium abbreviatum
Odontochitina porifera
Isabelidinium cretaceum
Trithyrodinium suspectum Nelsoniella aceras
Vozzhennikovia spinulosa
Isabelidinium korojonense
Isabelidinium pellucidum
Manumiella druggii/ Alterbidinium acutulum
(part)
Ev e c nt o mp level os s in ite c s ec alibr tio ate n d 10
20
30
40
50
60
70
80
90
100
uth
ara Wa
lat
ion Iso
ng
nM
ore Be
a Ko ra
na he Ow
ra go An
Ho B
irc h
llow Co
nw ay
Ha
um uri
ore
Up p Be er nM
ng
u Ke ke re
Nid
d
Ta ng
he aru
Wa it
aia
Age (Ma)
Figure 6. Inferred ages of event levels in the composite section, as derived from Figure 5, and New Zealand dinoflagellate zonal boundaries. Zone boundaries shown with a dashed line are those located using the lowest occurrences of three key or characteristic taxa (see text for further explanation). On the right are summary plots of the sections studied here; boxes indicate time that is represented in the sections, vertical lines mark time that is missing at unconformities. Black fill indicates mudstone-dominated strata; stippled fill marks intervals containing significant sandstone. Wavy lines denote unconformities that were recorded in the field (see supplementary material [see text footnote 1] for detailed section logs). Horizontal gray lines represent the positions of significant event clusters in the composite section (see Figs. 7–8).
90
85
80
75
70
(part)
Maastrichtian
Campanian
Sant.
Coniacian
Late
Early
Late
Middle
Early
L
M
65
G l ob (G al ge ra o d ste chro in et nolog a l . , 1 ical 99 sc 5) ale
L
ate
E
Geological Society of America Bulletin, July/August 2006
E
Ne w zo Zea ne la s/s nd u d b zo inof ne lag s ell
Trithyrodinium evittii
Composite section (percentage units)
idWa ip M
So ipa ra
984
Paleocene
Crampton et al.
Quantitative biostratigraphy and Cretaceous eustasy anomalously high numbers of apparently coeval events at a single composite event level. In addition, however, such clusters are also expected at intervals of anomalously high biotic turnover, either real or apparent (the latter resulting from monographic and/or taxonomic artifacts at key stratigraphic boundaries). Examination of Figure 7A reveals that there are, indeed, conspicuous event clusters in the composite. When interpreting this plot, however, it is necessary to account for event-level spacing in the composite—spacing that is determined by the spacing of sampled levels in the original sections (see Methods section). Clearly, if two event levels are widely spaced in the composite, one would expect that, on average, more events would be assigned to one or both of these levels than to closely spaced levels. Here, we have standardized for event-level spacing in two ways. First, we calculated the adjusted event clusters, E, as the number of events at a level divided by the average spacing between this level and those above and below (Fig. 7B):
E=
e1
( hl +1 − hl −1 ) / 2
.
(1)
Where el is the number of events at event-level l, and hl + 1 and hl – 1 are the heights of the adjacent event levels above and below (in percentage composite units). Secondly, we fitted a smooth curve—a sixth-order polynomial—to the cumulative distribution of events against time (Fig. 7E); the local slope of this curve represents the expected average accumulation rate of events at that point in time. From this, we calculated differences between the expected number of events at each level in the composite and the observed number of events; these differences are plotted in Figure 7C. Examination of Figures 7A–7C reveals that all the major event clusters are identified in all three plots, although the relative magnitudes of peaks vary somewhat. This correspondence suggests that the event clusters are more than just artifacts of event-level spacing in the composite. On Figure 7D, adjusted event clusters (E) are plotted in rank order on semilog axes. (Henceforth, we ignore the events below the 10% level and above 90% in the composite that may be, in part at least, edge effects or artifacts—see discussion of unreliable or poorly constrained events, and discussion of the Cretaceous-Paleocene boundary, below.) Whereas much of the ranked event-cluster distribution is log-linear, the curve below rank ≈36 is very markedly non-log-linear. This suggests that adjusted event clusters larger than E ≈10 may be qualitatively and genetically different from those smaller than
10. Thus, for the purposes of this discussion, we have adopted a value of E = 10 as the threshold between “significant” event clusters and “background” event densities. This threshold is marked on Figure 7B as a dashed line, and event clusters that exceed this threshold are indicated with arrows. These event clusters are interpreted in the next section. The largest event cluster shown on Figure 7A, 24 events at 93% in the composite section, requires special mention. According to our age model, this cluster corresponds to the Cretaceous-Paleocene boundary. It is likely to reflect four underlying causes. First, the CretaceousPaleocene boundary could mark an interval of truly enhanced microfloral turnover, although globally there is no evidence for a major turnover of dinoflagellates at this boundary (Galeotti et al., 2004). Secondly, apparent turnover may be exaggerated by monographic effects inasmuch as many studies have tended to extend up to, but not across, the boundary (e.g., Roncaglia et al., 1999; but see Willumsen, 2000). Thirdly, the Cretaceous-Paleocene boundary is known to correspond to an unconformity in many or most New Zealand sections (Hollis, 2003). Finally, event levels in the composite are relatively widely spaced around the CretaceousPaleocene boundary, and, thus, for the reasons outlined above, the size of the cluster is likely to be exaggerated (compare Figs. 7A and 7B). INTERPRETATION AND DISCUSSION The 15 major adjusted event clusters shown on Figures 6 and 8, we believe, most likely represent regional unconformities or condensed horizons that are present in all 15 of the studied sections. If this is correct, then their presence will affect the age model used here. In essence, regional unconformities will introduce nonlinearities into the simple, linear age model of Figure 5. Using the simple age model, events between each of the unconformities will be widely separated compared to their true spacing in time (Fig. 9). Discrepancies between apparent and true ages will vary in magnitude depending on the position of the event relative to the unconformities, the duration of the unconformities, their frequency and spacing in time, and the distribution of dated events used to calibrate the model. In particular, errors in the age calibration of individual events will be greatest close to the regional unconformities and smallest midway between the unconformities. In the absence of many more numerous and precise age-control datums, we are unable to recalibrate our age model to account for the presence of regional unconformities. Our age model remains, however, the best “average” model available to us.
Figure 8 shows the durations of hiatuses identified on the New Jersey coastal plain by Miller et al. (2003, 2004) and Van Sickel et al. (2004) that are inferred to mark eustatic falls. Using the age model developed here, 10 out of the 15 major event clusters fall within the durations of the New Jersey hiatuses (Fig. 8). One can estimate the probability of this correspondence occurring by chance alone by using the standard binomial probability equation:
P=
n!
( n − r )!.r !
( n−r )
× (1 − p )
× pr ,
(2)
where n is the number of trials, r is the number of successes, and p is the probability of a single success. The probability of a single success can be approximated by calculating the proportion of total age-calibrated composite event levels that lie within the age range of the hiatuses. In our case, n = 15, r = 10, p = 85⁄175 = 0.486, and P = 0.079, or ~8%. In other words, the probability that 10 of our 15 major clusters would lie within the time spans of the New Jersey hiatuses due to chance alone is ~8%. Also shown for reference on Figure 8 are the sequence boundaries identified by the Exxon Production Research team (EPR; Haq et al., 1987), as recalibrated to the time scale of Gradstein et al. (1995) by Miller et al. (2003). Miller et al. (2003) discussed correlations between their New Jersey hiatuses and the EPR sequence boundaries. In the present paper, we are concerned primarily with the relationships of our event clusters to the explicitly age-calibrated New Jersey hiatuses; the EPR sequence boundaries are not discussed here. Taken at face value, these results suggest that we can identify seven out of the eight Santonian to Maastrichtian hiatuses described by Miller et al. (2003, 2004) in the New Jersey coastal plain. We cannot exclude the possibility that this correspondence has arisen by chance alone, although the probability of this is small (~8%). The only New Jersey hiatus that is not observed in New Zealand defines the base of the Englishtown composite sequence (Fig. 8) and is in fact relatively poorly dated (K. Miller, December 2005, personal commun.). In addition, we can resolve other significant event clusters and inferred unconformities in the New Zealand sections that have not been observed in the New Jersey succession, at ca. 68.5 Ma, 72.1 Ma, 76.2 Ma, 78.2 Ma, and 78.8 Ma. Based on our preferred age model, therefore, we suggest that apparent unconformities resolved in New Zealand are approximately coeval with sequence boundaries identified in the Northern Hemisphere. This conclusion lends considerable weight to
Geological Society of America Bulletin, July/August 2006
985
Crampton et al.
A
B
C 100
65
24 16.7
90
80 70
60
75
50
80
40 98 255 45
30
85 20
10.0
10
14
90 0
2 4 6 8 Number of events clustered at composite section level
0
10 20 30 40 Adjusted event clusters, (E )
400
0 4 8 Differences between observed and expected based on polynomial fit (part E)
E
100
Rank ≈ 36, E ≈ 10 10
1
0
40
80 120 Rank order
160
Cumulative number of events in composite section
Adjusted event clusters (E )
D
-4
300
200
100
0 90
986
85
80
75 Age (Ma)
Geological Society of America Bulletin, July/August 2006
70
65
0
Composite section (percentage units)
Age (Ma)
70
Quantitative biostratigraphy and Cretaceous eustasy the suggestion that these sequence boundaries reflect eustatic changes in sea level. Given the uncertainties in our age model, we do not consider this to be the last word on the subject. In particular, the correspondence between our event clusters and the New Jersey hiatuses is sensitive to the position of the line of correlation shown on Figure 5. To test the most conservative scenario, we used a Monte Carlo approach to determine the distribution of matches between our event clusters and the New Jersey hiatuses, based on the uncertainty limits for the age-calibration points shown on Figure 5 and assuming uniform probability distributions within these limits. From 1000 Monte Carlo trials, we found that the average number of matches—6.8—is statistically indistinguishable from the random expectation (based on the relative durations of hiatus to nonhiatus in the New Jersey succession). This test represents the worst-case scenario—the uncertainty limits in the age model are highly conservative, the assumption of uniform distributions is unrealistic for at least some of the calibration points, and the regression lines of the Monte Carlo trials were not constrained to pass through point 7 at the Cretaceous-Paleocene boundary. This test serves, however, to highlight the difficulties inherent in long-distance biocorrelation, difficulties that will only be overcome by the refinement of our age model and incorporation of additional, age-calibrated data points. From palynofacies analysis of four of the sections used herein, Schiøler et al. (2002) inferred that sea-level cycles in the East Coast basin showed a poor correlation with the EPR cycle chart. As noted already, in the present study, we focused on the well-constrained hiatus intervals identified from the New Jersey coastal plain; however, we agree with Schiøler et al. that Late Cretaceous unconformities in eastern New Zealand appear not to correlate closely with the EPR sequence boundaries. Correlations between the New Jersey hiatuses and the EPR sequences are discussed in Miller et al. (2003).
Despite the different approaches, there is substantial congruence between apparent sealevel signatures detected here using high-resolution quantitative biostratigraphy and those inferred using palynofacies analysis (Schiøler et al., 2002). Although detailed interpretation of all sections will be undertaken elsewhere, here we use the example of the Ben More section (Fig. 4) to compare results from the present study with evidence from field data and results from Schiøler et al. (2002). 1. A possible bundling of composite events at 30 m in the section does not correspond to any recognized lithological or palynofacies signature of condensation. This bundling, however, lies immediately below a 5 m interval of no exposure and a gap in sampling, and may be nothing more than an artifact of sample spacing. 2. Sequence boundary 1 of Schiøler et al. (2002), at 45 m, was not detected in the CONOP composite section and, again, corresponds to a gap in exposure. 3. A pronounced bundling of composite events occurs at 82 m. Within sample spacing uncertainty, this corresponds to both a conspicuous interval of red and greenish sandstone interbeds and sequence boundary 2 of Schiøler et al. (2002). 4. Maximum flooding surface 2 of Schiøler et al. (2002), at ~115 m, is not detected in the composite section, but is placed at a shell bed of disarticulated and fragmented bivalves. 5. An event cluster occurs at 135 m. This level was not distinguished by Schiøler et al. (2002), but corresponds to a conspicuous shell bed of disarticulated and fragmented bivalves and a gradational increase in mudstone content upsection. It is possible that this level represents the maximum flooding surface 2 of Schiøler et al. (2002). 6. Two event clusters are included within an event bundle at 150 m in the section. Although not distinguished by Schiøler et al. (2002), this level corresponds to an abrupt increase in sandstone proportion and loss of siltstone interbeds.
7. A relatively minor, but almost certainly significant, bundling of composite events occurs at the Paton-Herring Formation contact, at 168.5 m. Schiøler et al. (2002) placed their sequence boundary 3 at this level, and noted that this represents a relatively short hiatus in deposition. 8. A minor bundling of composite events, which includes one event cluster, occurs at 179 m. This level was not distinguished by Schiøler et al. (2002) and is marked by only a small lithological change, namely the loss of minor, centimeter-bedded, glauconic, fine sandstone. 9. Maximum flooding surface 3 of Schiøler et al. (2002) was placed tentatively at 208 m. This level probably corresponds to an event cluster identified herein. 10. Sequence boundary 4 of Schiøler et al. (2002), at 229 m, corresponds to a major bundling of composite events, including two event clusters. From these results it is clear that three of the four sequence boundaries recognized in the Ben More section by Schiøler et al. (2002), using palynofacies analysis, were detected here using quantitative biostratigraphic analysis. Sequence boundary 1 of Schiøler et al. is not identified using the CONOP composite section, possibly a consequence of relatively poor resolution and wide spacing of events in the lower part of the composite. Finally, it is interesting to note just how little time is actually represented by rock in the sections examined here (Fig. 6). At the scale of the study, on average just 24% of time is recorded across all the sections (maximum 47%, minimum 14%). These values will be overestimates for two reasons. First, because there is a positive relationship between the total range of any section and its apparent stratigraphic completeness, the measured levels of completeness would decrease at finer scales of resolution (e.g., Sadler, 1981). Secondly, our estimates of completeness do not take into account the regional
Figure 7. Clustering of events in the composite section. (A) Plot of number of events clustered at discrete levels in the composite section against inferred age and position in the composite. (B) Adjusted event clusters, E: event clusters from A are scaled by the average spacing between the event level in question and levels immediately above and below (see Equation 1). Arrows to the right of the plot indicate values of E that exceed the threshold of 10 (see explanation in D). In some cases, pairs of closely spaced maxima in E are regarded as a single significant peak (indicated by brackets next to the arrows). (C) Differences between the expected numbers of events at each level, based on a sixth-order polynomial fit to the cumulative distribution of events against time (shown in E), and the observed numbers of events. (D) Rankorder distribution of adjusted event clusters, E (shown in B, arranged from largest to smallest) plotted on semilog axes (event levels below 10% and above 90% in the composite section are ignored). Much of this distribution is log-linear, as indicated by the thick gray line. The extremes of the distribution deviate from log-linearity—there are more large event clusters than predicted by the log-linear relationship, whereas the tail of small clusters is truncated. We interpret clusters larger than E ≈10 (rank 36) as “significant” (see text). This threshold is indicated on B with a dashed line. (E) Cumulative distribution of events against time and sixth-order polynomial curved fitted to these data. Differences between the expected number of events at each level, based on this curve, and the observed number are shown in C.
Geological Society of America Bulletin, July/August 2006
987
et 0 M al 3, aj ., 20 o 19 0 id r e 87 4; en v ) tif en ie t c d lu he s re ter in s
G
N
q
20
ew co Je m rs (M po ey ille sit re e t a se l., qu 20 en 0 3 ce EP ,2 s R 00 (M s 4) ille eq af r ue te et n rH a c a l., es
lo b (G al ra ge ds oc te hr in on et ol al og ., ica 19 l 95 sc ) ale
Crampton et al.
100 Paleocene (part)
65 Late
TA1.1
90
Navesink I UZA4.5
80
Early
Maastrichtian
70
Navesink II
UZA4.4 60
75
Englishtown
UZA4.2 Merchantville III
80 Early E L
L
30
UZA3.4 UZA3.3 UZA3.2
20
Magothy III E M
Sant. Coniacian
Merchantville II Merchantville I Cheesequake
Magothy II?
UZA3.1
10
(part)
0
10 20 30 40 Adjusted event clusters, ( E)
unconformities identified by the event clusters. Our results are comparable with other studies that have determined levels of completeness in continental-margin marine strata ranging from ~20% to 50% at similar scales of study (e.g., Barrell, 1917; Kamp and Turner, 1990; Callomon, 1995). CONCLUSIONS Distinguishing tectonic from supposedly eustatic Upper Cretaceous stratigraphic sequences is problematic and contentious. A key
988
UZA3.5
255 45
Turonian
90
UZA4.1 40
98
85
UZA4.3 50
Middle
Campanian
Age (Ma)
Late
Marshalltown
Composite section (percentage units)
70
Figure 8. Adjusted event clusters from Figure 7B and sequence stratigraphic interpretations from Miller et al. (2003, 2004). The thick horizontal gray bars correspond to the ages of hiatuses in the New Jersey coastal plain succession that were identified by Miller et al. (2003, 2004). “EPR sequences” are the Exxon Production Research sequence boundaries from Haq et al. (1987) as recalibrated to the Gradstein et al. (1995) time scale by Miller et al. (2003, 2004). The EPR sequence boundaries are placed at the inflexion points of the EPR sea-level curve as presented in Figure 1 of Miller et al. (2003). The major event clusters identified herein are explained in Figure 7.
question remains unresolved: were there highfrequency, eustatic fluctuations during the Late Cretaceous, and, if so, what was their cause? In order to answer these questions, it is necessary first to demonstrate that Upper Cretaceous sequence boundaries are synchronous across the globe. Given that the frequency of inferred sequences is at or below the finest resolution of most standard biostratigraphic subdivisions, demonstration of synchroneity is extremely difficult. This problem is exacerbated when considering stratigraphic successions in New Zealand, which are structurally and/or stratigraphically
0
complex and remote from well-studied Northern Hemisphere localities. In order to address this problem, we used a new quantitative biostratigraphic method to derive a comparatively very high-resolution correlation for 15 onshore sections, spanning the Coniacian to Maastrichtian stages in eastern New Zealand. The method, constrained optimization (program CONOP), can be thought of as automated multidimensional graphic correlation. Our analysis was based on 398 wellconstrained lowest- and highest-occurrence bioevent datums derived from 245 dinoflagellate
Geological Society of America Bulletin, July/August 2006
Section missing at unconformities
Age
Ag e un mo co de n Ag for l 2, e mi wi un mo tie th co de s nf l 1 or , m w iti ith es o u
species. These data yielded a composite section with 182 discrete event levels that can be correlated between sections. The age of the composite was calibrated using an age model based on six dated bioevents, two of which were tied to geochemical or paleomagnetic datums. Following age calibration, the 182 composite event levels yielded an average relative temporal resolution of 128 k.y. through the entire Coniacian to Maastrichtian and just over 100 k.y. through the latest Coniacian to earliest Maastrichtian; these figures do not reflect precision in absolute dating, which is comparatively poor due to uncertainties in our age model. CONOP cannot resolve unconformities or condensed intervals that occur across all studied sections. Here we argue, however, that such “global” features can be recognized by the placement of multiple events at a single composite level. Before assigning significance to such event clusters, it is first necessary to account for spurious effects related to variable spacing of event levels in the composite. This was done here in two ways: by dividing the number of events by the average distance between adjacent levels, and by calculating the difference between the observed and expected number of events at any level, the expected number being estimated using a polynomial fit. Following adjustment for event-level spacing, we identified 15 major event clusters in Coniacian to Maastrichtian strata of eastern New Zealand. Of these, 10 coincide closely with well-dated hiatuses recorded in Upper Cretaceous sequences of the New Jersey coastal plain by Miller et al. (2003, 2004). The probability of this coincidence arising by chance alone is ~8%. We conclude that, at a resolution of less than ~±130 k.y., 10 New Zealand Coniacian to Maastrichtian unconformities or condensed intervals are likely to be coeval with hiatuses identified from the New Jersey coastal plain. Although these results must be regarded as provisional, they suggest that high-resolution quantitative stratigraphy can provide a potent tool for the identification and correlation of stratigraphic sequences. Furthermore, our findings add significant support for recent studies that have argued for the presence of Late Cretaceous eustatic changes. Finally, our results support other high-resolution studies that have demonstrated that the stratigraphic record of continental-margin marine strata is substantially incomplete, with typically ~20% to 50% of time represented by rock when observed at the scale of ~1 to ~25 m.y. Given the scale-dependence of completeness estimates, it is clear that far more time is missing in most continental-margin geological records than is present.
t
Quantitative biostratigraphy and Cretaceous eustasy
Composite section Figure 9. Age calibration of a hypothetical composite section using a simple linear age model (thick gray line, model 1) and a “true” age model that takes account of global unconformities (thin black line, model 2) (compare with Fig. 5). Dotted arrows illustrate the way in which events in the composite are age calibrated by projecting through the age model line onto the y axis.
ACKNOWLEDGMENTS
Crampton acknowledges funding from the Foundation for Research, Science and Technology program entitled “Towards self sufficiency: the geological basis for increasing New Zealand’s petroleum reserves” (contract C05X0302). Schiøler acknowledges a research grant from the Carlsberg Foundation, Denmark. Roncaglia acknowledges funding from the Ministry of University, Scientific Research and Technology, Italy. For assistance in the field and discussion, we would like to thank Peter King. We are grateful to landowners for access to the measured section localities. Craig Jones helped with data preparation; palynological samples were processed by Roger Tremain. For extensive and helpful discussions, and for reviewing earlier drafts of this paper, we gratefully acknowledge Ken Miller, Peter Sadler, and Roger Cooper. For useful discussions, we also acknowledge Charles Mitchell and Raphael Benites. We thank Peter Sadler for his ongoing development of CONOP and his generosity in making this program freely available. This is Institute of Geological & Nuclear Sciences contribution no. 3287. REFERENCES CITED Abreu, V.S., Hardenbol, J., Haddad, G.A., Baum, G.R., Droxler, A.W., and Vail, P.R., 1998, Oxygen isotope synthesis: A Cretaceous ice-house?, in Graciansky, P.-C.d., Hardenbol, J., Jacquin, T., and Vail, P.R., eds., Mesozoic and Cenozoic sequence stratigraphy of European basins: SEPM (Society for Sedimentary Geology) Special Publication 60, p. 75–80. Agterberg, F.P., 1990, Automated stratigraphic correlation: Amsterdam, Elsevier, Developments in Palaeontology and Stratigraphy, v. 13, 424 p. Ballance, P.F., 1993, The paleo-Pacific, post-subduction, passive margin thermal relaxation sequence (Late Creta-
ceous–Paleogene) of the drifting New Zealand continent, in Ballance, P.F., ed., South Pacific sedimentary basins: Sedimentary basins of the world, Volume 2: Amsterdam, Elsevier Science Publishers B.V., p. 93–110. Barrell, J., 1917, Rhythms and the measurements of geologic time: Geological Society of America Bulletin, v. 28, p. 745–904. Barrera, E., and Savin, S.M., 1999, Evolution of late Campanian–Maastrichtian marine climates and oceans, in Barrera, E., and Johnson, C.C., eds., Evolution of the Cretaceous ocean-climate system: Geological Society of America Special Paper 332, p. 245–282. Batten, D.J., 1999, Small palynomorphs, in Jones, T.P., and Rowe, N.P., eds., Fossil plants and spores: Modern techniques: Geological Society [London], p. 15–19. Callomon, J.H., 1995, Time from fossils: S.S. Buckman and Jurassic high-resolution geochronology, in Le Bas, M.J., ed., Milestones in Geology: Geological Society [London] Memoir, v. 16, p. 127–150. Carter, R.M., 1988, Post-breakup stratigraphy of the Kaikoura Synthem (Cretaceous-Cenozoic), continental margin, southeastern New Zealand: New Zealand Journal of Geology and Geophysics, v. 31, p. 405–429. Cooper, R.A., 2004, The New Zealand geological timescale: Institute of Geological and Nuclear Sciences Monograph, v. 22, p. 1–284. Cooper, R.A., Crampton, J.S., Raine, J.I., Gradstein, F.M., Morgans, H.E.G., Sadler, P.M., Strong, C.P., Waghorn, D., and Wilson, G.J., 2001, Quantitative biostratigraphy of the Taranaki Basin, New Zealand: A deterministic and probabilistic approach: American Association of Petroleum Geologists (AAPG) Bulletin, v. 85, p. 1469–1498. Crame, J.A., Pirrie, D., Riding, J.B., and Thomson, M.R.A., 1991, Campanian and Maastrichtian stratigraphy of the James Ross Island area, Antarctica: Journal of the Geological Society of London, v. 148, p. 1125–1140. Crampton, J., Mumme, T., Raine, I., Roncaglia, L., Schiøler, P., Strong, P., Turner, G., and Wilson, G., 2000, Revision of the Piripauan and Haumurian local stages and correlation of the Santonian-Maastrichtian (Late Cretaceous) in New Zealand: New Zealand Journal of Geology and Geophysics, v. 43, p. 309–333.
Geological Society of America Bulletin, July/August 2006
989
Crampton et al. Crampton, J.S., Raine, I., Strong, P., and Wilson, G., 2001, Integrated biostratigraphy of the Raukumara Series stratotype (Cenomanian to Coniacian) at Mangaotane Stream, Raukumara Peninsula, New Zealand: New Zealand Journal of Geology and Geophysics, v. 44, p. 365–389. Crampton, J., Laird, M., Nicol, A., Townsend, D., and Van Dissen, R., 2003, Palinspastic reconstructions of southeastern Marlborough, New Zealand, for mid-Cretaceous–Eocene times: New Zealand Journal of Geology and Geophysics, v. 46, p. 153–175. Field, B.D., Browne, G.H., Davy, B., Herzer, R.H., Hoskins, R.H., Raine, J.I., Wilson, G.J., Sewell, R.J., Smale, D., and Watters, W.A., 1989, Cretaceous and Cenozoic sedimentary basins and geological evolution of the Canterbury region, South Island, New Zealand: New Zealand Geological Survey Basin Studies, v. 2, p. 1–94. Field, B.D., Uruski, C., Beu, A., Browne, G., Crampton, J., Funnell, R., Killops, S., Laird, M., Mazengarb, C., Morgans, H., Rait, G., Smale, D., and Strong, P., 1997, Cretaceous-Cenozoic geology and petroleum systems of the East Coast region, New Zealand: Institute of Geological and Nuclear Sciences Monograph, v. 13, p. 1–244 + enclosures. Galeotti, S., Brinkhuis, H., and Huber, M., 2004, Records of post-Cretaceous-Tertiary boundary millennial-scale cooling from the western Tethys: A smoking gun for the impact-winter hypothesis?: Geology, v. 32, p. 529– 532, doi: 10.1130/G20439.1. Gradstein, F.M., and Agterberg, F.P., 1998, Uncertainty in stratigraphic correlation, in Gradstein, F.M., and Sandvik, K.O., eds., Sequence stratigraphy—concepts and applications: Norwegian Petroleum Society Special Publication 8, p. 9–29. Gradstein, F.M., Agterberg, F.P., Ogg, J.G., Hardenbol, J., van Veen, P., Thierry, J., and Huang, Z., 1995, A Triassic, Jurassic and Cretaceous time scale, in Berggren, W.A., Kent, D.V., Aubry, M.P., and Hardenbol, J., eds., Geochronology, time scales and global stratigraphic correlation: SEPM (Society for Sedimentary Geology) Special Publication 54, p. 95–126. Gradstein, F.M., Ogg, J.G., and Smith, A.G., 2004, A geologic time scale: Cambridge, Cambridge University Press, 589 p. Haq, B.U., Hardenbol, J., and Vail, P.R., 1987, Chronology of fluctuating sea levels since the Triassic (250 million years ago to present): Science, v. 235, p. 1156–1167. Helby, R., Morgan, R., and Partridge, A.D., 1987, A palynological zonation of the Australian Mesozoic: Memoir of the Association of Australasian Palaeontologists, v. 4, p. 1–94. Hollis, C.J., 2003, The Cretaceous/Tertiary boundary event in New Zealand: Profiling mass extinction: New Zealand Journal of Geology and Geophysics, v. 46, p. 307–321. Huber, B.T., Norris, R.D., and MacLeod, K.G., 2002, Deep-sea paleotemperature record of extreme warmth during the Cretaceous: Geology, v. 30, p. 123–126, doi: 10.1130/00917613(2002)0302.0.CO;2. Kamp, P.J.J., and Turner, G.M., 1990, Pleistocene unconformity-bounded shelf sequences (Wanganui Basin, New Zealand) correlated with global isotope record: Sedimentary Geology, v. 68, p. 155–161, doi: 10.1016/0037-0738(90)90125-D. Kemple, W.G., Sadler, P.M., and Strauss, D.J., 1995, Extending graphic correlation to many dimensions: Stratigraphic correlation as constrained optimization, in Mann, K.O., and Lane, H.R., eds., Graphic Correlation: SEPM (Society for Sedimentary Geology) Special Publication 53, p. 65–82. King, P.R., Naish, T.R., Browne, G.H., Field, B.D., and Edbrooke, S.W., 1999, Cretaceous to Recent sedimentary patterns in New Zealand: Institute of Geological and Nuclear Sciences Folio Series, v. 1, p. 1–35. LeMasurier, W.E., and Landis, C.A., 1996, Mantle-plume activity recorded by low-relief erosion surfaces in West Antarctica and New Zealand: Geological Society of America Bulletin, v. 108, p. 1450–1466, doi: 10.1130/0016-7606(1996)1082.3.CO;2. MacLeod, N., 1998, Impacts and marine invertebrate extinctions, in Grady, M.M., Hutchinson, R., McCall, G.J.H., and Rothery, D.A., eds., Meteoritex: flux with time and impact effects: Geological Society [London] Special Publication 140, p. 217–246. MacPherson, E.O., 1948, The Upper Senonian transgression in New Zealand: New Zealand Journal of Science and Technology, v. 29 (section B), p. 280–296.
990
Marshall, N.G., 1990a, The dinoflagellate Canninginopsis Cookson & Eisenack 1962 from the Cretaceous of Perth and Gippsland basins, Australia: Alcheringa, v. 14, p. 77–87. Marshall, N.G., 1990b, Campanian dinoflagellates from southeastern Australia: Alcheringa, v. 14, p. 1–38. Miall, A.D., 1992, Exxon global cycle chart: An event for every occasion?: Geology, v. 20, p. 787–790, doi: 10.1130/00917613(1992)0202.3.CO;2. Miall, A.D., and Miall, C.E., 2004, Empiricism and modelbuilding in stratigraphy: Around the hermeneutic circle in the pursuit of stratigraphic correlation: Stratigraphy, v. 1, p. 27–46. Miller, K.G., Barrera, E., Olsson, R.K., Sugarman, P.J., and Savin, S.M., 1999, Does ice drive early Maastrichtian eustasy?: Geology, v. 27, p. 783–786, doi: 10.1130/00917613(1999)0272.3.CO;2. Miller, K.G., Sugarman, P.J., Browning, J.V., Kominz, M.A., Hernández, J.C., Olsson, R.K., Wright, J.D., Feigenson, M.D., and Van Sickel, W., 2003, Late Cretaceous chronology of large, rapid sea-level changes: Glacioeustasy during the greenhouse world: Geology, v. 31, p. 585– 588, doi: 10.1130/0091-7613(2003)0312.0.CO;2. Miller, K.G., Sugarman, P.J., Browning, J.V., Kominz, M.A., Olsson, R.K., Feigenson, M.D., and Hernández, J.C., 2004, Upper Cretaceous sequences and sea-level history, New Jersey coastal plain: Geological Society of America Bulletin, v. 116, p. 368–393, doi: 10.1130/ B25279.1. Miller, K.G., Wright, J.D., and Browning, J.V., 2005a, Visions of ice sheets in a greenhouse world: Marine Geology, v. 217, p. 215–231, doi: 10.1016/j.margeo.2005.02.007. Miller, K.G., Kominz, M.A., Browning, J.V., Wright, J.D., Mountain, G.S., Katz, M.E., Sugarman, P.J., Cramer, B.S., Christie-Blick, N., and Pekar, S.F., 2005b, The Phanerozoic record of global sea-level change: Science, v. 310, p. 1293–1298. Mohr, B.A.R., and Mao, S., 1997, Maastrichtian dinocyst floras from Maud Rise and Georgia Basin (Southern Ocean): Their stratigraphic and paleoenvironmental implications: Palynology, v. 21, p. 41–65. Moore, P.R., 1988a, Structural divisions of eastern North Island: New Zealand Geological Survey Record, v. 30, p. 1–24. Moore, P.R., 1988b, Stratigraphy, composition and environment of deposition of the Whangai Formation and associated Late Cretaceous–Paleocene rocks, eastern North Island, New Zealand: New Zealand Geological Survey Bulletin, v. 100, p. 1–82. Pedley, H.M., and Frostick, L.E., 1999, Introduction: Unravelling tectonic and climatic signals in sedimentary successions: Journal of the Geological Society of London, v. 156, p. 747. Pirrie, D., and Riding, J.B., 1988, Sedimentology, palynology and structure of Humps Island, northern Antarctic Peninsula: British Antarctic Survey Bulletin, v. 80, p. 1–19. Pirrie, D., Crame, J.A., Lomas, S.A., and Riding, J.B., 1997, Late Cretaceous stratigraphy of the Admiralty Sound region, James Ross Basin, Antarctica: Cretaceous Research, v. 18, p. 109–137. Price, G.D., 1999, The evidence and implications of polar ice during the Mesozoic: Earth-Science Reviews, v. 48, p. 183–210. Riding, J.B., Keating, J.M., Snape, M.G., Newham, S., and Pirrie, D., 1992, Preliminary Jurassic and Cretaceous dinoflagellate cyst stratigraphy of the James Ross Island area, Antarctic Peninsula: Newsletters on Stratigraphy, v. 26, p. 19–39. Roncaglia, L., Field, B.D., Raine, J.I., Schiøler, P., and Wilson, G.J., 1999, Dinoflagellate biostratigraphy of Piripauan-Haumurian (Upper Cretaceous) sections from northeast South Island, New Zealand: Cretaceous Research, v. 20, p. 271–314. Sadler, P.M., 1981, Sedimentation rates and the completeness of stratigraphic sections: The Journal of Geology, v. 89, p. 569–584. Sadler, P.M., 2001, Constrained optimization approaches to the paleobiologic correlation and seriation problems: A users’ guide and reference manual to the CONOP program family: Riverside, Peter M. Sadler, 165 p. Sadler, P.M., 2004, Quantitative biostratigraphy—Achieving finer resolution in global correlation: Annual Review of Earth and Planetary Sciences, v. 32, p. 187–213, doi: 10.1146/annurev.earth.32.101802.120428.
Sadler, P.M., and Cooper, R.A., 2003, Best-fit intervals and consensus sequences; comparisons of the resolving power of traditional biostratigraphy and computerassisted correlation, in Harries, P.J., ed., High-resolution approaches in stratigraphic paleontology: Topics in geobiology: Dordrecht, Kluwer Academic Publishers, p. 49–94. Sadler, P.M., and Cooper, R.A., 2004, Calibration of the Ordovician timescale, in Webby, B.D., and Droser, M.L., eds., The great Ordovician biodiversification event: New York, Columbia University Press, p. 48–51. Schiøler, P., and Wilson, G.J., 1998, Dinoflagellate biostratigraphy of the middle Coniacian–lower Campanian (Upper Cretaceous) in south Marlborough, New Zealand: Micropaleontology, v. 44, p. 313–349. Schiøler, P., Crampton, J., and King, P., 2000, Palynostratigraphic analysis of a measured section through the Karekare, Owhena and Whangai Formations (Upper Cretaceous) at Waitahaia River—“Owhena Stream” (Raukumara Peninsula): Institute of Geological and Nuclear Sciences Science Report, v. 2000/31, p. 1–24. Schiøler, P., Crampton, J.S., and King, P., 2001, Palynostratigraphic analysis of a measured section through the Karekare, Tahora and Whangai Formations (Upper Cretaceous) at Koranga Stream (Raukumara Peninsula): Institute of Geological and Nuclear Sciences Science Report, v. 2001/27, p. 1–20. Schiøler, P., Crampton, J.S., and Laird, M.G., 2002, Palynofacies and sea-level changes in the middle Coniacian–late Campanian (Late Cretaceous) of the East Coast basin, New Zealand: Palaeogeography, Palaeoclimatology, Palaeoecology, v. 188, p. 101–125, doi: 10.1016/S0031-0182(02)00548-5. Stoll, H.M., and Schrag, D.P., 2000, High-resolution stable isotope records from the Upper Cretaceous rocks of Italy and Spain: Glacial episodes in a greenhouse planet?: Geological Society of America Bulletin, v. 112, p. 308–319, doi: 10.1130/0016-7606(2000)1122.3.CO;2. Strong, C.P., Hollis, C.J., and Wilson, G.J., 1995, Foraminiferal, radiolarian, and dinoflagellate biostratigraphy of Late Cretaceous to middle Eocene pelagic sediments (Muzzle Group), Mead Stream, Marlborough, New Zealand: New Zealand Journal of Geology and Geophysics, v. 38, p. 171–212. Sutherland, R., King, P., and Wood, R., 2001, Tectonic evolution of Cretaceous rift basins in south-eastern Australia and New Zealand: Implications for exploration risk assessment, PESA Eastern Australasian Basins Symposium; Melbourne, Victoria, 25–28 November 2001: Melbourne, Petroleum Exploration Society of Australia, p. 3–13. Van Sickel, W.A., Kominz, M.A., Miller, K.G., and Browning, J.V., 2004, Late Cretaceous and Cenozoic sea-level estimates: Backstripping analysis of borehole data, onshore New Jersey: Basin Research, v. 16, p. 451– 465, doi: 10.1111/j.1365-2117.2004.00242.x. Williams, G.L., Brinkhuis, H., Pearce, M.A., Fensome, R.A., and Weegink, J.W., 2004, Southern Ocean and global dinoflagellate cyst events compared: Index events for the Late Cretaceous–Neogene, in Exon, N.F., Kennett, J.P., and Malone, M.J., eds., The Tasmanian Gateway: Cenozoic climatic and oceanographic development; sites 1168–1172; 11 March–11 May 2000: Proceedings of the Ocean Drilling Program, Scientific Results, Volume 189, p. 1–98 [CD-ROM]. Willumsen, P.S., 2000, Late Cretaceous to early Paleocene palynological changes in midlatitude Southern Hemisphere, New Zealand: GFF, v. 122, p. 180–181. Wilson, G.J., 1984, New Zealand Late Jurassic to Eocene dinoflagellate biostratigraphy—A summary: Newsletters on Stratigraphy, v. 13, p. 104–117. Wilson, G.J., 1987, Dinoflagellate biostratigraphy of the Cretaceous-Tertiary boundary, mid–Waipara River section, north Canterbury, New Zealand: New Zealand Geological Survey Record, v. 20, p. 8–15. Young, G.C., and Laurie, J.R., 1996, An Australian Phanerozoic timescale: Melbourne, Oxford University Press, 279 p. + charts. MANUSCRIPT RECEIVED 26 MARCH 2005 REVISED MANUSCRIPT RECEIVED 13 DECEMBER 2005 MANUSCRIPT ACCEPTED 10 MARCH 2006 Printed in the USA
Geological Society of America Bulletin, July/August 2006