SUPPLEMENTARY INFORMATION
Supplementary Figures
Supplementary Fig. 1: Bio-, and lithostratigraphy in the Anti-Atlas
Supplementary Fig. 2: Bio-, and litho-, and chemostratigraphy at Anticosti Island Supplementary Fig. 3: Stratigraphic sampling of shelf settings during glaciation
Supplementary Tables
Supplementary Table 1: 13C dataset from Katian to lowermost Silurian at the Anticosti Island
Supplementary Table 2: Reservoirs of the OrdovicianSilurian global carbon cycle Supplementary Table 3: Fluxes of the OrdovicianSilurian global carbon cycle
Supplementary Note: The essentials of the three Late Ordovician Glacial Cycles
Supplementary Discussion: Box model of the late Ordovician carbon cycle
Supplementary References
1
Supplementary Figure 1
Bio-, and lithostratigraphy in the Anti-Atlas.
a
d
) is
part of the up to 5 km thick Palaeozoic succession preserved in the central Anti-Atlas depositional trough
(e , modified from Destombes et al.1). It was folded during the Hercynian orogeny and currently offers a
world-class record of predominantly shallow-marine Cambrian to Carboniferous depositional
sequences1,2,3. There, the Ordovician succession is up to 2 km in thickness. The Lower Palaeozoic
succession is dominated by shallow-marine siliciclastics. During the latest Ordovician, offshore shales
prevailed in the basin4 (Bou Ingarf/Tazzarine area). They graded laterally at basin edge into shoreface to
tidal facies, which migrated basinward during regressive events. Only in the middle-late Hirnantian,
glaciomarine to fluvioglacial environments arose.
Lithofacies and interpreted depositional environments are shown for four sections distributed from the
latest Ordovician basin edge to the basin centre of the central Anti-Atlas depositional trough. The two
basin-edge sections on the left-hand side (
In southern Morocco, the end-Ordovician record (
-
a
b
) specifically
represents a synthetic section with superimposed glacial erosion surfaces, the depth of which is in fact
deeper than illustrated, in the 75-200 m range (
) specifically represent latest Ordovician stratigraphies in
locations that do not show a subsequent glacial overdeepening event. The third log (
-
c
). In this location, fully bioturbated shallow-marine
intervals noticeably occur in between glacial erosion surfaces. The stratigraphic distribution of a pre-
glacial, lower to middle Hirnantian,
Hirnantia or Hirnantia-
f
) depicts a basin-axis location, within which the ice sheet arrived later. Here, three individual
glacial erosion surfaces numbered 3 to 5 are documented7, yet the glacial erosion surfaces 1, 2 and 6 are
not observed (i.e., related ice fronts did not reach the basin centre, see
right (
d
related fauna is shaded1,5,6. The section to the
). The Ouzregui Beds4, at the
f
Katian/Hirnantian boundary, correlate to the Pernik Beds of the Prague Basin7,8.
The late Katian to Hirnantian chitinozoan biozones9 of the upper Ktaoua and lower and upper Second Bani
formations are illustrated for the two main sections (unpublished data at basin edge; at the basin
centre4,10). In the basin-edge section, Upper Ktaoua and Lower Ktaoua formation names are in
brackets to indicate that they are not coeval with formations in the basin-centre type section; chitinozoans
of the Upper Ktaoua Fm. in the basin-edge sections belong to the
elongata
those of the Upper Ktaoua Fm. at the basin axis indicate the
merga
Biozone (Hirnantian), while
Biozone (latest Katian). As a
consequence, the Upper Ktaoua Fm. at the basin-edge is coeval with the lower member of the Second
Bani formation.
The lower diagram (
) tentatively illustrates the spatial and temporal distribution of glacial erosion
surfaces numbered 1-6 that have been mapped out and correlated from the basin edge to the basin centre
in the last ten years. Glacial surfaces are essentially amalgamated toward the basin edge, with the
exception of deep downcuttings or palaeovalleys (mainly tunnel valleys11) and intervening successions
progressively open basinward. Glacial erosion surfaces 1 and 2 have the smallest extent, as they have not
been so far observed in the Bou Ingarf area4,12 (BI). The Tizi nTazougart palaeovalley (TzT) has been
illustrated before1,13. Glacial erosion surfaces 3 and 4 correspond to the Hirnantian glacial maxima and
expand northward at least into the High Atlas, and at least one of the two reaches the Meseta area14. Glacial erosion surfaces 5 and 6 mark ice-sheet front readvances occurring during the overall deglaciation
at the end of the Hirnantian, but may be only of regional significance. The glacial erosion surface 6 is
associated with the tunnel valley known as the Foum Larjam palaeovalley (FoL 1,14). The complexity generated by superimposed glacial erosion surfaces has been schematically accounted for in between the
two sedimentary logs illustrated in the figure 2.
Above glaciation-related deposits, a
Hirnantia
f
fauna is preserved in places, which represents a postglacial
recovery distinct from the earlier (shaded)
Hirnantia
fauna that predates the occurrence of glacial
surfaces. Stratigraphic relationships at the transition from latest Ordovician to early Silurian strata are to
some extent obscured by a Telychian unconformity. Rhuddanian strata (lowermost Silurian) have been
documented overlying the Ordovician sandstones1,15 in the area corresponding to the three most proximal
sections. Conversely, at the basin axis, an erosional surface (transgressive surface) truncates the
uppermost glaciation-related strata, which include cryogenic structures, and which are sealed by
Telychian sandstones and shales of the
spiralis
graptolite Biozone15,16.
2
Supplementary Figure 2
Bio-, litho-, and chemostratigraphy at Anticosti Island.
a
to Telychian Anticosti succession (
and
d
The upper 900 m of the >2 km thick Sandbian
) constitutes a comprehensive, latest Ordovician to early
Silurian, record of thick, storm-dominated depositional sequences17. Mid to outer ramp carbonate facies
that prevail in the western part of Anticosti Island grade eastward towards the basin margin into thinner,
more siliciclastic-rich inner to proximal mid ramp facies that include several local discontinuities18,19. The
first-order stratigraphic trends of this exposed succession display a long-term shelf aggradation from the
late Katian to the early Telychian, culminating into a shelf progradation and basin fill due to the reduced
post-Taconic tectonic subsidence during the Telychian. The shelf aggradation phase is, however,
interrupted in the late Hirnantian by the presence of atypical shallow water limestones bounded by two
regional disconformities extending far into the basin.
The late Katian to early Rhuddanian chitinozoan biozones20 of the upper Vaural, Ellis Bay, and lower
Becscie formations are illustrated for the western (
) and eastern (
) sections of Anticosti Island. From
the base of the Ellis Bay Formation to the base of the uppermost Laframboise Member at the west end of
the island, three chitinozoan zones are distinguished in ascending order: the
florentini-concinna
Zone19,21. These zones are all considered Hirnantian in age, based
on several concordant paleontological data related to the occurrence of pre- and post-extinction
Hirnantian biota. This statement harmonizes with those previously reached on the basis of acritarchs22, of
brachiopods23-29, of stromatoporoids30, and of graptolites31 (the black star in the West End section locates
a graptolite assemblage identifying the N.
persculptus
Biozone). Other faunal groups display similar
patterns including crinoids32, nautiloids33, and rugose corals (McLean and Copper, written commun. 2012).
In the western section, the member names of the Ellis Bay Formation are in brackets to indicate that they
are not coeval with their eponyms of the eastern section; as an example, chitinozoans of the Grindstone
Member in the eastern section belong to the H.
crickmayi
concinna
Zone. As for the Anti-Atlas, revised chitinozoan
biostratigraphy results in regional-scale correlations that are noticeably different from lithostratigraphic
schemes. The lithostratigraphic framework of the latest Ordovician strata exposed on Anticosti Island is
currently under revision (P. Copper, pers. commun. 2012).
Depositional facies at the highly subsiding western end of the island are dominated by mid- to outer-ramp,
storm-dominated carbonates with calcareous shales18,34. Storm-influenced siliciclastic shoreface to mid
ramp sediments prevail at the eastern end of the island. Oncolitic limestones associated with local
metazoan-calcimicrobial reef development are present along the entire outcrop belt in the Laframboise
Member of the uppermost Ellis Bay Formation23.
Our high-resolution 13C curve (n= 135 micrites; (
) and Supplementary Table 1) extends from the
Vaural Formation up to the lower Becscie Formation at the west end of the island35. For the first 100
meters, 13C values are relatively stable and contain values of approximately 0 to 1 , with a mean of
0.4 and a standard deviation of 0.25 ; the orange strip in (
) represents the 95 % confidence
interval. These represent the background values for 13C prior to the very latest Katian. The Hirnantian
age of the Ellis Bay Formation confirms that the Hirnantian isotopic carbon excursion (HICE) is not
restricted to the main peak in the Laframboise Member, but includes the smaller excursions in the lower
part of the formation and in the uppermost part of the Vaureal Formation. The 13C drops to pre-excursion
values in the A.
ellisbayensis
Zone19.
The Middle Ordovician to Llandovery portion of the surface and subsurface stratigraphy of the Anticosti
Basin is illustrated in (
) (modified from Long17). The subsidence curve of the Anticosti succession shows
periods of increased subsidence rates during the Sandbian-Katian and Aeronian that are related to
Taconian and Salinic tectonic events further south in the Humber Zone and Gasp Belt36. Coupled with a
sustained sediment supply within the basin, the Anticosti record is exceptionally thick (e.g. Sandbian to
Katian 1600 m, Hirnantian ~100 m, Rhuddanian to Telychian 500m), one or two orders higher than
present in age equivalent carbonate sections of other shallow epeiric or ramp settings37. This argues
against the proposition, based on chemostratigraphic analysis38-41 that the Hirnantian and its associated
HICE could be restricted to less than 10 m within the uppermost Ellis Bay Formation (see also Fig. 4 in the
main text).
b
c
Zone and the
taugourdeaui
gamachiana
Zone, the
the western section indicate the H.
florentini
-C.
Zone, while those of the Grindstone Member in
a
a
N.
persculptus
chitinozoan zone at the base of the Becscie Formation during the uppermost
e
4
Supplementary Figure 3
Stratigraphic sampling of shelf settings during glaciation.
Conceptual models, based on Jerveys
approach42, illustrate relationships between the rates of relative base-level change (glacio-eustasy and
subsidence), initial water depth and sediment supply for shelfal archives during icehouse conditions. The
resulting stratigraphic columns in siliciclastic shelf environments after three scenarios are discussed
below.
The glacio-eustatic forcing is based on the Quaternary glacial/interglacial model from Isotopic Stages 1 to
11. The relative sea-level in (
a
b
c
). The green time intervals are
times with corresponding depositional units, while the red segments represent hiatuses, the latter
subaerial erosion surfaces or time-transgressive sedimentary condensations43. In (
) includes a long-term subsidence-related component of 40 m for time
interval under consideration, while no subsidence is assumed in (
) and (
), subsidence
combined with moderate initial water depths (~ 100 m) and relatively high sediment supply result in
~42% stratigraphic sampling and 58% hiatuses. The initial glacio-eustatic oscillations (time intervals 11-
7) are well represented in the depositional succession, while hiatuses (time intervals 6 and 4-2)
correspond to most of the later lowstand events44. The resulting picture is that of high-frequency cycles,
particularly from intervals with the greater sea-level highstands.
a
b
a
but without subsidence and with an initial water depth reduced to < 100 m.
The outcome is a thin succession with a very low (<20%) stratigraphic sampling and multiple, severe,
amalgamated hiatuses. Only the record of the earliest sea-level evolution is preserved despite abundant
sediment supply. In contrast, for scenario (
The scenario in (
) is as in (
),
), with a significant initial water depth of > 100 m, in general a
deeper shelf beyond the shelf roll-over, the stratigraphic sampling is substantial, at ~60%; despite low
sediment supply in this case, only the later glacio-eustatic lowstand events correspond to erosion surfaces.
However, the resulting stratigraphic column is thin, with a poorly decipherable record.
For a given glacio-eustatic scenario, it is the rates of shelf subsidence, sediment supply, and initial
depositional depths at glaciation onset that control sampling and temporal extent of stratigraphic units45. At any rate, stratigraphic hiatuses account for 40% to 80% of the time span in shelfal domains. Basin
overfilling resulting in erosion and hiatuses is delayed when subsidence is active (
c
) or initial water depth
a
c
is significant (
). Active subsidence results in a great number of depositional units and a representative
record, providing the rate of sediment supply is adjusted to subsidence rates (
a
), or great initial water depths (
b
), result in a limited number of well differentiated
c
). Negligible or moderate
),
several low- to high-frequency cycles are amalgamated. Similar relationships are expected in carbonate
platforms and all three scenarios are potentially applicable to the latest Ordovician case studies. Applied
to the Hirnantian record, we suggest that the Anticosti Island stratigraphy resembles the (
b)
depositional units44: in (
subsidence rates (
, depositional units essentially superimpose a set of discrete cycles; in (
c
)-type diagram.
a
). The outer shelf
includes a pro- to inter-glacial stratigraphic record, generating sequences somewhat similar to the (
See also PermianCarboniferous case studies46,47.
For the glaciated shelf, contrasting records arise at outer and inner shelf settings (
d
-
a
c
)
scenarios (green parts) with hiatuses only at glacial maxima (in red). The inner shelf record, on the other
hand, samples major interglacials, and most of the time is represented by erosion surfaces48,49,50. Interestingly, four glacial time intervals and related glacial erosion surfaces are captured in both cases, but
they do not represent coeval glacial cycles. Clearly, the glaciated inner shelf record (e.g. Mauritania, Libya
or Niger in the end-Ordovician12,51,52) is not correlatable in a one-to-one manner to depositional units of
the far-field Anticosti Island stratigraphic sampler. The Anti-Atlas record that is understood as that of a
glaciated outer shelf12 is expected to correlate more closely with the Anticosti Island record, except at the
time of glacial climax.
6
Supplementary Table 1
13C dataset from Katian to lowermost Silurian at the Anticosti Island.
Stable isotope
measurements were made on either micritic matrixes or the micritic phase of pelloids in grainstone facies35. Some samples were analyzed more than once: "QCD" indicates a Quality Control Duplicate, which is used to test the accuracy of the mass spectrometer, "repeat" denotes a sample that was analyzed again due to it being an anomalous measurement, or having encountered a problem in the process of measurement. The data are listed relative to the PDB standard. Sampled section 1A, West End (Baie St. Claire Laframboise Section); its base is in the Homard Member of the Vaureal Formation; it spans the entire Ellis Bay Formation, and terminates in the lower Fox Point Member of the Becscie Formation.
Sample #
Height
(m) Description 13C () 18O ()
B-M-01 0 mudstone/packstone 0.3 -3.22
B-M-02 1.33 wackestone/packstone 0.24 -2.98
B-M-03 2.66 wackestone/packstone 0.18 -3.15
B-M-03 QCD 2.66 wackestone/packstone 0.16 -3.35
B-M-04 4 wackestone/packstone 0.32 -2.74
B-M-05 5.33 wackestone/packstone 0.05 -4.05
B-M-06 6.66 mudstone/packstone 0.26 -3
B-M-07 8 wackestone 0.2 -3.11
B-M-08 9.33 wackestone 0.07 -3.92
B-M-09 10.66 packstone 0.26 -3.68
B-M-10 12 mudstone/wackestone 0.08 -4.07
B-M-11 13.33 wackestone/packstone 0.4 -4.34
B-M-12 14.66 mudstone/packstone 0.41 -4.25
B-M-13 16 mudstone/wackestone 0.23 -4.66
B-M-13 QCD 16 mudstone/wackestone 0.26 -4.64
B-M-14 17.33 mudstone 0.59 -4.27
B-M-15 18.66 mudstone/wackestone 0.22 -4.32
B-M-16 20 mudstone/wackestone 0.49 -3.83
B-M-17 21.8 wackestone 0.66 -5.01
B-M-18 23.6 mudstone/wackestone -0.18 -4.74
B-M-19 25.4 wackestone 0.32 -4.83
B-M-20 27.3 wackestone/packstone 0.35 -4.13
B-M-21 29.1 wackestone/packstone 0.63 -4.02
B-M-22 29.9 packstone 0.18 -3.88
B-M-23 32.7 packstone -0.03 -4.13
B-M-24 34.6 wackestone 0.51 -4.45
B-M-24 QCD 34.6 wackestone 0.48 -4.47
B-M-25 36.4 wackestone 0.49 -4.13
B-M-26 38.2 mudstone/wackestone 0.52 -3.87
B-M-27 40 wackestone 0.43 -4.74
B-M-28 41.8 wackestone 1.53 -1.68
B-M-28
repeat 41.8 wackestone 1.75 -1.47
B-M-29 43.6 mudstone/wackestone 0.6 -3.42
B-M-30 45.4 mudstone/wackestone 0.3 -4.12
B-M-31 47.3 mudstone/wackestone 0.46 -4.48
B-M-32 49.1 wackestone 0.7 -4.26
8
B-M-33 49.9 wackestone/packstone 0.7 -3.64
B-M-33 QCD 49.9 wackestone/packstone 0.72 -3.96
B-M-34 52.7 mudstone/wackestone 0.68 -4.05
B-M-35 54.6 packstone 0.36 -4.23
B-M-36 56.4 mudstone/packstone 0.85 -3.58
B-M-37 58.2 wackestone 0.39 -3.83
B-M-38 60 wackestone/packstone 0.88 -3.48
B-M-39 61.8 mudstone/wackestone 0.58 -4.34
B-M-40 63.6 wackestone 0.15 -3.94
B-M-41 65.4 wackestone 0.37 -3.65
B-M-42 67.3 mudstone/wackestone 0.48 -3.53
B-M-43 69.1 wackestone/packstone 0.38 -4.18
B-M-43 QCD 69.1 wackestone/packstone 0.36 -4.18
B-M-44 69.9 mudstone/wackestone 0.52 -3.65
B-M-45 72.7 mudstone/packstone 0.62 -3.36
B-M-46 74.6 mudstone/packstone 0.32 -3.57
B-M-47 76.4 wackestone/packstone 0.71 -4.26
B-M-48 78.2 wackestone/packstone 0.01 -3.85
B-M-49 80 mudstone/wackestone -0.01 -4.23
B-M-50 81.66 wackestone/packstone 0.28 -3.72
B-M-51 83.33 mudstone/wackestone (stylolitic) -0.07 -4.62
B-M-52 85 mudstone/wackestone 0.29 -3.43
B-M-53 86.66 mudstone 0.44 -3.68
B-M-53 QCD 86.66 mudstone 0.43 -3.65
B-M-54 88.33 mudstone 0.57 -3.76
B-M-55 90 mudstone 0.3 -3.71
B-M-56 91.66 mudstone/wackestone 0.36 -3.62
B-M-57 93.33 mudstone/wackestone 0.33 -3.69
B-M-58 95 wackestone 0.4 -3.82
B-M-59 96.66 wackestone 0.45 -3.66
B-M-60 98.33 wackestone 0.16 -3.97
B-M-61 100 mudstone (stylolitic) 0.55 -3.45
B-M-62 101.66 mudstone 0.78 -3.46
B-M-63 103.33 mudstone/wackestone 0.97 -2.89
B-M-63 QCD 103.33 mudstone/wackestone 1 -2.9
B-M-64 105 wackestone/packstone 1.22 -3.02
B-M-65 106.66 mudstone/wackestone 0.93 -3.08
L-M-01 100 packstone (brachiopod) 0.8 -3.46
L-M-02 102.5 wackestone/packstone 0.71 -3.26
L-M-03 105 wackestone 1.06 -4.58
L-M-04 107.5 packstone 1.01 -3.21
L-M-05 110 packstone 0.84 -3.36
L-M-06 110.6 packstone 0.44 -3.76
L-M-07 111.1 packstone (brachiopod) 0.25 -3.13
L-M-08 111.7 wackestone/packstone 0.76 -3.04
L-M-08 QCD 111.7 wackestone/packstone 0.74 -2.99
L-M-09 112.2 packstone 0.61 -3.22
L-M-10 112.8 wackestone/packstone 0.98 -3.04
L-M-11 113.3 wackestone 0.26 -3.43
L-M-12 113.9 packstone 0.66 -3.33
L-M-13 114.4 packstone 1.16 -3.48
L-M-14 115 packstone/grainstone 0.85 -3.89
9
L-M-15 115.83 packstone 0.98 -3.67
L-M-16 116.66 packstone 1.61 -3.04
L-M-17 117.5 grainstone (peloidal) 1.94 -3.13
L-M-18 118.33 packstone 1.59 -3.44
L-M-18 QCD 118.33 packstone 1.61 -3.42
L-M-19 119.16 packstone 1.6 -3.24
L-M-20 120 packstone 0.28 -3.25
L-M-21 120.83 wackestone 2.26 -3.26
L-M-22 121.66 wackestone 1.47 -3.66
L-M-23 122.5 wackestone 2.36 -3.18
L-M-24 123.33 mudstone 2.35 -3.02
L-M-25 124.16 packstone 1.17 -3.55
L-M-26 125 packstone 1.1 -3.96
L-M-27 127.5 wackestone/packstone 1.69 -3.49
L-M-28 130 packstone 1.44 -3.5
L-M-28 QCD 130 packstone 1.51 -3.48
L-M-29 132.5 wackestone/packstone 1.32 -3.14
L-M-30 134.6 wackestone 1.9 -3.28
L-M-31 136.7 mudstone/wackestone 1.78 -3.32
L-M-32 138.8 mudstone/wackestone 1.56 -3.43
L-M-33 140.9 wackestone 1.29 -3.26
L-M-34 143 wackestone 0.98 -3.15
L-M-35 145.1 wackestone 1.15 -3.48
L-M-36 147.2 wackestone 1.19 -3.25
L-M-37 149.3 mudstone/wackestone 1.47 -3.15
L-M-38 151.4 wackestone 0.81 -3.28
L-M-38 QCD 151.4 wackestone 0.85 -3.29
L-M-39 153.5 mudstone 1.46 -3.36
L-M-39
repeat 153.5 mudstone 1.54 -3.75
L-M-40 155.6 mudstone/wackestone 0.95 -3.34
L-M-41 157.4 mudstone 0.95 -3.46
L-M-42 159.2 wackestone 1.04 -3.45
L-M-43 161 wackestone/packstone 1.08 -3.4
L-M-44 162.8 wackestone 0.9 -3.28
L-M-45 164.6 wackestone/packstone 0.54 -3.38
L-M-46 166.4 wackestone/packstone 0.66 -3.31
L-M-47 168.2 packstone 0.5 -3.73
L-M-48 170 wackestone 0.99 -3.51
L-M-49 171.8 wackestone 0.02 -3.18
L-M-50 173.6 mudstone/wackestone 0.72 -3.32
L-M-50 QCD 173.6 mudstone/wackestone 0.72 -3.34
L-M-51 175.4 mudstone/wackestone 0.82 -3.14
L-M-52 177.2 wackestone 0.84 -3.26
L-M-53 179 wackestone 1.11 -3.9
L-M-54 180.8 wackestone 1.04 -3.65
L-M-55 182.6 wackestone 1.02 -3.62
L-M-56 184.4 wackestone 0.62 -4.58
L-M-57 186.2 packstone 0.74 -3.88
L-M-58 188 packstone 0.49 -3.47
L-M-59 188.9 wackestone/packstone 0.47 -3.61
L-M-59 QCD 188.9 wackestone/packstone 0.55 -3.53
L-M-60 189.8 wackestone/packstone 0.29 -3.52
10
L-M-61 190.7 wackestone/packstone 0.2 -3.71
L-M-62 191.6 wackestone/packstone 0.26 -3.83
L-M-63 192.4 packstone (brachiopod) 0.47 -3.81
L-M-64 193.3 packstone 0.52 -4.06
L-M-65 194.2 packstone 0.69 -3.4
L-M-66 195.1 packstone 0.94 -3.1
L-M-67 196 wackestone/packstone 1.99 -3.03
L-M-68 197.5 mudstone 1.95 -3.71
L-M-69 199 mudstone/wackestone 1.93 -3.69
L-M-70 201 wackestone 1.92 -3.92
PL-4i 197.5 Peloidal grst/pakst 1.80 -2.75
PL-5i 198 Peloidal grst/packst 1.84 -3.12
PL-6i 198.5 Peloidal grst/packst 2.20 -3.53
PL-7i 199 Peloidal grst/packst 2.27 -3.42
PL-8i 200 Peloidal grst/packst (top Lousy Cove) 2.05 -3.87
PL-9i 200.5 Oncolitic packst/grst (base Laframboise) 2.69 -3.46
PL-10i 200.9 Oncolitic packst/grst 3.56 -2.73
PL-11i 201.1 Oncolitic packst/grst 3.85 -2.38
PL-12i 201.5 Oncolitic packst/grst 4.04 -2.74
PL-13i 201.8 Oncolitic packst/grst 3.78 -2.88
PL-14i 202.1 Oncolitic packst/grst 3.48 -2.69
PL-15i 202.8 Inter-reef wackst/packst 2.95 -2.96
PL-16i 203.2 Inter-reef wackst/packst 3.85 -2.75
PL-17Bi 203.8 Inter-reef wackst/packst 3.60 -2.85
PL-18i 204.1 Inter-reef wackst/packst (top Laframboise) 3.63 -2.28
PL-19Ai 204.3 packstone-wackestone (base Becscie) 2.18 -4.39
PL-20i 204.5 packstone-wackestone 2.31 -4.36
PL-21i 204.8 packstone-wackestone 1.49 -3.64
PL-22i 205.3 wackestone 1.16 -3.15
PL-23i 205.6 wackestone 1.18 -3.08
PL-24i 205.9 wackestone 0.58 -3.99
PL-25i 206.2 wackestone 0.56 -3.34
PL-26i 206.5 wackestone 0.00 -4.24
PL-27i 207.5 wackestone 0.42 -3.80
PL-28i 208.5 wackestone 0.34 -3.69
11
Supplementary Table 2
Reservoirs of the OrdovicianSilurian global carbon
Estimated quantities of the reservoirs of the OrdovicianSilurian global carbon cycle35. With flux quantities (Supplementary
Table 3), these are the basis for the box model presented in the Supplementary Discussion.
cycle.
Reservoir abbr. Gt C (today) Gt C (O-S) 13C ()
Lithosphere l 150,000,000 69,802,795 -6 Carbonates c 70,000,000 150,000,000 0 Fossil Fuels f 20,000 10,000 -28 Reactive Sediments r 3,000 18,000 0 Deep Ocean d 38,000 220,000 0 Surface Ocean s 1,000 6,000 +3 Phytomass p 500 5 -28 Soil 2500 0 Atmosphere a 800 9,000 -6
Supplementary Table 3
Fluxes of the OrdovicianSilurian global carbon cycle
. Estimated quantities for the fluxes of the OrdovicianSilurian global carbon cycle35. With reservoir quantities (Supplementary Table 2), these are the basis for the box model presented in the Supplementary Discussion. Flux abbr. Gt/yr (today) Gt/yr (O-S)
Terrestrial primary production a
a 63.1 0
Marine primary production a
a 50.5 80
Volatilization from soil a
a 62.5 0
CO2 dissolution & evasion a
a 96 1050
CaCO3 production & dissolution s
s 0.5 0.5
CO2 uptake by plants & humus a
p 0.6 0
CO2 used in weathering a
l 0.26 0.13
River input from silicates l
s 0.25 0.13
River input from carbonates c
s 0.13 0.06
River input from organic matter c
s 0.31 0
Ocean-atmosphere exchange s
a 0.48 0
CaCO3 storage in sediments s
c 0.38 0.38
Organic C storage in sediments s
c 0.1 0.1
Upwelling d
s 2.15 12.5
Volcanism & metamorphism l
a 0.12 0.18
Hydrothermal l
a 0.1 0.15
Uplift l
a 0.4 0.4
12
Supplementary Note
The essentials of the three Late Ordovician Glacial Cycles (LOGCs)
We here summarize the essentials that characterize the three Late Ordovician Glacial Cycles, as understood from the high-resolution sequence stratigraphic frameworks (Figure 2), in the near-field Anti-Atlas (Supplementary Fig. 1) and the far-field Anticosti Island (Supplementary Fig. 2).
Anti-Atlas
LOGC 1 includes a severe latest Katian sea-level fall reflected by a major facies shift at the basin edge and an ensuing important transgression with basin-wide sediment starvation and condensation in the very latest Katian. Maximum regressive and early transgressive facies together form the Ouzregui Beds4, coeval with a significant faunal turnover, corresponding to the replacement of the diversified Late Katian faunas by a poorly diversified
fauna, which is only present at basin edges, and not in deeper parts of the depocentre.
The lower to middle Hirnantian LOGC 2 commences with a highstand. Then, two high-order GSS with strikingly sharp-based, regressive depositional units, characterize its lower part. The older GSS is poorly developed in basinal position, while the younger one is best recognized in the basin centre. No subaerial exposure occurred at this time at the basin centre, but is suspected at the basin edge. Associated significant sedimentary aggradation suggests that this regressive succession is a lowstand wedge reflecting early time-transgressive conditions immediately following the glacial maximum of LOGC2. A relatively long-term transgressive trend followed that included well-defined higher-order oscillations capped by a major flooding surface with phosphogenesis.
The late Hirnantian LOGC 3 is essentially preserved at the basin axis, and/or within restricted glacially-related overdeepenings (Supplementary Fig. 1). Thin regressive nearshore facies ascribed to falling stage deposits are truncated by a glacial wedge (the glacial interval in Fig. 2) that includes several polyphased glacial erosion surfaces and related glaciomarine to fluvioglacial units. Within the glacial wedge, glacio-eustatic cycles are difficult to decipher because glacio-eustasy here is expected to have interfered with glacio-isostasy. The subsequent post-glacial transgression is associated with renewed deposition at the basin margin, re- colonisation by a
fauna, and a severe latest Hirnantian to Rhuddanian condensation1.
In the basin centre, an early Silurian unconformity of unknown origin associated with a ca. 7 myr long hiatus truncates the very latest deglacial Ordovician record16.
The first-order stratigraphic trends reveal a long-term shelf progradation through the latest Katian to the late Hirnantian, which was punctuated by multi-order regressive and/or glacial events. The Hirnantian glacial record included in LOGC 3 is only preserved in a lowstand, basinal position, with virtually no record (except in glacial overdeepenings) at the basin edge. Post- glacial flooding was non-accretionary53, suggesting high rates of sea-level rise in the very latest Hirnantian.
The glacial record (glacial erosion surface, glaciotectonic deformation, tunnel channels, ice- contact deposits) of LOGC 1 is known in Niger51 as re-interpreted in Loi et al.4, while that of LOGC 2 and 3 most likely correspond to the well known glacial successions in Libya, Algeria and Mauritania12,52,54-59.
13
-related
Hirnantia
Hirnantia
Anticosti Island
Several orders of depositional units, reflecting multi-order base-level changes, were identified in the Anticosti succession that display striking similarities with the time-equivalent Anti-Atlas sequence. During the latest Katian LOGC 1, a major sea-level fall was followed by a significant transgressive event. This event, together with the ensuing highstand in the latest part of the Katian, is associated with a faunal turnover during which Katian acritarchs, chitinozoans, conodonts, brachiopods, nautiloids, crinoids, stromatoporoids, and corals are replaced by taxa with either Hirnantian or Silurian affinities. The first perturbation in a long-lived relatively stable Katian 13C signal coincide with late regressive conditions in the latest Katian (Mill Bay), and not with the earlier sea-level fall (Joseph Point; Fig. 3 and Supplementary Fig. 2). Two sharp- based, regressive units representing sea-level drops are well expressed in the western distal basin sections during the early-middle Hirnantian LOGC 2. The older one, characterized by a greater facies offset than the younger one, is associated with a basal regressive surface of marine erosion resulting in a stratigraphic hiatus during the lowermost Hirnantian. Subaerial exposure did not occur at that time at the basin centre, but was present at the basin margin. The 13C
values are typically above the Katian background with a positive 2 excursion recorded above the first sharp-based surface. A well-expressed transgressive trend with higher-order oscillations is capped by a major flooding surface in the upper LOGC 2. This flooding event marks a return to typical Katian 13C values.
The middle-late Hirnantian LOGC 3 is composed of three distinct stratigraphic packages separated by two regional disconformities. The oldest package is a sharp-based regressive unit representing a major sea level drop. Its capping erosive surface recorded an emersion that was smoothed by ravinement during the ensuing transgression (see Fig. 4). This regressive unit coincides with a progressive increase in 13C values, up to +2. A second faunal turnover is recognized following the deposition and subsequent emersion of this initial package. This second turnover shows a more abrupt replacement of acritarchs, chitinozoans, conodonts, brachiopods, and corals than the first turnover, with the rapid disappearance of Ordovician taxa. The next package is composed of transgressive oncolitic calcirudites overlain locally by keep-up metazoan-calcimicrobial bioherms. The upper contact of the bioherms is erosional has local relief up to 10 m, and has a multi-phase origin including an initial emersion, a subsequent modification by a transgressive ravinement, and a final pyritic hardground development. The highest positive 13C values, up to 5 in places, are present in this middle package. The third package, locally onlaps and abuts against the exhumed massive bioherm cores of the underlying package. It displays a thin transgressive record at the more subsiding basin centre, but contains thicker, slightly older proximal ramp facies at the basin margin. This final package marks the return of pre-Hirnantian shelf aggradation architecture and displays a relatively rapid negative isotopic shift with return to 13C background values. The late Hirnantian LOGC-3 glacial far-field record is partially preserved at the basin centre, but reduced at the basin edge.
Within LOGCs 1 and 2, the 13Ccarb curve rises during the late and early regressive system tracts (lowstand and highstand conditions, respectively) and declines during transgressive and late regressive system tracts, respectively. Note that within LOGC 2, the excursion encompasses several higer-order stratigraphic cycles. The third and greatest excursion recorded in LOGC 3 amalgamates two signals, one predating and one postdating the LOGC 3 glacial maximum that is
14
represented by the unconformity at the base of the Laframboise Mb. These are time-regressive and time-transgressive, respectively (Fig. 4).
The two time intervals that correspond to the first and second faunal turnovers are not restricted to two short-term extinction events, supposedly glacial onset and termination, respectively. In fact, the first turnover is essentially coincident with the first interglacial that separates LOGCs 1 and 2. The ensuing lowermost Hirnantian stratigraphic hiatus in the Anticosti Island succession is likely responsible for its apparent sharpness (Fig. 3). The second turnover includes the entire glacial maximum of LOGC 3, commencing during the glacio-eustatic regression and terminating during the early deglaciation phase.
15
Supplementary Discussion
Box model of the late Ordovician carbon cycle
A variety of models have been proposed to explain the large isotopic excursion(s) in the Hirnantian; the productivity hypothesis, the weathering hypothesis and changes in oceanic circulation pattern. These are successively examined below, in the light of a box model for the Late Ordovician global carbon cycle35, with the conclusion that none of them can account for the amplitude of the observed anomaly at the global scale.
Presented here is a model designed to simulate the global carbon cycle for the Late Ordovician world. It is based mostly on the work of Mackenzie and Lerman60, a review of hundreds of scientific studies of the past and present carbon cycling. This overview quantifies carbon reservoirs and fluxes on global scales. The Ordovician-Silurian carbon cycle, based on Mackenzie and Lermans model, uses high Ordovician pCO2 values (~4000 ppm) and also takes into account minimal vascular land-plant cover. This model can be used to test some of the theoretical aspects of the hypotheses concerning the 13C excursions that occurred near the O-S Boundary. The estimated quantities of the reservoirs and the fluxes of the Ordovician-Silurian global carbon cycle are given in Supplementary Table 2 and Supplementary Table 3, respectively.
In this simplified model, the relationships between various reservoirs, fluxes and isotope values
of carbon are described by the conservation of mass Mx/t = F
i-x F
x-i (1) flux in flux out and a similar equation involving the enrichment of organic carbon. (Mx*x)/t = Fi-x * i F
x-i * x + F
i-x * (i+) F
x-i * (x+) (2) inorganic carbon organic carbon
Using the product rule and the two equations above we arrive at the equation of isotope continuity.
x/t = [ Fi-x * (i-x) + F
Ms = 6000 Gt = 28 Fl-s = 0.13 Gt/yr a = -6 Fc-s = 0.06 Gt/yr l = -6 Fd-s = 12.5 Gt/yr c = 0 Fs-a = 0 Gt/yr c = -28 Fs-c = 0.38 Gt/yr d = 0 Fs-c = 0.1 Gt/yr s = +3
Where: Mx represents the mass of C in a reservoir
i-x * (i+) F
x-i * (x+)] / Mx (3)
Fi-x is the flux of C from reservoir i into reservoir x Fx-i is the flux of organic C from reservoir x to reservoir i x is the isotopic value of a carbon reservoir
is the depletion factor for organic carbon
16
, a change in the rate of burial of organic carbon, Fs-c, brought the 13C of the surface ocean from 0 to 4 , in the timespan of approximately 100,000 to 500,000 years.
s * Ms /t = Fl-s * (l-s) + Fc-s * (c-s) + Fd-s * (d-s) - F
s-c * (s+) (4)
global productivity hypothesis
In the case of the
+4 * 6000 Gt/t = 0.13 Gt/yr * (-6 -3 ) + 0.06 Gt/yr * (0 -3 ) + 12.5 Gt/yr * (0 - 3 ) F
s-c * (3 -28 ) 24,000 Gt /t = (-39 Gt/yr + 25 * F
For the lower limit, t = 100,000 yr, F
s-c)
s-c = 1.564 Gt/yr
For the upper limit, t = 500,000 yr, F
s-c = 1.556 Gt/yr
Therefore, to produce a 13C increase of 4 in the surface ocean, the carbon burial rate has to increase to approximately 1.56 Gt/yr, 15 times the present day rate of carbon burial in the oceans (0.1 Gt/yr); an unsustainable proposition on a global scale.
Cramer and Saltzmans hypothesis61,62 for ocean state changes, the value for F
d-s, which represents upwelling of inorganic carbon from the deep ocean to the surface ocean, will have to change from 12.5 Gt/yr to 0, assuming stratified oceans with no active thermohaline circulation. This is difficult to conceive on a global scale but can be easily achieved on regional (basinal) scales.
+4 * 6000 Gt/t = 0.13 Gt/yr * (-6 -3 ) + 0.06 Gt/yr * (0 -3 ) + 0 Gt/yr * (0 -
3 ) F
s-c * (3 -28 ) 24,000 Gt /t = (-1.35 Gt/yr + 25 * F
t = 100,000 yr, F
s-c = 0.064 Gt/yr t = 500,000 yr, F
s-c = 0.056 Gt/yr
s-c)
Therefore, a 13C increase of 4 in the surface ocean is possible with modern day burial rates of organic carbon but only on regional scales and providing the upwelling of water from the deep ocean were to cease completely. Considering that the tide-related recirculation of deep, dense water masses to the surface ocean is enhanced during lowstand events63, such
circumstances are in fact unlikely.
Note that in the Late Ordovician world, organic carbon was produced exclusively in the Surface Ocean Reservoir. Taking the above qualifications into account, let us now consider the viability of these earlier advocated Late Ordovician scenarios in the context of geological framework.
1)
The
productivity hypothesis64
argues that phytoplankton blooms resulted in preferential removal of 13C from the water column, leading to a drawdown of atmospheric CO2 that initiated the Hirnantian glaciation, sea level drop, and generation of a widespread anoxia followed by the late Ordovician extinction event. The above model calculations35 show, however, that the rate of organic carbon burial would have to be 15 times that of its modern counterpart and sustained over 107108 years. This is an unrealistic proposition, even leaving aside the issue of the fate of the missing carbon-rich sediments in coeval sedimentary sections. This could have been a
17
viable scenario only if applied to localized basins within the broad epeiric seas of the Ordovician that may not have been strictly synchronous. 2)
The alternative
weathering hypothesis
and its modifications65-67 attributes the glaciation to CO2 drawdown initiated by enhanced silicate weathering related to the Taconic Orogeny. The erosion of platform carbonates subsequent to glacially-induced sea-level fall is then advocated as an explanation for the HICE. The "weathering" scenario requires as a starting assumption riverine flux of carbon that is significantly depleted in 13C. Such isotopically depleted carbon is presently derived from soil CO2 that originates from decomposition of land-based biomass. The positive carbon excursion in the ocean is then driven by diminution of the input from such sources. Yet, the land-based biosphere prior to Silurian was either absent or putative and the input from soil CO2 into the riverine systems has therefore been limited. At that time the dissolution of carbonates must have been dominated by carbonic acid derived mostly from ingassing of atmosperic CO2 and the isotope signal of the riverine carbon flux would have been around 0 %. Moreover, the presumed additional erosional source associated with low sea levels, the underlying Paleozoic rocks, have 13C depleted (01 )68 relative to HICE and thus cannot be the cause of the anomaly. The onset of 13C excursions during regressive time intervals thus cannot be the consequence of enhanced erosion of platform carbonates, unless a significant land cover (and related massive production of carbonic acid) can be demonstrated in the Ordovician. 3)
Another alternative argues that the sea level and/or climate triggered
, from upwelling dominated shelf circulations during highstands to downwelling during lowstand, resulted in redox stratification with a 13C-rich upper layer due to enhanced productivity and a 13C-depleted water body at depth. This hypothesis, in essence developed for Silurian 13C excursions, suffers the same limitations as (1) described above. It cannot be produced and sustained on the scale of global oceans. This scenario is feasible only for excursions developed on a basin scale during highstand conditions of a high-order GSS (see Fig. 4).
changes in circulation
patterns 61,62,69-71
18
Supplementary References
1. Destombes, J., Hollard, H. & Willefert, S. Lower Palaeozoic rocks of Morocco. In
Lower Palaeozoic of
(ed C.H. Holland) 91336 (John Wiley, New-York, 1985).
2. Burkhard, M., Caritg, S., Helg, U., Robert-Charrue, C. & Soulaimani, S. Tectonics of the Anti-Atlas of
Morocco.
C.R. Geoscience
North-western and West Central Africa
338, 1124 (2006).
3. Michard, A., Hoepffner, C., Soulaimani, A. & Baidder, L. The variscan belt. In
Continental Evolution:
the Geology of Morocco
(eds. A. Michard, O. Saddiqi, A. Chalouan & D. Frizon de Lamotte) 65132
(Lecture Notes in Earth Sciences Vol. 126, Springer, 2008).
4. Loi, A.
et al.
The Late Ordovician glacio-eustatic record from a highlatitude storm-dominated shelf
succession: the Bou Ingarf section (Anti-Atlas, Southern Morocco).
Palaeogeogr. Palaeoclim.
296, 332358 (2010).
5. Havlicek, V. Brachiopodes de l'Ordovicien du Maroc.
Notes Mm. Serv. Gol. Maroc
Palaeoeco.
230, 1135
(1971).
6. Sutcliffe, O.E., Harper, D.A.T., At Salem, A., Whittington, R.J. & Craig, J. The development of an
atypical Hirnantian brachiopod fauna and the onset of glaciation in the Late Ordovician of
Gondwana.
114 (2001).
7. Mergl, M. Faunal turnover near the Katian/Hirnantian boundary in the Prague Basin (Czech
Republic). In
Ordovician of the World
(eds J.C. Gutirrez-Marco, I. Rabano & D. Garcia-Bellido) 359
366 (Publ. Inst. Geol. Miner. Espana, Vol. 14, 2011).
8. Mitchell, C.E., torch, P., Holmden C., Melchin, M.J. & Gutirrrez-Marco, J.C. New stable isotope data
and fossils from the Hirnantian stage in Bohemia and Spain: Implications for correlation and
paleoclimate. In
Ordovician of the World
(eds J.C. Gutirrez-Marco, I. Rabano & D. Garcia-Bellido)
371378 (Publ. Inst. Geol. Miner. Espana, Vol. 14, 2011).
9. Paris, F. The Ordovician chitinozoan biozones of the northern Gondwana Domain.
Rev. Palaeobot.
Trans. Roy. Soc. Edinburgh: Earth Sci. 92,
66, 181209 (1990).
10. Bourahrouh, A., Paris, F. & Elaouad-Debbaj, Z. Biostratigraphy, bodiversity and palaeoenvironments
of the chitinozoans and associated palynomorphs from the Upper Ordovician of the Central Anti-
Atlas, Morocco.
Rev. Palaeobot. Palynol.
Palynol.
11. Clerc, S.
et al
. Subglacial to proglacial depositional environments in an Ordovician glacial tunnel
valley, Alnif, Morocco.
Palaeogeogr. Palaeoclim. Palaeoeco.
130, 1740 (2004).
370, 127144 (2013).
12. Ghienne, J.-F., Le Heron, D., Moreau, J., Denis, M. & Deynoux, M. The Late Ordovician glacial
sedimentary system of the North Gondwana platform. In
Glacial Sedimentary Processes and
Products
.) 295319 (International Association of Sedimentologists,
Special Publication, Blackwells, 2007).
13. Le Heron, D. Late Ordovician glacial record of the Anti-Atlas, Morocco.
Sediment. Geol.
Vol. 39 (eds M. Hambrey
et al
201, 93110
(2007).
14. Le Heron D., Ghienne J.-F., El Houicha M., Khoukhi Y. & Rubino J.-L. Maximum extent of ice sheets in
Morocco during the Late Ordovician glaciation.
Palaeogeogr. Palaeoclim. Palaeoeco.
245, 200226
49, 121200 (2000).
16. Nutz, A., Ghienne, J.-F. & torch, P. Circular, cryogenic structures from the Hirnantian deglaciation
sequence (Anti-Atlas,Morocco).
J. Sediment. Res.
(2007).
15. Lning, S., Craig, J., Loydell, D.K., torch, P. & Fitches, B. Lower Silurian 'hot shales' in North Africa
and Arabia: regional distribution and depositional model.
Earth-Sci. Rev.
83, 115131 (2013).
17. Long, D.G.F. Tempestite frequency curves: a key to Late Ordovician and Early Silurian subsidence,
sea-level change, and orbital forcing in the Anticosti foreland basin, Quebec, Canada.
Can. J. Earth
413431 (2007).
18. Desrochers, A., Farley. C., Achab, A., Asselin, E. & Riva, J.F. A far-field record of the end Ordovician
glaciation: The Ellis Bay Formation, Anticosti Island, Eastern Canada.
Palaeogeogr. Palaeoclim.
Sci. 44,
296, 248263 (2010).
Palaeoeco.
19
19. Achab, A., Asselin, E., Desrochers, A., Riva, J.F. & Farley, C. Chitinozoan biostratigraphy of a new
Upper Ordovician stratigraphic framework for Anticosti Island, Canada.
Geol. Soc. Am. Bull.
186205 (2011).
123,
109, 85-111 (2000).
21. Achab, A., Asselin, E., Desrochers, A. & Riva, J.F. The end-Ordovician chitinozoan zones of Anticosti
Island, Qubec: definition and stratigraphic position.
Rev. Palaeobot. Palynol.
20. Soufiane, A. & Achab, A. Chitinozoan zonation of the Late Ordovician and the Early Silurian of the
Island of Anticosti, Quebec, Canada:
Rev. Palaeobot. Palynol.
22. Delabroye, A.
et al
198, 92109 (2013).
. Phytoplankton dynamics across the Ordovician/Silurian boundary at low
palaeolatitudes: Correlations with carbon isotopic and glacial events.
Palaeogeogr. Palaeoclim.
7997 (2011).
23. Copper, P. Reefs during the multiple crises towards the Ordovician-Silurian boundary: Anticosti
Island, Eastern Canada, and worldwide.
Canad. J. Earth Sci.
Palaeoeco
.
312,
38, 153171 (2001).
24. Jin, J. & Copper, P. 2000. Late Ordovician and Early Silurian pentamerid brachiopods from Anticosti
Island, Qubec, Canada.
Palaeontograph. Canad.
, Vol. 18, 140pp.
25. Jin, J. & Copper, P. Origin and evolution of the Early Silurian (Rhuddanian) virgianid pentameride
brachiopods the extinction recovery fauna from Anticosti Island, eastern Canada.
Boll. Soc.
1 11 (2010).
26. Jin, J. & Copper, P. Parastrophinella (Brachiopoda): its paleogeographic significance at the
Ordovician/Silurian boundary.
J. Paleontol.
Paleontol. Ital. 49,
71, 369380 (1997).
27. Jin, J. & Copper, P. Response of brachiopod communities to environmental change during the Late
Ordovician mass extinction interval, Anticosti Island, eastern Canada.
Fossils and Strata
54, 4151
(2008).
28. Rongyu, L. & Copper, P. Early Silurian (Llandovery) orthide brachiopods from Anticosti Island,
eastern Canada: the O/S extinction recovery fauna.
Spec. Pap. Palaeont.
76, 171 (2006).
29. Jin, J. & Zhan, R. B. Late Ordovician orthide and billingsellide brachiopods from Anticosti Island,
Eastern Canada; diversity change through mass extinction.
National Research Council Research
Press
(Ottawa, Canada, 2008).
30. Nestor, H., Copper, P. & Stock, C.W. Late Ordovician and Early Silurian stromatoporoid sponges
from Anticosti Island, eastern Canada: crossing the O/S mass extinction boundary.
National
(Ottawa, 163 pp., 2010).
31. Melchin, M. J. Restudy of some OrdovicianSilurian boundary graptolites from Anticosti Island,
Canada, and their biostratigraphic significance.
Lethaia
Research Council Research Press
41, 5562 (2008).
32. Ausich, W.I. & Copper, P. The Crinoidea of Anticosti Island (Late Ordovician and Early Silurian).
Palaeontograph. Canad.
Vol. 29 (2010).
33. Holland, C.H. & Copper, P. Ordovician and Silurian nautiloid cephalopods from Anticosti island:
traject across the Ordovician-Silurian (O/S) mass extinction boundary.
Canad. J. Earth Sci.
45,
10151038 (2008).
34. Farley, C.
Sediment dynamics and stratigraphic architecture of a mixed carbonate-siliciclastic ramp:
the upper ordovician (Hirnantian) Ellis Bay Formation, Anticosti Island, Qubec, Canada
. Unpublished
M.Sc. thesis, University of Ottawa, Canada (2008).
35. Wickson, S.
High-Resolution Carbon Isotope Stratigraphy of the Ordovician-Silurian Boundary on
Anticosti Island, Quebec: Regional and Global Implications
. Unpublished M.Sc. thesis, University of
Ottawa, Canada (2010).
36. Pinet, N. Hinterland-directed transtensional faulting at an orogen structural front: the example of
the Cap-Chat Mlange, Quebec Appalachians.
Geol. Soc. Amer.
123, 22562265 (2011).
Bull.
37. Ainsaar, L.,
et al.
2010. Middle and Upper Ordovician carbon isotope chemostratigraphy in
Baltoscandia: A correlation standard and clues to environmental history.
Palaeogeogr. Palaeoclim.
294, 189201 (2010).
Palaeoeco.
38. Brenchley, P.J.
et al
. High-resolution stable isotope stratigraphy of Upper Ordovician sequences:
Constraints on the timing of bioevents and environmental changes associated with mass extinction
and glaciation.
Geol. Soc. Am. Bull.
115, 89104 (2003).
20
39. Bergstrm, S.M., Saltzman, M.R. & Schmitz, B. First record of the Hirnantian (Upper Ordovician) 13C excursion in the North American Midcontinent and its regional implications.
Geol. Mag. 143,
678 (2006).
40. Bergstrm, S.M., Kleffner, M., Schmitz, B. & Cramer, B. Revision of the position of the Ordovician
Silurian boundary in southern Ontario: regional chronostratigraphic implications of 13C chemostratigraphy of the Manitoulin Formation and associated strata.
Canad. J. Earth Sci.
657
48, 1447
1470 (2011).
41. Jones, D. S.
et al
Bull.
123, 16451664 (2011).
42. Jervey, M.T. Quantitative geological modeling of siliciclastic rock sequences and their seismic
expression. In
Sea Level Changes An Integrated Approach
. Terminal Ordovician carbon isotope stratigraphy and glacioeustatic sea-level
change across Anticosti Island (Qubec, Canada).
Geol. Soc. Amer.
(Society of Economic Paleontologists
and Mineralogists, Special Publication, 1988).
Vol. 42 (eds Wilgus, C.K.
et al
.) 4769
43. Posamentier, H.W. & Vail, P.R. Eustatic controls on clastic deposition II sequence and systems
tract models. In
Sea-Level Changes An Integrated Approach
(SEPM Special Publication, 1988).
44. Mountain, G.S.
et al.
Vol. 42 (eds C.K. Wilgus et al.) 125154
record. In
Continental Margin Sedimentation: From Sediment Transport to Sequence Stratigraph
The long-term stratigraphic record on continental margins - The long-term
37
(eds C.A. Nittrouer
et al.)
Vol.
381-458 (International Association of Sedimentologists, Special
Publication, Blackwells, 2007).
45. Meijer, X.D., Postma, G., Burrough, P.A. & De Boer P.L. Modelling the preservation of sedimentary
deposits on passive continental margins during glacial-interglacial cycles. In
Analogue and
Numerical Modelling of Sedimentary Systems: From Understanding to Prediction
. Vol. 40 (eds P. De
Boer et al.) 223-238 (International Association of Sedimentologists, Special Publication, Blackwells,
2008).
46. Grader, G.W., Isaacson, P.E., Daz-Martnez, E. & Pope, M.C. Pennsylvanian and Permian sequences
in Bolivia: Direct responses to Gondwana glaciation. In
Resolving the Late Paleozoic Ice Age in Time
and Space
Vol. 441 (eds C.R. Fielding, T.D. Frank & J.L. Isbell) 143159 (Geological Society of
America, Special Paper, 2008).
47. Heckel, P.H. Pennsylvanian cyclothems in Midcontinent North America as far-field effects of waxing
and waning of Gondwana ice sheets. In
Resolving the Late Paleozoic Ice Age in Time and Space
Vol.
441 (eds C.R. Fielding, T.D. Frank & J.L. Isbell) 275289 (Geological Society of America, Special
Paper, 2008).
48. ten Brink, U.S. & Schneider, C. Glacial morphology and depositional sequences of the Antarctic
continental shelf.
Geology
23, 580-584 (1995).
49. Pollard, D. & DeConto, R.M. A coupled ice-sheet/ice-shelf/sediment model applied to a marine-
margin flowline: forced and unforced variations. In
Glacial Sedimentary Processes and Products
39 (eds M. Hambrey
et al
Vol.
.) 3752 (International Association of Sedimentologists, Special
Publication, Blackwells, 2007).
50. Bart, P.J. & Anderson J.B. Seismic expression of depositional sequences associated with expansion
and contraction of ice sheets on the northwestern Antarctic Peninsula continental shelf. In
Geology
of Siliciclastic Shelf Seas
Vol. 171 (eds M. de Batist & P. Jacobs) 171186 (Geol. Soc. Spec. Publs,
London, 1996).
51. Denis, M., Guiraud, M., Konat, M. & Buoncristiani, J.-F. Subglacial deformation and water-pressure
cycles as a key for understanding ice stream dynamics: evidence from the Late Ordovician
succession of the Djado Basin (Niger).
Intern. J. Earth Sci.
223, 100125 (2010).
53. Helland-Hansen, W. & Hampson, G. Trajectory analysis: concepts and applications.
Basin Res.
99, 1399 1425 (2010).
52. Le Heron, D.P., Armstrong, H.A., Wilson, C., Howard, J.P. & Gindre, L. Glaciation and deglaciation of
the Libyan Desert: The Late Ordovician record.
Sediment. Geol.
21,
454483 (2009).
54. Le Heron, D.P. & Craig, J. First-order reconstructions of a Late Ordovician Saharan ice sheet.
J. Geol.
165, 1929 (2008).
55. Le Heron, D.P., Craig, J., Sutcliffe, O. & Whittington, R. Late Ordovician glaciogenic reservoir
heterogeneity: an example from the Murzuq Basin Libya.
Marine Petrol. Geol.
Soc.
23, 655677 (2006)
21
(2010).
57. Girard, F., Ghienne, J.-F. & Rubino, J.-L. Occurrence of hyperpycnal flows and hybrid event beds
related to glacial outburst vents in a Late Ordovician proglacial delta (Murzuq Basin, SW Libya).
J.
56. Videt. B. et al. (2010) Biostratigraphical calibration of the third-order Ordovician sequences on the
northern Gondwana platform.
Palaeogeogr. Palaeoclim. Palaeoeco. 298,
359375
82, 688708 (2012).
58. Ghienne, J.-F., Moreau, J., Degermann, L. & Rubino, J.-L. Lower Palaeozoic unconformities in an
intracratonic platform setting: glacial erosion versustectonics in the eastern Murzuq Basin
(southern Libya).
Sediment. Res.
455482 (2013).
59. Deschamps, R., Eschard, R. & Rousse, S. Architecture of Late Ordovician glacial valleys in the Tassili
NAjjer area (Algeria).
Sediment. Geol.
Intern. J. Earth Sci. 102,
60. Mackenzie, F.T. & Lerman, A.
Carbon in the geobiosphere: in Earths Outer Shell
289, 124147 (2013).
2006).
61. Cramer, B. & Saltzman, M.R. Sequestration of 12C in the deep ocean during the early Wenlock
(Silurian) positive carbon isotope excursion.
Palaeogeogr. Palaeoclim. Palaeoeco.
(Springer, London,
219, 333349
(2005).
62. Cramer, B.S., & Saltzman, M.R. Fluctuations in epeiric sea carbonate production during Silurian
positive carbon isotope excursions: A review of proposed paleoceanographic models.
Palaeogeogr.
245, 3745 (2007).
63. de Boer, P.L. & Alexandre, J.T. Orbitally forced sedimentary rhythms in the stratigraphic record: is
there room for tidal forcing?
Sedimentology
Palaeoclim. Palaeoeco.
64. Brenchley, P.J.
et al
59, 379392 (2012).
. Bathymetric and isotopic evidence for a shortlived Late Ordovician glaciation in
a greenhouse period.
Geology
65. Kump, L.R.
et al
. A weathering hypothesis for glaciation at high atmospheric
p
22, 295298 (1994).
Ordovician.
Palaeogeogr. Palaeoclim. Palaeoeco
CO2 during the Late
152, 173187 (1999).
66. Melchin, M.J. & Holmden, C. Carbon isotope chemostratigraphy in Arctic Canada: Sea-level forcing of
carbonate platform weathering and implications for Hirnantian global correlation.
Palaeogeogr.
234, 186200 (2006).
Palaeoclim. Palaeoeco.
67. LaPorte, D.F.
et al
. Local and global perspectives on carbon and nitrogen cycling during the
Hirnantian glaciation.
Palaeogeogr. Palaeoclim. Palaeoeco
68. Veizer, J.
et al.
87 Sr/86Sr, 13C and 18O evolution of Phanerozoic seawater.
276, 182195 (2009).
5988
Chem. Geol. 161,
(1999).
69. Bickert, T., Ptzold, J., Samtleben, C. & Munnecke, A. Paleoenvironmental changes in the Silurian
indicated by stable isotopes in brachiopod shells from Gotland, Sweden.
Geochim. Cosmochim. Acta
61, 27172730 (1997).
70. Jeppson, L. & Calner, M. The Silurian Mulde Event and a scenario for a secundo-secundo events.
135154 (2003).
71. Munnecke, A., Samtleben, C. & Bickert, T. The Ireveken event in the lower Silurian of Gotland,
Sweden relation to similar Palaeozoic and Proterozoic events.
Palaeogeogr. Palaeoclim.
Trans. Roy. Soc. Edinburgh, Earth Sci. 93,
195, 99124 (2003).
Palaeoeco.
22
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer
Copyright Nature Publishing Group Sep 2014
Abstract
The end-Ordovician was an enigmatic interval in the Phanerozoic, known for massive glaciation potentially at elevated CO2 levels, biogeochemical cycle disruptions recorded as large isotope anomalies and a devastating extinction event. Ice-sheet volumes claimed to be twice those of the Last Glacial Maximum paradoxically coincided with oceans as warm as today. Here we argue that some of these remarkable claims arise from undersampling of incomplete geological sections that led to apparent temporal correlations within the relatively coarse resolution capability of Palaeozoic biochronostratigraphy. We examine exceptionally complete sedimentary records from two, low and high, palaeolatitude settings. Their correlation framework reveals a Cenozoic-style scenario including three main glacial cycles and higher-order phenomena. This necessitates revision of mechanisms for the end-Ordovician events, as the first extinction is tied to an early phase of melting, not to initial cooling, and the largest δ13 C excursion occurs during final deglaciation, not at the glacial apex.
You have requested "on-the-fly" machine translation of selected content from our databases. This functionality is provided solely for your convenience and is in no way intended to replace human translation. Show full disclaimer
Neither ProQuest nor its licensors make any representations or warranties with respect to the translations. The translations are automatically generated "AS IS" and "AS AVAILABLE" and are not retained in our systems. PROQUEST AND ITS LICENSORS SPECIFICALLY DISCLAIM ANY AND ALL EXPRESS OR IMPLIED WARRANTIES, INCLUDING WITHOUT LIMITATION, ANY WARRANTIES FOR AVAILABILITY, ACCURACY, TIMELINESS, COMPLETENESS, NON-INFRINGMENT, MERCHANTABILITY OR FITNESS FOR A PARTICULAR PURPOSE. Your use of the translations is subject to all use restrictions contained in your Electronic Products License Agreement and by using the translation functionality you agree to forgo any and all claims against ProQuest or its licensors for your use of the translation functionality and any output derived there from. Hide full disclaimer