The structure, evolution and symmetry of the magma

Tectonophysics 468 (2009) 6–27
Contents lists available at ScienceDirect
Tectonophysics
j o u r n a l h o m e p a g e : w w w. e l s e v i e r. c o m / l o c a t e / t e c t o
The structure, evolution and symmetry of the magma-poor rifted margins of
the North and Central Atlantic: A synthesis
T.J. Reston
University of Birmingham, School of Geography, Earth and Environmental Sciences, Birmingham, United Kingdom
a r t i c l e
i n f o
Article history:
Received 31 August 2007
Accepted 3 September 2008
Available online 16 September 2008
Keywords:
Rifted margins
Rheological evolution
Symmetry of rifting
Serpentinization
Extension
a b s t r a c t
Magma-poor rifted margins consistently show extreme crustal thinning accompanied by normal faulting, the
serpentinization of the mantle beneath crust thinned to less than 8 ± 2 km, and the unroofing of a broad zone
of mantle within the continent–ocean transition, accompanied by the development of detachment and other
large-offset faults. These observations are the logical result of the progressive extension of cool lithosphere
away from thermal anomalies such as plumes. Although the paucity of magmatism may be explained by
depth-dependent extension of the lithosphere, and by pre-depleted sub-crustal lithosphere, rifting above
initially cool sub-lithospheric mantle also may explain the subsidence deficit observed at some margins:
synrift subsidence is buffered by the simultaneous influx of warmer oceanic asthenosphere. As it extends,
thins and cools, ductile creeping layers in the mid- and deep crust become progressively more brittle,
resulting in increased coupling between the upper and lower crust, and eventually the embrittlement of the
entire crust, faults cutting from the surface across the Moho, bringing water into the mantle and causing its
serpentinization. Increased coupling and the development of serpentine detachments predict the
development of late-stage asymmetry once the entire crust is brittle. Such detachments are imaged on
some margins and inferred on others; analysis of the crustal structure across conjugate margins shows that
these are approximately symmetric until this late stage when they become markedly asymmetric. Similar
analyses show that depth-dependent stretching of the crust is insufficient to explain the discrepancy
between the amount of the visible extension along faults and the amount of crustal thinning. Instead this
“extension discrepancy” may be related to the complex evolution of brittle deformation through multiple
phases and styles of faulting, related to the changes in the rheological character and strength of the
lithosphere as it is thinned. Complex polyphase faulting continues after complete crustal separation,
resulting in the exhumation of broad expanses of peridotitic basement, the top of which is everywhere
marked by an exhumed slip surface, similar to the corrugated surface observed at mid-ocean ridges. The
similarity in processes between mantle unroofing and seafloor spreading makes the distinction between the
COT and true oceanic crust difficult and possibly moot.
© 2008 Elsevier B.V. All rights reserved.
1. Introduction
The first order tectonic process of continental breakup initiates the
plate tectonic cycle of plate creation and destruction. Rifted margins are
the trailing edges of the continents that develop as the continents are
rifted apart (Fig. 1). In addition to their economic importance (hydrocarbons) and hazard potential (especially slope failure), they provide a
record of the processes that accompanied continental breakup.
In the past, rifted margins have been classified as volcanic or nonvolcanic margins. These terms are a bit misleading as even nonvolcanic margins exhibit magmatism, so in this paper I will prefer the
term “magma-poor” and “magma-dominated margin” (Sawyer et al.,
2007), depending whether less or more magmatism is observed than
might be expected during the rifting above normal asthenosphere.
E-mail address: [email protected].
0040-1951/$ – see front matter © 2008 Elsevier B.V. All rights reserved.
doi:10.1016/j.tecto.2008.09.002
This is undepleted mantle with a potential temperature of 1300 °C
(potential temperature Tp is the temperature it would have if brought
rapidly to the surface). During rifting, the upwelling of such mantle
leads to pressure release melting, producing 6–7 km for very rapid
rifting and hence upwelling, but perhaps 3–4 km for typical rift
durations (Minshull et al., 2001). Although some margins (e.g. West
Iberia) are so devoid of melt products that it is very likely that little
melting took place during rifting, at other margins the difficulties in
determining both the thickness of igneous addition to the lithosphere
(some may be trapped below the Moho) and precise rift durations
means that a practical definition of magma-poor margins may be
those where tectonic rather than magmatic processes dominate
during rifting. Judging by the magma-poor west Iberia margin,
tectonic rather than magmatic processes dominate when there has
been less than about 2–3 km of igneous addition prior to crustal
breakup. Note that by restricting the definition to what happens
during rifting, several margins and basins that had little magmatism
T.J. Reston / Tectonophysics 468 (2009) 6–27
7
Fig. 1. Top: Worldwide distribution of different margin types (modified after Boillot and Coulon, 1998) and the locations of the margins discussed in this paper (bold). Magma-poor, magma-dominated and transform are all approximately
equally numerous. Below: Cartoon sections summarising the architecture of magma-poor (based on observations made in this paper), magma-dominated (modified after Gernigon et al., 2004) and transform (based on Edwards et al., 1997)
margins. This paper investigates the first category.
8
T.J. Reston / Tectonophysics 468 (2009) 6–27
until breakup, although subsequently subjected to voluminous postrift
magmatism (e.g. Rockall Trough and the Norwegian Vøring (Gernigon et
al., 2004) and Møre Basins), have much in common with magma-poor
margins.
Other rifted margins (e.g. Woodlark Basin) appear intermediate
between the end-member “magma-poor” and “magma-dominated”
margins as they may exhibit approximately the amount of magmatism
expected for rapid mantle upwelling and do not show the extreme
tectonism of magma-poor margins described below. However in several
cases, the margin is simply not well enough characterised at synrift
levels to determine the mechanics of breakup. Finally, transform
margins (e.g., Edwards et al., 1997) form a completely separate class of
margin formed during continental breakup, generally exhibiting little
magmatism, but marked by abrupt change in crustal thickness.
This paper is concerned with the structure and tectonics of
magma-poor margins. By concentrating on magma-poor margins, it is
possible to avoid some of the geological and geophysical complexities
caused by magmatic addition both within the crust (affecting crustal
thickness, crustal rheology, and seismic velocity) and on top of the
crust as lava flows (causing imaging problems and thus obscuring the
tectonic structures beneath). As such it is intended that the paper
provide a counterpoint to the more magmatic theme of many papers
in this volume. In this paper, I combine observations of the structure of
the 22 best studied magma-poor margins of the North and Central
Atlantic with the predictions of the fundamental processes accompanying lithospheric extension and breakup.
2.
2. Magma-poor margins of the North and Central Atlantic and
neighbouring areas
Magma-poor margins are found in every ocean (Atlantic, Indian,
Southern, Pacific, Arctic), but those in the North and Central Atlantic
have been most intensively studied, are well constrained by
geophysical (both reflection and refraction) and ODP data, and for
the most part are unambiguously magma-poor (Figs. 2 and 3). Moving
from north to south, examples include parts of Labrador–West
Greenland (Chian et al., 1995; Chalmers and Pulvertaft, 2001), Flemish
Cap–Goban Spur (Keen et al., 1989; Bullock and Minshull, 2005), North
Biscay (Thinon et al., 2003), Flemish Cap–Galicia Bank (Funck et al.,
2003; Zelt et al., 2003; Hopper et al., 2004; Reston et al., 2007),
Newfoundland Basin–South Iberia Abyssal Plain (Krawczyk et al.,
1996; Pickup et al., 1996; Chian et al., 1999; Dean et al., 2001;
Whitmarsh et al., 2001; Van Avendonk et al., 2006), S Newfoundland
Basin–Tagus (Lau et al., 2006) and Nova Scotia–Morocco (Funck et al.,
2004; Contrucci et al., 2004; Wu et al., 2006; Maillard et al., 2006). In
addition are a number of failed rifts bordering the Atlantic where
rifting progressed so far that crustal separation may have taken place
or was at least imminent. These include Rockall Trough (O'Reilly et al.,
1996; Pérez-Gussinyé et al., 2001), Porcupine Basin (Reston et al., 2001,
2004; O'Reilly et al., 2006) and similar basins in the Mediterranean Sea:
the Ligurian Sea (Contrucci et al., 2001), the Tyrrhenian Sea, the
Alboran Sea, the South Balearic Basin and the Valencia Trough. Finally,
one example from further south is included: the Rio Muni (Equatorial
Guinea) margin (Wilson et al., 2003; Turner et al., 2003). Other margins
which appear very similar to the above either are not well enough
characterized geophysically (South Australia) or have thick sedimentary
sequences that obscure the deeper structure so that it is not always
totally clear whether they are magma-poor or not (e.g., Angola margin
in the South Atlantic).
The majority of magma-poor margins from the North and Central
Atlantic display a number of common features (Figs. 2 and 3; Table 1):
1. Extreme crustal thinning from ∼ 30 km to a few km in most cases
over a distance of 100–200 km. On most margins well-defined fault
blocks are imaged, but the amount of extension associated with
these is generally less than required to explain the crustal thinning.
3.
4.
5.
This is the so-called extension discrepancy, discussed below. The
thinning occurs most rapidly in a necking region between c 20 and
10 km, corresponding to β factors (McKenzie, 1978) between 1.5
and 3 and (1 − 1 / β) values of 0.33 to 0.67 (Fig. 4): most margins are
characterised by a broad shelf, a steep continental slope and an
extensive continental rise of very thin crust. Extension of the brittle
upper crust is accommodated by faulting, and summing the fault
heaves provides an estimate of the amount of extension. Other
parts of the lithosphere may initially deform by other means
(ductile flow or creep), although there is evidence that deformation
in parts of the lower crust and uppermost mantle is localised into
shear zones (Reston, 1988, 1993). However, as discussed below, the
deformation mechanics of the lithosphere evolve during extension
and thinning, as the initial rheological properties are modified by
strain, and by changes in temperature and pressure (PérezGussinyé and Reston, 2001).
A zone of anomalous basement between the last identifiable
continental crust and the first true oceanic crust. This unit typically
has a moderate velocity gradient (∼0.2 s− 1), passing down from
5 km/s to close to 8 km/s and weak magnetic anomalies. The unit
has been sampled by dredging, by submersible (Boillot et al., 1988)
and by drilling on the West Iberian margin (sites 637, 897, 899, 1068
and 1070 — Boillot and Winterer, 1988; Whitmarsh et al., 1996,
2001), by drilling off the Newfoundland margin (Tucholke et al.,
2007) and by dredging to the south of Australia (Nicholls et al.,
1981), and has been found to be partially serpentinized peridotites.
The anomalous crust is thus thought to be exhumed mantle, with
downward decreasing degrees of serpentinization explaining
the velocity gradient. Elsewhere, it has been identified by similar
geophysical characteristics (velocity), a lack of oceanic spreading
anomalies, and simply a non-continental appearance on the seismic reflection images.
Beneath the crust where it had thinned during rifting to less than
about 8 km, zones with velocity intermediate between those of the
crust and mantle (i.e. between 7 and 7.8 km/s) and a velocity
gradient of close to 0.1 s− 1. In most cases, the zones can be traced
laterally into regions of unroofed mantle (see above) and so are
interpreted as partially serpentinized peridotites, “undercrusting”
(Boillot et al., 1989) the crust. Where they cannot be traced laterally,
they might alternatively have been interpreted as mafic underplate
or a zone of mantle intruded with mafic intrusions, if it were not for
the paucity of magmatism at most of these margins. Furthermore,
the velocity gradient noted above is consistent with downward
decreasing degrees of serpentinization but less easily explained by
mafic underplate. It is however important to realise that nowhere
has serpentinized mantle actually been sampled beneath the
thinned continental crust.
At some margins, the unroofed mantle is covered with low velocity
and locally mounded units up to 3 km thick. As they are generally
observed atop the exhumed mantle and have seismic velocity
consistent with highly serpentinized peridotites with varying degrees of porosity, they are generally interpreted as such (Tucholke
et al., 2007): serpentine breccias and sedimentary serpentine have
been sampled at ODP sites 897 and 899 in the southern Iberia
Abyssal Plain (Whitmarsh et al., 1996). The mounded pattern may
reflect either lateral transport, such as the collapse of serpentinite
highs, or emplacement from below (Reston et al., 2001, 2004).
At some margins, the boundary between highly thinned continental crust and the underlying zone interpreted as serpentinized
mantle is marked by a sharp reflector, most clearly imaged on
depth images produced by prestack depth migration. Examples
include the S reflector west of Galicia (de Charpal et al., 1978;
Hoffmann and Reston, 1992), the H reflector in the Iberia Abyssal
Plain (Krawczyk et al., 1996), and the P reflection beneath the
Porcupine Basin (Reston et al., 2001, 2004). In each case, the
reflector does not only appear to separate crustal rocks from the
T.J. Reston / Tectonophysics 468 (2009) 6–27
Fig. 2. Sections through 14 conjugate or near-conjugate margins (including both sides of two basins) best constrained by combined seismic reflection and wide-angle data, all shown at the same scale, with no vertical exaggeration. Key velocity
boundaries and contours are marked; those in bold are used for the analyses in Figs. 8, 9 and 13. Also marked is the zone where the crust thins to the value at which it should have become brittle during rifting. This always occurs landward of the
first serpentinized peridotites that underlie the thinned crust of the deep margin.
9
10
T.J. Reston / Tectonophysics 468 (2009) 6–27
Fig. 3. Sections (no vertical exaggeration) through a further 8 margins, either non-conjugate or where full velocity control to the shelf is lacking. In all cases, the landward limit of the
serpentinites beneath the lies slightly oceanward of point at which the crust is thin enough to have become completely brittle during rifting.
underlying serpentinized mantle, but also to act as a detachment
on which the overlying fault blocks ride. Analysis of the physical
properties of the H and S reflection (Reston, 1996) have shown that
they are compatible with a relatively sharp boundary between
fractured crustal rocks and partially serpentinized peridotites,
perhaps supplemented by a thin zone of reduced velocity at the
detachment itself (Leythaeuser et al., 2005).
As many magma-poor margins display the same general characteristics, it seems likely that they also share a common mode of
formation. In the following sections I examine the processes that may
have led to the development of these characteristics at magma-poor
rifted margins and derive a generic model for the formation of
magma-poor margins, focussing on the large-scale structure of the
margins.
3. Large-scale lithospheric deformation — subsidence and melt
generation
Continental breakup occurs through lithospheric extension, thinning
and eventual rupture. Lithospheric stretching occurs at all levels, so that
the crust is highly extended, and separates during continental breakup.
Lithospheric thinning allows the underlying asthenospheric mantle to
well up towards the surface, undergoing partial pressure release melting
according to the amount and rate of upwelling, coupled with the
temperature and fertility of the asthenosphere (Bown and White, 1995).
The thinning of the crust causes subsidence, which is partly offset during
and shortly after rifting by the thermal uplift provided by the thinning of
the entire lithosphere (McKenzie, 1978) and the compression of the
lithospheric isotherms, and by the upwelling of the asthenosphere. As a
result, the distribution of both subsidence and magmatism in time and
space provides constraints on the amount of extension and thinning of
the different lithospheric levels.
A magma-poor margin is defined by a lack of significant magmatism. Excess melting at magma-dominated margins (Nielsen and
Hopper, 2002) is generally explained by a combination of high temperature asthenosphere (e.g. “hotspots”), secondary convection and
chemically enriched mantle (including the presence of water). Conversely, the melt deficiency of magma-poor margins may be explained
(Fig. 5) by a cool (Tp below 1300 °C — Reston and Phipps Morgan,
2004), or partially depleted sub-lithospheric mantle (Pérez-Gussinyé
et al., 2006), or one that does not rise completely due to incomplete
lithospheric separation, requiring a component of lithospheric-scale
depth-dependent stretching (DDS) if crustal separation has occurred
(Minshull et al., 2001). Higher viscosities resulting from cooler temperatures and lack of water are also likely to hinder secondary convection. These different explanations have different implications for
the subsidence history of the rifted margin.
Traditionally, the large-scale deformation of the lithosphere
beneath rifts has been determined from subsidence patterns (e.g.
White and McKenzie, 1988). In considering subsidence, lithospheric
thinning is usually considered in two parts: the thinning of the crust
and the thinning of the lithosphere. These generally occur at the same
time but may have a different spatial distribution. The thinning of the
whole lithosphere causes both a compression of the isotherms within
the lithosphere (note in general no part of the lithosphere actually
gets hotter, but hot rocks are brought towards the surface) and the
replacement of higher density mantle lithosphere by lower density
asthenosphere. These combine to cause a thermal uplift which decays with time after rifting as the isotherms revert to horizontal
(meaning that the uplifted and upwelled lithosphere and asthenosphere cool) giving a subsequent gradual thermal subsidence. The
isostatic response to crustal thinning is subsidence: as crustal thinning
is permanent it controls the eventual total amount of subsidence.
In principle, the subsidence at any time is that caused by crustal
thinning minus the effect of any residual thermal uplift. Thus the
Table 1
Geophysical characteristics of rifted margins characterized by mantle serpentinites and their conjugates
MCS Refraction DSDP, Final rift
ODP
duration
[m.y.]
References Underrusting — Thickness Gradient
vels [km/s]
[km]
[s− 1]
approximate
Extent Max
[km]
thickness
overlying
crust [km]
Beta factor
at landward
limit of
serpentinites
Unroofed
mantle max
thickness [km]
Velocity
[km/s]
Gradient
[s− 1]
Extent Highly
Max
[km]
serpentd
thickness
peridotites — [km]
vels km/s
Extent References
[km]
Y
OBS
30–50
6.4–7.6
b 4.5
0.27
50
4–6.5
5–7.5
7
6–7.6
0.23
30–66
∼ 4.8
2
35
W Greenland
Y
OBS
30–50
7.4–8.0
b4
0.15
54
9
3.5
7
6.2–7.4
0.17
54
∼ 4.5
3
35
Rockall
Y
OBS
40–70
7.5–7.8
3
0.1
30?
10–11.5
∼3
15
6–6.9
7.5–7.8
0.15.03
N150
∼ 4.5
2
N150?
Rockall II
(cont crust)
Porcupine
Y
OBS
40–70
7.5–7.8
3–10
0.1–0.03
200
10–11.5
∼3
N
–
–
–
N
–
–
Y
OBS
25–45
7.2–7.8
∼ 10
0.06
75
8–10
3.75–3
?
?
?
30
4.5
b3
30
Goban Spur
Y
OBS
30
7.0–7.5
b2
0.25
10
7
4.3
4.5
6–7
0.22
b39
4.5
2
b30
Armorican
Y
OBS
30
7.4–7.5
2
0.025
10
4
8
4
7.4–7.5
0.025
20
2
20
Ligurian Sea
Y
ESP
14
7.2–7.3
3
–
b50
6–9
5.3–3.6
N
–
–
–
4.6–5.0
(3C, 3D)
N
–
–
SE Flemish Cap
Y
OBS
15–25
7–6–8.0
4
0.1
24
4.5.6.5
6.5–4.5
?
?
?
–
Galicia Bank
IAP–ODP
149–173. LG12,
CAM144
Y
Y
OBS
OBS
103
149–
173
15–25
5–20
7.0–7.6
?
4
4
0.15
50
?
5–6
6.5
5–6
5–4.3
6
10
5.6–7.6
?
0.33
N10
∼100
–
∼4
–
b 3gradient
–
∼100
N Newf Basin
Y
OBS
210
5–20
–
–
–
.6
3
10
4
5–7.5
0.62
N20
N
–
–
S IAP–IAM9
S Newf Basin
Y
Y
OBS
OBS
(173)
5–20
10–25
N
7.6–8
–
6
–
0.07
–
95
–
8
–
3.4–4.3
b6.5
2.8
2.5–5
4–5
N 4.4
2
0.5gradient
100
30?
N
?
Y
OBS
OBS
OBS
10–25
b 55
7.2–7.6
7.2–7.6
N
5
4.5
0.08
0.09
25
70?
7
6.5
4–5
4–5.5
6.5
b6
N
0.3
1.33
0.5
.16–.08
0.08
N 0.07
–
N170
50
S Newf Basin
NovaScotia 1
Morocco
6–7.9
4.4–6.4
6.4–7.8
7.6–8.0
7.2–7.7
7.2–7.6
–
N20
60?
–
4.5–5
5.1
–
2–3
b2
100
60?
Nova Scotia 2
Rio Muni
?
Y
OBS
Grav
b 55
10–15
7.6–7.8
n/a
b 3.5
0.06
n/a
50
6.5
8
4–5.5
4
N
Y?
–
n/a
–
n/a
–
n/a
N? 5?
N
3?
–
35?
–
Chalmers +
Pulvertaft, 2001;
Chian et al., 1995,
Chian et al., 1995;
Reston + PérezGussinyé, 2007
Perez-Gussinye
et al., 2001;
O'Reilly et al.,
1996
O'Reilly et al.,
1996
O'Reilly et al.,
2006; Reston
et al., 2004;
Bullock and
Minshull (2005)
Thinon et al.
(2003)
Contrucci et al.
(2001)
Hopper et al.
(2004)
Zelt et al. (2003)
Krawczyk et al.,
1996; Wilson
et al., 2001;
Chian et al., 1999
Van Avendonk
et al. (2006)
Dean et al. (2001)
Lau et al.,
2006a; 2006b
T.J. Reston / Tectonophysics 468 (2009) 6–27
Labrador
Reid et al. (1994)
Funck et al.,
Maillard et al.,
2006;
Wu et al. (2006)
Wilson et al.
(2003)
Most margins are characterised by, beneath the feather edge of the continental crust, a layer of reduced mantle velocity, which passes laterally into a region of unroofed partially serpentinized mantle, in places overlain by a layer of highly
serpentinized peridotites.
11
12
T.J. Reston / Tectonophysics 468 (2009) 6–27
Fig. 4. Whole crustal thinning factors vs distance from the edge of the continental crust for the conjugate margins in Fig. 2. Some margins appear to thin rapidly and might be termed
“hard” margins (Davison, 1997), others thin more gradually and might be termed “soft”. However, there is no clear pattern of hard–soft asymmetry: some conjugate pairs are hard–
hard, some soft–soft and some hard–soft.
spatial distribution of crustal thinning controls the total tectonic
subsidence, whereas the spatial distribution of whole lithospheric
strain controls the distribution of that subsidence in time. One classic
example is the interpretation of postrift onlap onto prerift sediments
and basement as evidence that mantle extension and thinning is
distributed over a wider area than that of the crust (White and
McKenzie, 1988), consistent with the presence of outward dipping
mantle shear zones, both observed on deep seismic profiles (Reston,
1993) and predicted by numerical modelling (Harry and Sawyer, 1992;
Huismans and Beaumont, 2002). Such models of lithospheric depthdependent stretching are an effective way of suppressing melt
generation: where the crust thins more than does the lithospheric
mantle, the amount of melt generated will be less than expected for
the observed crustal thinning (Fig. 5).
Increasingly however it is being reported that rifted margins
exhibit a subsidence deficit: the deep margin subsided less during
rifting than might be expected from the observed crustal thinning (e.g.
Kusznir and Karner, 2007). This contradicts the DDS explanation for
the steer's head geometry and for the suppression of melting as a
synrift subsidence deficit would imply excessive thermal uplift and
hence more lithospheric than crustal thinning. Although, the
subsidence deficit can potentially be explained by major crustal
DDS, problems with such models are discussed below.
It may be possible to explain both the subsidence deficit and the
lack of voluminous magmatism at rifted margins if sub-lithospheric
mantle is cooler beneath the continents than the oceans, as implied by
several studies of mantle convection (e.g. Lenardic and Moresi, 2000;
Phipps Morgan et al., 1995). In support of these modelling studies,
independent evidence for cool sub-lithospheric mantle beneath
continents comes from a variety of sources, including cratonic
geotherms constructed from kimberlites and more general geotherms
deduced from seismic velocities (Reston and Phipps Morgan, 2004).
Even a slightly cool sub-lithospheric mantle (e.g. potential temperature Tp of 1200–1250 °C) will first of all lead to major melt suppression
in the same way that a “hotspot” might lead to an increase in melt
production (Fig. 5). However a ubiquitous cool sub-lithospheric mantle
cannot explain the observed 6–7 km thickness of the magmatic crust in
the oceans, which requires “normal” asthenospheric Tp of ∼ 1300 °C.
Reston and Phipps Morgan (2004) explained the switch from
suppressed to normal (oceanic) melt production through the influx
of oceanic asthenosphere (1300 °C Tp) beneath a cool continental rift
and pointed out that such an influx would also lead to thermal uplift
during the late stages of rifting leading to breakup. If the influx
Fig. 5. Melt thickness vs rift duration for mantle potential temperatures of 1300, 1250
and 1200 °C (after Minshull et al., 2001; Reston, 2007b). Magma-poor margins may be
defined as those that plot below the 1300 °C contour, due to the presence of cool or
depleted sub-lithospheric mantle, or by lithospheric depth-dependent stretching
(DDS), or as those where there is less than 2–3 km of igneous addition (shaded).
T.J. Reston / Tectonophysics 468 (2009) 6–27
occurred gradually as rifting progressed, the associated thermal uplift
would maintain the thinning crust at a higher level than might
otherwise be expected, appearing as a subsidence deficit. As the
McKenzie (1978) model assumes a constant Tp asthenosphere, this
explanation for the subsidence deficit may have been overlooked.
In general, detailed analysis of the subsidence patterns to determine
how the lithosphere was thinned is problematic. Where the margin is
well-supplied with sediment throughout its evolution, the deepest
synrift units are so deeply buried that they cannot be sampled; where
such margins are sediment-starved, most of the sediments are deposited
in deep water where the paleobathymetric controls are so poor that
meaningful estimates of tectonic subsidence after back-stripping are
impossible. Furthermore, the most reliable indicators of paleo-water
depth (e.g. erosional unconformities) may record not so much the
pattern of lithospheric deformation as the thermal structure of the
underlying mantle and specifically the passage of thermal anomalies
such as mantle plumes. Finally, it may in some circumstances be
misleading to relate paleobathymetry to global sea-level: during the
opening of the ocean basins (and the development of rifted margins)
isolated deep basins may develop in which the local sea-level may be far
below global. Such restricted basins may be characterized by lacustrine
to hypersaline conditions, leading to potentially important lacustrine
deposits (including potential source rocks) and to salt, depending on
climate and on the relative importance of influx from rivers and from the
global ocean. Thus the reported subsidence deficit at rifted margins may
also have a paleo-oceanographic explanation.
Although studies of active rifts could in principle help address these
problems (substituting heat flow measurement for postrift subsidence
as a proxy for the thinning of the lithosphere), there are few such
margins, and even fewer appropriate studies. Buck et al. (1988) did point
out that the structure and heat distribution of the northern Red Sea, once
invoked as an example of simple shear lithospheric extension
(Wernicke, 1985), could be best explained by a symmetric model in
which extension and thinning focused beneath the centre of the rift.
Thus for many ancient rifted margins, the measurable subsidence provides constraints on the pattern of whole crustal but not
lithospheric thinning. However, as discussed below, identical patterns
of whole crustal thinning can be produced by very different patterns
of lithospheric extension. The implication is that total subsidence (and
thus bathymetry), although invoked by some authors (e.g., Le Pichon
and Sibuet, 1981; Kusznir and Karner, 2007) as support for specific
models of lithospheric thinning, really only confirm the pattern of
whole crustal thinning (which also controls gravity) and not the
extension mechanism.
13
(Nova Scotia to 0.75 km/km (S Newfoundland Basin). The gradual
increase in the thinning gradient as the crust thins from ∼30 km
towards ∼ 25 km is reminiscent of a necking profile, except that the
thinning appears asymmetric relative to the edge of the continental
crust (Fig. 4). Using the terminology of Davison (1997), some margins
appear “hard” (abrupt transition from relatively normal thickness
crust) whereas at other (“soft”) margins, the change in crustal
thickness takes place far more gradually. Some conjugate pairs show
similar thinning profiles for both margins, for instance Morocco–Nova
Scotia (soft–soft), Porcupine (hard–hard), Rockall (hard–hard, with
the constraint that it is not known exactly where the COT occurs if at
all), whereas others appear markedly asymmetric (Labrador–W
Greenland; Flemish Cap–Galicia; south IAP–Newfoundland Basin).
The sort of asymmetry the last group displays might be interpreted as
implying extension by an asymmetric crustal-scale process, such as a
crustal-scale detachment, but asymmetry of the final margins does
not necessarily imply an asymmetric rifting process (Reston and
Pérez-Gussinyé, 2007). Rifted margins form through a long period of
extension and even if each rift phase is fundamentally symmetric, the
locus of rifting may change with time, resulting in asymmetric basins
and margins. The Galicia–Flemish Cap margin is a case in point: the
asymmetry here need not imply the presence of a lithospheric-scale
shear zone as has been suggested (Lau and Louden, 2007) but rather
may reflect the changing focus of rifting during prolonged extension
(Tucholke et al., 2007), resulting in an older abandoned rift basin
on one side (the Galicia Interior Basin –rifting) and a more abrupt
continental shelf on the other. This can be illustrated by modelling the
stretching across the Galicia margin as the superposition of two
symmetric rift phases (Fig. 6), each approximated by Gaussian necking
profiles: an older (Tithonian–Hauterivian) one centred on the Galicia
Interior Basin (Pérez-Gussinyé et al., 2003; Reston, 2005) but extending as far as the deep Galicia margin, and a later one (Hauterivian
or Barremian to Aptian) modelled on the shape of thinning across the
SE Flemish Cap margin, but applied to the Galicia margin to make this
second rift phase symmetric between the two margins. In applying
the curve from the Flemish Cap margin to the Galicia margin, the two
profiles are matched at the point of complete crustal embrittlement,
as discussed below. Summing the log of the stretching factors of these
two rift phases reproduces the log of the observed stretching factor
across the entire Galicia margin remarkably well, suggesting that the
main asymmetry between the Flemish Cap and Galicia margins may
4. Crustal thinning: distribution of strain
Subsidence and melt production are both related to the way the
lithosphere thins and the thermal structure of the lithosphere and the
underlying asthenosphere, with total subsidence being controlled in
particular by the way the crust thins. I now turn my attention to this
topic, and consider the distribution of crustal thinning across margins
and their conjugates, focussing in particular on symmetry.
Symmetry comes in all shapes and sizes. For instance, rifted margins
and rifts commonly are divided into domains or segments dominated by
faults dipping in one direction, but this is completely compatible with an
overall symmetric pattern of extension and thinning (e.g. Reston and
Pérez-Gussinyé, 2007). Here we are concerned instead with large-scale
asymmetry in margin structure and in the processes of lithospheric
extension. Most previous discussions (e.g. Davison, 1997) have been
rather qualitative, but the thinning profiles presented here provide a
more quantitative assessment, bearing in mind the limitations of the
seismic velocities for mapping out different crustal levels.
It is clear from both the crustal sections (Figs. 2 and 3) and thinning
factor profiles (Fig. 4) that the crust undergoes a relatively sharp
thinning from ∼ 25 to ∼10 km at a gradient of between 0.2 km/km
Fig. 6. Plot of the logarithm (base 10) of the stretching factors deduced from crustal
thinning across the Galicia and Flemish Cap margins against distance from the point
where the entire crust became brittle. In gray dash are simplified curves for the
stretching in the Galicia Interior Basin (GIB) and the Flemish Cap margin. The sum of the
two rift phases (solid grey line) fits the observed distribution of stretching across the
entire Galicia margin (solid black line), illustrating how an asymmetric margin pair can
be produced by two non-superimposed symmetric phases of rifting.
14
T.J. Reston / Tectonophysics 468 (2009) 6–27
be due to the presence of a major, older rift phase on the Galicia side.
Of course, these thinning profiles are only a crude representation of
the true distribution of thinning both in time and space and should
not be taken as a true representation of the rift history between Spain
and Flemish Cap (especially give the complexities associated with the
rotation of the Cap — Sibuet et al., 2007), but do demonstrate that a
strongly asymmetric rift can develop through the occurrence of two or
more symmetric phases of rifting at different locations. Such a pattern
might arise if the locus of rifting migrates during slow rifting due
to strain hardening beneath the rift axis as the stretched mantle
lithosphere and upwelling asthenosphere both cool and increase in
strength (Kusznir and Park, 1987; Bassi, 1995).
4.1. Determining thinning at different crustal levels
It is tempting to relate an asymmetric crustal thickness distribution to an asymmetric rift process, such as the simple shear model
(Wernicke, 1985). However, in its simplest form, this model predicts
symmetric thinning of the whole crust — and symmetric total
subsidence — as long as the master shear zone passes through the
crust at a constant dip (Voorhoeve and Houseman, 1988; Fig. 7).
Observations from mature rifted margins (e.g. Le Pichon and Sibuet,
1981) that crustal thinning and total subsidence give the same
stretching factor does not mean that extension occurred by symmetric, uniform pure shear (McKenzie, 1978) as the simple shear
model predicts exactly the same relationship. Although some
conceptual models for margin development with more complex
detachment geometries (e.g. “delamination” — Lister et al., 1991) may
produce “upper plate” and “lower plate” asymmetries, these are not a
fundamental prediction of lithospheric extension along a throughgoing master shear zone. As a result, the absence of the detailed
“predictions” of these models does not exclude a simple shear origin,
even if as discussed below, such a model is rheologically unlikely.
However, asymmetric rifting processes such as the simple shear
model and depth-dependent stretching models do predict an asymmetry in the distribution of upper and lower crust (Fig. 7). Consequently, it is useful to consider the distribution of upper crustal and
lower crustal as well as whole crustal thinning. Systematic differences
between the patterns of upper crustal and lower crustal thinning may
indicate an asymmetric rifting process; the upper plate margin and
lower plate margin for a simple shear model show clear and distinguishable divergence in upper crustal and lower crustal thinning.
Thus if the thickness of different crustal layers can be mapped, they
can be used to constrain the distribution of strain during rifting. High
quality, modern wide-angle data would seem to provide a means to do
just this, but raise new problems, discussed below.
Reston (2007a) discusses the main problems with using seismic
velocity structure to investigate changes in the thickness of various
crustal layers. These are: 1) the velocity of rocks need not stay constant
during rifting, 2) the probability that crustal layers are not of uniform
thickness to start with, and 3) the accuracy of the velocity models.
Overall crustal thickness is relatively well determined, especially
when PmP reflections are recorded, but the location of intracrustal
velocity boundaries is less well resolved. Typically, high quality
modern wide-angle data can determine the large-scale velocity
structure to about 0.1 km/s (e.g. Dean et al., 2001), but are less
capable of resolving small structures (Van Avendonk et al., 2006) and
have a vertical resolution of ∼ 1 km. The implication is that trends of
the thickness of velocity layers are likely to be correct (Reston, 2007a),
although the precise degree of such changes less so.
The effect of uncertainities in both seismic velocity and in the location
of intracrustal boundaries can be estimated, as shown here for two
profiles across the Newfoundland Basin (Fig. 8). Although the uncertainty
in the stretching factor β of the upper crust and of the lower crust
become very large towards the continent–ocean transition, that in the
reciprocal (1 /β) or the thinning factor (1 − 1 /β) remain more constant.
Similar results were found for the Iberia Abyssal Plain margin (Reston,
2007a). However the overlap of the thinning factor ranges towards the
COT imply that any depth-dependent variation in thinning (depthdependent stretching) is too small to be resolved using seismic velocities.
Less quantifiable are the changes in seismic velocity of any given
rock that are expected during rifting as a result of reductions in
temperature and pressure, fracturing and alteration. The effects of
temperature and pressure tend to cancel out, but both fracturing and
alteration can lead to a net reduction in velocity. As a result, lower
crustal rocks (initially with relatively high lower crustal velocity) can
appear after rifting to be upper crust on the basis of their reduced
seismic velocity. As a result, excess thinning of the lower crust should
be viewed with suspicion. However when comparing the symmetry of
a pair of conjugate margins, similar effects might be expected on both
Fig. 7. Crustal structure predicted by extension of the lithosphere along a single throughgoing, planar shear zone, modified after Voorhoeve and Houseman (1988). A: starting model.
Lithospheric shear zone dips at an angle ϕ, cutting through crust (thickness Zc, equally divided into upper crust and lower crust (shaded)) and whole lithosphere (thickness ZL).
B: Extension sufficient to exhume mantle rocks by movement along the shear zone. C: as B, but after isostatic compensation by vertical shear. Note simple symmetric thinning of the
crust, producing a symmetric basin. D: detail of crustal structure (V.E. of 2) resulting from complete crustal separation along a single shear zone. Note overall symmetry of crustal
thinning and thus of subsidence, but asymmetry in the distribution of upper crustal and lower crustal thinning. E: plot of 1 − 1 / β for the upper, lower and whole crust as a function of
distance from the edge of the continental crust. Asymmetric rifting on a crustal scale should leave a similar pattern at rifted margins.
T.J. Reston / Tectonophysics 468 (2009) 6–27
15
Fig. 8. A, B, C: Plot of whole crustal, upper crustal and lower crustal parameters across the north Newfoundland Basin margin. A) apparent stretching factor β for the whole crust (βc),
the upper crust (βuc), and the lower crust (βlc) showing error bars for a ∼0.1 km/s uncertainty in seismic velocity and a 1 km uncertainty on the depth to the boundary between the
mid and lower crust. B) the corresponding plot of 1 / β showing uncertainties; C) the corresponding plot of thinning factors (1 − 1 / β) for the whole crust, upper crust (with
uncertainity) and lower crust (with uncertainty). Note that uncertainty is generally less than 0.1 and that at thinning factors greater than ∼0.7, the thinning factors are the same
within error for upper crust, lower crust and whole crust. D) Plot of whole crustal thinning vs upper crustal thinning (including uncertainty) at the north Newfoundland Basin margin,
showing no discernible depth-dependent stretching. E) Plot of thinning factors (1 − 1 / β) across the south Newfoundland Basin margin. Note uncertainty is generally less than 0.1 and
how at high thinning factors upper crustal, whole crustal and lower crustal thinning factors are all approximately equal within error. F) upper crustal thinning vs whole crustal
thinning across the South Newfoundland Basin, showing no discernible crustal depth-dependent stretching.
16
T.J. Reston / Tectonophysics 468 (2009) 6–27
sides, so the relative (between the conjugate margins) thinning of the
different layers should still be relevant.
Finally, it is unlikely that the crust consists of horizontal layers of
uniform thickness at the start of rifting. Furthermore, rifts are likely to
initiate or focus on inhomogeneities, e.g. terrane boundaries and
sutures across which the velocity structures are likely to change. An
example is the Galicia Interior Basin which follows at least in part a
Variscan terrane boundary (Pérez-Gussinyé et al., 2003) and as a result
has a very different velocity structure on either side. As focussing of
rifting is most likely to occur at such boundaries, it is not surprising
that the velocity structures of conjugate margins do not match. As a
result, I do not directly compare these, but rather the way the different
crustal levels appear to thin on each margin. Furthermore, I only
consider those margins where velocity layering is constrained by good
quality wide-angle data from the COT until the crust is at least 20 km
thick (i.e. at or close to the continental shelf) and concentrate on
conjugate pairs of margins where any major asymmetries could be
readily identified. Of the 22 margins shown in Figs. 2 and 3, only the 14
in Fig. 2 meet these criteria and are used for analysis of the distribution
of crustal thinning across conjugate margins.
Bearing in mind these concerns, it is nevertheless instructive to
discuss the apparent thinning of the different crustal levels toward
the margin. To remove the effect of initial layer and crustal thickness it is
necessary to normalise the data by plotting the thinning factor (1
− 1 / β) rather than thickness, an approach also adopted by Lau and
Louden (2007). The uncertainty in the determination of seismic
velocity of ∼ 0.1 km/s in the crust where velocity typically varies by
∼ 1 km/s represents a 10% error in the location of any given velocity.
The vertical resolution of 1 km provides an additional 1 km uncertainty. Using these bounds, maximum and minimum estimates of
the location of the boundary between the upper and lower crust
have been computed for the SCREECH2 and SCREECH 3 examples
and used to estimate errors on the thinning factors shown (Fig. 8).
Fig. 9 shows the variation in (1− 1 /β) for the whole crust (βc — solid),
the upper crust (βuc — long dash) and the lower crust (βlc — short dash)
for all margins in Fig. 2 where the structure from the COT toward the
shelf is constrained by modern wide-angle data. In most cases a midcrustal boundary or velocity contour (marked in bold in Fig. 2) has been
used to split the crust into upper and lower halves. Margins to the west
are shown black; those to the east in gray, allowing conjugate or nearly
conjugate margins to be displayed on the same plots The likely errors for
each margin are shown once even if the margin is plotted twice to avoid
unnecessary overlap with their conjugates: error bars for Porcupine
Basin and Rockall Trough are similar to the margins shown but are not
shown as thinning appears so close to symmetric that the curves would
overlap and hence obscure each other.
4.2. Evidence for asymmetry of rifting processes
As described above, the variation in the relative thickness of
different crustal layers, particularly when the converse can be
observed on the conjugate margin, may constrain the asymmetry of
the rifting process, which I define as an asymmetric distribution of
strain associated with one episode of rifting, such as that caused by
movement along asymmetric structures such as lithospheric-scale
shear zones. I begin by plotting thinning factors against distance from
the feather edge of the continental crust, chosen as an easier point to
define in most cases than the edge of the unthinned crust landward.
Apparent divergence between upper and lower crustal thinning
can be observed for IAM9, for the Galicia Interior Basin (GIB) part of the
Galicia profile, for the landward portion of the Nova Scotia margin
(SMART2) and for the Moroccan margin. In each of these cases,
the pattern might indicate a degree of asymmetry during a phase of
rifting.
The structure of the Galicia Interior Basin (on Galicia profile at
∼ 200 km from the COT) is dominated by east-dipping faults (Pérez-
Gussinyé et al., 2003; Reston, 2005, 2007b), which although not
lithospheric or even crustal shear zones do result in a local offset of
upper crustal and lower crustal extension. Overall though the basin is
close to symmetric. The conjugate SE Flemish Cap margin shows no
significant divergence between lower and whole crust, and fluctuations in the thinning of the upper crust appear to be related to
individual fault blocks but do not form a consistent asymmetric
pattern.
Probably the most pronounced example of heterogeneous crustal
extension comes from the Moroccan margin. Here, the lowermost
crust pinches out well over 120 km from the COT. The distribution of
thinning of the crust of the Moroccan margin has been interpreted by
Maillard et al. (2006) as evidence for a lithospheric-scale shear zone
and the pattern does resemble the “upper plate” margin expected
above such a shear zone, although the shear zone would have to be
very low-angle (b15°) to explain the distance between the last lowermost crust and the COT. However the conjugate Nova Scotian margin
(SMART1 — Funck et al., 2004) switches from possible “lower plate”
characteristics (the upper crust thinned more than the lower crust
between 100 and 200 km from the COT) to more upper plate characteristics (lower crust apparently absent, but a thin upper crustal
remnant) over the 50 km leading to the edge of continental crust.
Furthermore, along strike the asymmetry on the Nova Scotian margin
(SMART2 — Wu et al., 2006) also resembles an upper plate margin.
Thus the structure of these margins does not provide compelling
evidence for a simple asymmetry in the rifting process, but rather for a
more complicated pattern of non-uniform extension.
IAM9 (Fig. 2) most closely resembles the predictions of the simple
shear model, in this case for a “lower plate” margin (Fig. 9) below a
shear zone dipping at ∼ 15° to the west. However, the conjugate
Newfoundland Basin margin (SCREECH 2) shows only limited evidence for asymmetry. The fluctuations between upper plate and lower
plate characteristics 100–200 km from the COT represent fluctuations
in the location of the velocity contours: fine, largely unresolved
variations in velocity produce large apparent variations in the
thickness of the upper and lower crust. The last 100 km of continental
crust (0–100 km from the edge of continental crust) appear to show
more lower crustal than upper crustal thinning, consistent with an
“upper plate” margin, but these variations may be beyond the
resolution of the seismic velocities. A more prosaic explanation for
the structure of the IAM9 margin may be a combination of poorly
resolved velocity of the less extended crust landward (mainly
constrained by gravity data — Dean et al., 2001), leading to incorrect
estimates of the relative degree of thinning of the upper and lower
crusts, and the loss (or reworking as sediment) of some upper crust
through mass-wasting at the steep continental slope. It should also be
noticed that the heterogeneity of the IAM9 crust appears to disappear
as the crust thins to below ∼ 8 km (β N 4; 1 − 1 / β N 0.75).
The eastern margin of Porcupine also could be interpreted as a
typical upperplate margin, with apparently more rapid thinning of the
lower crust. However, as discussed below, the detachment imaged at
this margin cuts to the west, making the eastern margin a lower plate
margin if anything. Furthermore, the western margin displays similar
characteristics, although less pronounced. The thinning across the
Rockall Trough shows almost identical patterns of lower crustal
thinning on both sides, somewhat similar whole crustal thinning, but
rather different upper crustal thinning. This does not resemble a largescale simple shear but does perhaps indicate a degree of heterogeneity
of extension with the locus of upper crust thinning and separation
being shifted somewhat.
In summary, the relative distributions of upper crustal and lower
crustal thinning across these margins may indicate that extension is
heterogeneous and not uniform, symmetric pure shear, but does not
provide any evidence for asymmetry on the scale predicted by the
Wernicke model, particularly when the problems in using seismic
velocity as a proxy for crustal structure are borne in mind. The error
T.J. Reston / Tectonophysics 468 (2009) 6–27
17
Fig. 9. Plots of whole crustal, upper crustal and lower crustal thinning deduced from the velocity structure of 14 margins in Fig. 2. (S)NB — (South)Newfoundland Basin; SIAP — South
Iberia Abyssal Plain. Conjugate and near conjugates are plotted on same graph: black for western side, grey for eastern. Thinning of whole crust is solid line; upper crust long dash line
with vertical stripe error bounds, lower crust (whole crust–upper crust generally) as short dash line with grey error bounds. Lowermost crust for Moroccan margin plotted as short
bars. All are plotted against distance from the oceanward limit of continental crust except for Porcupine Basin where they are plotted from the centre of the basin. Bottom right: at
same scale thinning factors expected for shear zones cutting through 30 km thick crust at angles of 15, 30 and 45° for the upper plate margin and 15° for the lower plate margin. See
text for discussion.
bars for both lower crustal and upper crustal thinning are likely to
have a value of (1 − 1 / β) ∼0.1 (Fig. 8), sufficient to cast doubts on any
perceived asymmetry in the distribution of upper crustal and lower
crustal thinning. Interestingly, the margins tend to resemble upper
plate margins more than lower plate margins. This might represent a
systematic depth-dependent stretching (Driscoll and Karner, 1998),
but may also be influenced by the decrease in the seismic velocity as a
result of continuing extension (fracturing and alteration as described
above and by Reston, 2007a). The upper crustal velocities noted for the
Porcupine, Newfoundland Basin, Labrador, southern Iberia Abyssal
Plain (drilling transect) and northern Nova Scotia margin are
completely compatible with the presence of highly fractured lower
crustal rocks: basement has only been sampled on one of these
margins (the southern IAP), recovering dominantly lower crustal rocks
formed during the Variscan orogeny and uplifted through the Ar–Ar
closing temperature of feldspar of 200 °C at 137 Ma (Whitmarsh and
18
T.J. Reston / Tectonophysics 468 (2009) 6–27
Fig. 10. Illustration of the rheological evolution of the upper lithosphere during pure shear stretching. Initial weak zones within the mid-crust and lowermost crust gradually
disappear as the crust thins and creeping rocks become brittle (Pérez-Gussinyé and Reston, 2001). Once the entire crust becomes brittle, water can reach the mantle, causing
serpentinization. UC upper crust; LC lower crust; CMB crust–mantle boundary; WQ wet quartz; AN anorthite; OL olivine.
Wallace, 2001), that is during rifting, but with a seismic velocity of
∼ 5.5 km/s (Krawczyk et al., 1996).
The sections in Fig. 9 do however appear to show a large-scale
asymmetry between conjugate margin pairs of Labrador–West Greenland, Galicia–SE Flemish Cap, Newfoundland Basin–Iberia Abyssal
Plain, and a weaker asymmetry between Nova Scotian and Moroccan
margins, with in each case the latter margin being thinned more
abruptly that the former. Such an asymmetry may indicate the way
rifting focussed prior to breakup. However it is in part a function of
plotting the thinning curves vs distance from the edge of the continental crust. One possible influence on the development of
asymmetry may be the changing rheology of the lithosphere, and in
particular when the crust becomes brittle (Bassi, 1995).
5. Rheological evolution: embrittlement, serpentinization and
coupling
Apart from the thinning of the crust to zero and the lack of
magmatism at magma-poor margins, one of their other key characteristics (Figs. 1–3) is the presence of serpentinized mantle, both
within the continent–ocean transition (COT) oceanward of the last
continental crust, and beneath the feather edge of the crust. The
formation of these serpentinites results from the reaction of mantle
peridotites with water and thus requires that sufficient volumes of
water can be brought into the mantle. As Pérez-Gussinyé and Reston
(2001) pointed out, this can only be along faults linking the peridotites
with the surface, in turn requiring the whole crust to be brittle. As an
entirely brittle crust does not have decoupling ductile zones, whereas
serpentinites are relatively weak, serpentinization may be also be
linked to changes in tectonic style.
The development of an entirely brittle crust comes about through
changes in rheology accompanying lithospheric extension, changes
which have implications for the coupling between crust and mantle as
well as for crustal hydrogeology. It is generally well understood that the
interplay between downward increasing pressure and temperature and
changes in the dominant mineralogy of the lithosphere from granitic
(quartz-dominated), through mafic-gabbroic (anorthosite dominated)
to ultramafic (olivine dominated) results in a rheological stratification
and the development of weak zones where creep probably dominates in
the mid-crust and at the base of the crust, and intervening strong zones
including the seismogenic brittle upper crust which is 8–10 km thick
(Kusznir and Park, 1987; Jackson and White, 1989). As a result, initial
extension in the uppermost crust is likely to be by brittle faulting, in the
weak middle and lowermost crust by ductile creep, with perhaps
boudinage of any intervening strong layer. However these zones do not
remain fixed: as the crust thins and cools, the reduction in overburden
pressure and temperature means that rocks which originally deformed
by plastic creep gradually become brittle and fracture. The result is that
the initial weak zones in the mid-crust and deep crust disappear and the
entire crust becomes brittle (Fig. 10).
Fig. 11. Summary of numerical modelling results (Huismans and Beaumont, 2002), showing how deformation of a decoupled lithosphere (weak lower crust) results in overall
symmetric extension. However for a coupled lithosphere, extension is asymmetric as throughgoing structures develop.
T.J. Reston / Tectonophysics 468 (2009) 6–27
The embrittlement of the crust is perhaps counter-intuitive as
there is a common misconception that lithospheric extension brings
hot rocks up towards the surface and as a result, the lithosphere is
heated and rocks move from being brittle to ductile. However, even if
the vertical distance between the surface and the base of the
lithosphere decreases so that the geothermal gradient increases, the
temperatures of individual pieces of the lithosphere (rocks) do not
(see Reston, 2007b for a complete discussion).
The rheological evolution culminating in crustal embrittlement
has two important consequences. First, complete crustal embrittlement implies that brittle faults can for the first time cut down from the
surface into the mantle, bringing aqueous fluids that can serpentinize
substantial volumes of the mantle (Pérez-Gussinyé and Reston, 2001).
Second, Huismans and Beaumont (2002) show that the large-scale
symmetry of the rifting process is largely controlled by whether the
different lithospheric levels are coupled or decoupled. In the simplest
coupled models, through-going shear zones may develop, but in more
realistic models for lithospheric strength, the presence of weak
decoupling layers in the mid-crust and just above the Moho make
single through-going shear zones unlikely. As a result, deformation is
initially distributed among several structures and overall is approximately symmetric (Fig. 11). However, as extension proceeds and the
weak zones disappear, the strong layers become ever more tightly
coupled to another, and in particular the crust becomes coupled to the
upper mantle, allowing deformation at a late stage to become strongly
asymmetric (Reston and Pérez-Gussinyé, 2007).
The stretching factor at which the entire crust should become
brittle (βb) is to some extent controlled by rift duration and the
original lithospheric configuration (Pérez-Gussinyé et al., 2001 —
Fig. 12); changing the mineralogy and rheology of the lower crust in
particular may shift the curve up and down respectively within the
shaded area. The end-member lower crustal rheologies (dry quartz
and anorthite) probably underestimate and overestimate lower
crustal strength respectively, so also shown is an intermediate rheology calculated using polymineralic flow laws (Tullis et al., 1991) for a
50:50 mix of the two.
19
Although young hot orogens (e.g. Woodlark Basin) behave quite
differently (Pérez-Gussinyé et al., 2001), the results for both cooled
post-orogenic and cratonic lithosphere are quite similar (Fig. 12).
Increasing the rift duration leads to crustal embrittlement at lower
stretching factors as the lower crust has time to cool; moving from a
strong (anorthitic) to weak (dry quartz) lower crust increases the
stretching factor at which embrittlement occurs. It can be seen that
the stretching factor at which the entire crust becomes brittle is
dependent on both rift duration and the rheology of the lower crust,
itself dependent on lower crustal mineralogy and hence composition.
In general, the modelling predicts that the entire crust should
become brittle at stretching factors between about 3 and 5, depending
on rift duration. As noted by Pérez-Gussinyé and Reston (2001),
this corresponds reasonably well with the crustal thickness at the
landward limit of serpentinized mantle for several margins. The range,
based on uncertainties in rift duration and lower crustal composition, at
which the crust should have become brittle during rifting (Figs. 2 and 3)
tends to occur just landward of the landward limit of serpentinized
mantle, emphasising the correlation between crustal embrittlement and
serpentinization. At some margins however serpentinization only starts
a long way further oceanward, if at all. This probably reflects other
influences in the hydrogeology and permeability of the crust, including
the spacing of faults and the presence of thick sequences of sediments
especially evaporites (e.g. on the Moroccan margin).
5.1. Change from symmetry to asymmetry — the effect of embrittlement
It has been shown that conjugate margins appeared quite asymmetric when thinning profiles were plotted as a function of distance
from the edge of the continental crust (Figs. 4 and 9). Bearing in mind
the possible importance of coupling and decoupling in controlling rift
(a)symmetry, a more useful way of aligning the curves may be at the
point at which the crust becomes brittle as this might be expected to
mark a transition from symmetric (decoupled lithosphere) to a later
asymmetric extension (Whitmarsh et al., 2000; Reston and PérezGussinyé, 2007) when the development of large throughgoing shear
Fig. 12. Plot showing the stretching factor at which the entire crust should become brittle (βb) as a function of rift duration for two different lithospheric models (Pérez-Gussinyé et al.,
2001). Broad grey band — the width represents variations depending on the relative importance of feldspar (lower limit) and dry quartz (upper limit) in the lower crust. Grey dashed
line is for a 50:50 dry quartz/feldspar aggregate. Estimates of rift duration for margins shown (Nova Scotia — Funck et al. (2004), Wu et al. (2006); Morocco — Maillard et al. (2006);
Porcupine Basin — O'Reilly et al. (2006); Goban Spur — Bullock and Minshull (2005); Armorican — Thinon et al. (2003); South Newfoundland Basin — Reid (1994),Lau et al. (2006);
Galicia Bank — Mauffret and Montardet (1988), Reston (2007a,b); SE Flemish Cap — Tucholke et al. (2007); Iberia Abyssal Plain — Wilson et al. (2003); Minshull et al. (2001);Tucholke
et al. (2007); Rockall Trough — Pérez-Gussinyé et al. (2001); West Greenland, Labrador — Chalmers and Pulvertaft (2001); Reston and Pérez-Gussinyé (2007). As βb varies more with
lower crustal rheology than with uncertainty in rift duration at any given margin, the former is used to give first estimate of βb uncertainty (Galicia is illustrated) used to plot
predicted embrittlement in Figs. 2 and 3.
20
T.J. Reston / Tectonophysics 468 (2009) 6–27
zones is possible (coupled lithosphere — Huismans and Beaumont,
2002).
To separate out such late-stage asymmetry from any early
asymmetry, Fig. 13 shows the thinning curves of individual margins
shifted laterally (arrows) so that they do not coincide at a thinning
factor of 1 (the edge of the continental crust), but rather at a thinning
factor of between 0.67 and 0.8, corresponding to the stretching factor
βb of 3–5 at which the margin in question should become brittle
(Fig. 12 — Pérez-Gussinyé and Reston, 2001). Much of the apparent
asymmetry between margin pairs (hard vs soft — Fig. 4) is removed:
the margins are symmetric (curves overlap) from thinning factors
of zero (little or no extension) until βb is reached (thinning factor of
1 − 1 / βb). Only at thinning factors above this do the conjugate margin
pairs become asymmetric (Fig. 13). This implies that the asymmetry
for example between Labrador–W Greenland, between Newfoundland Basin (north and south)–Iberia Abyssal Plain, and between Nova
Scotian and Moroccan margins only developed once the entire crust
had become brittle and coupled, exactly as predicted by Huismans
and Beaumont (2002). The asymmetry between SE Flemish Cap and
Galicia Bank is reduced as the thin crust above the S reflector is
ignored in aligning the thinning curves, leaving an asymmetry that
can easily be explained by the overlapping of two rift phases (Fig. 6).
In each case where an asymmetry at a stretching factor greater
than βb can be detected, the margin not shifted to the right (i.e. the
margin with a broader expanse of very thin crust) appears to have the
seismic structure of an upper plate margin. However, as discussed
above, these characteristics are equally compatible with the reduction
of velocity (and hence the transfer of seismically defined lower crust
into upper crust) accompanying faulting, fracturing and pressure
reduction. Thus the velocities may simply be confirming that highly
thinned crust is also highly fractured and faulted.
5.2. Development of serpentine detachments
As discussed above, crustal embrittlement leads to both increased
coupling between the crust and the mantle and the serpentinization
of the mantle beneath the thinned and fractured crust. As serpentinites have a friction coefficient considerably below that of the most
rocks, they might be expected to decouple deformation between the
crust and the mantle. However, the serpentinites develop where fluids
can penetrate the mantle, that is along faults and fractures and so
instead of forming a general decoupling zone may tend to form
local detachments where major faults cut across the crust–mantle
boundary. As detachment faults are fundamentally asymmetric
Fig. 13. Thinning factor curves plotted against distance from point where entire crust should have become brittle (βb, filled black circle — based on rift duration and starting model).
Undercrusting serpentinites occur in the region marked by bold lines. Most of the asymmetries in Fig. 9 are removed, leaving variable widths of highly extended, completely brittle crust to
the left of the embrittlement point. This implies that most of the apparent asymmetry in crustal structure across conjugate margin pairs may develop after crustal embrittlement and
increased coupling between crust and mantle. The implication is that most of the extension is approximately symmetric, but that breakup and crustal separation are asymmetric.
T.J. Reston / Tectonophysics 468 (2009) 6–27
21
Fig. 14. Asymmetric P detachment beneath the Porcupine Basin (SW Ireland — A) and the S detachment at the west Galicia margin (depth image — B). Note that both detachments
developed when the crust was already thin, cut across the CMB and follow the boundary between crust and serpentinized mantle. Postrift sediments are partially blanked out and
crustal basement lightened to help emphasise fault block structure. C–F: illustration of the possible sequential evolution of the S detachment and overlying fault blocks during
polyphase faulting at the west Galicia margin. S may have started as a steep fault only becoming a low-angle detachment once the entire crust was brittle and mantle serpentinites
had formed. Note that S forms between blocks III and II and that block III is dismembered by 2 generations of later faults detaching onto S (D, E and F). Block III is the footwall to TBF
(“top basement fault”) which is cut by later faults and forms top basement to several small fault blocks above S (D–F).
structures, their presence (Fig. 14) is direct evidence for the
development of late-stage asymmetry once the crust had been prethinned to less than 8 ± 2 km (stretching factors between 3 and 5 —
Fig. 12) and thus become entirely brittle. The best known example is
the S reflector observed on the Galicia Bank margin; this structure cuts
to structurally deeper levels to the west. Its continuation may be
interpreted on the Flemish Cap margin (SCREECH 1 — Hopper et al.,
2004) where it follows the CMB, cutting to depth to the west (Reston
et al., 2007), but no east-cutting mirror-image is seen. A similar
detachment (P) imaged beneath the Porcupine Basin (Reston et al.,
2001) is also an asymmetric structure, cutting from the east flank of
the basin across the crust before following the CMB in the centre of the
basin and beneath its western flank (Fig. 14).
Reconstructions (Reston, 2005; Reston et al., 2007) suggest that
S developed as a high angle fault and only became a low-angle
detachment once the crust had thinned to about 7 km (Fig. 14),
compatible with the prior onset of serpentinization. Relationships
between wedges of synrift sediment and the block-bounding faults that
detach onto S indicate that S was active at angles down to or below 15°
(Reston et al., 2007). This can be explained by the extremely low friction
coefficients associated with serpentine and related minerals, or by the
transient development of high fluid pressures (Reston et al., 2007). The
development of detachments thus follows the onset of serpentinization,
itself a result of crustal embrittlement and the hydrogeology of the
fractured crust and is one stage in a long and complex polyphase faulting
history. This appears to have resulted in a focussing of extension towards
22
T.J. Reston / Tectonophysics 468 (2009) 6–27
the site of eventual crustal separation: the faulting above the S reflector
occurred later than the development of the larger fault blocks further to
the east (Reston, 2005) and the fault blocks overlying S developed
through the dismemberment of such larger, earlier fault blocks (Reston
et al., 2007). Similar polyphase faulting and focussing of extension
towards the COT is likely to have occurred at other rifted margins and
represents a simple focussing of strain at the weakest point of the rift,
that is the rift axis.
6. The extension discrepancy and depth-dependent stretching
Rifted margins form by the extension and eventual separation of
the crust and lithosphere. Extension results in thinning, apparent on
Figs. 2 and 3, and is accommodated by a combination of brittle faulting
on scales ranging from detachment faults through to sub-seismic
fractures, and ductile shear, flow and creep. As all lithospheric levels
are pulled apart, albeit possibly not in a uniform manner, the total
amount of extension going from the continental interior to the onset
of seafloor spreading should be equal at all levels. However, it is
commonly reported that the extension of the upper crust, as determined from fault geometries (βf), in much less than that required to
explain the crustal thinning (βc) either observed from wide-angle
velocity structure (Figs. 2 and 3) or deduced from subsidence (Fig. 15 —
Sibuet, 1992; Ziegler, 1983; Davis and Kusznir, 2004). This is the
extension discrepancy, expressed by rifted margins plotting in the
lower right quadrant of a (1 − 1 / βc) vs (1 − 1 / βf) plot, whereas
less extended rift basins plot along the diagonal (no extension
discrepancy — Kusznir and Karner, 2007).
There are two end-member explanations for the extension
discrepancy. First that crustal, as well as lithospheric, extension is
strongly depth-dependent and that at rifted margins the extension of
the upper crust is far less than the extension and thinning of the whole
crust. However, as noted above, the total amount of extension across
the whole margin must be equal at all lithospheric and crustal levels,
so that somewhere extension and thinning of the upper crust must
exceed that of the lower crust, although this has not yet been
observed. Alternatively, the other end-member explanation is that not
all the brittle extension has been measured, discussed later.
One possibility is that the margin exhibiting the extension discrepancy forms above a major shear zone, which transfers the deeper
crust out from beneath a virtually unfaulted upper plate margin.
However, as extension discrepancies are observed at virtually all
margins, this would imply that all are upper plate margins, the socalled upper plate paradox (Driscoll and Karner, 1998). Instead, more
symmetric patterns of depth-dependent stretching have been proposed, in which the lower crust on both sides of the incipient ocean
has somehow been displaced, leaving behind thin but apparently little
extended upper crustal section.
The sections in Fig. 2 and the thinning profiles in Fig. 9 can be used
to assess whether crustal depth-dependent stretching (as opposed to
lithospheric-scale depth-dependent stretching described above) is a
major process. The velocity structure and resulting thinning profiles in
Figs. 9 and 13 do show that many margins appear to exhibit upper
plate characteristics, although not consistent with movement along a
major crustal-scale detachment. However, such upper plate characteristics may simply represent the tendency for rifting to reduce seismic velocity and so give lower crustal rocks upper crustal velocities, as
implied by the presence of prerift lower crustal rocks with upper
crustal velocities on the southern IAP margin.
A degree of crustal depth-dependent stretching is to be expected as
it unlikely that the crust would extend perfectly uniformly, but the
problem with invoking DDS to explain the extension discrepancy is
the extreme amount of DDS required. By assuming that βf = βuc, and
that the extension discrepancy is caused by DDS, it is possible to invert
the βf (= βuc) and βc values and reconstruct a section, constrained
further by the condition that the upper crust can nowhere be thicker
than the whole crust (Fig. 15). The result is a section where the entire
crust thins from ~ 25 km (right) to at most 5 km (left), while the upper
crust maintains a maximum 5 km thickness. Given that the thickness
of the seismogenic brittle upper layer is 8–10 km (Jackson and White,
1989), this seems problematic.
The variation in the apparent thinning of the different crustal
layers shown in Figs. 9 and 14 does suggest a degree of heterogeneous
crustal extension (DDS), but does not show that this is sufficient to
explain the extension discrepancy. On a crossplot of upper crustal vs
whole crustal thinning (Reston, 2007a), significant crustal DDS should
Fig. 15. The extension discrepancy. A: The stretching factor deduced from fault geometries (βf) is considerably less than that deduced from crustal thinning (βc) across the Goban Spur
(A) and West Galicia Bank (B) margins. This means that these and other rifted margins plot below and to the right of the diagonal on a plot of thinning factors deduced from crustal
thickness and fault geometry data (C). This is the extension discrepancy. If fault geometries genuinely represent the amount of extension of the brittle upper crust, the data can be
inverted and thickness of the brittle upper crust and of the entire crust reconstructed and compared (D). As the upper crust can nowhere be thicker than the whole crust, the two
curves cannot cross, constraining the brittle upper crustal thickness to ∼3 km, which is unreasonably low. Even by maximising the thickness of the brittle upper crust (E), it only
reaches ∼ 5 km, still lower than evidenced by earthquake and rheological data. An alternative is that βf does not record all the extension of the brittle upper crust, meaning there no
need for depth-dependent thinning of the crust and that the reconstructions in D and E are incorrect.
T.J. Reston / Tectonophysics 468 (2009) 6–27
be represented as a systematic deviation away from the diagonal: if
crustal DDS is the cause of the observed extension discrepancy, the
data should plot in the lower right quadrant. However, all 14 margins
that provide sufficient control from the shelf to the COT plot close to
the diagonal (Fig. 16). Although several show some evidence for minor
crustal DDS (e.g., SCREECH 3 or Labrador where upper crust appears to
continue out over serpentinized mantle), the DDS is insufficient to
explain the extension discrepancy. Even where apparently observed
(e.g. SCREECH 3 and Labrador), the upper crust is also very strongly
thinned.
23
The problems in using velocity structure to infer the thinning of
different crustal levels do not help, as the general effect of rifting
would be to reduce velocity through fracturing and alteration, thus
giving lower crust rocks an upper crustal velocity and meaning that in
places the thin upper crustal layer such as observed on SCREECH 3
may actually be lower crustal rocks. This is borne out by the result of
ODP legs 149 and 173 across the west Iberia margin: although the
crustal blocks typically have velocities well below 6.5 km/s, the basement rocks sampled at Site 1067 and 900 were tonalites, amphibolites,
anorthosites and mafic granulites, i.e. mid-crustal to lower crustal
Fig. 16. Plots of upper crustal (1 − 1 / βuc) vs whole crustal thinning factors (1 − 1 / βc) for conjugate margin pairs from the North and North-Central Atlantic magma-poor margins. In
each case, the western half of the transect is shown black and the eastern half bold gray. In each case, the data plot close to the diagonal, indicating no significant depth-dependent
thinning resolvable from the velocity structure. Local deviations correspond to fault block topography, e.g. the lettered points which correspond to the letters on the corresponding
profiles (Fig. 2). A plot of the whole crust vs upper crustal thinning from all the margins (bottom right — including error bounds) emphasises that there DDS does not cause the
extension discrepancy.
24
T.J. Reston / Tectonophysics 468 (2009) 6–27
rocks which formed during the Variscan orogeny and which were
exhumed through 200 °C at 135 Ma, during early rifting (Whitmarsh
et al., 2000). Thus the geophysical “upper plate” characteristics of the
highly thinned margin near the COT are misleading and cast doubts on
the interpretation of other highly thinned crust as being exclusively of
upper crustal origin and evidence for crustal DDS on the scale required
to explain the extension discrepancy.
The plot in Fig. 16 only covers the 14 conjugate margins where the
velocity structure can be traced far enough to estimate the relative
thinning of the crustal layers. However, the sort of DDS required to
explain the extension discrepancy (Fig. 15) is also unlikely to apply at
several other deep margins as these exhibit velocity structures that are
comparable with those were crustal DDS can be ruled out (Fig. 16). For
instance, both the thin Goban Spur section and the deep Armorican
margin (Fig. 3) have upper crustal velocities around or below 6 km/s
underlain by lower crustal velocities close to 7 km/s, very similar to
the other margins where such layers can be traced landwards. At both
margins, the velocity structure suggests the presence of lower crust
right out near the COT. Although several margins (Labrador,
Porcupine, S. Newfoundland Basin, do show a thin layer of upper
crustal velocities extending tens of km further oceanward than the
higher velocity deeper layers, in each case the upper crustal layer is far
thinner than its equivalent further landward, requiring far more
extension than apparently observed: the data plot in the top right not
lower right quadrant of the thinning factor crossplot (compare Figs. 15
and 16). Furthermore, the only margin where this anomalous layer
with only upper crustal velocities has been sampled (the IAP margin
along the drilling transect) has been found to contain both upper
crustal and lower crustal rocks shuffled together by extensional
tectonics; the seismic velocity of the lower crustal lithologies has been
reduced by intense fracturing and alteration.
In summary, although it is quite possible that significant crustal
DDS does occur at other, warmer margins where lower crustal flow is
possible (Hopper and Buck, 1996), there is no evidence that the
extension discrepancy is caused by crustal DDS at any of the magmapoor margins of the North and Central Atlantic. The extension
discrepancy thus needs another explanation. One possibility (Reston,
2005, 2007a) is unrecognised faulting, such as polyphase faulting and
large unrecognised top basement faults. Put another way, the
extension discrepancy can be summarised as βf ≪ βc. In the DDS
explanation, it was assumed that βf ∼ βuc and that as a result βuc ≪ βc.
However, above we have shown that βuc ∼ βc, and thus by implication
that βf ≪ βuc. This means that the seismic images are somehow not
measuring all the extension and thinning of the upper crust.
The amount of brittle extension may be underestimated if either
distributed deformation accompanies the observed faulting, if the
extent of the observed faults has been substantially underestimated, or if
the observed faults are not the only phase of faulting to have occurred.
Davis and Kusznir (2004) show that distributed deformation is
insufficient to explain the extension discrepancy and also argue that
the discrepancy cannot be caused by a phase of faulting that post-dates
the one imaged. However, there is strong evidence at some margins (e.g.
Galicia) for faulting and pre-thinning of the crust prior to the
development of the latest fault blocks (Reston, 2005). There are good
reasons why complex polyphase faulting (Proffett, 1977; Jackson and
White, 1989) may be difficult to image, especially on time sections
(Reston, 2007a), and the same is true for faults that have flexurally
rotated to follow top basement (Reston, 2007a) so that their extent can
be underestimated. In short, there are ample reasons why earlier phases
of faulting are not always recognized at rifted margins (Reston, 2005,
2007a). However, a long and complex rifting and faulting history can be
deduced from field studies (Manatschal, 2004), is predicted by
numerical modelling (Lavier and Manatschal, 2006) and is completely
compatible with the structure and rift history of rifted margins (Reston,
2005, 2007a). A single generation of faults should not accommodate
more than about 100% extension (stretching factor of 2) before rotating
to low-angle and locking up, so that the highly extended deep margins
should have undergone extension along several generations of faults
(Reston, 2005). Indeed, Lavier and Manatschal explicitly recognise three
phases of extension during the formation of a rifted margin, each with
distinctive modes of brittle extension. Fig. 14 shows how simple
polyphase faulting can both closely reproduce the sort of geometries
observed on the seismic images, and be difficult to recognize. It should
not be considered a restoration as there is too little calibration by well
data to confirm the interpretation.
7. Unroofed mantle within the COT
A final feature of many of the margins discussed here is the
presence of unroofed mantle peridotites within the COT (Figs. 2, 3
and 17). These appear to be the geophysical continuation of the
reduced velocity mantle undercrusting the feather edge of the
Fig. 17. A: cartoon after Pickup et al. (1996) showing exhumation of a broad expanse of mantle peridotites within the continent–ocean transition. However this illustration fails to explain
how the peridotites were unroofed, i.e. how the overlying rocks were removed B: detail of IAM 9 (approximate position marked by box in A) showing mounded layer (serpentine breccias —
Reston et al., 2004) above serpentinized basement, interpreted as being topped by series of top basement faults that can be traced to depth as landward-dipping reflections.
T.J. Reston / Tectonophysics 468 (2009) 6–27
continental crust. However geophysics cannot discriminate between
sub-continental lithospheric mantle that has been exhumed within
the COT, and exhumed asthenospheric mantle. These distinctions are
important to choose between the different explanations for magmapoor margins and to constrain the large-scale kinematics of the
stretching process. If the unroofed mantle is new lithosphere resulting
from the cooling of asthenosphere, it must have either been cool or
depleted to explain the lack of melting (Reston and Phipps Morgan,
2004; Pérez-Gussinyé et al., 2006). If it is old sub-crustal lithosphere
pulled out from beneath the continent, then lithospheric-scale
extension must have been strongly heterogeneous, helping to explain
some if perhaps not all, of the melt-deficit at magma-poor margins
(Reston, 2007b). Although the prevailing view is that the exhumed
mantle was sub-continental lithosphere in nature (e.g., Müntener and
Manatschal, 2006), others dispute this (e.g., Abe, 2001). Further
sampling work is needed to distinguish between these models.
In either case, the mantle must have been exhumed through the
removal of overlying portions of the lithosphere. As final exhumation
to the surface must have occurred in the brittle regime, everywhere
this unroofed mantle must be topped by faults which may predate
the structures that offset the top of the mantle basement and which
define the fault block topography within the COT (Fig. 17). The
process of mantle exhumation thus is the result of continued tectonic extension as the crust thins through zero and faults exhume
and expose at the seafloor the underlying mantle (Manatscha et al.,
2001).
Finally, the distinction between unroofing mantle within the
continent–ocean transition and seafloor spreading is partly a matter of
degree. Slow-spread crust is characterised by mantle unroofing,
suppressed magmatism, and detachment faulting (e.g., Tucholke and
Lin, 1994; Reston et al., 2002), features also apparent within the COT.
25
Although magmatism may be more important during the formation of
oceanic crust, probably due to increased focussing of mantle
upwelling and hence melting beneath the spreading centre (Minshull
et al., 2001), perhaps the most meaningful definition of the boundary
between continent and ocean is the boundary between old and new
lithosphere. This has implications for the nature of the COT: the thin
oceanic crust interpreted by Hopper et al. (2004) on profile SCREECH 1
may certainly be thin magmatic crust, but if it is floored by unroofed
sub-continental lithosphere may not be truly oceanic as it would then
not be underlain by oceanic lithosphere.
8. Conclusions
This paper has focussed on the large-scale structure and symmetry
of magma-poor rifted margins whereas the structural evolution has
been discussed more elsewhere (Reston, 2007b). Their overall
structure is the logical result of progressive extension, summarised
in Fig. 18. A number of key elements and stages in the evolution of
such margins can be identified
1. The extreme thinning of the crust (Fig. 18A–D) is the result of
multiple phases of extensional tectonics: that not all of these are
recognised is to be expected, leading to the extension discrepancy.
Substantial depth-dependent stretching of the crust is neither
allowed by the velocity structure of the margins nor expected for a
cool lithosphere with an initially substantial brittle lid.
2. Rifting appears to be a largely symmetric process on a crustal-scale
until the entire crust has become brittle (Figs. 15 and 18C,D),
although the distribution of extension with time may lead to
asymmetric rifts. When this occurs, the complete coupling
between the crustal and mantle allows the development of single
Fig. 18. Schematic model for lithospheric extension resulting in the formation of conjugate magma-poor margins modified from Reston (2007b). A: Initial lithospheric extension by
∼100% (β of 2). Rift focusses above axis of mantle thinning: polyphase faulting in the crust thins this to less than a quarter of its original thickness, until the entire crust is brittle.
(B). Weak serpentinites develop at the base of the thinned crust (B,C), allowing the development of serpentine detachment (C) as observed west of Iberia (S, H) and in the Porcupine
Basin (P). D: Movement along the serpentine detachment and on subsequent detachments results in complete crustal separation and the unroofing of mantle rocks within the COT.
26
T.J. Reston / Tectonophysics 468 (2009) 6–27
large asymmetric faults and shear zones which leave much of the
highly thinned crust on one side of the conjugate margin pair,
producing a late-stage asymmetry. Early asymmetry, such as
between Flemish Cap and Galicia is probably due to the distribution
of different rift phases rather than due to a large-scale asymmetry
in the rifting process.
3. The undercrusting of the crust by serpentinized mantle is the
logical result of the complete embrittlement of the crust once a
stretching factor of about 4 has been reached (Fig. 18B–C).
Serpentinites are only observed where the overlying crust is so
thin that should have become brittle during rifting.
4. The serpentinized mantle beneath the crust is laterally continuous
with serpentinized mantle that forms the basement between the
last continental crust and the first unambiguous oceanic crust.
Mantle unroofing is the logical result of continued extension once
the entire crust has become brittle and detachment faults develop
in the weak mantle serpentinites.
5. The transition to seafloor spreading (Fig. 18D) is marked by a
transition to more magma-dominated and less tectonic divergence,
but may be most rigorously defined as a transition from old to new
lithosphere, marked by the transition from exhumed mantle
peridotites were originally part of the continental lithosphere to
those that represent younger lithosphere formed by the cooling
and depletion (through melt extraction) of asthenosphere.
One important control on the development of such margins
appears to be the rheological evolution of the lithosphere, with crustal
embrittlement both allowing mantle serpentinization and controlling
the transition from a decoupled to a coupled lithosphere with a
related change from symmetric to late-stage asymmetric extension.
Acknowledgements
This paper has benefited from numerous discussions with
colleagues, in particular Marta Pérez-Gussinyé, Jonathan Turner and
Gianreto Manatschal. It is the result of a long programme of support by
funding agencies, in particular the DFG.
References
Abe, N., 2001. Petrochemistry of serpentinized peridotite from the Iberia Abyssal Plain
(ODP Leg 173): its character intermediate between sub-oceanic and subcontinental upper mantle. Geological Society of London, Special Publication 187,
143–159.
Bassi, G., 1995. Relative importance of strain rate and rheology for the mode of
continental extension. Geophysical Journal International 122, 195–210.
Boillot, G., Winterer, E., 1988. Drilling on the Galicia margin: retrospect and prospect.
Proceedings Ocean Drilling Program. Scientific Results, vol. 103, pp. 809–828.
Boillot, G., Coulon, C., 1998. La déchirure continentale et l'ouverture océanique. Gordon
and Breach Science publishers, Amsterdam. 208 pp.
Boillot, G., Comas, M.C., Girardeau, J., Kornprobst, J., Loreau, J.P., Malod, J., Mougenot, D.,
Moullade, M., 1988. Preliminary results of the Galinaute cruise: dives of the
submersive Nautile on the western Galicia margin. Proceedings Ocean Drilling
Program. Scientific Results, vol. 103, pp. 37–53.
Boillot, G., Feraud, G., Recq, M., Girardeau, J., 1989. Undercrusting by serpentinite
beneath rifted margins. Nature 341, 523–525.
Bown, J.W., White, R.S., 1995. Effect of finite extension rate on melt generation at rifted
continental margins. Journal of Geophysical Research 100, 18,011–18,029.
Buck, W.R., Martinez, F., Steckler, M., Cochran, J., 1988. Thermal consequences of
lithospheric extension: pure and simple. Tectonics 7, 213–234.
Bullock, A.D., Minshull, T.A., 2005. From continental extension to seafloor spreading:
crustal structure of the Goban Spur rifted margin, SW of the UK. Geophysical
Journal International 163, 527–546.
Chalmers, J.A., Pulvertaft, T.C.R., 2001. Development of the continental margins of the
Labrador Sea: a review. In: Wilson, R.C.L., Whitmarsh, R.B., Taylor, B., Froitzheim, N.
(Eds.), Non-Volcanic Rifting of Continental Margins: A Comparison of Evidence
from Land and Sea. Geological Society of London, Special Publication, vol. 187,
pp. 77–105.
Chian, D., Louden, K., Reid, I., 1995. Crustal structure of the Labrador Sea conjugate
margin and implications for the formation of nonvolcanic continental margins.
Journal of Geophysical Research 100, 24,239–24,253.
Chian, D., Louden, K.E., Minshull, T.A., Whitmarsh, R.B., 1999. Deep structure of the
ocean–continent transition in the southern Iberia abyssal plain from seismic
refraction profiles; Ocean Drilling Program (legs 149 and 173) transect. Journal of
Geophysical Research 104, 7443–7462.
Contrucci, I., Nercessian, A., Béthoux, N., Mauffret, A., Pascal, G., 2001. A Ligurian
(Western Mediterranean Sea) geophysical transect revisited. Geophysical Journal
International 146, 74–97.
Contrucci, I., Klingelhofer, F., Perrot, J., Bartolome, R., Gutscher, M.A., Sahabi, M., Malod,
J., Rehault, J.P., 2004. The crustal structure of the NW Moroccan continental margin
from wide-angle and reflection seismic data. Geophysical Journal International 159,
117–128.
Davis, M., Kusznir, N.J., 2004. Depth-dependent lithospheric stretching at rifted
continental margins. In: Karner, G.D. (Ed.), Proceedings of NSF Rifted Margins
Theoretical Institute. Columbia University Press, pp. 92–136.
Davison, I., 1997. Wide and narrow margins of the Brazilian South Atlantic. Journal of
Geological Society of London 154, 471–476.
Dean, S.M., Minshull, T.A., Whitmarsh, R.B., Louden, K.E., 2001. Deep structure of the
ocean–continent transition in the southern Iberia abyssal plain from seismic
refraction profiles; the IAM-9 transect at 40 degrees 20′N. Journal of Geophysical
Research 105, 5859–5885.
De Charpal, O., Guennoc, P., Montadert, L., Roberts, D.G., 1978. Rifting and subsidence of
the northern continental margin of the Bay of Biscay. Nature 275, 706–711.
Driscoll, N.W., Karner, G.D., 1998. Lower crustal extension across the Northern
Carnarvon basin, Australia: evidence for an eastward dipping detachment. Journal
of Geophysical Research 103, 4975–4991.
Edwards, R.A., Whitmarsh, R.B., Scrutton, R.A., 1997. The crustal structure across the
transform continental margin of Ghana, eastern equatorial Atlantic. Journal of
Geophysical Research 102, 747–772.
Funck, T., Hopper, J.R., Larsen, H.C., Louden, K.E., Tucholke, B.E., Holbrook, W.S., 2003.
Crustal structure of the ocean–continent transition at Flemish Cap: seismic
refraction results. Journal of Geophysical Research 108 (B11) Art No 2531.
Funck, T., Jackson, H.R., Louden, K.E., Dehler, S.A., Wu, Y., 2004. Crustal structure of the
northern Nova Scotia rifted continental margin (Eastern Canada). Journal of
Geophysical Research 109, B09102. doi:10.1029/2004JB003008.
Gernigon, L., Ringenbach, J.C., Planke, S., Le Gall, B., 2004. Deep structures and breakup
along volcanic rifted margins: insights from integrated studies along the outer
Vøring Basin (Norway). Marine and Petroleum Geology 21, 363–372.
Harry, D.L., Sawyer, D., 1992. A dynamic model of extension in the Baltimore Canyon
Trough region. Tectonics 11, 420–436.
Hoffmann, H.J., Reston, T.J., 1992. The nature of the S reflector beneath the Galicia Bank rifted
margin: preliminary results from pre-stack depth migration. Geology 20, 1091–1094.
Hopper, J., Buck, W.R., 1996. The effect of lower crustal flow on continental extension
and passive margin formation. Journal of Geophysical Research 101, 20175–20194.
Hopper, J.R., Funck, T., Tucholke, B.E., Larsen, H.C., Holbrook, W.S., Louden, K.E.,
Shillington, D., Lau, H., 2004. Continental breakup and the onset of ultraslow
seafloor spreading off Flemish Cap on the Newfoundland rifted margin. Geology 32,
93–96.
Huismans, R.S., Beaumont, C., 2002. Asymmetric lithospheric extension; the role of
frictional plastic strain softening inferred from numerical experiments. Geology 30,
211–214.
Jackson, J.A., White, N.J., 1989. Normal faulting in the upper continental crust: observation
from regions of active extension. Journal of Structural Geology 11, 15–36.
Keen, C., Peddy, C., de Voogd, B., Matthews, D., 1989. Conjugate margins of Canada and
Europe; results from deep reflection profiling. Geology 17, 173–176.
Krawczyk, C.M., Reston, T.J., Beslier, M.O., Boillot, G., 1996. In: Whitmarsh, R., Sawyer, D.,
Klaus, A. (Eds.), Evidence for detachment tectonics on the Iberia Abyssal Plain
margin. Proc ODP, Scientific Results, vol. 149. Ocean Drilling Program, College
Station (TX), pp. 603–615.
Kusznir, N., Park, R.G., 1987. The extensional strength of continental lithosphere: its
dependence on geothermal gradient, and crustal composition and thickness. In:
Coward, M.P., Dewey, J.F., Hancock, P.L. (Eds.), Continental Extensional tectonics.
Geological Society of London, Special Publication, vol. 28, pp. 35–52.
Kusznir, N., Karner, G., 2007. Continental lithospheric thinning and breakup in response
to upwelling divergent mantleflow: application to the Woodlark, Newfoundland
and Iberia margins. Geological Society of London, Special Publication 282, 389–419.
Lau, K.W.H., Louden, K.E., 2007. Systematic patterns of asymmetric extension and breakup
based on seismic models of the Iberia–Newfoundland and other North Atlantic nonvolcanic conjugate margins. Geophysical Research Abstracts 9, 04527 2007.
Lau, K.W.H., Louden, K.E., Funck, T., Tucholke, B.E., Holbrook, W.S., Hopper, J.R., Larsen,
H.C., 2006. Crustal structure across the Grand Banks–Newfoundland Basin
continental margin (part I)—results from a seismic refraction profile. Geophysical
Journal International 167, 127–156.
Lavier, L.L., Manatschal, G., 2006. A mechanism to thin the continental lithosphere at
magma-poor margins. Nature 440, 324–328.
Lenardic, A., Moresi, L., 2000. A new class of equilibrium geotherms in the deep thermal
lithosphere of the continents. Earth and Planetary Science Letters 176, 331–338.
Le Pichon, X., Sibuet, J.C., 1981. Passive margins; a model of formation. Journal of
Geophysical Research 86, 3708–3720.
Leythaeuser, T., Reston, T.J., Minshull, T.A., 2005. Waveform inversion of the S reflector
west of Spain; fine structure of a detachment fault. Geophysical Research Letters 32,
L22304. doi:10.1029/2005GL024026.
Lister, G.S., Etheridge, M.A., Symonds, P.A., 1991. Detachment models for the formation
of passive continental margins. Tectonics 10, 1038–1063.
Maillard, A., Malod, J., Thiébot, E., Klingelhoefer, F., Réhault, J.-P., 2006. Imaging a
lithospheric detachment at the continent–ocean crustal transition off Morocco.
Earth and Planetary Science Letters 241, 686–698.
Manatscha, l.G., Froitzheim, N., Rubenach, M., Turrin, B.D., 2001. The role of detachment
faulting in the formation of an ocean–continent transition: insights from the Iberia
T.J. Reston / Tectonophysics 468 (2009) 6–27
Abyssal Plain. In: Wilson, R.C.L., Whitmarsh, R.B., Taylor, B., Froitzheim, N. (Eds.),
Non-Volcanic Rifting of Continental Margins: A Comparison of Evidence from Land
and Sea. Geological Society of London, Special Publication, vol. 187, pp. 405–428.
Manatschal, G., 2004. New models for evolution of magma-poor rifted margin based on
a review of data and concepts from West Iberia and the Alps. International Journal
of Earth Sciences 93, 432–466.
Mauffret, A., Montardet, L., 1998. Seismic stratigraphy off Galicia (Spain). Proc. Ocean
Drill. Program Sci. Results 103, 13–20.
McKenzie, D.P., 1978. Some remarks on the development of sedimentary basins. Earth
and Planetary Science Letters 40, 25–32.
Minshull, T.A., Dean, S.M., White, R.S., Whitmarsh, R.B., 2001. Anomalous melt
production after continental break-up in the southern Iberia Abyssal Plain. In:
Wilson, R.C.L. (Ed.), Non-volcanic rifting of continental margins: A comparison of
evidence from land and sea. Geological Society of London, Special Publication, vol. 187,
pp. 537–550.
Müntener, O., Manatschal, G., 2006. High degrees of melt extraction recorded by spinel
harzburgite of the Newfoundland margin: the role of inheritance and consequences
for the evolution of the southern North Atlantic. Earth and Planetary Science Letters
252, 437–452.
Nicholls, I.A., Ferguson, J., Jones, H., Marks, G.P., Mutter, J.C., 1981. Ultramafic blocks from
the seafloor southwest of Australia. Earth and Planetary Science Letters 56,
362–374.
Nielsen, T.K., Hopper, J.R., 2002. Formation of volcanic rifted margins: are temperature
anomalies required? Geophysical Research Letters 29 Article 2022.
O'Reilly, B.M., Hauser, F., Jacob, A.W.B., Shannon, P.M., 1996. The lithosphere below the
Rockall Trough: wide-angle seismic evidence for extensive serpentinization.
Tectonophysics 255, 1–23.
O'Reilly, B., Hauser, F., Ravaut, C., Shannon, P., Readman, P., 2006. Crustal thinning,
mantle exhumation and serpentinization in the Porcupine Basin, offshore Ireland:
evidence from wide-angle seismic data. Journal of the Geological Society, London
163, 775–787.
Pérez-Gussinyé, M., Reston, T.J., 2001. Rheological evolution during extension at passive
non-volcanic margins: onset of serpentinization and development of detachments
to continental break-up. Journal of Geophysical Research 106, 3691–3975.
Pérez-Gussinyé, M., Reston, T.J., Phipps Morgan, J., 2001. Rheological and magmatic
evolution during extension at passive non-volcanic margins: the effect of initial
lithospheric structure. In: Wilson, R.C.L., Whitmarsh, R.B., Taylor, B., Froitzheim, N.
(Eds.), Non-Volcanic Rifting of Continental Margins: A Comparison of Evidence
from Land and Sea. Geological Society of London, Special Publication, vol. 187,
pp. 551–576.
Pérez-Gussinyé, M., Ranero, C.R., Reston, T.J., Sawyer, D., 2003. Structure and
mechanisms of extension at the Galicia Interior Basin, west of Iberia. Journal of
Geophysical Research 108. doi:101029/2001JB000901.
Pérez-Gussinyé, M., Phipps Morgan, J., Reston, T.J., Ranero, C.R., 2006. The rift to drift
transition at non-volcanic margins: insights from numerical modelling. Earth and
Planetary Science Letters 244, 458–473.
Phipps Morgan, J., Morgan, W.J., Zhang, Y.-S., Smith, W.H.F., 1995. Observational hints for
a plume-fed, suboceanic asthenosphere and its role in mantle convection. Journal of
Geophysical Research 100, 12,753–12,767.
Pickup, S., Whitmarsh, R.B., Fowler, C.M.R., Reston, T.J., 1996. Insight into the nature of
the ocean–continent transition off West Iberia from a deep multichannel seismic
reflection profile. Geology 24, 1079–1082.
Proffett, J.M., 1977. Cenozoic geology of the Yerington district, Nevada, and implications
for the nature of Basin and Range faulting. Geological Society of America Bulletin
88, 247–266.
Reid, I.D., 1994. Crustal structure of a nonvolcanic rifted margin east of Newfoundland.
Journal of Geophysical Research 99, 15,161–15,180.
Reston, T.J., 1988. Evidence for shear zones in the lower crust offshore Britain. Tectonics
7, 929–945.
Reston, T.J., 1993. Evidence for extensional shear zones in the mantle offshore Britain,
and their implications for the extension of the continental lithosphere. Tectonics 12,
492–506.
Reston, T.J., 1996. The S reflector west of Galicia: The seismic signature of a detachment
fault. Geophysical Journal International 127, 230–244.
Reston, T.J., 2005. Polyphase faulting during the development of the west Galicia rifted
margin. Earth and Planetary Science Letters 237, 561–576.
Reston, T.J., 2007a. The extension discrepancy at North Atlantic non-volcanic rifted
margins: depth-dependent stretching or unrecognised faulting? Geology 35,
367–370.
Reston, T.J., 2007b. The formation of non-volcanic rifted margins by the progressive
extension of the lithosphere. In: Karner, G.D., Manatschal, G., Pinheiro, L.M. (Eds.),
Imaging, Mapping and Modelling Continental Lithosphere Extension and Breakup.
Geological Society of London, Special Publication, vol. 282, pp. 77–110. doi:10.1144/
SP282.5.
Reston, T.J., Phipps Morgan, J., 2004. The continental geotherm and the evolution of
rifted margins. Geology 32, 133–136.
Reston, T.J., Pérez-Gussinyé, M., 2007. Lithospheric extension from rifting to continental
breakup at magma-poor margins: rheology, serpentinisation and symmetry.
International Journal of Earth Sciences. doi:10.1007/s00531-006-0161-z.
27
Reston, T.J., Pennell, J., Stubenrauch, A., Walker, I., Pérez-Gussinyé, M., 2001. Detachment
faulting, mantle serpentinization and serpentinite mud volcanism beneath the
Porcupine Basin SW Ireland. Geology 29, 587–590.
Reston, T.J., Weinrebe, W., Grevemeyer, I., Flueh, E.R., Mitchell, N.C., Kirstein, L., Kopp, C.,
Kopp, H., 2002. A rifted inside corner massif on the Mid-Atlantic ridge at 5°S. Earth
and Planetary Science Letters 200, 255–269.
Reston, T.J., Gaw, V., Pennell, J., Klaeschen, D., Stubenrauch, A., Walker, I., 2004. Extreme
crustal thinning in the south Porcupine Basin and the nature of the Porcupine
Median High: implications for the formation of non-volcanic rifted margins. Journal
of Geological Society of London 161, 783–798.
Reston, T.J., Leythaeuser, T., Booth-Rea, G., Sawyer, D., Klaeschen, D., Long, C., 2007.
Movement along a low-angle normal fault. The S reflector west of Spain.
Geochemistry Geophysics Geosystems 8, Q06002. doi:10.1029/2006GC001437.
Sawyer, D.S., Coffin, M.F., Reston, T.J., Stock, J.M., Hopper, J.R., 2007. COBBOOM: The
Continental Breakup and Birth of Oceans Mission. Scientific Drilling, No. 5,
September 2007.
Sibuet, J.C., 1992. New constraints on the formation of the nonvolcanic continental
Galicia–Flemish Cap conjugate margins. Journal of Geological Society of London
149, 829–840.
Sibuet, J.C., Srivastava, S.P., Enachescu, M., Karner, G.D., 2007. Early Cretaceous motion of
Flemish Cap with respect to North America: implications on the formation of
Orphan Basin and SE Flemish Cap–Galicia Bank conjugate margins. Geological
Society of London, Special Publication 282, 63–76.
Thinon, I., Matias, L., Rehault, J.P., Hirn, A., Fidalgo-Gonzalez, L., Avedik, F., 2003. Deep
structure of the Armorican Basin (Bay of Biscay): a review of Norgasis seismic
reflection and refraction data. Journal of the Geological Society of London 160,
99–116.
Tucholke, B., Lin, J., 1994. A geological model for the structure of ridge segments in slow
spreading oceanic crust. Journal of Geophysical Research 99, 11937–11958.
Tucholke, B.E., Sawyer, D.S., Sibuet, J.-C., 2007. Breakup of the Newfoundland–Iberia Rift.
Geological Society of London, Special Publication 282, 9–46.
Tullis, T.E., Horowitz, F.G., Tullis, J., 1991. Flow law of polyphase aggregates from end–
member flow laws. Journal of Geophysical Research 96, 8081–8096.
Turner, J.P., Rosendahl, B.R., Wilson, P.G., 2003. Structure and evolution of an obliquely
sheared continental margin: Rio Muni, West Africa. Tectonophysics 374, 41–55.
Van Avendonk, H.J.A., Holbrook, W.S., Nunes, G.T., Shillington, D.J., Tucholke, B.E.,
Louden, K.E., Larsen, H.C., Hopper, J.R., 2006. Seismic velocity structure of the rifted
margin of the eastern Grand Banks of Newfoundland, Canada. Journal of
Geophysical Research 111, B11404. doi:10.1029/2005JB004156, 2006.
Voorhoeve, H., Houseman, G., 1988. The thermal evolution of lithosphere extending on a
low-angle detachment zone. Basin Research 1, 1–9.
Wernicke, B., 1985. Uniform-sense normal simple shear of the continental lithosphere.
Canadian Journal of Earth Sciences 22, 331–339.
White, N., McKenzie, D.P., 1988. Formation of the “steer's head” geometry of
sedimentary basins by differential stretching of the crust and mantle. Geology 16,
250–253.
Whitmarsh, R.B., Wallace, P.J., 2001. In: Beslier, M.-O., Whitmarsh, R.B., Wallace, P.J.,
Girardeau, J. (Eds.), The rift-to-drift development of the west Iberia nonvolcanic
continental margin: a summary and review of the contribution of Ocean Drilling
Program Leg 173. Proc. ODP, Sci. Results, vol. 173, pp. 1–36.
Whitmarsh, R.B., Sawyer, D.S., Klaus, A., Masson, D.G., 1996. Editors: Proceedings of the
Ocean Drilling Program. Scientifc Results, vol. 149.
Whitmarsh, R.B., Dean, S.M., Minshull, T.A., Tompkins, M., 2000. Tectonic implications of
exposure of lower continental crust beneath the Iberia Abyssal Plain, Northeast
Atlantic Ocean: Geophysical evidence. Tectonics 19, 919–942.
Whitmarsh, R.B., Manatschal, G., Minshull, T.A., 2001. Evolution of magma-poor
continental margins from rifting to seafloor spreading. Nature 413, 150–154.
Wilson, R.C.L., Wise, S., Manatschal, G., 2001. Non-Volcanic Rifting of Continental
Margins: A Comparison of Evidence from Land and Sea. In: Wilson, R.C.L.,
Whitmarsh, R.B., Taylor, B., Froitzheim, N. (Eds.), Non-Volcanic Rifting of
Continental Margins: A Comparison of Evidence from Land and Sea. Geological
Society of London, Special Publication, vol. 187, pp. 429–452.
Wilson, P.G., Turner, J.P., Westbrook, G.K., 2003. Structural architecture of the ocean–
continent boundary at an oblique transform margin through deep-imaging seismic
interpretation and gravity modelling: Equatorial Guinea, West Africa. Tectonophysics 374, 41–55.
Wu, Y., Louden, K.E., Funck, T., Jackson, H.R., Dehler, S.A., 2006. Crustal structure of the
central Nova Scotia margin off Eastern Canada. Geophysical Journal International
166, 878–906.
Zelt, C.A., Sain, K., Naumenko, J.V., Sawyer, D.S., 2003. Assessment of crustal velocity
models using seismic refraction and reflection tomography. Geophysical Journal
International 153, 609–626.
Ziegler, P.A., 1983. Comment on: Crustal thinning and subsidence in the North Sea.
Nature 304, 561.