The Rapoport effect and the climatic variability ... - Wiley Online Library

[Palaeontology, 2018, pp. 1–18]

THE RAPOPORT EFFECT AND THE CLIMATIC VARIABILITY HYPOTHESIS IN EARLY JURASSIC AMMONITES  LAFFONT 1 , by AXELLE ZACA€I1 , 2 , ARNAUD BRAYARD 1 , R EMI JEAN-LOUIS DOMMERGUES 1 , CHRISTIAN MEISTER 3 and EMMANUEL FARA 1 1

Biogeosciences, UMR 6282, CNRS, Universite Bourgogne Franche-Comte, 6 Boulevard Gabriel, 21000, Dijon, France; [email protected] Current address: UMR 8198, CNRS, Universite de Lille, Avenue Paul Langevin, 59655, Villeneuve d’Ascq, France 3 Departement de Geologie et de Paleontologie, Museum d’Histoire Naturelle de Geneve, 1 Route de Malagnou, 1211, Geneve 6, Switzerland 2

Typescript received 8 February 2018; accepted in revised form 30 May 2018

alternative climatic variability hypothesis that better captures latitudinal seasonal variations and outlines the influence of the intertropical zone, characterized by stable and homogeneous climate that allows species to disperse over very large areas, regardless of their thermal tolerance. Accordingly, the NWE province probably displayed a gradient of seasonal climatic variations which caused the emergence of a Rapoport effect, whereas the MED province was probably located in the intertropical zone where no gradient in species range size is expected. Our multiscale approach further shows that the Rapoport effect is scale-dependent and may be labile through time. This probably explains the conflicting results of previous studies carried out at various spatiotemporal scales.

Abstract: The increase in species range size towards high latitudes, known as Rapoport’s rule, remains one of the most debated and poorly understood macroecological patterns. Numerous studies have challenged both its universality and the main mechanism originally proposed to explain it: the climatic variability hypothesis. Here, we study this pattern using a group of fossil marine organisms: the early Pliensbachian ammonites of the western Tethys. We further take into account the influence of the marked provincialism prevailing at that time, with a Mediterranean province (MED) and a north-west European province (NWE) located on each side of a latitudinally oriented palaeobiogeographical barrier. We find that only species from the NWE province display a Rapoport effect, whereas species from the more tropical MED province show a boundary effect and have larger range sizes on average. This dual pattern can be explained by an

Key words: Rapoport effect, range size, latitude, ammonites, Jurassic, climatic variability hypothesis.

T H E most pervasive patterns in biogeography are those related to latitude (Wallace 1878; Fischer 1960; Hubbell 2001; Blackburn & Gaston 2003; Willig et al. 2003). Among those, the positive correlation between latitude and species geographical range size, originally observed in mammalian subspecies by Rapoport (1982), is one of the most intriguing and poorly understood patterns. Since its establishment as ‘Rapoport’s rule’ by Stevens (1989), this trend of increasing geographical range size towards high latitudes has received much attention. It has been reported for numerous organisms including trees, mammals, birds, marine molluscs, reptiles, freshwater amphipods and fishes, but mostly in the northern hemisphere and for terrestrial organisms (Gaston et al. 1998; Gaston 2003; Table 1). Stevens (1989) initially related this pattern to the existence of latitudinal diversity gradients and proposed an explanatory mechanism called the ‘climatic variability hypothesis’. Based on Janzen’s (1967) work on

mountain passes, it assumes that greater seasonal climatic variability at high latitudes favours species with broad environmental tolerances and thus broad range sizes. Conversely, more stable climates observed between the tropics are expected to be favourable to species with narrower environmental tolerances and thus small range sizes, thereby allowing more species to coexist at low latitudes. Intuitively attractive and supported by deep historical roots (Hutchins 1947; Dobzhansky 1950; Janzen 1967; Brattstrom 1968), the climatic variability hypothesis has largely contributed to the popularization of Rapoport’s rule. However, many studies have challenged the generality of this ‘rule’ (Rohde et al. 1993; Roy et al. 1994; Rohde 1996, 1999; Gaston et al. 1998; Gaston & Chown 1999). Indeed, the increase in species range size with latitude seems inconsistent under tropical latitudes (Stevens 1989; Rohde et al. 1993; Letcher & Harvey 1994; Ruggiero 1994; Gaston et al. 1998; Sax 2001) and

doi: 10.1111/pala.12389 © 2018 The Authors. Palaeontology published by John Wiley & Sons Ltd on behalf of The Palaeontological Association. This is an open access article under the terms of the Creative Commons Attribution License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

1

2

PALAEONTOLOGY

equivocal in the southern hemisphere (Rohde 1996; Gaston et al. 1998; Gaston & Chown 1999; Cardillo 2002; Chown et al. 2004; Hernandez Fernandez & Vrba 2005). Some authors have thus suggested that it could be a ‘local’ phenomenon that essentially applies to the northern hemisphere, out of the tropics (Rohde et al. 1993; Rohde 1996; Blackburn & Gaston 1996; Gaston et al. 1998; Gaston & Chown 1999). For this reason, Blackburn & Gaston (1996) advocated replacing the term ‘Rapoport’s rule’ with ‘Rapoport effect’; a terminology that we follow in this work. The debate is still open, but as noted by Gaston & Chown (1999), it is hindered by the non-randomness of the published case studies, which have mainly focused on extant terrestrial vertebrates from the Nearctic, Palaearctic and Australasian biogeographical realms, and have barely investigated low-latitude and deep-time patterns (Table 1). In addition, some confusion has arisen from the diversity of approaches used to investigate the Rapoport effect: (1) empirical tests of the climatic variability hypothesis with known thermal latitudinal gradients (e.g. Morin & Lechowicz 2011; Pintor et al. 2015; Tomasovych et al. 2015; Stuart-Smith et al. 2017); (2) empirical tests with both thermal gradient and species thermal tolerances known (e.g. Khaliq et al. 2017; Gr€ unig et al. 2017); and (3) predictions of a Rapoport effect tacitly assuming (or not) the presence of a thermal latitudinal gradient (e.g. Stevens 1989; Rohde 1996; Fortes & Absal~ao 2004; Weiser et al. 2007; Veter et al. 2013). Here, we test for the existence of a Rapoport effect among fossil representatives of a marine group of invertebrates: the ammonites (Cephalopoda). Our study focuses on a c. three million-year-long interval in the Early Jurassic (the early Pliensbachian), and covers what was a vast tropical epicontinental sea at that time (present-day Europe, the Middle East and North Africa; Fig. 1). Such a case study presents several advantages. First, it uses one of the most highly resolved fossil datasets available, because the taxonomy of early Pliensbachian ammonites has been thoroughly revised up to the species level by Dommergues et al. (2009) and species occurrences are documented at the scale of the chronozone (i.e. c. 1 million-year-long intervals). As such, it provides a high-resolution palaeontological framework in which to investigate the Rapoport effect through time, hence offering more possibilities than the binary response that can be obtained when studying extant organisms (i.e. existence/absence of a Rapoport effect). This is especially interesting because the end of the early Pliensbachian records climatic changes that affected the phylogenetic conservatism of ammonite range sizes (Zaca€ı et al. 2017) and could also have modified their latitudinal distribution. Second, the singular configuration of the area under study, divided into two contiguous marine biogeographical provinces (Dommergues 1982; Dommergues & Meister 1991;

Meister & Stampfli 2000; Dommergues et al. 2009), allows us to test for the influence of a latitudinally oriented biogeographical barrier on the geographical distribution of species range sizes and on the existence of a Rapoport effect. Third, assuming the validity of the climatic variability hypothesis for extant organisms, the dataset seems appropriate to test for the applicability of this hypothesis and the existence of a Rapoport effect in deep time, because the environmental and supposedly climatic contrast that characterized the two provinces is known to have influenced both ammonite dispersal (Zaca€ı et al. 2016) and phylogenetic conservatism of species range size (Zaca€ı et al. 2017). Fourth, it extends the study of the Rapoport effect to the marine realm and to the fossil record, both of which have been largely neglected in previous studies (Table 1). Our goal is thus to test for the existence of a Rapoport effect in ammonites and for the applicability of the climatic variability hypothesis in a fossil marine environment, and to estimate the influence of the singular geographical and environmental configuration of the early Pliensbachian western Tethys on the latitudinal distribution of species range size.

MATERIAL AND METHOD Dataset We used Dommergues et al.’s (2009) dataset of early Pliensbachian ammonites (Early Jurassic, c. 191.36
1 Ma to 187.6
1 Ma; Gradstein et al. 2012; Ogg et al. 2016; Fig. 1A). It consists of 214 revised species distributed among 104 fossil localities of the western Tethys and adjacent areas (present-day Europe, the Middle East and North Africa; Fig. 1B). Palaeocoordinates for each fossil locality were calculated following the procedure thoroughly described by Zaca€ı et al. (2017). The three successive chronozones of the early Pliensbachian, that is Jamesoni, Ibex and Davoei (Page 2003; Fig. 1A), are represented, respectively, by 88, 94 and 80 localities, at a total of 104 localities. Species occurrences are documented at the scale of the chronozone. The early Pliensbachian is a time of marked faunal provincialism between two marine areas: the Mediterranean Tethys (MED) in the south, located approximately between 25°N and 45°N, and north-west Europe (NWE) in the north, located approximately between 40°N and 65°N (Fig. 1B). These two provinces have been identified for ammonites (Donovan 1967; Hallam 1969; Dommergues 1982; Geczy 1984; Dommergues & Meister 1991; Meister & Stampfli 2000; Dommergues et al. 2009), bivalves (Hallam 1977; Liu et al. 1998), brachiopods (Ager 1973; V€ or€ os

Z A C A €I E T A L . : R A P O P O R T E F F E C T I N E A R L Y J U R A S S I C A M M O N I T E S TABLE 1.

3

Synthesis of a survey of published studies treating the relationship between range size and latitude.

Taxon

Region

R*

M†

P‡

Lat. range§

Support¶

References

Birds Birds Birds Birds

North America Global Global New World

a l l,a l

m m s,m s

x x U –

25–60°N ? c. 90°S–90°N c. 60°S–80°N

Brown (1995, p. 113) Cardillo (2002) Orme et al. (2006) Tomasovych et al. (2016)

Birds Non-migratory birds Birds

Soviet Union Soviet Union

l l

s s

– –

40–80°N 40–80°N

25–40°N U c. 0–90°N 0–60°S (genus level) – 40–80°N

New World

a

s

–

55°S–70°N

17–70°N

Blackburn & Gaston (1996)

Wildfowl Phylloscopus leaf warblers Mammals and birds Mammals Mammals

Global Old World

a a l

m m m

– U U

55°S–70°N 50°S–75°N 15–65°N

U – 15–65°N

Price et al. (1997)

North America

l

s

–

c. 10–70°N

c. 30–65°N

Sax (2001)

North America North America Palaearctic

Mammals

Australia

Mammals Mammals Mammals

South America North America Africa

s s m s m s m s s m

– – U – U – x – – U

35–80° N 5–80°N 5–80°N 20–80°N 25–70°N 10–40°S 10–40°S c. 0–60°S ? 40°S–40°N

35–80°N 5–80°N U 30–80°N 30–60°N – – c. 40–50°S U 0–40°N

Stevens (1989, 1992) Pagel et al. (1991)

Mammals

l a a l,a l,a l a l l l,a

Mammals Mammals (weak support) Primates

North America Australia

l l

s,m m

– x

c. 0°–70°N 11–39°N

– 11–39°N

Africa

l

m

x

c. 37°S–37°N

c. 0–37°N

Bats

New World

Forest mammals Fossil mammals (Cenozoic) Carnivora, primates Artiodactyla, Hystricognathi, marsupials Marsupials

Australia North America

l l l l

s m m m

– x x –

60°S–65°N 45°S–45°N 11–39°S c. 25°–60°N

60°S–65°N 15°N–45°N 11–39°S Pleistocene

South America South America

l l

s s

– –

55°S–10°N 55°S–10°N

55°S–10°N 55°S–10°N**

Ruggiero (1994)

New World Europe

s m s

– x –

60°S–45°N 55°S–25°N c. 35–75°N

0°–45°N ? – c. 40–65°N

Lyons & Willig (1997)

Mammals and vascular plants Reptiles and amphibians Lizards Snakes Brachinus beetles Moths

l l l

Sax (2001)

North America

l

s

–

30–65°N

30–60°N

Stevens (1989, 1992)

Australia Global North America Switzerland

s,m s/m m m

– x x x

40°S–10°S c. 90°S–90°N ? 600–2400 m

2 out of 3 clades c. 90°S–90°N U –

Pintor et al. (2015) B€ ohm et al. (2017) Juliano (1983) Beck et al. (2016)

Human pathogens

Global

l a l,a Elevational range l

m

–

c. 35°S–60°N

c. 20°S–30°N

Guernier & Guegan (2009)

Stevens (1989)

Letcher & Harvey (1994) Smith et al. (1994) Sax (2001) Arita et al. (2005) Hernandez Fernandez & Vrba (2005) Lyons & Smith (2010) Johnson (1998) Cowlishaw & Hacker (1997) Lyons & Willig (1997) Johnson (1998) Veter et al. (2013)

(continued)

4

PALAEONTOLOGY

TABLE 1.

(Continued)

Taxon

Region

R*

M†

P‡

Lat. range§

Support¶

References

Plants

New World

a

–

–

?

Trees Trees

North America North America and Europe North America North America

l l

s m

– U

25–70°N 15°–71°N

North America 25–70°N 15°–71°N

Morueta-Holme et al. (2013) Stevens (1989, 1992) Morin & Chuine (2006)

l l,a

m m

– x

30°–55°N 6.4°N–74.3°N

30°–55°N U

Temperate rainforest trees Woody plants Acacia

South America

l

s

x

38–54°S

38–54°S

Lane (2007) Morin & Lechowicz (2011) Arroyo et al. (1996)

New World Australia

l l,a

s,m m

(U)

c. 60°S–75°N ?

c. 40°–60°N –

Eucalyptus Eucalypts Pinus Columnar cacti

Australia Australia North America Argentina

l,a l,a l l,a

m m s m

(U) x – x

? 10–45°S 5–70°N 23–33°S

– – 15–70°N –

Columnar cacti

Mexico

Fishes Freshwater fishes

North America North America

Freshwater fishes

Northern Europe

Freshwater fishes

Australia

Cyprinella minnows Freshwater amphipods Freshwater crayfish Freshwater and coastal fishes Marine fishes Marine teleost fishes

North North North North

l a l l l l l l l a l l l

m m s s m s m s m m s s s

x x – – x – x – x (U) – – –

15–28°N 15–28°N c. 30–60°N 0°–85°N 0°–70°N 20–75°N 40–65°N 45–10°S 45–10°S ? 30–70°N 25–65°N 25–80°N

15–28°N U c. 30–60°N 35–85°N 25–60°N 45–75°N 45–55°N – – U 30–70°N 25–65°N 35–80°N

Marine teleost fishes

Indo-Pacific

Marine teleost fishes

Various

Teleosts

East Atlantic

l l l l l l l l

s s m s m s m s

– – x – x – x –

30–65°N 55°S–90°N 45°S–75°N 55°S–75°N 40°S–50°N 55°S–75°N 45°S–55°N 35°S–70°N

– – – – – – –†† 35°S–70°N

Corals and fishes Fishes and invertebrates Elasmobranchs

Indo-Pacific Atlantic

l l

m s

– –

c. 35°S–35°N 70°S–80°N

– –

East Atlantic

l

s

–

35°S–70°N

–‡‡

Marine molluscs Marine molluscs

North America Western Atlantic and eastern Pacific Western & eastern Pacific and western Atlantic

l l

s s,m

– –

25–70°N c. 90°S–90°N

25–70°N U

MacPherson & Duarte (1994) Stevens (1989, 1992) Fortes & Absal~ao (2004)

l

s

–

77°S–83°N

–

Tomasovych et al. (2015)

Trees Trees

Marines bivalves

America America America America

Pacific Ocean Atlantic

Weiser et al. (2007) Edwards & Westoby (1996) Hughes et al. (1996) Stevens & Enquist (1998) Mourelle & Ezcurra (1997)

Sax (2001) Rohde et al. (1993)

Taylor & Gotelli (1994) France (1992) Stevens (1989, 1992) Stevens (1996) Rohde & Heap (1996)

Rohde et al. (1993) MacPherson & Duarte (1994) Connolly et al. (2003) MacPherson (2003)

(continued)

Z A C A €I E T A L . : R A P O P O R T E F F E C T I N E A R L Y J U R A S S I C A M M O N I T E S TABLE 1.

5

(Continued)

Taxon

Region

R*

M†

P‡

Lat. range§

Support¶

References

Marine bivalves

Western & eastern Pacific and western Atlantic Western Atlantic and eastern Pacific Western Atlantic and eastern Pacific Atlantic

l

s

–

c. 80°S–80°N

species-rich genera

l

m

–

5°S–65°N ?

–

Tomasovych et al. (2016) and Tomasovych & Jablonski (2017) Roy et al. (1998)

l

s

–

0–70°N

–

Roy et al. (1994)

l& depth

s

–

c. 60°S–70°N

Western Atlantic

Rosa et al. (2008)

Marine prosobranch gastropods Prosobranchs Coastal cephalopods

*Range size measure: l, latitudinal extent; a, range area. †Method used to test for Rapoport effect: s, Stevens method; m, midpoint method. ‡If the midpoint method was used, whether or not the phylogenetic non-independence of species was accounted for; U indicates that a phylogenetic effect was considered but not explicitly controlled for. §The approximate latitudinal range over which Rapoport effect was tested. ¶The approximate latitudinal range over which there was evidence for a Rapoport effect (there is an element of subjectivity in these judgments in several cases); U, there was evidence for the pattern but it is not clear over what latitudes; –, there was no evidence for the pattern. Where this range embraces the equator, there was evidence for a Rapoport effect in both hemispheres. **The pattern is found across all species, but across species in only three of seven orders when analysed separately. ††Separate analyses for families and genera did not find evidence for a Rapoport effect. ‡‡Of separate analyses for elasmobranchs, and for deep-living, littoral and pelagic teleosts, only the last found evidence for a Rapoport effect.

1977), belemnites (Doyle 1987, 1994) and ostracods (Lord 1988; Arias & Whatley 2005). It seems that those benthic and pelagic organisms were unable to cross the latitudinally oriented biogeographical barrier (Fig. 1B), except in two areas located in the MED province where ‘transprovince species’ occur: the Austro-Alps and the Pontides (Meister & Stampfli 2000; Dommergues et al. 2009). The proportion of transprovince species (i.e. species that occur in both provinces) tends to increase through time: 16% during the Jamesoni chronozone, 23% during the Ibex chronozone and 31% during the Davoei chronozone. This increasing trend may be related to: (1) the extinction of endemic species during the Davoei chronozone (Dommergues & El Hariri 2002; Hardy et al. 2012), which could have allowed the establishment of cosmopolitan taxa in vacated niches (Aberhan & F€ ursich 1996, 2000; Dera et al. 2011a); and (2) ingressions of MED species into the NWE province (Meister & Stampfli 2000; Dommergues et al. 2001), probably facilitated by the incursion of eastern Tethyan waters into the NWE province (Dera et al. 2009a, 2015). This concordance between organisms which are physiologically and phylogenetically unrelated and have markedly differing modes of life suggests common regional palaeobiogeographical constraints (Dera et al. 2011a). In addition to the faunal differentiation, the difference in prevailing facies on both sides of the boundary, which indicates that both environmental and tectonic constraints were clearly distinct in the NWE and MED provinces (Dommergues et al. 2009), also supports the hypothesis of a physical boundary

between the two provinces. The NWE province corresponds to an epicontinental sea not exceeding 200 m in depth, surrounded by continents and punctuated by numerous archipelagos (Fig. 1B). This singular configuration led to a strong influence of numerous perturbations, such as freshwater influxes, nutrient inputs, anoxia, incursions of cold and warm waters and changes in water currents (Rosales et al. 2004a; Wignall et al. 2005; Arias 2007; McArthur et al. 2008; Dera et al. 2009a, b, 2011a; Dommergues et al. 2009). Conversely, the MED province corresponds to a series of deep interconnected basins influenced by warm and nutrient-rich surface waters from the Tethys Ocean and therefore represents a physically more stable open-marine palaeoenvironment (Hallam 1972; Smith & Tipper 1986; Dommergues et al. 2009; Dera et al. 2011a). For all these reasons, the most largely accepted hypothesis regarding the nature of the biogeographical barrier separating the NWE and MED province is that it comprised a narrow chain of emerged or shallows areas, acting as a full physical barrier or a faunal filter, depending on sea-level fluctuations (Dommergues 1982; Dercourt et al. 1985; Dommergues & Meister 1991; Meister & Stampfli 2000; Dommergues et al. 2009; Dera et al. 2011b).

Relationship between range size and latitude Stevens (1989) initially documented the positive relationship between species latitudinal range size and latitude by

PALAEONTOLOGY

186

188

190

Stage

Ammonite chronozones

B

60°N

Spinatum

Late

Eurasia

50°N

Iberia

Margaritatus

Ad

Anatolia 40°N

Davoei

Early

Age (Ma)

A

PLIENSBACHIAN

6

Ibex

Western Tethys Africa 30°N

Jamesoni

190.8 Deep basins Epicontinental seas

Emerged landmasses Biogeographical boundary

NWE locality MED locality

Biostratigraphical and palaeogeographical context of the early Pliensbachian. A, ammonite chronozones (Page 2003; Gradstein et al. 2012) and their absolute ages according to Gradstein et al. (2012) and Ogg et al. (2016). B, palaeogeographical reconstruction of the western Tethys and adjacent areas for the early Pliensbachian, and palaeolocation of the 104 fossil localities used in this study. NWE localities are represented by cyan circles and MED localities are represented by orange squares. Abbreviations: Ad, Adria bloc; Africa, African bloc; Anatolia, Anatolian bloc; Eurasia, Eurasian bloc; Iberia, Iberian bloc; MED, Mediterranean; NWE, north-west European. Colour online.

FIG. 1.

comparing the average range size of species occurring in 5° latitudinal bands. This approach has been criticized because it suffers from several statistical problems (Rohde et al. 1993; Colwell & Hurtt 1994; Letcher & Harvey 1994; Gaston et al. 1998), the most salient being the non-independence of mean range sizes since large-range species contribute to mean range size at multiple latitudes (Rohde et al. 1993; Letcher & Harvey 1994; Connolly 2009). To circumvent this problem, Rohde et al. (1993) proposed a different approach called the ‘midpoint method’, which associates species latitudinal range sizes with their latitudinal midpoints and treats them as independent data. However, this method is affected by triangular geometric constraints, which implies that: (1) small variations in range size of species at the boundaries of the area under study will have higher weight in the analyses than substantial variations in range size of species distributed more centrally within the study area; and (2) large-range sizes are more susceptible to being located at the centre of the study area (Blackburn & Gaston 1996; Lyons & Willig 1997; Lyons & Smith 2010; Pintor et al. 2015). As noted by Rohde et al. (1993), Fortes & Absal~ao (2004) and Connolly (2009), the midpoint approach can be misleading in some particular cases, especially when there are many species with very large ranges. A

survey of the literature reveals that the Stevens and midpoint methods have had the same appeal to researchers despite their respective advantages and drawbacks (Table 1). About half of the published studies failed to consider the issue of the non-independence of statistical observations. When taken into account, this dependence can be treated either spatially (e.g. Tomasovych et al. 2015, 2016), or phylogenetically (e.g. Morin & Lechowicz 2011, Beck et al. 2016, B€ ohm et al. 2017), or by the mean of null model simulations through bootstrapping approaches (e.g. Arita et al. 2005, Pintor et al. 2015). Overall, at least 17 different statistical approaches have been used to investigate the relationship between range size and latitude, making the published studies hardly comparable (Table 1). Here, we adopt the midpoint method and we control for the phylogenetic dependence of statistical observations. To allow for future investigations, we also provide data matrices corresponding to the Stevens method and to the study of the subsequent spatial non-independence of observations (Zaca€ı et al. 2018, table S1). We analysed the relationship between species range size and latitude using generalized estimating equations (GEE; see Morin & Chuine 2006 and Morin & Lechowicz 2011). GEE fit regression models by taking into account the

Z A C A €I E T A L . : R A P O P O R T E F F E C T I N E A R L Y J U R A S S I C A M M O N I T E S

relatedness among observations (Liang & Zeger 1986). Paradis & Claude (2002) proposed using GEE in the context of comparative analyses, the relatedness among observations (i.e. species) being computed from the phylogeny. These authors also underlined the ability of GEE to deal with phylogenies containing polychotomies, which is the case in our dataset (see Zaca€ı et al. 2018, fig. S1, and Hardy et al. 2012). Here, we fit a GEE model to explain variation in species latitudinal range size (response variable) depending on variation of species latitudinal midpoints (explanatory variable), while removing phylogenetic autocorrelation. Considering the continuous nature of the response variable, we used a Gaussian error distribution for this variable, with an identity link function in the GEE model. This model has the following form (Paradis & Claude 2002): E½yi  ¼ b0 þ b1 xi1 ; where E[yi] is the expectation of the response variable (latitudinal range) for the ith species; b0 and b1 are the coefficients of the model (i.e. intercept and slope, respectively); and xi1 is the explanatory variable (species latitudinal midpoint). We computed this model using the compar.gee function (Paradis & Claude 2002) as implemented in the R package ape (Paradis et al. 2004; R Core Team 2016), with a significance level set at a = 0.05. We determined palaeolatitudinal range size (i.e. the difference, in degrees of latitude, between the northernmost and southernmost points of a geographical range) and midpoint (i.e. mean of these two extremes) of each species for the early Pliensbachian and for the three chronozones of this substage. These range sizes correspond to ‘maximal’ species range sizes as they are calculated based on the sum of species occurrences for the early Pliensbachian and for each chronozone. Thus, they do not represent instantaneous ranges but rather the whole range explored by a species during its lifetime or during a specific time interval (Zaca€ı et al. 2017). We conducted the analyses: (1) at the scale of the western Tethys, with no distinction between species occurring exclusively in one province (hereafter referred as ‘MED species’ and ‘NWE species’) and those occurring in both provinces (transprovince species); and (2) at the province scale, that is, considering either MED or NWE species only. For the temporal scale, we conducted the analyses for the entire early Pliensbachian but also for each of its chronozone in order to track potential changes in Rapoport effect through time. We paid special attention to the Davoei chronozone because it records several climatic and environmental changes, including an increase in sea-surface temperatures (Rosales et al. 2004a; Dera et al. 2009a, b, 2011b; Price et al. 2016), a sea-level rise (Hallam 2001) and more humid climatic conditions (Dera et al.

7

2009b). For analyses on chronozone/province subdatasets, interspecies phylogenetic distances were calculated on subtrees containing only species occurring in the chronozone/ area of interest (see Zaca€ı et al. 2018). Analyses were conducted on ‘complete’ dataset and subdatasets, meaning that all species were conserved, even those showing only one occurrence and thus null latitudinal range sizes. We also conducted analyses on datasets containing only species showing more than one occurrence in order to evaluate the influence of null range sizes on the relationship between species range size and latitude.

RESULTS At the scale of the early Pliensbachian, latitudinal range sizes of western Tethys species increase significantly with latitude (Fig. 2A; Table 2) and thus show a Rapoport effect. However, when species are separated according to their province of occurrence (i.e. MED, NWE and transprovince species), only NWE species exhibit such a pattern (see Fig. 2C, Table 1 and Zaca€ı et al. (2018, table S2) for analyses of transprovince species). This positive relationship between NWE species range size and latitude could be driven by an outlier, that is, the species showing the largest range at the highest latitude (Beaniceras luridum Simpson, 1855). However, it remains significant even when this species is removed from the dataset analysed (Zaca€ı et al. 2018, table S2). The relationship observed at the scale of the western Tethys is therefore probably driven by NWE species. MED species have their largest latitudinal range sizes near the centre of the province, and those ranges decrease towards the northern and southern province boundaries (Fig. 2B). This distribution pattern certainly results from a triangular geometric constraint (see Material and Method). For the Jamesoni chronozone, only the analysis at the spatial scale of the western Tethys shows a significant Rapoport effect (Fig. 2J; Table 2). When considered separately, neither the NWE and MED species nor the transprovince species display an increase in range size with latitude (Fig. 2K–L; Zaca€ı et al. 2018, table S2). For the Ibex chronozone, only NWE species show a significant Rapoport effect (Fig. 2I; Table 2). However, this relationship is no more significant when the outlier represented by Beaniceras luridum is removed from the analysis (see Zaca€ı et al. 2018, table S2). For the Davoei chronozone, NWE species show a clear increase in their range sizes with latitude but this relationship is only almost significant (p = 0.059; Fig. 2F), probably because of the small sample size (8 species). Conversely, analyses at the scale of the western Tethys and MED province exhibit a significant decrease in species range size with latitude (Fig. 2D– E and Table 2). Results of analyses conducted on datasets

8

PALAEONTOLOGY

Western Tethys

NWE province

MED province

40

Latitudinal range

A

B

C

30

20

10

0 25

30

35

40

45

50

55

25

30

Latitudinal midpoint

35

40

44

45

46

Latitudinal midpoint

48

50

52

54

Latitudinal midpoint

Early Pliensbachian D

Latitudinal range

25

E

F

20 15 10 5 0 35

40

45

50

34

36

38

40

42

46

48

50

Davoei 40

Latitudinal range

G

H

I

30

20

10

0 25

30

35

40

45

50

25

30

35

40

44

46

48

50

52

54

Ibex Latitudinal range

25

J

K

L

20 15 10 5 0 25

30

35

40

45

Latitudinal midpoint

50

55 25

30

35

40

Latitudinal midpoint

Jamesoni

45

46

48

50

52

Latitudinal midpoint

54

Z A C A €I E T A L . : R A P O P O R T E F F E C T I N E A R L Y J U R A S S I C A M M O N I T E S TABLE 2.

9

Results of GEE models testing the relationship between species latitudinal midpoint and range size.

Temporal interval

Early Pliensbachian

Davoei

Ibex

Jamesoni

Spatial scale

n Estimate SE n Estimate SE n Estimate SE n Estimate SE

Western Tethys

MED

NWE

214 b1 = 1.29 9 101** 3.86e 9 102 45 b1 = 9.89 9 101*** 1.92 9 101 103 b1 = 7.56 9 102 7.78 9 102 132 b1 = 2.53 9 101*** 4.33 9 102

110 b1 = 9.42 9 102 6.33 9 102 23 b1 = 2.13*** 2.63 9 101 56 b1 = 1.75 9 102 1.07 9 101 81 b1 = 9.34 9 102 8.08 9 102

52 b1 = 7.64 9 101*** 1.39 9 101 8 b1 = 1.12 3.36 9 101 23 b1 = 1.06** 2.16 9 101 29 b1 = 3.39 9 101 3.99 9 101

Phylogenetic linear regression of species latitudinal range size against latitudinal midpoint for early Pliensbachian ammonites of the western Tethys calculated with GEE. n, number of species considered in the analysis; SE, standard error; b1, slope of the regression; MED, Mediterranean province; NWE, north-west European province. ***p < 0.001; **p < 0.01; *p < 0.05; no symbol = not significant.

containing only species with more than one occurrence are consistent with those conducted on complete datasets (Zaca€ı et al. 2018, table S3), indicating that endemic species do not have a noticeable influence on the overall latitudinal distribution of species range sizes and on the emergence of a Rapoport effect.

DISCUSSION Rapoport effect during the early Pliensbachian and the modified climatic variability hypothesis In the climatic variability hypothesis, Stevens (1989) regarded seasonal climatic variations as the main factor generating a Rapoport effect: lower seasonal climatic variability observed between the tropics is favourable to stenotopic species, while greater seasonal climatic variability at high latitudes favours eurytopic species, thereby resulting in an increase of species range size with latitude. Gaston & Chown (1999) later added the mean annual temperature as an explicative agent of this pattern and also revised its form. Today, mean annual temperatures increase from the poles to a latitude of about 20–25° north and south, after which they level off (Terborgh 1973; Gaston & Chown 1999; Colwell et al. 2008). Thus, environmental/climatic conditions are more homogeneous across the tropics (i.e.

over a very large latitudinal extent). For this reason, the difference of thermal tolerances between tropical and temperate species anticipated under Janzen’s (1967) hypothesis should not produce a Rapoport pattern within the tropics. Indeed, a tropical species with a narrow temperature tolerance may theoretically have a broader latitudinal range than a temperate species with a greater temperature tolerance, given the thermal plateau observed between 20– 25°S and 20–25°N (Terborgh 1973; Gaston & Chown 1999; Colwell 2011; Jablonski et al. 2013). This modified climatic variability hypothesis implies that large latitudinal ranges are also expected at low latitudes and that the decreasing trend of species range size from high to low latitudes should vanish in the intertropical zone. A broad dispersal has been observed in many extant groups of tropical species such as neotropical birds, mammals and reptiles (Terborgh 1973), American wasps (Richards 1978), insects (Gaston 1991), Costa Rican hymenoptera (Gaston & Gauld 1993; Gaston 1996) and marine bivalves (Jablonski et al. 2013; Tomasovych et al. 2015, 2016; Tomasovych & Jablonski 2017). Some groups even show the pattern expected under Gaston & Chown (1999) hypothesis, that is broader latitudinal extents towards lower latitudes, or at least a breakdown of the Rapoport effect towards lower latitudes (e.g. teleost fishes, Rohde et al. 1993; molluscs from the eastern Pacific, Roy et al. 1994; North American birds and mammals, Sax 2001).

Relationship between species latitudinal midpoint and latitudinal range size for early Pliensbachian ammonites of the western Tethys and adjacent areas, estimated with GEE. A–C, analyses at the scale of the substage (early Pliensbachian). D–L, analyses at the scale of the chronozone (Jamesoni, Ibex and Davoei). Trend lines show significant relationships (see Table 2 for details). Colour online.

FIG. 2.

10

PALAEONTOLOGY

However, it is important to note that the majority of these groups are terrestrial, and more generally, that the climatic variability hypothesis has been established for the terrestrial realm (Stevens 1989). The latitudinal distribution of species range sizes for early Pliensbachian ammonites of the western Tethys and adjacent areas is consistent with the pattern corresponding to the modified climatic variability hypothesis of Gaston & Chown (1999). At the scale of the early Pliensbachian, a Rapoport effect is observed for NWE species, that is for species with range midpoints located at mid- to high latitudes (between c. 43°N and 55°N), whereas MED species, that is species with midpoints located at lower latitudes (between c. 25°N and 45°N), do not show any increase in their range size with latitude (Fig. 2; Table 2), but they disperse over larger distances than the northern ones on average (Fig. 3). In the marine realm, the influence of climatic variability is: (1) essentially limited to the epipelagic zone (i.e. between 0 and about 200 m); (2) of lower amplitude than in terrestrial environments; and (3) stronger at midlatitudes corresponding to temperate zones (c. 35–60°N and S) than in polar regions (Helmuth et al. 2002; Parmesan et al. 2005; Clarke 2009). The influence of seasonal climatic variations on water temperature during the early Pliensbachian was thus probably very different between the NWE and MED provinces, because of their latitudinal position but also of their palaeogeographical configuration. The shallowness of the epicontinental sea and the proximity of emerged lands (Fig. 1B) probably made the NWE province much more sensitive to seasonal climate changes (Arias 2007) than the MED province. The latter was probably much less subject to seasonal variations in temperature and much more stable and homogeneous in terms of climate and environments, due to its more tropical position (Fig. 1B), its deep basins and its large opening to the Tethys Ocean (Hallam 1972; Smith & Tipper 1986; Dommergues et al. 2009; Dera et al. 2011a). Hence, the Rapoport effect observed in the NWE province at the scale of the early Pliensbachian (Fig. 2C) probably results from the influence of climatic variability on species latitudinal distribution. On the other hand, the distribution of latitudinal range sizes of MED species, which shows a boundary effect (Fig. 2B), indicates that their range limits were probably mainly constrained not by temperature but by the geomorphological setting of the province, with emerged and nearemerged landmasses delimiting the northern and southern parts of the province (Fig. 1B). Colwell & Hurtt (1994) and Colwell & Lees (2000) have described a similar pattern using null models and show that boundary constraints force a steady decrease in range size as range midpoints become closer to a boundary, regardless of the spatial distribution and size of ranges. Since the latitudinal gradient of temperature during the Jurassic was weaker than today (e.g. Bjerrum et al. 2001; van de Schootbrugge et al. 2005;

Brigaud et al. 2008; Dera et al. 2009b), the northern and southern limits of the intertropical zone were undoubtedly located higher in latitude during the early Pliensbachian. Regarding the studied area, the loss of a Rapoport effect in the intertropical zone expected under Gaston & Chown (1999) hypothesis should thus have been shifted to the north during this interval. Additionally, in our case, the presence of a latitudinally oriented biogeographical barrier largely limiting faunal exchanges in the western Tethys (Fig. 1B) may blur the existence of such a break in the graph resulting from the analysis at the scale of the early Pliensbachian (Fig. 2A). However, considering the distinct patterns of species range size distribution observed between the two provinces (Fig. 2B, C), we assume that the limit of the intertropical zone where the Rapoport effect should have vanished was located approximately at the same latitude as the biogeographical barrier (i.e. between 40°N and 45°N; Fig. 1B). Accordingly, the NWE province would have been characterized by a latitudinal gradient of temperatures and marked seasonal variations, and the MED province would have corresponded to a much more stable and homogeneous environment in terms of temperature fluctuations. Differences in environmental stability during the Pliensbachian have been documented for these two provinces based on sedimentary, geochemical and palaeontological data (Hallam 1972; Smith & Tipper 1986; Dommergues & El Hariri 2002; Wignall et al. 2005; Arias 2007; Dera et al. 2009a, b; Dommergues et al. 2009; Dera et al. 2011a). Although several studies have estimated seawater palaeotemperatures for the NWE province (Saelen et al. 1996; Rosales et al. 2004a, b; van de Schootbrugge et al. 2005; Arias 2007; Dera et al. 2009a, 2011a), they did not formally establish the existence of a temperature gradient in this province during the early Pliensbachian. As for the MED province, samples (oyster shells and belemnite rostra) are apparently not well enough preserved to conduct geochemical analyses to estimate palaeotemperatures (Dera 2009). However, data based on palaeobotany, sedimentary indices, clay mineralogy, modelling experiments, and oxygen isotope analyses reveal a marked division of climatic zones for the longer interval of the Early Jurassic (Chandler et al. 1992; McArthur et al. 2000; Rees et al. 2000; Bailey et al. 2003; Wignall et al. 2005; Dera et al. 2009a, b, 2011a). The NWE province was characterized by a subtropical humid climate affected by seasonal mega monsoons leading to frequent northerly freshwater influxes, whereas warm and semiarid conditions prevailed in the MED province. Furthermore, the northern latitudinal position of the NWE province and the presence of connections with cold waters via the Viking Corridor to the north (e.g. Westermann 1993; Arias 2006; Dera et al. 2009a, 2015) and with warm waters of the Tethys Ocean via the Austro-Alpine and

Z A C A €I E T A L . : R A P O P O R T E F F E C T I N E A R L Y J U R A S S I C A M M O N I T E S

Jamesoni

11

Davoei

Ibex

Number of species

60

40

40

20

20

10

0

0 0

10

20

0

30

0

10

20

30

40

0

10

20

30

Latitudinal range (°) NWE species

MED species

Trans-province species

Distribution of species latitudinal range sizes of ammonites from the western Tethys and adjacent areas for the three chronozones of the early Pliensbachian. Colour online.

FIG. 3.

Pontic provinces to the south (Dommergues & Meister 1991; Dera et al. 2009a; Dommergues et al. 2009) strongly suggest the presence of a latitudinal gradient of temperatures in the NWE province. As for the MED province, palaeoclimatic simulations have shown that: (1) the warmer sea-surface temperatures prevailing during the Jurassic implied larger heat transfers to the poles (Chandler et al. 1992); and (2) the peculiar geomorphology of the western Tethys margin (i.e. a north-westward bottleneck) was prone to tropical cyclone development (Marsaglia & Klein 1983; Krencker et al. 2015), which also favoured north-western heat transport and thus a homogenization of seawater temperature in this area. This supports the hypothesis of a Rapoport effect driven by climatic variability in the NWE province; the MED province belonging to the climatically more stable and homogeneous intertropical zone during the early Pliensbachian.

Rapoport effect through the early Pliensbachian The early Pliensbachian corresponds to a climatically stable interval, except for its last chronozone during which several palaeoclimatic and palaeoenvironmental changes are recorded in the NWE province. The Davoei chronozone is indeed characterized by an increase in seasurface temperatures of about 4°C (Rosales et al. 2004a;

Dera et al. 2009a, b, 2011b; Price et al. 2016), a sea-level rise (Hallam 2001), more humid climatic conditions (Dera et al. 2009b) and a potential incursion of warm water masses from the Tethys Ocean (Dera et al. 2009a, 2015). These modifications probably obliterated the phylogenetic conservatism of ammonite range sizes (Zaca€ı et al. 2017) and could also potentially have modified their latitudinal distribution (this work). Latitudinal range sizes of NWE species do not increase with latitude during the Jamesoni chronozone but show a significant Rapoport effect during the Ibex chronozone and obviously increase with latitude during the Davoei chronozone, although the relationship is not significant (Fig. 2; Table 2). This lack of significance is probably due to the low species richness characterizing the Davoei chronozone (Dommergues et al. 2009; Hardy et al. 2012; Table 2) and to a slightly less even spatial distribution of fossil localities during this chronozone compared to the two previous ones (Zaca€ı et al. 2018, fig. S2). This decrease in species richness could be related to the increase in sea-surface temperatures previously mentioned and to a fall in origination rates (Dera et al. 2010). Hence, contrary to what could have been expected based on the analysis at the scale of the early Pliensbachian (Fig. 2C) and on the known palaeoclimatic changes, the only chronozone that does not show a Rapoport effect for the NWE province is not Davoei but Jamesoni. For the MED province, each chronozone shows a different pattern of latitudinal distribution of species range

12

PALAEONTOLOGY

size: no particular trend during Jamesoni, a boundary effect similar to the one observed at the scale of the early Pliensbachian during Ibex, and a significant decrease in range size with latitude (i.e. inverse of a Rapoport effect) during Davoei (Fig. 2). Latitudinal midpoints of MED species being exclusively located in the northern half of the province during the Davoei chronozone, this decrease in species range size with latitude could simply have resulted from a geographical bias: if no occurrences are recorded in the southernmost localities during this interval, only the northern half of the boundary effect recognized during the Ibex chronozone and at the scale of the early Pliensbachian (Fig. 2H and B, respectively) would be visible (Fig. 2E). However, the southernmost localities of the MED province do record several occurrences, meaning that this pattern does not result from a geographical bias. It is probably related to the absence of endemic species (i.e. species showing only one occurrence and thus null range sizes) in those southernmost localities. Those that occurred during Ibex went extinct, and no other endemic species replaced them. The climatic change is likely to have played a role in these extinctions and in the absence of new endemic species in these southern localities. At the higher scale of the entire western Tethys, we observe a Rapoport effect during Jamesoni, no noticeable pattern during Ibex and an inverse Rapoport effect during Davoei (Fig. 2). These patterns are hardly interpretable as such because of the existence of the biogeographical barrier, but they can be compared with those found at the smaller scale of the province and are not predictable based on the latter (Fig. 2). For instance, the distribution of species latitudinal range sizes at the scale of the western Tethys shows a Rapoport effect, whereas it does not at the scale of the MED and NWE provinces (Fig. 2J–L). We may have concluded the existence of a genuine Rapoport effect and of a latitudinal gradient of climatic variability at the scale of the western Tethys if the biogeographical structuring had not been documented by previous studies (e.g. Dommergues 1982; Dommergues & Meister 1991; Meister & Stampfli 2000; Dommergues et al. 2009). Similarly, without knowledge of the biogeographical barrier, we could have ignored the existence of a Rapoport effect at the scale of the NWE province. This suggests that the influence of climatic variability on species latitudinal distribution is confined to ‘consistent’ biogeographical areas where species dispersal is not limited by physical barriers, which can explain previous findings stating that the Rapoport effect is a local phenomenon that essentially applies to the northern hemisphere, out of the tropics (Rohde et al. 1993; Blackburn & Gaston 1996; Rohde 1996; Gaston et al. 1998; Gaston & Chown 1999). Temporal scale effects are also observed; patterns observed at the scale of the early Pliensbachian (Fig. 2A–

C) do not result from the addition of patterns observed at the scale of the chronozone, regardless the geographical area under study (Fig. 2D–L). Thus, conclusions drawn at a given temporal scale can be inconsistent at another temporal scale. This calls for a multiscale approach when studying spatial distribution of extant species range size, in order to distinguish scale effects from genuine macroecological signal.

Climatic variability hypothesis and species thermal tolerance Acknowledging the fact that climatic variability influences species range size means recognizing that the size and limits of their geographical ranges relate to their temperature tolerance. In the case of marine invertebrates, this relationship remains poorly known (Witman & Roy 2009). Because ammonites are extinct organisms, one can only rely on: (1) indirect evidence of such a relationship as an empirical tests are impossible; and (2) current knowledge of modern coleoids such as squids and octopods, to which ammonoids are more closely related than to nautiloids (e.g. Ritterbush et al. 2014 and references therein). Zaca€ı et al. (2016) showed that long-distance dispersal of early Pliensbachian ammonites of the western Tethys was not related to morphological features (i.e. shell size and shape) but rather to environmental characteristics of the province to which they belong, and thus indirectly to temperature. In a series of studies on Early Triassic ammonoid biogeography, Brayard et al. (2006, 2007, 2009) argued that temperature directly influences species dispersal and distribution, meaning that these organisms were temperature-dependent. Numerous analyses of ammonoid shell d18O show variation through ontogenetic stages, which probably means that these were exposed to different temperatures during their life (e.g. Lukeneder 2015 and references therein). Studies on the relationship between ocean temperature and cephalopod species range size are very scarce, although temperature might play an important role in the spatial distribution of these poikilothermic organisms (Lee et al. 1998; Boyle & Rodhouse 2005). In their study of Atlantic cephalopod biogeography, Rosa et al. (2008) have shown that the combined latitudinal and bathymetric range size of western neritic cephalopods increases with latitude and is strongly associated with the sea-surface temperature gradient and the extension of coral reef regions. According to Voss (1983) and Rathjen & Voss (1987), temperature is the prime factor controlling aspects of squid and octopod migration; and studies in the English Channel and the North Sea have shown that temperature and climatic variability play an important role in squid distribution and abundance (Robin & Denis 1999; Challier et al. 2005; van der Kooij et al. 2016).

Z A C A €I E T A L . : R A P O P O R T E F F E C T I N E A R L Y J U R A S S I C A M M O N I T E S

Expanding this relationship to marine invertebrates shows that a consensus is far from being reached. For Briggs (1995), even if distribution patterns in local situations may be affected by factors such as food, salinity and oxygen, it is temperature that primarily controls species distribution on an oceanic or global scale. Sunday et al. (2012) studied the relationship between species temperature tolerance and range size for more than a hundred marine ectotherms and found that: (1) their observed geographical range boundaries closely match their potential latitudinal ranges (estimated on the basis of their thermal tolerance and extreme temperatures across latitudes); and (2) this relationship may lie as much with temperature as with factors closely correlated to temperature, such as dissolved oxygen availability. The analysis of Belanger et al. (2012) indicated that the large-scale geographical distribution of shallow marine bivalves (i.e. molluscs) can be predicted with 89–100% accuracy by temperature plus salinity and productivity (the last two having lesser but significant effects). Tomasovy ch et al. (2015) showed that bivalve species inhabiting seasonally variable mid-latitude environments have broader thermal ranges than species inhabiting less seasonal tropical environments, supporting Stevens’ (1989) climatic variability hypothesis. By contrast, in a study based on shallow marine invertebrates of the northern hemisphere, Pappalardo et al. (2015) found that range limits predicted from water temperature and species temperature tolerance did not match well with real observations but rather corresponded to patterns of ocean circulation. In addition, using the marine mollusc database of Roy et al. (1998), Parmesan et al. (2005) showed the absence of correlation between seasonal variations in temperature at a given latitude and the median range of species living there, probably because of the broad ranges of tropical species associated with low seasonal variations. These contrasting results show that the relationship among ocean temperature, species temperature tolerance and species range size is not straightforward and may be taxon-specific and scale-dependent. For the early Pliensbachian ammonites of the western Tethys and adjacent areas, temperature seems to be an important controlling factor of their latitudinal distribution.

CONCLUSIONS Known to biogeographers for more than a century, belatedly popularized by Stevens (1989) and extensively studied since then, Rapoport’s rule remains one of the most controversial patterns in macroecology. Numerous studies, mainly focusing on extant terrestrial vertebrates, have cast doubts on both its ubiquity and on the reliability of Stevens (1989) climatic variability hypothesis as an explanatory mechanism. The study of early Pliensbachian ammonites

13

range sizes reveals patterns that are not as straightforward as those expected according to Stevens’ (1989) hypothesis. Early Pliensbachian ammonites have a dual pattern of distribution of latitudinal range sizes: in the northern biogeographic province (NWE), species range sizes show a Rapoport effect, while in the more tropical province (MED), they show a boundary effect and are in average larger than those of the NWE species. Considering the palaeoclimate, palaeoenvironment and the palaeogeographical configuration of the two provinces, this dual pattern can be explained by the alternative climatic variability hypothesis proposed by Gaston & Chown (1999). The latter better captures latitudinal seasonal variations and outlines the existence of a large intertropical zone with stable and homogeneous climate. Since climate was globally warmer in the Early Jurassic than in the Recent, the northern limit of this intertropical zone was certainly located at a higher latitude than today and probably coincided with the biogeographical barrier separating the NWE and MED provinces. This could explain both the larger dispersal of MED species, essentially limited by the presence of emerged lands in the northern and southern margins of the province, and the increase in species range sizes with latitude in the NWE province. This increase probably corresponds to a gradient of seasonal variations of temperatures from mid to high latitudes within this province. These results support the hypothesis that the Rapoport effect is a phenomenon that occurs mainly out of the tropics (Rohde et al. 1993; Blackburn & Gaston 1996; Rohde 1996; Gaston et al. 1998; Gaston & Chown 1999). Our multiscale approach further shows that this macroecological pattern is scale-dependent and may be labile through time. This probably explains the conflicting results of previous studies carried out at various spatiotemporal scales. Acknowledgements. We are grateful to Kenneth De Baets, Sally Thomas and to an anonymous reviewer for their constructive comments, which helped us to improve the manuscript. We thank Bruno Vrielynck (UPMC) for the reconstruction of the early Pliensbachian palaeogeographic map of the western Tethys and for the calculation of fossil locality palaeocoordinates that boosted several studies, including this work. Romain Vaucher is thanked for insightful comments regarding the relationship between western Tethys geomorphology and climate. This is a contribution by the team ‘SAMBA’ of the laboratory Biogeosciences, Dijon (France).

DATA ARCHIVING STATEMENT Data for this study (including phylogenetic hypothesis, locality distribution, dataset and GEE model results) are available in the Dryad Digital Repository: https://doi.org/10.5061/dryad.dh46139

Editor. Dieter Korn

14

PALAEONTOLOGY

REFERENCES € S I C H , F. T. 1996. Diversity analysis A B E R H A N , M. and F UR of Lower Jurassic bivalves of the Andean Basin and the Pliensbachian-Toarcian mass extinction. Lethaia, 29, 181–195. -- 2000. Mass origination versus mass extinction: the biological contribution to the Pliensbachian–Toarcian extinction event. Journal of the Geological Society, 157, 55–60. A G E R , D. V. 1973. Mesozoic brachiopoda. 431–436. In H A L L A M , A. (ed.) Atlas of palaeobiogeography. Elsevier. A R I A S , C. 2006. Northern and Southern Hemispheres ostracod palaeobiogeography during the Early Jurassic: possible migration routes. Palaeogeography, Palaeoclimatology, Palaeoecology, 233, 63–95. -2007. Pliensbachian–Toarcian ostracod biogeography in NW Europe: evidence for water mass structure evolution. Palaeogeography, Palaeoclimatology, Palaeoecology, 251, 398–421. - and W H A T L E Y , R. 2005. Palaeozoogeography of Western European Lower Jurassic (Pliensbachian and Toarcian) Ostracoda. Geobios, 38, 697–724.  QUEZA R I T A , H. T., R O D R IG U E Z , P. and V AZ D O MIN G U E Z , E. 2005. Continental and regional ranges of North American mammals: Rapoport’s rule in real and null worlds. Journal of Biogeography, 32, 961–971. ~ L O Z A , A., A R R O Y O , M. T. K., R I V E R O S , M., P E NA C A V I E R E S , L. and F A G G I , A. M. 1996. Phytogeographic relationships and regional richness patterns of the cool temperate rainforest flora of southern South America. 134–172. In L A W F O R D , R. G., A L A B A C K , P. and F U E N T E S , E. (eds). High-latitude rainforests and associated ecosystems of the West Coast of the Americas. Springer, New York. B A I L E Y , T. R., R O S E N T H A L , Y., M c A R T H U R , J. M., S C H O O T B R U G G E , B. VAN DE and T H I R L W A L L , M. F. 2003. Paleoceanographic changes of the Late Pliensbachian– Early Toarcian interval: a possible link to the genesis of an Oceanic Anoxic Event. Earth & Planetary Science Letters, 212, 307–320. B E C K , J., L I E D T K E , H. C., W I D L E R , S., A L T E R M A T T , F., L O A D E R , S. P., H A G M A N N , R., L A N G , S. and F I E D L E R , K. 2016. Patterns or mechanisms? Bergmann’s and Rapoport’s rule in moths along an elevational gradient. Community Ecology, 17, 137–148. B E L A N G E R , C. L., J A B L O N S K I , D., R O Y , K., B E R K E , S. K., K R U G , A. Z. and V A L E N T I N E , J. W. 2012. Global environmental predictors of benthic marine biogeographic structure. Proceedings of the National Academy of Sciences, 109, 14046–14051. B J E R R U M , C. J., S U R L Y K , F., C A L L O M O N , J. H. and S L I N G E R L A N D , R. L. 2001. Numerical paleoceanographic study of the Early Jurassic transcontinental Laurasian Seaway. Paleoceanography, 16, 390–404. B L A C K B U R N , T. M. and G A S T O N , K. J. 1996. Spatial patterns in the geographic range sizes of bird species in the New World. Philosophical Transactions of the Royal Society B, 351, 897–912. -2003. Macroecology: Concepts and consequences. 43rd Symposium of the British Ecological Society. Cambridge University Press.

€ M , M., K E M P , R., W I L L I A M S , R., D A V I D S O N , A. B OH D., G A R C I A , A., M c M I L L A N , K. M., B R A M H A L L , H. R. and C O L L E N , B. 2017. Rapoport’s rule and determinants of species range size in snakes. Diversity & Distributions, 23, 1472–1481. B O Y L E , P. R. and R O D H O U S E , P. 2005. Cephalopods: ecology and fisheries. Blackwell Science. B R A T T S T R O M , B. H. 1968. Thermal acclimation in Anuran amphibians as a function of latitude and altitude. Comparative Biochemistry & Physiology, 24, 93–111. B R A Y A R D , A., B U C H E R , H., E S C A R G U E L , G., F L U T E A U , F., B O U R Q U I N , S. and G A L F E T T I , T. 2006. The Early Triassic ammonoid recovery: paleoclimatic significance of diversity gradients. Palaeogeography, Palaeoclimatology, Palaeoecology, 239, 374–395. - E S C A R G U E L , G. and B U C H E R , H. 2007. The biogeography of Early Triassic ammonoid faunas: clusters, gradients, and networks. Geobios, 40, 749–765. € W I L E R T. 2009. Smithian and --- and B R UH Spathian (Early Triassic) ammonoid assemblages from terranes: paleoceanographic and paleogeographic implications. Journal of Asian Earth Sciences, 36, 420–433.  T , E., P E L L E N A R D , P., V I N B R I G A U D , B., P U C EA C E N T , B. and J O A C H I M S K I , M. M. 2008. Climatic fluctuations and seasonality during the Late Jurassic (Oxfordian– Early Kimmeridgian) inferred from d18O of Paris Basin oyster shells. Earth & Planetary Science Letters, 273, 58–67. B R I G G S , J. C. 1995. Global biogeography. Developments in Palaeontology and Stratigraphy, 14. Elsevier. B R O W N , J. H. 1995. Macroecology. University of Chicago Press. C A R D I L L O , M. 2002. The life-history basis of latitudinal diversity gradients: how do species traits vary from the poles to the equator? Journal of Animal Ecology, 71, 79–87. C H A L L I E R , L., R O Y E R , J., P I E R C E , G. J., B A I L E Y , N., R O E L , B. and R O B I N , J.-P. 2005. Environmental and stock effects on recruitment variability in the English Channel squid Loligo forbesi. Aquatic Living Resources, 18, 353– 360. C H A N D L E R , M. A., R I N D , D. and R U E D Y , R. 1992. Pangaean climate during the Early Jurassic: GCM simulations and the sedimentary record of paleoclimate. Geological Society of America Bulletin, 104, 543–559. C H O W N , S. L., S I N C L A I R , B. J., L E I N A A S , H. P. and G A S T O N , K. J. 2004. Hemispheric asymmetries in biodiversity—a serious matter for ecology. PLoS Biology, 2, e406. C L A R K E , A. 2009. Temperature and marine macroecology. 250–278. In W I T M A N , J. D. and R O Y , K. (eds). Marine macroecology. University of Chicago Press. C O L W E L L , R. K. 2011. Biogeographical gradient theory. 391– 421. In S C H E I N E R , S. M. and W I L L I G , M. R. (eds). The theory of ecology. University of Chicago Press. -and H U R T T , G. C. 1994. Nonbiological gradients in species richness and a spurious Rapoport effect. The American Naturalist, 144, 570–595. -and L E E S , D. C. 2000. The mid-domain effect: geometric constraints on the geography of species richness. Trends in Ecology & Evolution, 15, 70–76.

Z A C A €I E T A L . : R A P O P O R T E F F E C T I N E A R L Y J U R A S S I C A M M O N I T E S

 , C. L., G I L M A N , A. C. and G., C A R D E L US L O N G I N O , J. T. 2008. Global warming, elevational range shifts, and lowland biotic attrition in the wet tropics. Science, 322, 258–261. C O N N O L L Y , S. R. 2009. Macroecological theory and the analysis of species richness gradients. 279–309. In W I T M A N , J. D. and R O Y , K. (eds). Marine macroecology. University of Chicago Press. -B E L L W O O D , D. R. and H U G H E S , T. P. 2003. IndoPacific biodiversity of coral reefs: deviations from a middomain model. Ecology, 84, 2178–2190. C O W L I S H A W , G. and H A C K E R , J. E. 1997. Distribution, diversity, and latitude in African primates. The American Naturalist, 150, 505–512. D E R A , G. 2009. Le r^ ole des changements paleoclimatiques sur l’evolution de la biodiversite au Pliensbachien et au Toarcien. PhD thesis, Universite de Bourgogne, Dijon, 352 pp.  T , E., P E L L E N A R D , P., N E I G E , P., D E L - P U C EA S A T E , D., J O A C H I M S K I , M. M., R E I S B E R G , L. and M A R T I N E Z , M. 2009a. Water mass exchange and variations in seawater temperature in the NW Tethys during the Early Jurassic: evidence from neodymium and oxygen isotopes of fish teeth and belemnites. Earth & Planetary Science Letters, 286, 198–207. - P E L L E N A R D , P., N E I G E , P., D E C O N I N C K , J.-F.,  T , E. and D O M M E R G U E S , J.-L. 2009b. DistriP U C EA bution of clay minerals in Early Jurassic Peritethyan seas: palaeoclimatic significance inferred from multiproxy comparisons. Palaeogeography, Palaeoclimatology, Palaeoecology, 271, 39–51. -N E I G E , P., D O M M E R G U E S , J.-L., F A R A , E., L A F F O N T , R. and P E L L E N A R D , P. 2010. High-resolution dynamics of Early Jurassic marine extinctions: the case of Pliensbachian-Toarcian ammonites (Cephalopoda). Journal of the Geological Society, 167, 21–33. - and B R A Y A R D A. 2011a. Ammonite paleobiogeography during the Pliensbachian–Toarcian crisis (Early Jurassic) reflecting paleoclimate, eustasy, and extinctions. Global & Planetary Change, 78, 92–105. - B R I G A U D , B., M O N N A , F., L A F F O N T , R.,  T , E., D E C O N I N C K , J.-F., P E L L E N A R D , P., P U C EA J O A C H I M S K I , M. M. and D U R L E T , C. 2011b. Climatic ups and downs in a disturbed Jurassic world. Geology, 39, 215–218. - P R U N I E R , J., S M I T H , P. L., H A G G A R T , J. W., P O P O V , E., G U Z H O V , A., R O G O V , M., D E L S A T E ,  T , E., C H A R B O N D., T H I E S , D., C U N Y , G., P U C EA N I E R , G. and B A Y O N , G. 2015. Nd isotope constraints on ocean circulation, paleoclimate, and continental drainage during the Jurassic breakup of Pangea. Gondwana Research, 27, 1599–1615. D E R C O U R T , J., Z O N E N S H A I N , L. P., R I C O U , L. E., K A Z M I N , V. G., P I C H O N , X. LE, K N I P P E R , A. L., G R A N D J A C Q U E T , C., S B O R S H C H I K O V , I. M., B O U L I N , J. and S O R O K H T I N , O. 1985. Presentation de 9 cartes paleogeographiques au 1/20 000 000e s’etendant de l’Atlantique au Pamir pour la periode du Lias a l’Actuel. Bulletin de la Societe geologique de France, 1, 637–652. -B R E H M ,

15

D O B Z H A N S K Y , T. 1950. Evolution in the tropics. American Scientist, 38, 209–221. D O M M E R G U E S , J.-L. 1982. Le provincialisme des ammonites nord-ouest europeennes au Lias moyen. Une crise faunique sous contr^ ole paleogeographique. Bulletin de la Societe geologique de France, 24, 1047–1051. -and E L H A R I R I , K. 2002. Endemism as a palaeobiogeographic parameter of basin history illustrated by early- and mid-Liassic peri-Tethyan ammonite faunas. Palaeogeography, Palaeoclimatology, Palaeoecology, 184, 407–418. - and M E I S T E R , C. 1991. Area of mixed marine faunas between two major paleogeographical realms, exemplified by the Early Jurassic (Late Sinemurian and Pliensbachian) ammonites in the Alps. Palaeogeography, Palaeoclimatology, Palaeoecology, 86, 265–282. -L A U R I N , B. and M E I S T E R , C. 2001. The recovery and radiation of Early Jurassic ammonoids: morphologic versus palaeobiogeographical patterns. Palaeogeography, Palaeoclimatology, Palaeoecology, 165, 195–213. -F A R A , E. and M E I S T E R , C. 2009. Ammonite diversity and its palaeobiogeographical structure during the early Pliensbachian (Jurassic) in the western Tethys and adjacent areas. Palaeogeography, Palaeoclimatology, Palaeoecology, 280, 64–77. D O N O V A N , D. T. 1967. The geographical distribution of Lower Jurassic ammonites in Europe and adjacent areas. Aspects of Tethyan Biogeography, 7, 111–134. D O Y L E , P. 1987. Lower Jurassic-Lower Cretaceous belemnite biogeography and the development of the Mesozoic Boreal Realm. Palaeogeography, Palaeoclimatology, Palaeoecology, 61, 237–254. -1994. Aspects of the distribution of Early Jurassic belemnites. Palaeopelagos Special Publication, 1, 109–120. E D W A R D S , W. and W E S T O B Y , M. 1996. Reserve mass and dispersal investment in relation to geographic range of plant species: phylogenetically independent contrasts. Journal of Biogeography, 23, 329–338. F I S C H E R , A. G. 1960. Latitudinal variations in organic diversity. Evolution, 14, 64. ~ , R. S. 2004. The applicability F O R T E S , R. R. and A B S A L AO of Rapoport’s rule to the marine molluscs of the Americas. Journal of Biogeography, 31, 1909–1916. F R A N C E , R. 1992. The North American latitudinal gradient in species richness and geographical range of freshwater crayfish and amphipods. The American Naturalist, 139, 342–354. G A S T O N , K. J. 1991. The magnitude of global insect species richness. Conservation Biology, 5, 283–296. -1996. Species-range-size distributions: patterns, mechanisms and implications. Trends in Ecology & Evolution, 11, 197–201. - 2003. The structure and dynamics of geographic ranges. Oxford University Press. -and C H O W N , S. L. 1999. Why Rapoport’s rule does not generalise. Oikos, 84, 309. -and G A U L D , I. D. 1993. How many species of pimplines (Hymenoptera: Ichneumonidae) are there in Costa Rica? Journal of Tropical Ecology, 9, 491–499. -B L A C K B U R N , T. M. and S P I C E R , J. I. 1998. Rapoport’s rule: time for an epitaph? Trends in Ecology & Evolution, 13, 70–74.

16

PALAEONTOLOGY

G E C Z Y , B. 1984. Jurassic ammonite provinces of Europe. Acta Geologica Hungarica, 27, 67–71. G R A D S T E I N , F. M., S C H M I T Z , M. D., O G G , G. M. and O G G , J. G. 2012. The geologic time scale. Elsevier. € I G , M., B E E R L I , N., B A L L E S T E R O S - M E J I A , L., G R UN K I T C H I N G , I. J. and B E C K , J. 2017. How climatic variability is linked to the spatial distribution of range sizes: seasonality versus climate change velocity in sphingid moths. Journal of Biogeography, 44, 2441–2450.  A N , J.-F. 2009. May Rapoport’s G U E R N I E R , V. and G U EG rule apply to human associated pathogens? EcoHealth, 6, 509–521. H A L L A M , A. 1969. Faunal realms and facies in the Jurassic. Palaeontology, 12, 1–18. - 1972. Diversity and density characteristics of Pliensbachian-Toarcian molluscan and brachiopod faunas of the North Atlantic margins. Lethaia, 5, 389–412. -1977. Jurassic bivalve biogeography. Paleobiology, 3, 58–73. - 2001. A review of the broad pattern of Jurassic sea-level changes and their possible causes in the light of current knowledge. Palaeogeography, Palaeoclimatology, Palaeoecology, 167, 23–37. H A R D Y , C., F A R A , E., L A F F O N T , R., D O M M E R G U E S , J.-L., M E I S T E R , C. and N E I G E , P. 2012. Deep-time phylogenetic clustering of extinctions in an evolutionarily dynamic clade (Early Jurassic ammonites). PLoS One, 7, e37977. H E L M U T H , B., H A R L E Y , C. D. G., H A L P I N , P. M., O ’ D O N N E L L , M., H O F M A N N , G. E. and B L A N C H E T T E , C. A. 2002. Climate change and latitudinal patterns of intertidal thermal stress. Science, 298, 1015–1017.  D E Z F E R N AN  D E Z , M. and V R B A , E. S. 2005. H E R N AN Rapoport effect and biomic specialization in African mammals: revisiting the climatic variability hypothesis. Journal of Biogeography, 32, 903–918. H U B B E L L , S. P. 2001. The unified neutral theory of biodiversity and biogeography. Monographs in Population Biology, 32. Princeton University Press. H U G H E S , L., C A W S E Y , E. M. and W E S T O B Y , M. 1996. Geographic and climatic range sizes of Australian eucalypts and a test of Rapoport’s rule. Global Ecology & Biogeography Letters, 5, 128. HUTCHINS, L. W. 1947. The bases for temperature zonation in geographical distribution. Ecological Monographs, 17, 325–335. J A B L O N S K I , D., B E L A N G E R , C. L., B E R K E , S. K.,  H, H U A N G , S., K R U G , A. Z., R O Y , K., T O M A SO V YC A. and V A L E N T I N E , J. W. 2013. Out of the tropics, but how? Fossils, bridge species, and thermal ranges in the dynamics of the marine latitudinal diversity gradient. Proceedings of the National Academy of Sciences, 110, 10487–10494. J A N Z E N , D. H. 1967. Why mountain passes are higher in the tropics. The American Naturalist, 101, 233–249. J O H N S O N , C. N. 1998. Rarity in the tropics: latitudinal gradients in distribution and abundance in Australian mammals. Journal of Animal Ecology, 67, 689–698. J U L I A N O , S. A. 1983. Body size, dispersal ability, and range size in North American species of Brachinus (Coleoptera: Carabidae). The Coleopterists’ Bulletin, 37, 232–238. € N I N G - G A E S E , K., P R I N Z I N G E R , R., K H A L I Q , I., B OH P F E N N I N G E R , M. and H O F , C. 2017. The influence of

thermal tolerances on geographical ranges of endotherms. Global Ecology & Biogeography, 26, 650–668. K O O I J , J. VA N D E R, E N G E L H A R D , G. H. and RIGHTON, D. A. 2016. Climate change and squid range expansion in the North Sea. Journal of Biogeography, 43, 2285–2298. K R E N C K E R , F.-N., B O D I N , S., S U A N , G., H E I M H O F E R , U., K A B I R I , L. and I M M E N H A U S E R , A. 2015. Toarcian extreme warmth led to tropical cyclone intensification. Earth & Planetary Science Letters, 425, 120–130. L A N E , C. S. 2007. Latitudinal range variation of trees in the United States: a reanalysis of the applicability of Rapoport’s rule. The Professional Geographer, 59, 115–130. L E E , P. G., T U R K , P. E., F O R S Y T H E , J. W. and D I M A R C O , F. P. 1998. Cephalopod culture: physiological, behavioral and environmental requirements. Aquaculture Science, 46, 417–422. L E T C H E R , A. J. and H A R V E Y , P. H. 1994. Variation in geographical range size among mammals of the Palearctic. The American Naturalist, 144, 30–42. L I A N G , K.-Y. and Z E G E R , S. L. 1986. Longitudinal data analysis using generalized linear models. Biometrika, 73, 13–22. € S I C H , F. T. 1998. Bivalve L I U , C., H E I N Z E , M. and F UR provinces in the Proto-Atlantic and along the southern margin of the Tethys in the Jurassic. Palaeogeography, Palaeoclimatology, Palaeoecology, 137, 127–151. L O R D , A. R. 1988. Ostracoda of the Early Jurassic Tethyan ocean. 855–868. In H A N A I , T., I K E Y A , N. and I S H I Z A K I , K. (eds). Evolutionary biology of Ostracoda. Developments in Palaeontology & Stratigraphy, 11. Elsevier. L U K E N E D E R , A. 2015. Ammonoid habitats and life history. 689–791. In K L U G , C., K O R N , D., D E B A E T S , K., K R U T A , I. and M A P E S , H. R. (eds). Ammonoid paleobiology: From anatomy to ecology. Springer. L Y O N S , S. K. and S M I T H , F. A. 2010. Using macroecological approach to study geographic range, abundance and body size in the fossil record. The Paleontological Society Papers, 16 (Quantitative Methods in Paleobiology), 117–141. -and W I L L I G , M. R. 1997. Latitudinal patterns of range size: methodological concerns and empirical evaluations for New World bats and marsupials. Oikos, 79, 568–580. M A C P H E R S O N , E. 2003. Species range size distributions for some marine taxa in the Atlantic Ocean. Effect of latitude and depth. Biological Journal of the Linnean Society, 80, 437–455. -and D U A R T E , C. M. 1994. Patterns in species richness, size, and latitudinal range of East Atlantic fishes. Ecography, 17, 242–248. M A R S A G L I A , K. M. and K L E I N , G. D. 1983. The paleogeography of Paleozoic and Mesozoic storm depositional systems. The Journal of Geology, 91, 117–142. M c A R T H U R , J. M., D O N O V A N , D. T., T H I R L W A L L , M. F., F O U K E , B. W. and M A T T E Y , D. 2000. Strontium isotope profile of the early Toarcian (Jurassic) oceanic anoxic event, the duration of ammonite biozones, and belemnite palaeotemperatures. Earth & Planetary Science Letters, 179, 269–285. -A L G E O , T. J., S C H O O T B R U G G E , B. VAN DE, L I , Q. and H O W A R T H , R. J. 2008. Basinal restriction, black shales, Re-Os dating, and the Early Toarcian (Jurassic) oceanic anoxic event. Paleoceanography & Paleoclimatology, 23, PA4217.

Z A C A €I E T A L . : R A P O P O R T E F F E C T I N E A R L Y J U R A S S I C A M M O N I T E S

M E I S T E R , C. and S T A M P F L I , G. 2000. Les ammonites du Lias moyen (Pliensbachien) de la Neotethys et de ses confins; compositions fauniques, affinites paleogeographiques et biodiversite. Revue de Paleobiologie, 19, 227–292. M O R I N , X. and C H U I N E , I. 2006. Niche breadth, competitive strength and range size of tree species: a trade-off based framework to understand species distribution: Species range patterns and processes. Ecology Letters, 9, 185–195. - and L E C H O W I C Z , M. J. 2011. Geographical and ecological patterns of range size in North American trees. Ecography, 34, 738–750. M O R U E T A - H O L M E , N., E N Q U I S T , B. J., M c G I L L , B. J., B O Y L E , B., J Ø R G E N S E N , P. M., O T T , J. E., P E E T ,  I., S L O A T , L. L., T H I E R S , B., V I O L L E , R. K., SIM O V A, C., W I S E R , S. K., D O L I N S , S., D O N O G H U E , J. C., K R A F T , N. J. B., R E G E T Z , J., S C H I L D H A U E R , M., S P E N C E R , N. and S V E N N I N G , J.-C. 2013. Habitat area and climate stability determine geographical variation in plant species range sizes. Ecology Letters, 16, 1446–1454. M O U R E L L E , C. and E Z C U R R A , E. 1997. Rapoport’s rule: a comparative analysis between South and North American columnar cacti. The American Naturalist, 150, 131–142. O G G , J. G., O G G , G. M. and G R A D S T E I N , F. M. 2016. A concise geologic time scale: 2016. Elsevier. O R M E , C. D. L., D A V I E S , R. G., O L S O N , V. A., T H O M A S , G. H., D I N G , T.-S., R A S M U S S E N , P. C., R I D G E L Y , R. S., S T A T T E R S F I E L D , A. J., B E N N E T T , P. M. and O W E N S , I. P. 2006. Global patterns of geographic range size in birds. PLoS Biology, 4, e208. P A G E , K. N. 2003. The Lower Jurassic of Europe: its subdivision and correlation. Geological Survey of Denmark & Greenland Bulletin, 1, 23–59. P A G E L , M. D., M A Y , R. M. and C O L L I E , A. R. 1991. Ecological aspects of the geographical distribution and diversity of mammalian species. The American Naturalist, 137, 791–815. P A P P A L A R D O , P., P R I N G L E , J. M., W A R E S , J. P. and B Y E R S , J. E. 2015. The location, strength, and mechanisms behind marine biogeographic boundaries of the east coast of North America. Ecography, 38, 722–731. P A R A D I S , E. and C L A U D E , J. 2002. Analysis of comparative data using generalized estimating equations. Journal of Theoretical Biology, 218, 175–185. -and S T R I M M E R K. 2004. APE: analyses of phylogenetics and evolution in R language. Bioinformatics, 20, 289–290. P A R M E S A N , C., G A I N E S , S., G O N Z A L E Z , L., K A U F M A N , D. M., K I N G S O L V E R , J., T O W N S E N D P E T E R S O N , A. and S A G A R I N , R. 2005. Empirical perspectives on species borders: from traditional biogeography to global change. Oikos, 108, 58–75. P I N T O R , A. F. V., S C H W A R Z K O P F , L. and K R O C K E N B E R G E R , A. K. 2015. Rapoport’s Rule: do climatic variability gradients shape range extent? Ecological Monographs, 85, 643–659. P R I C E , T. D., H E L B I G , A. J. and R I C H M A N , A. D. 1997. Evolution of breeding distributions in the Old World leaf warblers (genus Phylloscopus). Evolution, 51, 552–561.  P R I C E , G. D., B A K E R , S. J., V A N D E V E L D E , J. and C L EM E N C E , M.-E. 2016. High-resolution carbon cycle and

17

seawater temperature evolution during the Early Jurassic (Sinemurian-Early Pliensbachian). Geochemistry, Geophysics, Geosystems, 17, 3917–3928. R C O R E T E A M (2016). R: A language and environment for statistical computing. Version 3.3.1. R Foundation for Statistical Computing, Vienna, Austria. https://www.R-project.org/ R A P O P O R T , E. H. 1982. Areography: Geographical strategies of species, Vol. 1. Pergamon Press, Oxford. R A T H J E N , W. F. and V O S S , G. L. 1987. The cephalopod fisheries: a review. Cephalopod Life Cycles, 2, 253–275. R E E S , P. M., Z I E G L E R , A. M. and V A L D E S , P. J. 2000. Jurassic phytogeography and climates: new data and model comparisons. 297–318. In H U B E R , B. T., M A C L E O D , K. G. and W I N G , S. L. (eds). Warm climates in earth history. Cambridge University Press. RICHARDS, O. W. 1978. Social wasps of the Americas excluding the Vespinae. British Museum (Natural History), London. R I T T E R B U S H , K. A., H O F F M A N N , R., L U K E D E N E R , A. and D E B A E T S , K. 2014. Pelagic palaeoecology: the importance of recent constraints on ammonoid palaeobiology and life history. Journal of Zoology, 292, 229–241. R O B I N , J.-P. and D E N I S , V. 1999. Squid stock fluctuations and water temperature: temporal analysis of English Channel Loliginidae. Journal of Applied Ecology, 36, 101–110. R O H D E , K. 1996. Rapoport’s rule is a local phenomenon and cannot explain latitudinal gradients in species diversity. Biodiversity Letters, 3, 10–13. -1999. Latitudinal gradients in species diversity and Rapoport’s rule revisited: a review of recent work and what can parasites teach us about the causes of the gradients? Ecography, 22, 593–613. -and H E A P , M. 1996. Latitudinal ranges of teleost fish in the Atlantic and Indo-Pacific oceans. The American Naturalist, 147, 659–665. -and H E A P D. 1993. Rapoport’s rule does not apply to marine teleosts and cannot explain latitudinal gradients in species richness. The American Naturalist, 142, 1–16. R O S A , R., D I E R S S E N , H. M., G O N Z A L E Z , L. and S E I B E L , B. A. 2008. Ecological biogeography of cephalopod molluscs in the Atlantic Ocean: historical and contemporary causes of coastal diversity patterns. Global Ecology & Biogeography, 17, 600–610. R O S A L E S , I., Q U E S A D A , S. and R O B L E S , S. 2004a. Paleotemperature variations of Early Jurassic seawater recorded in geochemical trends of belemnites from the Basque–Cantabrian basin, northern Spain. Palaeogeography, Palaeoclimatology, Palaeoecology, 203, 253–275. -R O B L E S , S. and Q U E S A D A , S. 2004b. Elemental and oxygen isotope composition of Early Jurassic belemnites: salinity vs. temperature signals. Journal of Sedimentary Research, 74, 342–354. R O Y , K., J A B L O N S K I , D. and V A L E N T I N E , J. W. 1994. Eastern Pacific molluscan provinces and latitudinal diversity gradient: no evidence for ‘Rapoport’s rule’. Proceedings of the National Academy of Sciences, 91, 8871–8874. -and R O S E N B E R G G. 1998. Marine latitudinal diversity gradients: tests of causal hypotheses. Proceedings of the National Academy of Sciences, 95, 3699–3702.

18

PALAEONTOLOGY

R U G G I E R O , A. 1994. Latitudinal correlates of the sizes of mammalian geographical ranges in South America. Journal of Biogeography, 21, 545–559. S A E L E N , G., D O Y L E , P. and T A L B O T , M. R. 1996. Stableisotope analyses of belemnite rostra from the Whitby Mudstone Formation, England; surface water conditions during deposition of a marine black shale. Palaios, 11, 97–117. S A X , D. F. 2001. Latitudinal gradients and geographic ranges of exotic species: implications for biogeography. Journal of Biogeography, 28, 139–150. S C H O O T B R U G G E , B. VAN DE, B A I L E Y , T. R., R O S E N T H A L , Y., K A T Z , M. E., W R I G H T , J. D., M I L L E R , K. G., F E I S T - B U R K H A R D T , S. and F A L K O W S K I , P. G. 2005. Early Jurassic climate change and the radiation of organic-walled phytoplankton in the Tethys Ocean. Paleobiology, 31, 73–97. S I M P S O N , M. 1855. The fossils of the Yorkshire Lias, described from nature, with a short outline of the Yorkshire coast. Whittaker & Co., London & Whitby, 149 pp. S M I T H , P. L. and T I P P E R , H. W. 1986. Plate tectonics and paleobiogeography; Early Jurassic (Pliensbachian) endemism and diversity. Palaios, 1, 399–412. S M I T H , F. D. M., M A Y , R. M. and H A R V E Y , P. H. 1994. Geographical ranges of Australian mammals. The Journal of Animal Ecology, 63, 441–450. S T E V E N S , G. C. 1989. The latitudinal gradient in geographical range: how so many species coexist in the tropics. The American Naturalist, 133, 240–256. -1992. Spilling over the competitive limits to species coexistence. 40–58. In E L D R E D G E , N. (ed.) Systematics, ecology, and the biodiversity crisis. Columbia University Press. - 1996. Extending Rapoport’s rule to Pacific marine fishes. Journal of Biogeography, 23, 149–154. -and E N Q U I S T , B. J. 1998. Macroecological limits to the abundance and distribution of Pinus. 183–190. In R I C H A R D S O N , D. M. (ed.) Ecology and biogeography of Pinus. Cambridge University Press. S T U A R T - S M I T H , R. D., E D G A R , G. J. and B A T E S , A. E. 2017. Thermal limits to the geographic distributions of shallow-water marine species. Nature Ecology & Evolution, 1, 1846–1852. S U N D A Y , J. M., B A T E S , A. E. and D U L V Y , N. K. 2012. Thermal tolerance and the global redistribution of animals. Nature Climate Change, 2, 686–690. T A Y L O R , C. M. and G O T E L L I , N. J. 1994. The macroecology of Cyprinella: correlates of phylogeny, body size, and geographical range. The American Naturalist, 144, 549–569. T E R B O R G H , J. 1973. On the notion of favorableness in plant ecology. The American Naturalist, 107, 481–501.  H , A. and J A B L O N S K I , D. 2017. Decoupling T O M A SO V YC of latitudinal gradients in species and genus geographic range size: a signature of clade range expansion. Global Ecology & Biogeography, 26, 288–303. -B E R K E , S. K., K R U G , A. Z. and V A L E N T I N E , J. W. 2015. Nonlinear thermal gradients shape broad-scale patterns in geographic range size and can reverse Rapoport’s rule: latitudinal gradients in geographic range size. Global Ecology & Biogeography, 24, 157–167.

-K E N N E D Y ,

J. D., B E T Z N E R , T. J., K U E H N L E , N. B., E D I E , S., K I M , S., S U P R I Y A , K., W H I T E , A. E., R A H B E K , C., H U A N G , S., P R I C E , T. D. and J A B L O N S K I , D. 2016. Unifying latitudinal gradients in range size and richness across marine and terrestrial systems. Proceedings of the Royal Society B, 283, 20153027. V E T E R , N. M., D E S A N T I S , L. R. G., Y A N N , L. T., D O N O H U E , S. L., H A U P T , R. J., C O R A P I , S. E., F A T H E L , S. L., G O O T E E , E. K., L O F F R E D O , L. F., R O M E R , J. L. and V E L K O V S K Y , S. M. 2013. Is Rapoport’s rule a recent phenomenon? A deep time perspective on potential causal mechanisms. Biology Letters, 9, 20130398. € OS, € A. 1977. Provinciality of the Mediterranean Lower JurasV OR sic brachiopod fauna: causes and plate-tectonic implications. Palaeogeography, Palaeoclimatology, Palaeoecology, 21, 1–16. V O S S , G. L. 1983. A review of cephalopod fisheries biology. Memoirs of the National Museum of Victoria, 44, 229–241. W A L L A C E , A. R. 1878. Tropical nature and other essays. Macmillan & Company, London. W E I S E R , M. D., E N Q U I S T , B. J., B O Y L E , B., K I L L E E N , T. J., J Ø R G E N S E N , P. M., F O N S E C A , G., J E N N I N G S , M. D., K E R K H O F F , A. J., L A C H E R J R , T. E., M O N T E A G U D O , A., V A R G A S , M. P. N., P H I L L I P S , O. L., S W E N S O N , N. G. and M A R T IN E Z , R. V. 2007. Latitudinal patterns of range size and species richness of New World woody plants. Global Ecology & Biogeography, 16, 679–688. W E S T E R M A N N , G. E. G. 1993. Global bio-events in midJurassic ammonites controlled by seaways. 187–226. In H O U S E , M. R. (ed.) The Ammonoidea: Environment, ecology and evolutionary change. Systematics Association Special Volume, 47. Oxford University Press. W I G N A L L , P. B., N E W T O N , R. J. and L I T T L E , C. T. 2005. The timing of paleoenvironmental change and causeand-effect relationships during the Early Jurassic mass extinction in Europe. American Journal of Science, 305, 1014–1032. W I L L I G , M. R., K A U F M A N , D. M. and S T E V E N S , R. D. 2003. Latitudinal gradients of biodiversity: pattern, process, scale, and synthesis. Annual Review of Ecology, Evolution, & Systematics, 34, 273–309. W I T M A N , J. D. and R O Y , K. 2009. Marine macroecology. University of Chicago Press. Z A C A €I, A., B R A Y A R D , A., D O M M E R G U E S , J.-L., M E I S T E R , C., E S C A R G U E L , G., L A F F O N T , R., V R I E L Y N C K , B. and F A R A , E. 2016. Gauging scale effects and biogeographical signals in similarity distance decay analyses: an Early Jurassic ammonite case study. Palaeontology, 59, 671–687. - F A R A , E., B R A Y A R D , A., L A F F O N T , R., D O M M E R G U E S , J.-L. and M E I S T E R , C. 2017. Phylogenetic conservatism of species range size is the combined outcome of phylogeny and environmental stability. Journal of Biogeography, 44, 2451–2462. - B R A Y A R D , A., L A F F O N T , R., D O M M E R G U E S , J.-L., M E I S T E R , C. and F A R A , E. 2018. Data from: The Rapoport effect and the climatic variability hypothesis in Early Jurassic ammonites. Dryad Digital Repository. https://doi.org/ 10.5061/dryad.dh46139