
==== Front
Ecol Evol
Ecol Evol
10.1002/(ISSN)2045-7758
ECE3
Ecology and Evolution
2045-7758
John Wiley and Sons Inc. Hoboken

10.1002/ece3.70133
ECE370133
ECE-2023-10-01836.R2
Biogeography
Macroecology
Paleoecology
Research Article
Research Article
Do modern climatic niches distinguish extinct and extant plant genera in New Zealand?
Schlenker et al.
Schlenker Nora https://orcid.org/0000-0002-3693-5946
1 nschlenker@wisc.edu

Lee William G. https://orcid.org/0000-0001-7717-0807
2
Reichgelt Tammo https://orcid.org/0000-0001-8652-5489
3
Ohlemüller Ralf https://orcid.org/0000-0001-9102-6481
4
1 Department of Geography University of Wisconsin‐Madison Madison Wisconsin USA
2 Manaaki Whenua – Landcare Research Dunedin New Zealand
3 Department of Earth Sciences University of Connecticut Storrs Connecticut USA
4 School of Geography University of Otago Dunedin New Zealand
* Correspondence
Nora Schlenker, Department of Geography, University of Wisconsin‐Madison, Madison, WI, USA.
Email: nschlenker@wisc.edu

03 9 2024
9 2024
14 9 10.1002/ece3.v14.9 e7013325 6 2024
19 10 2023
24 7 2024
© 2024 The Author(s). Ecology and Evolution published by John Wiley & Sons Ltd.
https://creativecommons.org/licenses/by/4.0/ This is an open access article under the terms of the http://creativecommons.org/licenses/by/4.0/ License, which permits use, distribution and reproduction in any medium, provided the original work is properly cited.

Abstract

Past climate changes have had large impacts on modern ecological patterns. Understanding if legacies are distinguishable in the climatic niches of extant and locally extinct taxa can provide insight into the importance of climate in extinction events. To better understand mid‐ to late‐Cenozoic New Zealand plant extinctions, which are often attributed to Cenozoic climate cooling, we identify 13 con‐familial extinct and extant New Zealand genus pairs, which have modern distributions in Australia. Using climatic niches derived from current geographic distributions in Australia, we compared (i) total niche breadth, (ii) niche overlap, and (iii) individual climate parameters, to investigate potential climate drivers of intrafamilial extinction and persistence patterns in New Zealand. A majority of New Zealand extinct genera (9 out of 13 pairs) do not indicate climate niche legacies consistent with susceptibility to extinction from changing climates, while the remaining four extinct/extant pairs show slight climatic niche legacies. Three extinct genera have warmer niches than their extant counterpart, which is consistent with extinction reflecting intolerance of cooling Cenozoic climates. The other genus pair with a climatic niche legacy has an extinct genus that is distinguished by a niche with smaller precipitation seasonality than its extant counterpart, suggesting that climate metrics other than temperature may also be important extinction drivers in some taxa. Our results show that the mechanisms of Cenozoic extinctions of New Zealand genera are likely more complex than taxa reaching environmental tolerances due to cooling climates. Comparisons of current climatic niches between extant and extinct sister taxa can provide useful insights into large‐scale, long‐term climatic legacies but more analyses, including trait and phylogeographic patterns, would lead to additional insights into alternative pathways of extinction.

We aim to better understand New Zealand plant extinctions, often attributed to Cenozoic climate cooling, through determining if legacies of past extinction events are detectable in the climatic niches of con‐familial genera in the modern Australian floral assemblage. We find that the mechanisms of Cenozoic extinctions of New Zealand genera may be more complex than taxa reaching environmental tolerances from cooling climates. Comparisons of current climate niches between extant and extinct sister taxa can provide useful insights into large‐scale, long‐term climatic legacies.

Australia
Cenozoic
climate
cooling
extinction
flora
New Zealand
niche
Schoener's D
source-schema-version-number2.0
cover-dateSeptember 2024
details-of-publishers-convertorConverter:WILEY_ML3GV2_TO_JATSPMC version:6.4.8 mode:remove_FC converted:04.09.2024
Schlenker, N. , Lee, W. G. , Reichgelt, T. , & Ohlemüller, R. (2024). Do modern climatic niches distinguish extinct and extant plant genera in New Zealand? Ecology and Evolution, 14 , e70133. 10.1002/ece3.70133
==== Body
pmc1 INTRODUCTION

Climate is a key determinant of the large‐scale spatial distribution of plant species and ecosystems (McGill, 2010; Willis & Whittaker, 2002). Multiple lines of evidence show changes in climate as a cause of recent and paleontological floral extinctions (Cahill et al., 2013; McElwain & Punyasena, 2007) and as a predicted cause for future extinctions under anthropogenic climate change (Thomas et al., 2004; Urban, 2015). It is important to understand how species niches have responded to past climate changes to understand the susceptibility of flora to climate‐mediated extinctions. Plant niches can be described as the set of environmental conditions in which a taxa's population can persist (Hutchinson, 1957). The fundamental niche includes all possible factors required for a species to survive, while the realized climatic niche is the estimation of the niche based on the current climate‐distribution relationship and is limited by the breadth of current climate space (Jackson & Overpeck, 2000). Therefore, the realized climatic niche of a species can change through time and space as areas of the fundamental niche appear or disappear based on the available climate within geographic space (Jackson & Overpeck, 2000), while at some taxonomic levels (i.e. within genera), the fundamental climatic niche can evolve through speciation (Emery et al., 2012; Lu et al., 2023).

Plants can respond to climate changes in three ways: move, adapt, or go extinct. Plant movement in response to climate change has been extensively studied. Plant species may respond to climate changes through geographic migration to track suitable climates thus retaining their climatic niches (i.e. niche conservatism). For example, in response to deglacial warming, many taxa in Eastern North America migrated northward to follow cooler climates (Davis, 1981; Davis & Shaw, 2001; Prentice et al., 1991) and have maintained their climatic niches through the last 18,000 years (Wang et al., 2023). A similar pattern can be seen in response to cooling climates between the last interglacial and the last glacial maximum in flowering stonecrops (Rhodiola spp.) on the Qinghai‐Tibetan Plateau which migrated to lower elevations and latitudes to track warmer climates (You et al., 2018). On longer time scales, the genus Laurus has had a stable climatic niche over the last 3 million years, as the reconstructed paleoclimate of Laurus' paleogeographic range matches its modern climatic niche (Rodríguez‐Sánchez & Arroyo, 2008). This demonstrates that through the three climate shifts analyzed – the middle Pliocene, the last glacial maximum, and modern – Laurus climatic tolerances remained stable, and individuals moved to track suitable climates (Rodríguez‐Sánchez & Arroyo, 2008). Niche conservation has been invoked for many palaeoecological studies that use modern climatic tolerances of taxa as analogs to estimate past climates based on fossil records (Chevalier et al., 2020; Prebble et al., 2017; Thompson et al., 2012).

Adaptation in response to climate change can also occur in plant taxa where plants evolve to survive in new climates as observed in changes to the taxon's climatic niche, i.e. niche divergence. Paired phylogenetic and niche studies can track niche divergence through speciation events (Lu et al., 2023). For example, the niches of species within the Chrysanthemum zawadskii complex have diverged since the Pliocene in response to range expansion (Lu et al., 2023). There is evidence for family and genus level niche conservatism over large spatial, temporal, and taxonomic scales but with niche divergence at the species level. For example, most genera within the family Cyatheaceae, scaly tree ferns, show global niche conservatism, generally occupying areas that are warm, frost‐free, and have consistent rainfall (Bystriakova et al., 2011). However, geographic isolation and phylogenetic dynamics contribute to Cyatheaceae species‐level niche divergence (Bystriakova et al., 2011). Modeling efforts have also shown that the degree of speciation and niche divergence can impact species ability to survive climatic changes (Aguilée et al., 2016; Rangel et al., 2018). Over shorter time periods niche divergence has also been witnessed within species that have been introduced to new geographic areas, likely due to competitive release or expanded availability of favorable climates (Guisan et al., 2014). However, these niche changes exhibit interspecies variability, with niche conservatism observed in some and niche divergence in others (Pearman, Guisan, et al., 2008). Niche divergence is also sensitive to spatial scale, being more common at smaller, habitat scales, while niche conservatism appears more prominent at larger climatic scales (Emery et al., 2012; Silvertown et al., 2005). This suggests niche divergence may result from competition and local adaptation, but that these dynamics are not distinguishable at larger spatial scales (Emery et al., 2012, Silvertown et al., 2005).

Climate mediated extinction events in plant taxa are less well understood than migrations and adaptations. For example, of the five global mass extinction events, none coincide with mass extinctions of plant taxa (Wing, 2004). However, regional extinctions are more prevalent, and are often associated with major abiotic events (e.g. Carvalho et al., 2021; Jaramillo et al., 2010). A sister genera analysis of Pliocene‐Pleistocene floral extinctions of European temperate tree taxa with high levels of niche conservation shows that extinct genera typically occurred in warmer areas than widespread genera, suggesting that glacial cycles in combination with migration barriers caused warm‐adapted genera to go extinct while more cold‐adapted taxa persisted (Svenning, 2003).

Past climate changes have resulted in many legacies in the species composition and spatial distributions of ecosystems today (Svenning et al., 2015). To further understand how past climate changes may have interacted with extinction events we determine if legacies of past extinction events can be distinguished within taxa climatic niches. We define niche legacies as a substantial difference in the climatic niches of closely related extinct and extant taxa in the modern assemblage which could be hypothesized to have played a determining role in an extinction event. Large climate legacies would suggest that differences in climate tolerances made the extinct taxa more susceptible to extinction from past climate changes than their extant counterparts. The flora of New Zealand gives us a rare opportunity to test if taxa retain climatic niche legacies of extinction due to its exemplary Cenozoic floral extinction record, and because many genera that went extinct in New Zealand are still extant in the modern Australian flora (Lee et al., 2001).

The Cenozoic was a time of substantial climatic, geologic, and ecological change in New Zealand. Consistent with global long‐term trends of cooling climates (e.g. Westerhold et al., 2020), New Zealand lost much of the tropical and subtropical climates that existed during early periods of the Cenozoic (Conran, Bannister, et al., 2016; Crouch et al., 2020; Kennedy et al., 2002; Pancost et al., 2013), transforming into the largely temperate climates seen today (Nelson & Cooke, 2001; Prebble et al., 2017). In addition to climate changes, New Zealand underwent geological changes as terrain expanded into the alpine bioclimatic zone (Prebble et al., 2021; Winkworth et al., 2005). Although the current New Zealand flora has a high level of species endemism, there is low endemism at the genus and family level due to common Gondwanan ancestry, either through persistence since separation or later dispersal from Australia and elsewhere (McGlone et al., 2001; Winkworth et al., 2002). Therefore, the current genus‐level composition of the New Zealand flora is a result of a shared Gondwanan history and subsequent genus‐level extinctions throughout the Cenozoic. The many genera that went locally extinct in New Zealand during the Cenozoic has largely been attributed to climatic cooling but few comparative studies of the environmental tolerance of extinct and extant genera have been undertaken (Lee et al., 2001).

In New Zealand, Early Cenozoic climates were tropical to subtropical, with mean annual temperatures ranging from 20 to 23°C in the Eocene (56–34 mya, Conran, Bannister, et al., 2016; Pancost et al., 2013) but cooled to approximately 18°C in the Oligocene (34–23 mya), before increasing in the middle of the Miocene (23–5 mya), to between 18 and 20°C, accompanied by increasing temperature seasonality (Prebble et al., 2017; Reichgelt et al., 2019). During the Oligocene and Miocene, marine transgression submerged as much as 80% of the New Zealand land area (Cooper & Cooper, 1995; Lee et al., 2001). Mountain building, relative plate motion changes, and global sea level drop then increased land area resulting in the creation of cooler high‐elevation climates by the late Miocene (Barth et al., 2014; Cooper & Cooper, 1995; Lee et al., 2001; McGlone et al., 2001). The Pliocene (5–2.6 mya) and Pleistocene (2.6 mya–10,000 years ago) are marked by a sharp decrease in mean annual temperature, to 15°C and below, and increasing seasonality (Prebble et al., 2017). Average annual precipitation also increased through the Pliocene and Pleistocene, exceeding 2500 mm/year (Prebble et al., 2017). As Cenozoic climates changed in New Zealand, climate tracking may have become impossible for species as tropical and subtropical climates became geographically limited and eventually disappeared (Lee et al., 2001). Since there was no northward landmass for New Zealand genera favoring subtropical climates to migrate to, the available options were to either adapt to changing climates or go extinct.

The New Zealand Cenozoic fossil record is exemplary, including numerous sedimentary types (coal measures, lacustrine deposits) that preserve pollen, leaves, flowers, and wood of many taxa (e.g., Lee et al., 2001, 2016). A notable feature of this record is the late Cenozoic extinction of many groups that currently remain in Australia (Lee et al., 2001; Prebble et al., 2021). New Zealand and Australia share an early Cenozoic biota related to their common Gondwanan ancestry but have different climatic and extinction histories since separation in the late Cretaceous (Strogen et al., 2023), ~80 million years ago (MYA). Since separating from Gondwana, Australia maintained a larger subaerial landmass (>7 million km2) spanning a broad latitudinal range (currently 10–43° S), while New Zealand, in contrast, has less landmass (<0.3 million km2) and encompasses a narrower latitudinal range (currently 35–47° S). The persistence of many extinct New Zealand genera in Australia has been interpreted to reflect the consistent representation of tropical and subtropical environments throughout the Cenozoic (Lee et al., 2001).

The temperatures and precipitation experienced in New Zealand throughout the Cenozoic (Table 1, subtropical to temperate) also correspond with a substantial portion of the current Australian climate space, where many genera once in New Zealand now reside, mostly concentrated along the Australian east coast (Figure 1). This provides an opportunity to understand the climatic niche legacies of past New Zealand extinctions because we can directly compare climatic niches of extinct and extant taxa controlling for the available climate space. Understanding these legacies can help determine why and how some genera remained in New Zealand while others became extinct.

TABLE 1 Past climate estimates of New Zealand for mean annual temperature and annual precipitation from Conran, Bannister, et al. (2016), Mildenhall (1980), Prebble et al. (2017), and Reichgelt et al. (2019).

Cenozoic period	Temperature	Precipitation	
Eocene	20–23°C	1200–8000 mm	
Oligocene to Early Miocene	18–20°C	1000–2500 mm	
Middle Miocene to Pliocene	15–18°C	1000–3000 mm	
Pleistocene	<15°C	1600–5000 mm	

FIGURE 1 Mean annual temperature and annual precipitation of modern Australia from WorldClim Version 2 (Fick & Hijmans, 2017). Climatic niches are computed from the temperate to tropical portions of Australia (hashed area), determined by the Köppen–Geiger climate zones (Beck et al., 2018).

To test whether New Zealand climate changes influenced the extinction of genera in New Zealand, we compare the climatic niches of extant and extinct New Zealand genera based on their modern Australian distribution (Figure 2). We use the current climate distributions of genera to understand if modern realized niches show a legacy of climatic constraints that may have played a role in the extinction of genera from New Zealand. Specifically, this study looks to see if New Zealand extinct plant genera occupy different climatic niches than closely related New Zealand extant plant genera in the temperate and tropical portions of Australia which represent the full breadth of climates experienced by New Zealand during the Cenozoic (Figure 2). Niche conservatism over long time periods has been observed in other taxa (Huntley et al., 1989; Svenning, 2003) and using climate niches of modern taxa to estimate paleoclimates has been employed in many studies (Chevalier et al., 2020; Prebble et al., 2017; Reichgelt et al., 2015; Utescher et al., 2000).

FIGURE 2 Conceptual model of the link between the changing New Zealand climate space through time, the fundamental niches of extinct and extant New Zealand taxa, and the modern Australian climate space. As the New Zealand climate space shifts (orange dashed) it may result in New Zealand no longer containing climate within the fundamental niche of some taxa resulting in extinction of that taxa (NZ extinct taxa). Other taxa fundamental niches may still overlap with the modern New Zealand climate allowing for the species to survive the New Zealand climate shift (NZ extant taxa). To understand differences in past niche dynamics these species can be observed in the modern Australian climate.

We hypothesize that as the New Zealand climate cooled during the Miocene, Pliocene, and Pleistocene, the available climates in New Zealand moved outside of the fundamental niches of those species that went extinct in New Zealand (Figure 2). Because Cenozoic New Zealand climates are contained within modern‐day tropical and temperate zones of Australia (Table 1, Figure 1), we expect that the modern Australian realized niches of extinct New Zealand genera will be warmer than extant New Zealand genera, indicating that the loss of warm climates contributed to the extinction of those genera in New Zealand. If climate legacies of Cenozoic extinction events are present in the modern assemblage, then paired extinct and extant taxa will have substantially different climatic niches. To identify potential climatic correlates of extinctions of taxa in New Zealand we (i) quantify differences in total niche breadth, (ii) analyze similarities of climatic niches, and (iii) compare individual climate parameters between extinct and extant New Zealand genera in Australian climate space.

2 METHODS

2.1 Taxa selection and data

Con‐familial focal genera were selected to represent families that are present in Australia and have at least one genus that is extant in New Zealand and one that is only represented in the fossil record in New Zealand (i.e. Figure 2). This process resulted in the selection of 13 genus pairs across nine families. Global occurrence data was downloaded from the Global Biodiversity Information Facility (GBIF, 2023) for all nine families and filtered to include focal genera and distribution points in Australia (Appendix 1). Only data based on living or herbarium specimens (fossil specimens were not used) with a geographic location with uncertainty of less than half the width of a grid cell from the climate data (2250 m) were used.

TABLE 2 Climatic niche similarity between a New Zealand‐extinct genus and a related genus from the same family that is extant in New Zealand and Australia for nine families.

Family	Pair	Schoener's D	Niche similarity test	Niche difference test	
Araucariaceae	Araucaria vs. Agathis	0.33	–	–	
Argophyllaceae	Argophyllum vs. Corokia	0.37	‡	–	
Elaeocarpaceae	Sloanea vs. Elaeocarpus	0.75	‡	–	
Sloanea vs. Aristotelia	0.09	–	–	
Euphorbiaceae	Mallotus vs. Euphorbia	0.44	–	–	
Fabaceae	Acacia vs. Sophora	0.32	†	–	
Caesalpinia vs. Sophora	0.72	‡	–	
Lauraceae	Cryptocarya vs. Beilschmiedia	0.83	‡	–	
Cryptocarya vs. Litsea	0.85	‡	–	
Onagraceae	Ludwigia vs. Epilobium	0.06	–	–	
Podocarpaceae	Microcachrys vs. Podocarpus	0.24	†	–	
Sapindaceae	Cupaniopsis vs. Alectryon	0.64	‡	–	
Mischocarpus vs. Alectryon	0.52	†	–	
Note: Climatic niche similarity is based on the realized climatic niche of the two genera derived from their current Australian temperate and tropical distributions. Schoener's D values range from 0 (no overlap) to, 1 = full overlap = identical climatic niches and niche similarity and niche difference tests calculates a significance value for one genus while using a null model for the other genus, therefore each pair has two results (p‐values adjusted for multiple comparisons, – = both directions >.05, † = one direction <.05, ‡ = both directions <.05). Genus pairs that are non‐similar in at least one test are indicated in bold. Pairs are shown as New Zealand extinct vs. New Zealand extant genera.

2.2 Environmental space and climatic niche data

Using con‐familial genera of regionally extinct species attempts to control for the evolutionary changes to climatic niches that may occur over millions of years due to niche conservatism within many plant families that preserve genus‐climate relationships through time (Bystriakova et al., 2011; Svenning, 2003). Due to the millions of years involved, it is possible that these genera have adapted through time to novel climates in Australia, especially dryland ecosystems. We therefore restricted this study only to the temperate, sub‐tropical, and tropical portions of the Australian continent (Figure 1, hashed area). Here we compare the climatic niches of con‐familial pairs with one genus that is extinct in New Zealand and one genus that is extant in New Zealand that are both present in the tropical and temperate portions of Australia (Figure 2), which is a proxy for Cenozoic New Zealand climates (Table 1). Climatic niches for each taxon were quantified using their current distribution in temperate and tropical Australia. The areas that qualify as temperate and tropical Australia were determined by the Köppen‐Geiger climate zones (Beck et al., 2018) that are broadly consistent with past New Zealand climates (Table 1).

To account for geographic bias in GBIF data due to the clustering of points in cities or along major roads, the data was downscaled to match the resolution of the climate grid, 2.5 arc‐min, so that each grid cell associated with genera presence has the same weight in niche calculations. Some GBIF species occurrence points (less than 1 percent) did not overlap with the climate grid due to the climate grid's lower precision along coastlines and waterways and were excluded. After downscaling over 457,000 occurrence points remained for all nine plant families (Appendix 1).

The climate grid contains all 19 climate variables from WorldClim Version 2 at a spatial resolution of 2.5 arc‐min (approximately 4.5 km2) which are commonly used for niche modeling (Fick & Hijmans, 2017; Warren et al., 2008). Bioclimate variables included are annual mean temperature (bio 1), mean diurnal range (bio 2), isothermality (bio 3), temperature seasonality (bio 4), the maximum temperature of the warmest month (bio 5), minimum temperature of the coldest month (bio 6), temperature annual range (bio 7), mean temperature of the wettest and driest quarter (bio 8 and 9), mean temperature of warmest and coldest quarter (bio 10 and 11), annual precipitation (bio 12), precipitation of wettest and driest month (bio 13 and 14), precipitation seasonality (bio 15), precipitation of wettest and driest quarter (bio 16 and 17), and precipitation of warmest and coldest quarter (bio 18 and 19).

2.3 Niche overlap statistics

Niche overlap was calculated between genus pairs based on climatic niches in Australia with the aim of finding differences that could identify climate parameters responsible for the extinction of genera in New Zealand. The climatic niches of genera were determined using climate data associated with the distribution of each genus and performing a principal component analysis (PCA) to analyze climatic niches in two dimensions. The PCA was completed in R version 4.2.2 (R Core Team, 2022) using the package ade4 version 1.7‐20 (Dray & Dufour, 2007) and the package ecospat (Di Cola et al., 2017). The method is consistent with calculating the most accurate measure of climatic niche overlap as described in Broennimann et al. (2012).

2.3.1 Schoener's D

The Schoener's D metric (Schoener, 1968) was used to quantify the similarities and differences of the climatic niche of paired taxa. The metric calculates the amount of overlap between niche accounting for the frequency of occurrence for each input variable for both species, respectively. The Schoener's D metric is a value between 0 and 1 representing no niche overlap to complete niche overlap, respectively (Schoener, 1968). Schoener's D is typically used for niche comparison of microhabitat or diet and may not be appropriate for all environmental niche models (Warren et al., 2008). However, when compared against results from the Hellinger distance, no significant difference was observed (Warren et al., 2008). Calculation of Schoener's D first calculates an occurrence density grid for each genus using inputs from the first two PCA axes for each genus and the full available environmental space. Schoener's D overlap is then calculated based on the genera density grids while correcting for the density of the availability of climates within the study area (Di Cola et al., 2017). Schoener's D was completed in R version 4.2.2 (R Core Team, 2019) using the package ecospat version 3.4 (Di Cola et al., 2017).

2.3.2 Niche similarity

The statistical significance of Schoener's D was determined following the methods of Warren et al. (2008) which compare the climatic niche of one taxon to null model that incorporates differences in the available climate space between ranges to determine if taxa are more similar than random. This is achieved by calculating the niche overlap using the actual, observed niche of one genus and comparing a niche that is randomly generated from the available climate space of this first genus, in this case the available climate space is the same for both species (Broennimann et al., 2012; Warren et al., 2008). A histogram of the 1000 null model generated Schoener's D values is then compared to the actual Schoener's D value using a two‐tailed test on the null distribution model to determine if the actual Schoener's D is greater than the 95% confidence interval's upper value, which tests for niche similarity, or lower value, which tests for niche difference test (Broennimann et al., 2012; Warren et al., 2008). This calculation is directional, with each pair having two estimates for niche similarity and niche difference. Genus pairs with at least one direction with p < .05 in the similarity test were considered to have similar climatic niches and those with at least one direction with p < .05 in the difference test were considered to have different climatic niches. Calculation of niche similarity was competed in R version 4.2.2 (R Core Team, 2019) using the package ecospat version 3.4 (Di Cola et al., 2017).

2.3.3 Kernel density

In addition to analyses of overlap in multidimensional climatic niche space, climate overlap of each genus pair was analyzed separately for all climate variables using kernel density. Climatic distributions of genera are not normally distributed; therefore, the non‐parametric kernel density estimation (KDE) was used. To estimate overlap, KDE estimations were calculated for each genus and climate variable individually, these KDE estimates were then compared within the genera‐pairs and the percent overlap was computed. The KDE calculations were performed in R (R Core Team, 2019) using the package overlapping version 2.1 (Pastore, 2018).

2.4 Niche volume

Niche volume was calculated as the overlap between each individual genus climatic niche and the temperate and tropical Australian climate space. This was completed in two ways: (i) using Schoener's D as described above, and (ii) using the percent overlap in a multidimensional kernel density estimation. The Schoener's D niche volume and overlap value was generated for each genus with values between 0 and 1 representing none to full overlap between each genus and the temperate and tropical Australian climate space. Schoener's D was completed in R version 4.2.2 (R Core Team, 2019) using the package ecospat version 3.4 (Di Cola et al., 2017). The kernel density estimation of niche volume creates a two‐dimensional grid (PCA1 and PCA2) of kernel density and delimits the 95th quantile area for the temperate and tropical Australian climate space and each individual genus niche. The niche volume is then calculated as the percent of the temperate and tropical Australian climate space that is occupied by each genus. Kernel density niche volume analysis was completed in R version 4.2.2 (R Core Team, 2019) using the package ks version 1.14.1 (Duong, 2023).

3 RESULTS

The extinction dates of genera in New Zealand are based on the youngest known fossils. All extinctions occurred within the Cenozoic (Figure 3). Caesalpinia had the earliest extinction, during the Oligocene approximately 28 million years ago (MYA), and Acacia and Microcachrys were the most recent extinctions, occurring approximately 1 MYA in the Pleistocene (Figure 3, Appendix 2).

FIGURE 3 Presumed date of extinction from New Zealand for each focal genus and full list of extinct and extant genera. Extinction dates are based on the youngest known fossil record in New Zealand. Dates displayed in millions of years ago (MYA) and genera are color coded by family; for references see Appendix 1.

Australia currently contains climate regimes ranging from subtropical/tropical climates in the north, to temperate climates along portions western and eastern coasts, to dry desert climates in the interior (Figure 1). The temperate and tropical zones of Australia average annual temperature ranges from 5.0 to 29.1°C and annual precipitation from 318 to 3880 mm (Figure 1) which covers the entire range of New Zealand Cenozoic temperatures (>15 to 23°C, Table 1) and most of the New Zealand Cenozoic precipitation with the exception of extreme high precipitation in the Eocene and Pleistocene (1000 to 8000 mm, Table 1). The first axis in the temperate and tropical Australian climate PCA (PCA 1) explains 63.7% of the variation in the data and is strongly negatively correlated with temperature and precipitation seasonality variables (Figure 4). The second axis in the PCA (PCA 2) explains 20.3% of the variation in the data and is strongly positively correlated with annual precipitation and both daily and yearly temperature seasonality (Figure 4). Cenozoic New Zealand temperatures fall within the range of climates currently experienced in temperate and tropical Australia (Table 1). Cenozoic New Zealand precipitation is also well represented in the modern Australian climate space, however, precipitation extremes that are seen in Cenozoic New Zealand climate are missing from modern Australian climates (Figure 1 and Table 1). The modern New Zealand climate space largely overlaps with that of modern Australian climate space but is densely concentrated in the temperate and cool portions of Australia (Figure 4).

FIGURE 4 PCA results showing (a) the 99.99% density of tropical and temperate Australian climate points included in PCA construction (solid black outline) and the 99.95% density of climate space of different ecological areas: Tropical Australia (dashed green) and temperate Australia (dashed purple). Also shown are the distribution of areas not included in construction of the PCA: Desert Australia (dashed light gray) and modern New Zealand (dashed dark gray). (b) Description of main climate axes and where general climate zones occur. (c) Contribution of each bioclimatic variable along the first two PCA axes. Bioclimate variables that are included in the PCA are described in Section 2.

3.1 Niche volume

We anticipated that if extinct New Zealand genera were extirpated by cooling Cenozoic climates in New Zealand, they would be restricted to the warmer climates of Australia and would have smaller climatic niches than extant genera in Australia. However, no significant difference was found in climatic niche volume between genus pairs. The climatic niche volume amongst genera varied considerably from low niche volume in Australian temperate and tropical climate space (Corokia, Schoener's D = 0.06, kernel density = 0.01) to moderate (Euphorbia, Schoener's D = 0.62, kernel density = 0.79). Although individual taxa vary widely, no difference is observed between the niche volume of New Zealand extant and extinct genera, with median niche volumes of Schoener's D values 0.37 and 0.41, respectively (p = .33 two‐tailed paired t‐test, Figure 5) and kernel density values of 0.38 and 0.38, respectively (p = .94 two‐tailed paired t‐test, Figure 5). Additionally, no correlation in niche volume is seen when comparing individual genus pairs (Schoener's D p = .72, r 2 = .01, kernel density p = .76, r 2 = .01; Figure 5).

FIGURE 5 (a) Niche volume for extant (blue) and extinct (orange) genera in temperate and tropical Australia. Niche volume is calculated as the niche overlap (Schoener's D) and kernel density overlap of each genus' niche with the temperate and tropical Australian climate space. Comparison shows no significant difference in climatic niche volume between New Zealand extant and extinct genera in temperate and tropical Australia (Schoener's D p = .32, Kernel density p = .94, paired two‐tailed t‐test). (b) Comparison of niche volume between genus pairs, shows no correlation between individual genus pairs (Schoener's D, black points and black line, p = .72, r 2 = .01; Kernel density, open points and dashed line, p = .76, r 2 = .01).

3.2 Niche overlap and similarity

A climatic niche legacy of extirpation by Cenozoic climate changes in New Zealand would be indicated by extinct climatic niches that are either not significantly similar to their extant niche counterpart (similarity test p > .05 in both directions) or are significantly different from their extant niche counterpart (difference test p < .05 in both directions). It is expected that if cooling Cenozoic climates were the extinction driver, extinct niches would occupy warmer portions of the modern Australian climate space than their extant counterparts. Comparison of the shape and location of extinct and extant genera's climatic niche within the Australian climate space shows that a majority of the genus pairs have significantly similar niches (69% or 9 out of 13 pairs; Figure 6, Table 2) and no niche pairs have significantly different niches (Table 2).

FIGURE 6 Realized climatic niche overlap for New Zealand extant (blue) and extinct (orange) genus pairs in Australian climate space for (a) climatically non‐similar pairs (similarity test p > .05 in both directions) and (b) climatically similar pairs. Outline of the tropical and temperate Australian climate space is depicted by the gray line.

Of the four pairs that are not statistically similar (similarity test p > .05 in both directions), three pairs show the expected warmer extinct niche pattern, Sloanea vs. Aristotelia, Ludwigia vs. Epilobium and Mallotus vs. Euphorbia (Figure 6). Both Sloanea vs. Aristotelia and Ludwigia vs. Epilobium also have very low niche overlap (Schoener's D = 0.09 and 0.06, respectively) while Mallotus vs. Euphorbia has comparatively high overlap (Schoener's D = 0.44; Figure 6, Table 2). This can be seen when comparing niche differences along individual climate gradients where both Sloanea vs. Aristotelia and Ludwigia vs. Epilobium have markedly different kernel density overlap along the mean annual temperature gradient (12% and 7%, respectively) with extinct genera occupying warmer areas. The kernel density overlap of precipitation seasonality for the two pairs is slightly more similar (29% and 19%, respectively) with extinct genera in areas with larger precipitation seasonality (Figure 7). This suggests that the extinct genera Sloanea and Ludwigia occupy climates that are warmer with more distinct wet and dry seasons than their extant counterparts Aristotelia and Epilobium. These patterns are also seen with Mallotus vs. Euphorbia but with greater kernel density overlap values than the other two pairs (Figure 7). Kernel density overlap in temperature annual range and annual precipitation are large and likely not driving the differences between these genera (Figure 7).

FIGURE 7 Kernel density overlap between taxa with significantly different climatic niches for key climate variables based on climatic niche of temperate and tropical Australia. Showing extinct niches in orange (Sloenea, Ludwigia, Mallotus, and Araucaria) and extant niches in blue (Aristotelia, Epilobium, Euphorbia, and Agathis). Percent kernel density overlap between each pair is shown within each graph.

The final genus pair that is not similar, Araucaria vs. Agathis, shows the opposite trend with the extinct genus occupying areas that are cooler with less precipitation seasonality (Schoener's D = 0.33; Figure 6 and Table 2). Comparison of individual climate gradients show low overlap with the extinct genus found in cooler areas that have consistent year‐round precipitation but also have slightly more seasonal temperature variation suggesting that the extinct Araucaria occupies more temperate climates than its extant counterpart Agathis (Figures 6 and 7).

4 DISCUSSION

4.1 Potential climate legacies of extinction events

Through analysis of climatic niche similarities and differences within genus pairs, we found that four extinct taxa show potential climate legacies with non‐similar niches to their extant counterpart in the similarity test; however, this did not translate to the difference test which showed no pairs had significantly different niches. Three extinct taxa have warmer niches than their closely related extant taxa pair (Sloanea, Ludwigia, and Mallotus; Figures 6 and 7) which provides support for the importance of climate cooling events in the extinction of genera from New Zealand during the Cenozoic. The last occurrence of two taxa with warmer niches (Sloanea and Mallotus) coincided with two periods of pronounced climate cooling in New Zealand during the last 30 million years (Prebble et al., 2017). Although the last occurrence of Sloanea (15–18 MYA) occurs before the Middle Miocene Climatic Optimum (MMCO, ~14 MYA), the extinction date of Sloanea is particularly susceptible to the Signor–Lipps effect, which is the low likelihood that the fossil of the last individual is recovered therefore the true extinction date would occur after the date of the last fossil recovered (Signor III & Lipps, 1982). This is due to macrofossils recovered from abundant early‐middle Miocene lacustrine sediments of the Manuherikia Group (Pole, 1993, 2007), whereas fossil‐bearing strata of the late Miocene are much rarer. Therefore, the extinction date of Sloanea may have occurred in the late Miocene during climate cooling (Prebble et al., 2017) but be undetectable in the fossil record. The late Miocene cooling after the MMCO was particularly pronounced, with minimum temperatures cooling by approximately 6°C, even while maximum temperatures remained stable at approximately 20°C (Prebble et al., 2017). This change in temperature may have had a pronounced impact on Sloanea whose average mean annual temperature is 6.4°C warmer than its extant counterpart Aristotelia (Figure 7). The last occurrence of Mallotus (2–3 MYA) coincides with a second period of pronounced climatic cooling at the Pliocene‐Pleistocene transition when both maximum and minimum temperatures in New Zealand decreased by 1.5 and 6°C, respectively (Prebble et al., 2017). This is consistent with the warmer climatic niche of Mallotus, which occupies areas that are on average 2.2°C warmer than its extant counterpart Euphorbia.

The third taxa with a warmer niche, Ludwigia, has an extinction in New Zealand (23–25 MYA) that coincides with the Oligocene‐Miocene transition (~23.7–22.7 MYA), which is a period of regional cooling associated expansion of Antarctic glaciation (Liebrand et al., 2017). During a period of 1 million years, ice sheets expanded to approximately the size of the modern East Antarctic Ice Sheet and subsequently disintegrated to the pre‐transition volume (Liebrand et al., 2017). This suggests a period of rapid regional cooling that is consistent with extinction of Ludwigia, which occupies areas that are on average 10.8°C warmer than its extant counterpart Epilobium (Figure 7).

Adaptations to climate cooling involve multiple physiological processes and, in thermophilic species, may limit growth and disrupt species interactions (Cahill et al., 2013). In Europe, most widespread temperate taxa have a higher degree of cold‐tolerance than regionally extinct or relic species due to environmental filtering during the cool Pliocene and Pleistocene (Svenning, 2003). This is also seen in Japan and Tasmania where many cold‐intolerant plant species went extinct in response to early Pleistocene cooling (Huang et al., 2018; Jordan, 1997). However, in New Zealand many temperate tree taxa still share functional traits typically associated with tropical species, such as thin bark and gender dimorphism (McGlone et al., 2016). Additionally, warm adapted taxa, such as Ludwigia, have recently been introduced to New Zealand and in other introduced regions Ludwigia is predicted to expand its range in response to increased temperatures from anthropogenic climate change (Gillard et al., 2017).

Araucaria occupies areas in Australia that are cooler and with less seasonal variation in precipitation and temperature than its New Zealand extant counterpart Agathis but has similar overall precipitation (Figures 6 and 7). The last occurrence of Araucaria (2–3 MYA) occurs during a period of climate cooling during the Pliocene‐Pleistocene transition. This was also a period of increased total and more variable precipitation in New Zealand (Prebble et al., 2017) which may have been a determinant of Araucaria extinction in New Zealand due to Araucaria occupying areas with lower precipitation seasonality than its extant counterpart Agathis (Figure 7). Araucaria is also no longer present in southern Australia although fossil evidence indicates it was common in the past (Hill, 1990) suggesting that Araucaria migrated northward to more favorable temperate and subtropical climates in Australia. However, migration may not have been possible in New Zealand, due to the limited northern extent of the landmass. The temperature trends during the Pliocene‐Pleistocene transition would suggest that Araucaria would outcompete its extant counterpart during this period, with an average mean annual temperature of 2.5°C cooler than the extant Agathis. This effect is evident in South America where Araucaria forests have decreased their range since the LGM in response to warming and are predicted to suffer further range restrictions due to anthropogenic climate change (Bergamin et al., 2019). However, pronounced glacial cycles that developed during the Pliocene‐Pleistocene transition (Mildenhall et al., 2004; Naish et al., 1997) may negatively impact taxa that occupy areas with low seasonal variation such as Araucaria. For example, the disappearance of Miocene seasonal climates has been suggested to have led to the loss of other flora species in New Zealand, such as the deciduous Nothofagaceae species which still occupy highly seasonal climates in South America (Reichgelt et al., 2022).

4.2 Extinction events lacking climatic niche legacies

Although significant climatic niche legacies are not evident in the results of most taxa investigated here, many went extinct during the Miocene, Pliocene, and Pleistocene epochs (Acacia, Araucaria, Argophyllum, Cryptocarya, Cupaniopsis, Ludwigia, Mallotus, Microcachrys, Mischocarpus, and Sloanea) associated with pronounced long‐term, and sometimes more accelerated climate cooling or rapid climate fluctuations. Only one taxa extinction occurred during the relatively climatically stable Oligocene (Caesalpinia). The lack of a consistent discriminating climate indicator for these genus pairs suggests that extinction in New Zealand may be explained by two mechanisms: (i) niches are conserved, suggesting that the mechanism of extinction was not climate change directly but either an indirect effect impacting taxa environmental tolerance or the result of a non‐climatic factor, such as habitat loss, or (ii) the climatic niches have evolved to occupy overlapping climates in Australia and any climate legacies that may have existed at the time of extinction have been lost.

The degree to which climatic niches are conserved through time may impact our ability to detect climatic niche legacies, with extinct taxa that underwent high levels of niche evolution potentially obscuring detection of climatic niche legacies. However, niche conservation during periods of large climatic changes has been reported in various taxa. A study comparing sister taxa in Europe, North America, and East Asia found that cool‐temperate tree genera have conserved their niches since the Pliocene (Svenning, 2003). Niche conservation has also been witnessed in Fagus between the North American and European species which were geographically separated during the breakup of Laurasia (Huntley et al., 1989). The climatic niche of Laurus (Lauraceae) in Europe has also remained conserved since the middle Pliocene (3 MYA), consistently occupying areas that are warm and wet with little seasonal variation (Rodríguez‐Sánchez & Arroyo, 2008). The climatic niche of the tree species Carpinus betulus and Picea abies in Europe have been conserved since the middle Holocene despite large changes in environmental conditions (Pearman, Randin, et al., 2008). Assuming niche conservatism within taxa suggests that the lack of climate niche legacy is the result of the first proposed mechanism – that niches are conserved and that either non‐climatic or indirect climatic factors were the drivers of extinction.

Non‐climatic factors such as habitat loss may have played a role in the Caesalpinia whose extinction age coincides with the period of maximum marine transgression in New Zealand (28 MYA) and associated habitat loss from rising sea levels (Cooper & Cooper, 1995). Extinction caused habitat loss from sea‐level rise was also seen during the Pleistocene in Bermuda when a large reduction in land area was followed by the extinction of many bird species (Olson & Wingate, 2000). Marine transgression also caused New Zealand to split into a series of small islands (Cooper & Cooper, 1995) which may have increased the extinction rate during this period due to decreased land area, genetic isolation, or increased competition for limited resources (Rosenzweig, 1995). However, there is little evidence that changes in land area is a contributor to extinctions in New Zealand, as the maximum marine transgression of the Oligocene witnessed markedly few extinctions (Lee et al., 2001). Other potential non‐climatic extinction mechanisms could include geologic changes in the mid‐ to late‐Cenozoic resulting in the loss of soil types, such as deep, highly weathered nutrient poor soils, that developed during the early Cenozoic (McGlone et al., 2016).

Climate changes may have also had indirect impacts on genera's ability to survive through changes in species interactions or disturbance regimes that are not reflected in the current realized niche (Cahill et al., 2013). Families' dependent on insect pollination, for example, may have been susceptible to changes in abundance of sub‐tropical insects. Conran et al. (2014) found that of the 30 genera in 23 families identified at the 23 million‐year‐old Foulden Maar deposit in New Zealand, all had retained their reproductive niches, suggesting co‐extinction could play a role in plant extinctions. Additionally, it is suggested that after the maximum marine transgression of the late Oligocene, New Zealand may have contained many relatively vacant niches resulting in a high level of long‐distance dispersal and colonization from Australia, including putative evidence for Eucalyptus (Lee et al., 2012; Pole, 1993), which appeared in New Zealand during the Miocene (Lee et al., 2016; McGlone et al., 2001). These colonization events may have introduced competitors, leading to the extinction of genera in the early to middle Miocene (extinction of Sloanea, Cryptocarya, and Argophyllum). However, further studies are needed to determine if these events resulted in extinction. Additionally, changes in disturbance regimes may have influenced taxa survival. There is evidence that fire was more prevalent in the mid‐Miocene when fire adapted genera, such as Acacia, were prevalent in New Zealand, but subsequently went extinct (McGlone et al., 2016). A gradual change to less fire‐prone vegetation, as rainforests expanded, may have resulted in the extinction of pyrophilic genera, as is observed in southeast Australia today (Read & Hill, 1985).

Although many taxa conserve their climatic niches through time, some taxa respond to environmental changes through niche evolution. Studies that aim to determine niche conservatism and divergence within families show evidence of general niche conservatism but with occasional taxa undergoing niche divergence at the species level (Bystriakova et al., 2011). Tracing niche evolution within genera shows that speciation can cause niche divergence. This is seen in the Chrysanthemum zawadskii species complex in which niche evolution occurred with speciation as they migrated into new geographic areas (Lu et al., 2023). This can also be seen within a single species, such as the European Juniperus communis, whose niche has diverged from its middle Holocene niche, likely due to changes in habitat and competition dynamics rather than climate tolerance (Pearman, Randin, et al., 2008). Enhanced evolutionary pressure from competition can also be seen in other aspects of species' traits, such as stature, with more trait divergence amongst species with competitors than those without (Burns, 2022). Studies inspecting evolutionary responses to competition and habitat availability have largely been focused on at the species level (e.g. Ackerly et al., 2006; Emery et al., 2012; Lu et al., 2023) and it remains unclear if the same dynamics occur at higher taxonomic levels.

5 CONCLUSION

Using climatic niche comparisons, we show no universal climatic niche legacy in the Australian modern assemblage of New Zealand extinct and extant genera. We do provide some evidence that cold intolerance was potentially responsible for three extinctions during Cenozoic cooling events. This supports previous hypotheses about the importance of climate filtering in the composition of modern flora (Lee et al., 2001). In addition to temperature, differences in precipitation seasonality were found to be important in one investigated genus pair. However, the majority of the taxa examined here did not show detectable climatic niche legacies of extinction. This suggests that either climate had an indirect role in extinction event, such as through changes in species interactions, or that changes in non‐climatic mechanisms, such as geologic or tectonic changes, prevented persistence in New Zealand. Our results suggest that the mechanisms of Cenozoic extinctions of New Zealand plant genera may be more complex than taxa reaching environmental tolerance limits from cooling climates and more research is needed to determine New Zealand Cenozoic extinction drivers. Recommendations for future work include: (i) analyzing the relationship between phylogenetic similarity and niche similarity at various taxonomic levels to better understand how taxa evolution impacts current environmental niches; (ii) conducting niche comparison studies within genera to determine the level of niche conservatism, especially for those with distributions spanning the global temperate to tropical climate space, such as South America; (iii) analyzing individual genus‐specific species interactions to determine how prevalent obligate mutualisms are amongst extinct genera and how susceptible taxa with symbiosis may have been to cooling climates.

AUTHOR CONTRIBUTIONS

Nora Schlenker: Data curation (lead); formal analysis (lead); investigation (lead); methodology (equal); writing – original draft (lead); writing – review and editing (equal). William G. Lee: Conceptualization (equal); resources (equal); writing – review and editing (equal). Tammo Reichgelt: Conceptualization (equal); resources (equal); writing – review and editing (equal). Ralf Ohlemüller: Conceptualization (equal); methodology (equal); resources (equal); writing – original draft (supporting); writing – review and editing (equal).

CONFLICT OF INTEREST STATEMENT

The authors declare no known competing interests relating to the work presented.

ACKNOWLEDGEMENTS

We acknowledge the efforts of volunteer recorders and nature observers who populate databases that are used for studies of the type presented here. Without their contributions, large‐scale and long‐term insights into global biodiversity patterns would not be possible.

DATA AVAILABILITY STATEMENT

The data that support the findings of this study are openly available from GBIF (2023) with specific data references included in Appendix 1. Code for all analyses is openly available at https://zenodo.org/doi/10.5281/zenodo.12534920.

APPENDIX 1 Full references for species distribution data downloads and summary of data points

1.1 GBIF.org 29 April 2019a. GBIF Occurrence Download https://doi.org/10.15468/dl.3hry6r.

GBIF.org 29 April 2019b. GBIF Occurrence Download https://doi.org/10.15468/dl.m71cud.

GBIF.org 29 April 2019c. GBIF Occurrence Download https://doi.org/10.15468/dl.mrehyy.

GBIF.org 29 April 2019d. GBIF Occurrence Download https://doi.org/10.15468/dl.nlwct2.

GBIF.org 29 April 2019e. GBIF Occurrence Download https://doi.org/10.15468/dl.pljxpc.

GBIF.org 29 April 2019f. GBIF Occurrence Download https://doi.org/10.15468/dl.pqi1zj.

GBIF.org 29 April 2019g. GBIF Occurrence Download https://doi.org/10.15468/dl.t472lw.

GBIF.org 29 April 2019h. GBIF Occurrence Download https://doi.org/10.15468/dl.ukoj92.

GBIF.org 29 April 2019i. GBIF Occurrence Download https://doi.org/10.15468/dl.xhh2os.Family	Genera	GBIF occurrence points	Climate grid cells	Number of species	
Araucariaceae	Araucaria	1645	644	6	
Agathis	808	154	4	
Argophyllaceae	Argophyllum	634	154	4	
Corokia	252	28	1	
Elaeocarpaceae	Sloanea	2766	762	4	
Elaeocarpus	16,789	3408	26	
Aristotelia	347	217	3	
Euphorbiaceae	Mallotus	5969	1419	15	
Euphorbia	13,175	3706	79	
Fabaceae	Acacia	351,665	22,421	854	
Caesalpinia	1047	203	16	
Sophora	588	62	7	
Lauraceae	Cryptocarya	23,078	2644	52	
Beilschmiedia	3244	505	12	
Litsea	5406	1471	13	
Onagraceae	Ludwigia	3646	1603	10	
Epilobium	9599	4256	22	
Podocarpaceae	Microcachrys	137	43	1	
Podocarpus	2534	502	12	
Sapindaceae	Cupaniopsis	5545	1326	15	
Mischocarpus	2588	559	11	
Alectryon	5806	1303	14	

APPENDIX 2 Presumed extinction dates with associated references.

2.1 Family	Genera	Last occurrence (MYA)	References	
Araucariaceae	Araucaria	2–3	Pole (2008a)	
Argophyllaceae	Argophyllum	15–18	Pole (2008b)	
Elaeocarpaceae	Sloanea	15–18	Pole (1993b), Conran et al. (2014)	
Euphorbiaceae	Mallotus	2–3	Lee et al. (2010), Conran, Lee, and Reichgelt (2016)	
Fabaceae	Acacia	1	Mildenhall (1980), Mildenhall (1989)	
Caesalpinia	28	Pocknall (1982)	
Lauraceae	Cryptocarya	15–18	Pole (2007)	
Onagraceae	Ludwigia	23–25	Pocknall (1991), Bannister et al. (2005)	
Podocarpaceae	Microcachrys	1	Jordan et al. (2010), Mildenhall (1980)	
Sapindaceae	Cupaniopsis	2–3	Mildenhall (1980), Lee et al. (2016)	
Mischocarpus	2–3	Mildenhall et al. (1992)
==== Refs
REFERENCES

Ackerly, D. D. , Schwilk, D. W. , & Webb, C. O. (2006). Niche evolution and adaptive radiation: Testing the order of trait divergence. Ecology, 87 , S50–S61.16922302
Aguilée, R. , Raoul, G. , Rousset, F. , & Ronce, O. (2016). Pollen dispersal slows geographical range shift and accelerates ecological niche shift under climate change. Proceedings of the National Academy of Sciences of the United States of America, 113 , E5741–E5748. 10.1073/pnas.1607612113 27621443
Bannister, J. M. , Lee, D. E. , & Raine, J. I. (2005). Morphology and palaeoenvironmental context of Fouldenia staminosa, a fossil flower with associated pollen from the Early Miocene of Otago, New Zealand. New Zealand Journal of Botany, 43 (2 ), 515–525.
Barth, N. C. , Kulhanek, D. K. , Beu, A. G. , Murray‐Wallace, C. V. , Hayward, B. W. , Mildenhall, D. C. , & Lee, D. E. (2014). New c. 270 kyr strike‐slip and uplift rates for the southern alpine fault and implications for the New Zealand plate boundary. Journal of Structural Geology, Continental Transpressive Fault Zones, 64 , 39–52. 10.1016/j.jsg.2013.08.009
Beck, H. E. , Zimmermann, N. E. , McVicar, T. R. , Vergopolan, N. , Berg, A. , & Wood, E. F. (2018). Present and future Köppen‐Geiger climate classification maps at 1‐km resolution. Scientific Data, 5 , 180214. 10.1038/sdata.2018.214 30375988
Bergamin, R. S. , Debastiani, V. , Joner, D. C. , Lemes, P. , Guimarães, T. , Loyola, R. D. , & Müller, S. C. (2019). Loss of suitable climatic areas for araucaria forests over time. Plant Ecology and Diversity, 12 , 115–126. 10.1080/17550874.2019.1618408
Broennimann, O. , Fitzpatrick, M. C. , Pearman, P. B. , Petitpierre, B. , Pellissier, L. , Yoccoz, N. G. , Thuiller, W. , Fortin, M.‐J. , Randin, C. , Zimmermann, N. E. , Graham, C. H. , & Guisan, A. (2012). Measuring ecological niche overlap from occurrence and spatial environmental data. Global Ecology and Biogeography, 21 , 481–497. 10.1111/j.1466-8238.2011.00698.x
Burns, K. C. (2022). The paradox of Island evolution. Journal of Biogeography, 49 , 248–253. 10.1111/jbi.14303
Bystriakova, N. , Schneider, H. , & Coomes, D. (2011). Evolution of the climatic niche in scaly tree ferns (Cyatheaceae, Polypodiopsida). Botanical Journal of the Linnean Society, 165 , 1–19. 10.1111/j.1095-8339.2010.01092.x
Cahill, A. E. , Aiello‐Lammens, M. E. , Fisher‐Reid, M. C. , Hua, X. , Karanewsky, C. J. , Yeong Ryu, H. , Sbeglia, G. C. , Spagnolo, F. , Waldron, J. B. , Warsi, O. , & Wiens, J. J. (2013). How does climate change cause extinction? Proceedings of the Royal Society B: Biological Sciences, 280 , 20121890. 10.1098/rspb.2012.1890
Carvalho, M. R. , Jaramillo, C. , de la Parra, F. , Caballero‐Rodríguez, D. , Herrera, F. , Wing, S. , Turner, B. L. , D'Apolito, C. , Romero‐Báez, M. , Narváez, P. , Martínez, C. , Gutierrez, M. , Labandeira, C. , Bayona, G. , Rueda, M. , Paez‐Reyes, M. , Cárdenas, D. , Duque, Á. , Crowley, J. L. , … Silvestro, D. (2021). Extinction at the end‐cretaceous and the origin of modern Neotropical rainforests. Science, 372 , 63–68. 10.1126/science.abf1969 33795451
Chevalier, M. , Davis, B. A. S. , Heiri, O. , Seppä, H. , Chase, B. M. , Gajewski, K. , Lacourse, T. , Telford, R. J. , Finsinger, W. , Guiot, J. , Kühl, N. , Maezumi, S. Y. , Tipton, J. R. , Carter, V. A. , Brussel, T. , Phelps, L. N. , Dawson, A. , Zanon, M. , Vallé, F. , … Kupriyanov, D. (2020). Pollen‐based climate reconstruction techniques for late quaternary studies. Earth‐Science Reviews, 210 , 103384. 10.1016/j.earscirev.2020.103384
Conran, J. G. , Bannister, J. M. , Reichgelt, T. , & Lee, D. E. (2016). Epiphyllous fungi and leaf physiognomy indicate an ever‐wet humid mesothermal (subtropical) climate in the late Eocene of southern New Zealand. Palaeogeography, Palaeoclimatology, Palaeoecology, 452 , 1–10. 10.1016/j.palaeo.2016.03.032
Conran, J. G. , Lee, D. E. , & Reichgelt, T. (2016). Malloranga dentata (Euphorbiaceae: Acalyphoideae): A new fossil species from the Miocene of New Zealand. Review of Palaeobotany and Palynology, 226 , 58–64.
Conran, J. G. , Lee, W. G. , Lee, D. E. , Bannister, J. M. , & Kaulfuss, U. (2014). Reproductive niche conservatism in the isolated New Zealand flora over 23 million years. Biology Letters, 10 , 20140647. 10.1098/rsbl.2014.0647 25319820
Cooper, A. , & Cooper, R. A. (1995). The Oligocene bottleneck and New Zealand biota: Genetic record of a past environmental crisis. Proceedings of the Biological Sciences, 261 , 293–302.
Crouch, E. M. , Shepherd, C. L. , Morgans, H. E. G. , Naafs, B. D. A. , Dallanave, E. , Phillips, A. , Hollis, C. J. , & Pancost, R. D. (2020). Climatic and environmental changes across the early Eocene climatic optimum at mid‐Waipara River, Canterbury Basin, New Zealand. Earth‐Science Reviews, 200 , 102961. 10.1016/j.earscirev.2019.102961
Davis, M. B. (1981). Quaternary history and the stability of Forest communities. In D. C. West , H. H. Shugart , & D. B. Botkin (Eds.), Forest succession: Concepts and application (pp. 132–153). Springer Advanced Texts in Life Sciences. Springer. 10.1007/978-1-4612-5950-3_10
Davis, M. B. , & Shaw, R. G. (2001). Range shifts and adaptive responses to quaternary climate change. Science, 292 , 673–679. 10.1126/science.292.5517.673 11326089
Di Cola, V. , Broennimann, O. , Petitpierre, B. , Breiner, F. T. , D'Amen, M. , Randin, C. , Engler, R. , Pottier, J. , Pio, D. , Dubuis, A. , Pellissier, L. , Mateo, R. G. , Hordijk, W. , Salamin, N. , & Guisan, A. (2017). Ecospat: An R package to support spatial analyses and modeling of species niches and distributions. Ecography, 40 , 774–787. 10.1111/ecog.02671
Dray, S. , & Dufour, A.‐B. (2007). The ade4 package: Implementing the duality diagram for ecologists. Journal of Statistical Software, 22 , 1–20. 10.18637/jss.v022.i04
Duong, T. (2023). Ks: Kernel smoothing . R package version 1.14.1, https://CRAN.R‐project.org/package=ks
Emery, N. C. , Forrestel, E. J. , Jui, G. , Park, M. S. , Baldwin, B. G. , & Ackerly, D. D. (2012). Niche evolution across spatial scales: Climate and habitat specialization in California Lasthenia (Asteraceae). Ecology, 93 , S151–S166. 10.1890/11-0504.1
Fick, S. E. , & Hijmans, R. J. (2017). WorldClim 2: New 1‐km spatial resolution climate surfaces for global land areas. International Journal of Climatology, 37 , 4302–4315. 10.1002/joc.5086
Gillard, M. , Grewell, B. J. , Deleu, C. , & Thiébaut, G. (2017). Climate warming and water primroses: Germination responses of populations from two invaded ranges. Aquatic Botany, 136 , 155–163. 10.1016/j.aquabot.2016.10.001
Guisan, A. , Petitpierre, B. , Broennimann, O. , Daehler, C. , & Kueffer, C. (2014). Unifying niche shift studies: Insights from biological invasions. Trends in Ecology & Evolution, 29 , 260–269. 10.1016/j.tree.2014.02.009 24656621
Hill, R. S. (1990). Araucaria (Araucariaceae) species from Australian tertiary sediments—A micromorphological study. Australian Systematic Botany, 3 , 203–220. 10.1071/sb9900203
Huang, Y. , Momohara, A. , & Wang, Y. (2018). Selective extinction within a tertiary relict genus in the Japanese Pleistocene explained by climate cooling and species‐specific cold tolerance. Review of Palaeobotany and Palynology, 258 , 1–12. 10.1016/j.revpalbo.2018.06.009
Huntley, B. , Bartlein, P. J. , & Prentice, I. C. (1989). Climatic control of the distribution and abundance of beech (Fagus L.) in Europe and North America. Journal of Biogeography, 16 , 551–560. 10.2307/2845210
Hutchinson, G. E. (1957). Concluding remarks. Cold Spring Harbor Symposia on Quantitative Biology, 22 , 415–427. 10.1101/sqb.1957.022.01.039
Jackson, S. T. , & Overpeck, J. T. (2000). Responses of plant populations and communities to environmental changes of the late Quaternary. Paleobiology, 26 (S4 ), 194–220. 10.1017/s0094837300026932
Jaramillo, C. , Ochoa, D. , Contreras, L. , Pagani, M. , Carvajal‐Ortiz, H. , Pratt, L. M. , Krishnan, S. , Cardona, A. , Romero, M. , Quiroz, L. , Rodriguez, G. , Rueda, M. J. , de la Parra, F. , Morón, S. , Green, W. , Bayona, G. , Montes, C. , Quintero, O. , Ramirez, R. , … Vervoort, J. (2010). Effects of rapid global warming at the Paleocene‐Eocene boundary on Neotropical vegetation. Science, 330 , 957–961. 10.1126/science.1193833 21071667
Jordan, G. J. (1997). Evidence of Pleistocene plant extinction and diversity from regatta point, western Tasmania, Australia. Botanical Journal of the Linnean Society, 123 , 45–71. 10.1111/j.1095-8339.1997.tb01404.x
Jordan, G. J. , Bannister, J. M. , Mildenhall, D. C. , Zetter, R. , & Lee, D. E. (2010). Fossil Ericaceae from New Zealand: Deconstructing the use of fossil evidence in historical biogeography. American Journal of Botany, 97 , 59–70.21622367
Kennedy, E. M. , Spicer, R. A. , & Rees, P. M. (2002). Quantitative palaeoclimate estimates from late cretaceous and Paleocene leaf floras in the northwest of the South Island, New Zealand. Palaeogeography, Palaeoclimatology, Palaeoecology, 184 , 321–345. 10.1016/S0031-0182(02)00261-4
Lee, D. , Lee, W. , Jordan, G. , & Barreda, V. (2016). The Cenozoic history of New Zealand temperate rainforests: Comparisons with southern Australia and South America. New Zealand Journal of Botany, 54 , 100–127. 10.1080/0028825X.2016.1144623
Lee, D. E. , Bannister, J. M. , Raine, J. I. , & Conran, J. G. (2010). Euphorbiaceae: Acalyphoideae fossils from early Miocene New Zealand: Mallotus‐Macaranga leaves, fruits, and inflorescence with in situ Nyssapollenites endobalteus pollen. Review of Palaeobotany and Palynology, 163 , 127–138.
Lee, D. E. , Conran, J. G. , Lindqvist, J. K. , Bannister, J. M. , & Mildenhall, D. C. (2012). New Zealand Eocene, Oligocene and Miocene macrofossil and pollen records and modern plant distributions in the southern hemisphere. Botanical Review, 78 , 235–260. 10.1007/s12229-012-9102-7
Lee, D. E. , Lee, W. G. , & Mortimer, N. (2001). Where and why have all the flowers gone? Depletion and turnover in the New Zealand Cenozoic angiosperm flora in relation to palaeogeography and climate. Australian Journal of Botany, 49 , 341–356. 10.1071/bt00031
Liebrand, D. , de Bakker, A. T. M. , Beddow, H. M. , Wilson, P. A. , Bohaty, S. M. , Ruessink, G. , Pälike, H. , Batenburg, S. J. , Hilgen, F. J. , Hodell, D. A. , Huck, C. E. , Kroon, D. , Raffi, I. , Saes, M. J. M. , van Dijk, A. E. , & Lourens, L. J. (2017). Evolution of the early Antarctic ice ages. Proceedings of the National Academy of Sciences of the United States of America, 114 , 3867–3872. 10.1073/pnas.1615440114 28348211
Lu, W.‐X. , Wang, Z.‐Z. , Hu, X.‐Y. , & Rao, G.‐Y. (2023). Echoes of the past: Niche evolution, range dynamics, and their coupling shape the distribution of species in the Chrysanthemum zawadskii species complex. Frontiers in Ecology and Evolution, 11 , 1–14.
McElwain, J. C. , & Punyasena, S. W. (2007). Mass extinction events and the plant fossil record. Trends in Ecology & Evolution, 22 , 548–557. 10.1016/j.tree.2007.09.003 17919771
McGill, B. J. (2010). Matters of scale. Science, 328 , 575–576. 10.1126/science.1188528 20431001
McGlone, M. , Buitenwerf, R. , & Richardson, S. (2016). The formation of the oceanic temperate forests of New Zealand. New Zealand Journal of Botany, 54 , 128–155. 10.1080/0028825X.2016.1158196
McGlone, M. S. , Duncan, R. P. , & Heenan, P. B. (2001). Endemism, species selection and the origin and distribution of the vascular plant Flora of New Zealand. Journal of Biogeography, 28 , 199–216.
Mildenhall, D. C. (1980). New Zealand late cretaceous and cenozoic plant biogeography: A contribution. Palaeogeography, Palaeoclimatology, Palaeoecology, 31 , 197–233. 10.1016/0031-0182(80)90019-X
Mildenhall, D. C. (1989). Summary of the age and paleoecology of the Miocene Manuherikia Group, Central‐Otago, New‐Zealand. Journal of the Royal Society of New Zealand, 19 , 19–29.
Mildenhall, D. C. , Hollis, C. J. , & Naish, T. R. (2004). Orbitally‐influenced vegetation record of the mid‐Pleistocene climate transition, offshore eastern New Zealand (ODP leg 181, site 1123). Marine Geology, Cenozoic Oceanographic Evolution of the Southwest Pacific Gateway, ODP Leg, 181 (205 ), 87–111. 10.1016/S0025-3227(04)00019-2
Mildenhall, D. C. , Stokes, S. , & Nelson, C. S. (1992). Palynology, age, and paleoenvironments of Carbonaceous facies in the Kaihu Group: (late Pliocene‐Pleistocene), northern North Island. DSIR Geology & Geophysics.
Naish, T. , Kamp, P. J. J. , & Pillans, B. (1997). Recurring global sea‐level changes recorded in shelf deposits near the G/M polarity transition, Wanganui Basin, New Zealand: Implications for redefining the Pliocene‐Pleistocene boundary. Quaternary International, 40 , 61–71. 10.1016/S1040-6182(96)00062-6
Nelson, C. S. , & Cooke, P. J. (2001). History of oceanic front development in the New Zealand sector of the Southern Ocean during the Cenozoic—A synthesis. New Zealand Journal of Geology and Geophysics, 44 , 535–553. 10.1080/00288306.2001.9514954
Olson, S. L. , & Wingate, D. B. (2000). Two new species of flightless rails (Aves: Rallidae) from the Middle Pleistocene “crane fauna” of Bermuda .
Pancost, R. D. , Taylor, K. W. R. , Inglis, G. N. , Kennedy, E. M. , Handley, L. , Hollis, C. J. , Crouch, E. M. , Pross, J. , Huber, M. , Schouten, S. , Pearson, P. N. , Morgans, H. E. G. , & Raine, J. I. (2013). Early Paleogene evolution of terrestrial climate in the SW Pacific, southern New Zealand. Geochemistry, Geophysics, Geosystems, 14 , 5413–5429. 10.1002/2013GC004935
Pastore, M. (2018). Overlapping: A R package for estimating overlapping in empirical distributions. Journal of Open Source Software, 3 , 1023. 10.21105/joss.01023
Pearman, P. B. , Guisan, A. , Broennimann, O. , & Randin, C. F. (2008). Niche dynamics in space and time. Trends in Ecology & Evolution, 23 , 149–158. 10.1016/j.tree.2007.11.005 18289716
Pearman, P. B. , Randin, C. F. , Broennimann, O. , Vittoz, P. , van der Knaap, W. O. , Engler, R. , Lay, G. L. , Zimmermann, N. E. , & Guisan, A. (2008). Prediction of plant species distributions across six millennia. Ecology Letters, 11 , 357–369. 10.1111/j.1461-0248.2007.01150.x 18279357
Pocknall, D. T. (1982). Palynology of late Oligocene Pomahaka Estuarine Bed sediments, Waikoikoi, Southland, New Zealand. New Zealand Journal of Botany, 20 (3 ), 263–287.
Pocknall, D. T. (1991). Palynostratigraphy of the Te Kuiti Group (late Eocene–Oligocene), Waikato Basin, New Zealand. New Zealand Journal of Geology and Geophysics, 34 (4 ), 407–417.
Pole, M. (1993a). Early Miocene flora of the Manuherikia group, New Zealand. 6. Lauraceae. Journal of the Royal Society of New Zealand, 23 , 303–312. 10.1080/03036758.1993.10721228
Pole, M. (1993b). Early Miocene Flora of The Manuherikia Group, New‐Zealand. 2. Conifers. Journal of the Royal Society of New Zealand, 22 , 287–302.
Pole, M. (2007). Lauraceae macrofossils and dispersed cuticle from the Miocene of southern New Zealand. Palaeontologia Electronica, 10 , 1–38.
Pole, M. (2008a). The record of Araucariaceae macrofossils in New Zealand. Alcheringa, 32 , 405–426.
Pole, M. (2008b). Dispersed leaf cuticle from the Early Miocene of southern New Zealand. Palaeontologia Electronica, 11 , 1–117.
Prebble, J. G. , Reichgelt, T. , Mildenhall, D. C. , Greenwood, D. R. , Raine, J. I. , Kennedy, E. M. , & Seebeck, H. C. (2017). Terrestrial climate evolution in the Southwest Pacific over the past 30 million years. Earth and Planetary Science Letters, 459 , 136–144. 10.1016/j.epsl.2016.11.006
Prebble, J. G. , Kennedy, E. M. , Reichgelt, T. , Clowes, C. , Womack, T. , Mildenhall, D. C. , Raine, J. I. , & Crouch, E. M. (2021). A 100 million year composite pollen record from New Zealand shows maximum angiosperm abundance delayed until Eocene. Palaeogeography, Palaeoclimatology, Palaeoecology, 566 , 110207. 10.1016/j.palaeo.2020.110207
Prentice, I. C. , Bartlein, P. J. , & Webb, T. (1991). Vegetation and climate change in eastern North America since the last glacial maximum. Ecology, 72 , 2038–2056. 10.2307/1941558
R Core Team . (2022). R: A language and environment for statistical computing. R Foundation for Statistical Computing.
Rangel, T. F. , Edwards, N. R. , Holden, P. B. , Diniz‐Filho, J. A. F. , Gosling, W. D. , Coelho, M. T. P. , Cassemiro, F. A. S. , Rahbek, C. , & Colwell, R. K. (2018). Modeling the ecology and evolution of biodiversity: Biogeographical cradles, museums, and graves. Science, 361 , eaar5452. 10.1126/science.aar5452 30026200
Read, J. , & Hill, R. S. (1985). Dynamics of Nothofagus‐dominated rainforest on mainland Australia and lowland Tasmania. Vegetatio, 63 , 67–78. 10.1007/BF00032607
Reichgelt, T. , Kennedy, E. M. , Conran, J. G. , Mildenhall, D. C. , & Lee, D. E. (2015). The early Miocene paleolake Manuherikia: vegetation heterogeneity and warm‐temperate to subtropical climate in southern New Zealand. Journal of Paleolimnology, 53 (4 ), 349–365. 10.1007/s10933-015-9827-5
Reichgelt, T. , Kennedy, E. M. , Conran, J. G. , Lee, W. G. , & Lee, D. E. (2019). The presence of moisture deficits in Miocene New Zealand. Global and Planetary Change, 172 , 268–277. 10.1016/j.gloplacha.2018.10.013
Reichgelt, T. , Lee, W. G. , & Lee, D. E. (2022). The extinction of Miocene broad‐leaved deciduous Nothofagaceae and loss of seasonal forest biomes in New Zealand. Review of Palaeobotany and Palynology, 307 , 104779. 10.1016/j.revpalbo.2022.104779
Rodríguez‐Sánchez, F. , & Arroyo, J. (2008). Reconstructing the demise of Tethyan plants: Climate‐driven range dynamics of Laurus since the Pliocene. Global Ecology and Biogeography, 17 , 685–695. 10.1111/j.1466-8238.2008.00410.x
Rosenzweig, M. L. (1995). Species diversity in space and time .
Schoener, T. W. (1968). The Anolis lizards of Bimini: Resource partitioning in a complex Fauna. Ecology, 49 , 704–726. 10.2307/1935534
Signor, P. W., III , & Lipps, J. H. (1982). Sampling bias, gradual extinction patterns and catastrophes in the fossil record. In L. T. Silver & P. H. Schultz (Eds.), Geological implications of impacts of large asteroids and comets on the earth. Geological Society of America. 10.1130/SPE190-p291
Silvertown, J. , McConway, K. , Gowing, D. , Dodd, M. , Fay, M. F. , Joseph, J. A. , & Dolphin, K. (2005). Absence of phylogenetic signal in the niche structure of meadow plant communities. Proceedings of the Royal Society B: Biological Sciences, 273 , 39–44. 10.1098/rspb.2005.3288
Strogen, D. P. , Seebeck, H. , Hines, B. R. , Bland, K. J. , & Crampton, J. S. (2023). Palaeogeographic evolution of Zealandia: Mid‐cretaceous to present. New Zealand Journal of Geology and Geophysics, 66 , 528–557. 10.1080/00288306.2022.2115520
Svenning, J.‐C. (2003). Deterministic Plio‐Pleistocene extinctions in the European cool‐temperate tree flora. Ecology Letters, 6 , 646–653. 10.1046/j.1461-0248.2003.00477.x
Svenning, J.‐C. , Eiserhardt, W. L. , Normand, S. , Ordonez, A. , & Sandel, B. (2015). The influence of paleoclimate on present‐day patterns in biodiversity and ecosystems. Annual Review of Ecology, Evolution, and Systematics, 46 , 551–572. 10.1146/annurev-ecolsys-112414-054314
Thomas, C. D. , Cameron, A. , Green, R. E. , Bakkenes, M. , Beaumont, L. J. , Collingham, Y. C. , Erasmus, B. F. N. , de Siqueira, M. F. , Grainger, A. , Hannah, L. , Hughes, L. , Huntley, B. , van Jaarsveld, A. S. , Midgley, G. F. , Miles, L. , Ortega‐Huerta, M. A. , Townsend Peterson, A. , Phillips, O. L. , & Williams, S. E. (2004). Extinction risk from climate change. Nature, 427 (6970 ), 145–148. 10.1038/nature02121 14712274
Thompson, R. S. , Anderson, K. H. , Pelltier, R. T. , Strickland, L. E. , Bartlein, P. J. , & Shafer, S. L. (2012). Quantitative estimation of climatic parameters from vegetation data in North America by the mutual climatic range technique. Quaternary Science Reviews, 51 , 18–39. 10.1016/j.quascirev.2012.07.003
Urban, M. C. (2015). Accelerating extinction risk from climate change. Science, 348 , 571–573. 10.1126/science.aaa4984 25931559
Utescher, T. , Mosbrugger, V. , & Ashraf, A. R. (2000). Terrestrial climate evolution in Northwest Germany over the last 25 million years. PALAIOS, 15 , 430–449. 10.2307/3515514
Wang, Y. , Pineda‐Munoz, S. , & McGuire, J. L. (2023). Plants maintain climate fidelity in the face of dynamic climate change. Proceedings of the National Academy of Sciences of the United States of America, 120 , e2201946119. 10.1073/pnas.2201946119 36745797
Warren, D. L. , Glor, R. E. , & Turelli, M. (2008). Environmental niche equivalency versus conservatism: Quantitative approaches to niche evolution. Evolution, 62 , 2868–2883. 10.1111/j.1558-5646.2008.00482.x 18752605
Winkworth, R. C. , Wagstaff, S. J. , Glenny, D. , & Lockhart, P. J. (2002). Plant dispersal N.E.W.S from New Zealand. Trends in Ecology & Evolution, 17 (11 ), 514–520. 10.1016/s0169-5347(02)02590-9
Winkworth, R. C. , Wagstaff, S. J. , Glenny, D. , & Lockhart, P. J. (2005). Evolution of the New Zealand mountain flora: Origins, diversification and dispersal. Organisms Diversity & Evolution, 5 (3 ), 237–247. 10.1016/j.ode.2004.12.001
Westerhold, T. , Marwan, N. , Drury, A. J. , Liebrand, D. , Agnini, C. , Anagnostou, E. , Barnet, J. S. K. , Bohaty, S. M. , De Vleeschouwer, D. , Florindo, F. , Frederichs, T. , Hodell, D. A. , Holbourn, A. E. , Kroon, D. , Lauretano, V. , Littler, K. , Lourens, L. J. , Lyle, M. , Pälike, H. , … Zachos, J. C. (2020). An astronomically dated record of Earth's climate and its predictability over the last 66 million years. Science, 369 , 1383–1387. 10.1126/science.aba6853 32913105
Willis, K. J. , & Whittaker, R. J. (2002). Species diversity—Scale matters. Science, 295 , 1245–1248. 10.1126/science.1067335 11847328
Wing, S. L. (2004). Mass extinctions in plant evolution. In P. D. Taylor (Ed.), Extinctions in the history of life (pp. 61–98). Cambridge University Press. 10.1017/CBO9780511607370.004
You, J. , Qin, X. , Ranjitkar, S. , Lougheed, S. C. , Wang, M. , Zhou, W. , Ouyang, D. , Zhou, Y. , Xu, J. , Zhang, W. , Wang, Y. , Yang, J. , & Song, Z. (2018). Response to climate change of montane herbaceous plants in the genus Rhodiola predicted by ecological niche modelling. Scientific Reports, 8 , 5879. 10.1038/s41598-018-24360-9 29651147
