• Energy Research
• SE close
• AU close
• FI close
• English close

Although GPS coordinates for current populations are not included due to the potential threat of poaching, the climate variables for each species are provided. The records for extant gecko and skinks mainly came from the New Zealand's Department of Conervation Herpetofauna Database. After updating the taxonomy and cleaning the data to reflect the taxonomy as at 2019 of 43 geckos speceis recognised across seven genera and 61 species in genus, we then thinned the occurrence records at a 1 km resolution for all species then predicted distributions for those with > 15 records using species distribution models. The climate variables for each species were selected among annual mean temperature (bio1), maximum temperature of the warmest month (bio5), minimum temperature of the coldest month (bio6), mean temperature of driest quarter (bio9), mean temperature of wettest quarter (bio10), and precipitation of the driest quarter (bio17). To reduce multicollinearity in species distribution models for each species, we only retained climate variables with a variable inflation factor < 10. The climate variables were from the CHELSA database (https://chelsa-climate.org/), which can be freely downloaded for current and future scenarios. We also provide MCC tree files for the geckos and skinks. The phylogenetic trees have been constructed for NZ geckos by (Nielsen et al., 2011) and for NZ skinks by (Chapple et al., 2009). For geckos we used a subset of the sequences used by Nielsen et al. (2011) for four genes, two nuclear (RAG 1, PDC) and two mitochondrial (16S, ND2 along with flanking tRNA sequences). For skinks, we used sequences from Chapple et al. (2009) for one nuclear (RAG 1) and five mitochondrial (ND2, ND4, Cyt b, 12S and 16S) genes, and additional ND2 sequences for taxa not included in the original phylogeny (Chapple et al., 2011, p. 201). In total we used sequences for all recognised extant taxa (Hitchmough et al., 2016) as at 2019 except for three species of skink (O. aff. inconspicuum “Okuru”, O. robinsoni, and O. aff. inconspicuum “North Otago”) and two species of gecko (M. “Cupola” and W. “Kaikouras”) for which genetic data were not available. Aim: The primary drivers of species and population extirpations have been habitat loss, overexploitation, and invasive species, but human-mediated climate change is expected to be a major driver in future. To minimise biodiversity loss, conservation managers should identify species vulnerable to climate change and prioritise their protection. Here, we estimate climatic suitability for two speciose taxonomic groups, then use phylogenetic analyses to assess vulnerability to climate change. Location: Aotearoa New Zealand (NZ) Taxa: NZ lizards: diplodactylid geckos and eugongylinae skinks Methods: We built correlative species distribution models (SDMs) for NZ geckos and skinks to estimate climatic suitability under current climate and 2070 future-climate scenarios. We then used Bayesian phylogenetic mixed models (BPMMs) to assess vulnerability for both groups with predictor variables for life history traits (body size and activity phase) and current distribution (elevation and latitude). We explored two scenarios: an unlimited dispersal scenario, where projections track climate, and a no-dispersal scenario, where projections are restricted to areas currently identified as suitable. Results: SDMs projected vulnerability to climate change for most modelled lizards. For species’ ranges projected to decline in climatically suitable areas, average decreases were between 42–45% for geckos and 33–91% for skinks, although area did increase or remain stable for a minority of species. For the no-dispersal scenario, the average decrease for geckos was 37–52% and for skinks was 33–52%. Our BPMMs showed phylogenetic signal in climate change vulnerability for both groups, with elevation increasing vulnerability for geckos, and body size reducing vulnerability for skinks. Main conclusions: NZ lizards showed variable vulnerability to climate change, with most species’ ranges predicted to decrease. For species whose suitable climatic space is projected to disappear from within their current range, managed relocation could be considered to establish populations in regions that will be suitable under future climates.

This database includes more than 100 decarbonisation innovations in Paper, Plastic, Steel and Meat & Dairy sectors, across their value chains, as well as in Finance. For each innovation there is a description, information about its contribution to decarbonisation, actors and collaborators involved, sources of funding, drivers, (co)benefits and disadvantages. More information on the method for selecting innovations for the database is available here. The database was created as part of REINVENT – a Horizon 2020 research project funded by the European Commission (grant agreement 730053). REINVENT involves five research institutions from four countries: Lund University (Sweden), Durham University (United Kingdom), Wuppertal Institute (Germany), PBL Netherlands Environmental Assessment Agency (the Netherlands) and Utrecht University (the Netherlands). More information can be found on our website: www.reinvent-project.eu.

The atmosphere concentration of CO2 is steadily increasing and causing climate change. To achieve the Paris 1.5 or 2 oC target, negative emissions technologies must be deployed in addition to reducing carbon emissions. The ocean is a large carbon sink but the potential of marine primary producers to contribute to carbon neutrality remains unclear. Here we review the alterations to carbon capture and sequestration of marine primary producers (including traditional ‘blue carbon’ plants, microalgae, and macroalgae) in the Anthropocene, and, for the first time, assess and compare the potential of various marine primary producers to carbon neutrality and climate change mitigation via biogeoengineering approaches. The contributions of marine primary producers to carbon sequestration have been decreasing in the Anthropocene due to the decrease in biomass driven by direct anthropogenic activities and climate change. The potential of blue carbon plants (mangroves, saltmarshes, and seagrasses) is limited by the available areas for their revegetation. Microalgae appear to have a large potential due to their ubiquity but how to enhance their carbon sequestration efficiency is very complex and uncertain. On the other hand, macroalgae can play an essential role in mitigating climate change through extensive offshore cultivation due to higher carbon sequestration capacity and substantial available areas. This approach seems both technically and economically feasible due to the development of offshore aquaculture and a well-established market for macroalgal products. Synthesis and applications: This paper provides new insights and suggests promising directions for utilizing marine primary producers to achieve the Paris temperature target. We propose that macroalgae cultivation can play an essential role in attaining carbon neutrality and climate change mitigation, although its ecological impacts need to be assessed further. To calculate the parameters presented in Table 1, the relevant keywords "mangroves, salt marshes, macroalgae, microalgae, global area, net primary productivity, CO2 sequestration" were searched through the ISI Web of Science and Google Scholar in July 2021. Recent data published after 2010 were collected and used since area and productivity of plants change with decade. For data with limited availability, such as net primary productivity (NPP) of seagrasses and global area and NPP of wild macroalgae, data collection was extended back to 1980. Total NPP and CO2 sequestration for mangroves, salt marshes, seagrasses and wild macroalgae were obtained by the multiplication of area and NPP/CO2 sequestration density and subjected to error propagation analysis. Data were expressed as means ± standard error.

Heavy trucks contribute significantly to climate change, and in 2020 were responsible for 7% of total Swedish GHG emissions and 5% of total global CO2 emissions. Here we study the full lifecycle of cargo trucks powered by different energy pathways, comparing their biomass feedstock use, primary energy use, net biogenic and fossil CO2 emission, and cumulative radiative forcing. We analyse battery electric trucks with bioelectricity from standalone or combined heat and power (CHP) plants, and pathways where bioelectricity is integrated with wind and solar electricity. We analyse trucks operated on fossil diesel fuel and on dimethyl ether (DME). All energy pathways are analysed with and without carbon capture and storage (CCS). Bioelectricity and DME are produced from forest harvest residues. Forest biomass is a limited resource, so in a scenario analysis we allocate a fixed amount of biomass to power Swedish truck transport. Battery lifespan and chemistry, the technology level of energy supply, and the biomass source and transport distance are all varied to understand how sensitive the results are to these parameters. The scenario spans 100 years into the future. We find that pathways using electricity to power battery electric trucks have much lower climate impacts and primary energy use, compared to diesel and DME based pathways. The pathways using bioelectricity with CCS result in negative emissions leading to global cooling of the earth. The pathways using diesel and DME have significant and very similar climate impact, even with CCS. The robust results show that truck electrification and increased renewable electricity production is a much better strategy to reduce the climate impact of cargo transport and much more primary energy efficient than the adoption of DME trucks. This climate impact analysis includes all fossil and net biogenic CO2 emissions as well as the timing of these emissions. Considering only fossil emissions is incomplete and could be misleading. This dataset contains data on 4 metrics (primary energy use, biomass feedstock use, cumulative CO2 emissions, and cumulative radiative forcing) resulting from scenario modeling of cargo truck use in Sweden powered by different energy pathways. The energy pathways include battery electric trucks powered by bioelectricity, solar photovoltaic electricity and wind electricity, and internal combustion trucks powered by fossil diesel and dimethyl ether. The scenario spans 100 years into the future. The Excel sheet "tables" contains input data for the scenario modeling, with sources listed where applicable. The remaining sheets contains the modeled results and generated figures that are also a published in the associated article Sathre & Gustavsson (2023). Refer to the method description and reference list in the included documentation files for details. Tunga lastbilar bidrar kraftigt till klimatförändringarna och stod 2020 för 7% av de totala svenska växthusgasutsläppen och 5% av de totala globala CO2-utsläppen. Här studerar vi hela livscykeln för lastbilar som drivs av olika energivägar, jämför deras användning av biomassaråvaror, primär energianvändning, biogena och fossila CO2-utsläpp netto och kumulativ strålningstvingning. Vi analyserar batterielektriska lastbilar med bioel från fristående eller kraftvärmeverk och vägar där bioel integreras med vind- och solkraft. Vi analyserar lastbilar som drivs med fossilt dieselbränsle och med dimetyleter (DME). Alla energivägar analyseras med och utan avskiljning och lagring av koldioxid (CCS). Bioelektricitet och DME produceras av skogsavverkningsrester. Skogsbiomassa är en begränsad resurs, så i en scenarioanalys avsätter vi en fast mängd biomassa för att driva svenska lastbilstransporter. Batteriets livslängd och kemi, tekniknivån för energiförsörjning och biomassakällan och transportavståndet varierar alla för att förstå hur känsliga resultaten är för dessa parametrar. Scenariot sträcker sig 100 år in i framtiden. Vi finner att vägar som använder el för att driva batterielektriska lastbilar har mycket lägre klimatpåverkan och primär energianvändning, jämfört med diesel- och DME-baserade vägar. De vägar som använder bioelektricitet med CCS resulterar i negativa utsläpp som leder till global kylning av jorden. Vägarna med diesel och DME har betydande och mycket liknande klimatpåverkan, även med CCS. De robusta resultaten visar att elektrifiering av lastbilar och ökad förnybar elproduktion är en mycket bättre strategi för att minska godstransporternas klimatpåverkan än införandet av DME-lastbilar, och mycket mer primärenergieffektiv. Denna klimatkonsekvensanalys omfattar alla fossila och biogena CO2-utsläpp samt tidpunkten för dessa utsläpp. Att bara ta hänsyn till fossila utsläpp är ofullständigt och kan vara missvisande. Detta dataset innehåller data om 4 mätvärden (primär energianvändning, biomassaråvara, kumulativa CO2-utsläpp och kumulativ strålkraftspåverkan) som härrör från scenariomodellering av lastbilsanvändning i Sverige som drivs av olika energivägar. Energivägarna inkluderar batterielektriska lastbilar som drivs av bioelektricitet, solcellselektricitet och vindkraft samt förbränningsbilar som drivs av fossil diesel och dimetyleter. Scenariot sträcker sig 100 år in i framtiden. På arket "tables" i Excelfilen återfinns den indata som använts i modelleringen med angivna källor där detta är tillämpligt. Övriga ark innehåller resultat samt figurer som också publiceras i den samhörande artikeln Sathre & Gustavsson (2023). Se metodbeskrivning samt referenslista i tillhörande dokumentationsfiler för detaljer.

# Data from: Eocene Shark Teeth from Peninsular Antarctica: Windows to Habitat Use and Paleoceanography. [https://doi.org/10.5061/dryad.qz612jmq2](https://doi.org/10.5061/dryad.qz612jmq2) The repository folder includes scripts and spreadsheets for phosphate oxygen stable isotope (δ18Op) analysis measured from shark tooth biogenic apatite collected from the Eocene deposits of the La Meseta and Submeseta formations (West Antarctica, Seymour Island). It also contains Fourier-Transform Infrared Spectroscopy (FTIR) analysis, a Bayesian model for temperature estimates, and model output extraction scripts from the iCESM simulation for the Early Eocene (Zhu et al., 2020). Scripts and data are stored in specific folders on the type of analysis. All scripts are in R or Python language. **Usage notes** **1 "iCESM modeling scripts" directory** The folder includes scripts in Jupiter Notebook format for extracting and plotting iCESM seawater outputs for the Eocene. The folder includes two files: 1) “d18Ow Analysis Script.ipynb” - This is a Python script primarily using the XArray library, to import iCESM output from Zhu et al. (2020), calculating δ18Ow, and reorganizing the output into monthly time intervals along 25 m and 115 m depth slices, while also averaging output down to these depths; 2) “NetCDF Plotting.ipynb” - this is a Python script primarily using the XArray, Matplotlib, and Cartopy libraries. The script writes a single callable function that creates Matplotlib contour plots from iCESM history output. Variables include temperature, salinity, ideal age, oxygen isotopes, and neodymium isotopes, and map projections include Plate Carree, Mollweide, and orthographic (centering on the Drake Passage). Options are built to enable scale normalization or to set maximum and minimum values for data and select colormaps from a predefined selection of Matplotlib’s “Spectral”, “Viridis”, “Coolwarm”, “GNUplot2”, “PiYG”, “RdYlBu”, and “RdYlGn”. For further questions on model output scripts, please email Adam Aleksinski at [aaleksin@purdue.edu](https://datadryad.org/stash/dataset/doi:10.5061/aaleksin@purdue.edu). **2 "d18O data and maps" directory** The folder includes δ18Op of shark tooth bioapatite and other datasets to interpret shark paleoecology. These datasets include: · δ18Op of shark tooth bioapatite (“shark FEST d18Op.csv”). Isotope measurements were run at the Stable Isotope Ecosystem Laboratory of (SIELO) University of California, Merced (California, USA). · Reference silver phosphate material δ18Op for analytical accuracy and precision (“TCEA reference materials.csv"). Isotope measurements were run at the Stable Isotope Ecosystem Laboratory of (SIELO) University of California, Merced (California, USA). · Bulk and serially sampled δ18Oc data of co-occurring bivalves (Ivany et al., 2008; Judd et al., 2019) (“Ivany et al. 2008_bulk.csv” and “Judd et al., 2019_serial sampling.csv"). · iCESM model temperature and δ18Ow outputs at 3x and 6x pre-industrial CO2 levels for the Early Eocene (Zhu et al., 2020) (“SpinupX3_25m_Mean_Monthly.nc”, “SpinupX6_25m_Mean_Monthly.nc.”, and “CA_x3CO2.csv”). Simulations are integrated from the surface to 25 m. · δ18O values of invertebrate species published in Longinelli (1965) and Longinelli & Nuti (1973), used to convert bulk δ18Oc (V-SMOW) data of bivalves into δ18Op (V-SMOW) values after δ18Oc (V-PDB) - δ18Oc (V-SMOW) conversion found in Kim et al. (2015) (“d18O carbonate and phosphate references.csv”). · R script for data analysis ("d18O data and maps.Rmd”). The script provides annotation through libraries, instrumental accuracy and precision tests, tables, statistical analysis, figures, and model output extractions. . ("TELM_diversity.csv") displays diversity trends of chondrichthyans across TELMs in one of the main figures of the manuscript. **2.1 Dataset description** **shark FEST d18Op.csv** · *Sample_ID*: Identification number of tooth specimens. · *Other_ID*: Temporary identification number of tooth specimens. · *Taxon*: Species assigned to shark tooth specimens. · *TELM*: Stratigraphic units of La Meseta (TELM 2-5; ~45 to ~37 Ma) and Submeseta formations (TELMs 6 and 7; ~37 to ~34 Ma) (Amenábar et al., 2020; Douglas et al., 2014; Montes et al., 2013). · *d18Op*: Mean δ18Op values of silver phosphate crystals precipitated from shark tooth bioapatite. Specimens were run in triplicates, corrected, and standardized on the V-SMOW scale. · *sd*: Standard deviation of silver phosphate triplicate samples per specimen. · *Protocol*: Silver phosphate protocols used to precipitate crystals from shark tooth bioapatite. We adopted the Rapid UC (“UC_Rapid”) and the SPORA (“SPORA”) protocols after Mine et al. and (2017) Larocca Conte et al. (2024) based on the tooth specimen size and sampling strategy. Descriptions of the methods are included in the main manuscript. · *Environment*: Inferred shark habitat based on taxonomy classified as benthic or pelagic environment. · *Collection*: Institutional abbreviations of museum collections from which shark tooth specimens are housed. NRM-PZ is the abbreviation for the Swedish Natural History Museum (Stockholm, Sweden), PRI is the abbreviation for the Paleontological Research Institute (Ithaca, New York, United States), and UCMP is the University of California Museum of Paleontology (Berkeley, California, United States). **TCEA reference materials.csv** · *Identifier_1*: unique identifier number per sample. · *sample*: reference silver phosphate materials (USGS 80 and USGS 81). · *amount*: weight of samples in mg. · *Area 28*: peak area of mass 28 (12C16O). · *Area 30*: peak area of mass 30 (12C18O). · *d18O_corrected*: corrected δ18Op value of reference materials following drift correction, linearity correction, and 2-point calibration to report values on the V-SMOW scale. **Ivany et al. 2008_bulk.csv** · *Telm*: Stratigraphic units of La Meseta (TELM 2-5; ~45 to ~37 Ma) and Submeseta formations (TELMs 6 and 7; ~37 to ~34 Ma) (Amenábar et al., 2020; Douglas et al., 2014; Montes et al., 2013). · *Locality*: Locality code from which bivalves were collected. · *Genus*: Genera of bivalves. Specimens are assigned to *Cucullaea* and *Eurhomalea* genera. · *Line*: Sampling areas of specimens. The sampling strategy is described in Ivany et al. (2008). · *d13C*: δ13C values of specimens from sampled lines. Values are reported in the V-PDB scale. · *d18Oc_PDB*: δ18Oc values of specimens from sampled lines. Values are reported in the V-PDB scale. **Judd et al., 2019_serial sampling.csv** · *Horizon:* horizons of the TELM 5 unit (La Meseta Formation) from which bivalves were collected. Horizon 1 is stratigraphically the lowest, while horizon 4 is the highest (Judd et al., 2019). · *ID*: Identification number of specimens. · *Latitude*: Geographic coordinate where bivalve specimens were collected. · *Longitude*: Geographic coordinate where bivalve specimens were collected. · *Surface sampled*: Specific sampling area, indicating whether sampling occurred in the interior or exterior portion of shells. · *distance*: The distance from the umbo in mm from which sampling occurred along a single shell. · *d18Oc_PDB*: δ18Oc values of specimens from sampled areas of shells. Values are reported on the V-PDB scale. **SpinupX3_25m_Mean_Monthly.nc** See section 1 ("iCESM modeling scripts" directory, “d18Ow Analysis Script.ipynb” script) for a full description of the iCESM model output extraction. **SpinupX6_25m_Mean_Monthly.nc** See section 1 ("iCESM modeling scripts" directory, “d18Ow Analysis Script.ipynb” script) for a full description of the iCESM model output extraction. **CA_x3CO2.csv** · *lat*: Geographic coordinate where temperature and δ18Ow model values are extracted from the iCESM simulation scaled at 3x preindustrial CO2 levels (values averaged within a seawater column depth of 25 m). · *long*: Geographic coordinate where temperature and δ18Ow model values are extracted from the iCESM simulation scaled at 3x preindustrial CO2 levels (values averaged within a seawater column depth of 25 m). · *T_mean*: Simulated seawater temperature values in °C. · *d18Ow*: Simulated seawater δ18Ow values (V-SMOW). · *d18Op*: Simulated seawater δ18Op values (V-SMOW). Values were calculated by using seawater temperature and δ18Ow arrays following the paleothermometer equation after Lécuyer et al. (2013). **d18O carbonate and phosphate references.csv** · *species*: Species of invertebrate taxa. · *type*: Specimen type, including barnacles, brachiopods, crabs, and mollusks. · *depth*: Depth of seawater column where specimens were collected, reported in meters below sea level when specified. · *d18Op*: δ18Op values of invertebrate specimens (V-SMOW). · *d18Oc_PDB*: δ18Oc values of invertebrate specimens (V-PDB). · *Reference*: Citations from which data were taken to build the dataset (Longinelli, 1965; Longinelli & Nuti, 1973). **TELM diversity.csv** · *genus:* genera of sharks and rays compiled from literature (Engelbrecht et al., 2016a, 2016b, 2017a, 2017b, 2019; Kriwet, 2005; Kriwet et al., 2016; Long, 1992; Marramá et al., 2018). · *species*: species of sharks and rays compiled from literature (Engelbrecht et al., 2016a, 2016b, 2017a, 2017b, 2019; Kriwet, 2005; Kriwet et al., 2016; Long, 1992; Marramá et al., 2018). · *Environment*: Inferred shark habitat based on taxonomy classified as benthic or pelagic environment. · *TELM*: Stratigraphic units of La Meseta (TELM 1-5; ~44 to ~37 Ma) and Submeseta formations (TELMs 6 and 7; ~37 to ~34 Ma) (Amenábar et al., 2020; Douglas et al., 2014; Montes et al., 2013). **3 “FTIR data” directory** The folder includes FTIR acquisitions and data analysis scripts on reference materials and shark tooth bioapatite for quality checks to test diagenesis effects on δ18Op of sharks. The folder includes: · The R project file “apatite_ftir.Rproj”. This project file navigates through scripts for raw data processing and data analysis. The background of the raw data was processed following custom R functions from Trayler et al. (2023; [https://github.com/robintrayler/collagen_demineralization](https://github.com/robintrayler/collagen_demineralization)). · The “.Rproj.user” folder includes project-specific temporary files (e.g. auto-saved source documents, window-state, etc.) stored by the R project file “apatite_ftir.Rproj”. The folder may be hidden depending on directory view options. · The “raw data” directory stores spectra acquisitions as .dpt files. Spectra files are stored in the folders “apatite” and “calcite” based on the material type. Spectra were obtained in the 400 – 4000 cm⁻¹ range using a Bruker Vertex 70 Far-Infrared in ATR located at the Nuclear Magnetic Resonance Facility at the University of California Merced (California, USA). · The “processed” directory includes processed spectra stored as .csv files (“apatite_data.csv” and “calcite_data.csv”) following the background correction (Trayler et al., 2023) and processed infrared data from Larocca Conte et al. (2024) (“Larocca Conte et al._SPORA_apatite_data.csv”) from which the NIST SRM 120c spectrum was filtered. Infrared spectra data in “Larocca Conte et al._SPORA_apatite_data.csv” were obtained and corrected following the same methodologies mentioned above. · The “R” directory includes R scripts of customized source functions for background correction (Trayler et al., 2023; inspect the "functions" directory and the R script "0_process_data.R") and data analysis (“data_analysis.R”). The scripts provide annotation through libraries and functions used for data processing and analysis. · Additional datasets. The “data_FTIR_d18O.csv” includes infrared data and δ18Op values of specimens, while the “Grunenwald et al., 2014_CO3.csv” is the dataset after Grunenwald et al. (2014) used to predict carbonate content from the materials featured in this work. **3.1 Dataset description** Spreadsheets included in the “processed” directory The datasets “apatite_data.csv”, “calcite_data.csv”, and “Larocca Conte et al._SPORA_apatite_data.csv” are structured with the following variables: · *wavenumber*: infrared wavenumber in cm-1. · *absorbance*: infrared absorbance value. · *file_name:* .dpt file name from which infrared wavenumber and absorbance values were obtained following the background correction. **data_FTIR_d18O.csv** · *file_name:* .dpt file name from which infrared wavenumber and absorbance values were obtained following the background correction. · *v4PO4_565_wavenumber*: Wavenumber of maximum infrared absorbance around the first νPO4 band, usually at 565 cm-1. · *v4PO4_565*: Peak absorbance value of the first ν4PO4 band (~565 cm-1). · *v4PO4_valley_wavenumber*: Wavenumber of valley between ν4PO4 bands. · *v4PO4_valley*: Absorbance value of the valley between ν4PO4 bands. · *v4PO4_603_wavenumber*: Wavenumber of maximum infrared absorbance around the second ν4PO4 band, usually at 603 cm-1. · *v4PO4_603*: Peak absorbance value of the second ν4PO4 band (~603 cm-1). · *CI*: Crystallinity index calculated after equation provided in (Shemesh, 1990) as (*v4PO4_565* + *v4PO4_603* / *v4PO4_valley*) (i.e., the sum of peak absorbance of νPO4 bands divided by the absorbance value of the valley between peaks). · *material*: Material type of samples (i.e., standard material, enameloid, dentin sampled from the crown or root area of shark teeth, and enameloid mixed with dentin). · *AUC_v3PO4*: Area under the curve of the ν3PO4 and ν1PO4 bands where maximum absorbance is at ~1025 cm-1 and ~960 cm-1, respectively. · *AUC_v3CO3*: Area under the curves of Type-A and Type-B carbonate bands having maximum infrared absorbance at ~1410 (Type-B), ~1456 (Type-B), and ~1545 cm-1 (Type-A). · *v3CO3_v3PO4_ratio*: Ratio between area under the curves of carbonate and phosphate bands (i.e., *AUC_v3CO3* / *AUC_v3PO4*). · *CO3_wt*: Estimated mean carbonate content following the equation in Grunenwald et al. (2014) (i.e. *CO3_wt* = 28.4793 (±1.4803) *v3CO3_v3PO4_ratio* + 0.1808(±0.2710); R2 = 0.985). · *CO3_wt_sd*: Standard deviation of estimated carbonate content calculated by propagating the error around coefficients provided in the Grunenwald et al. (2014) equation (see full equation in *CO3_wt*). · *Taxon*: Species assigned to shark tooth specimens. · *TELM*: Stratigraphic units of La Meseta (TELM 2-5; ~45 to ~37 Ma) and Submeseta formations (TELMs 6 and 7; ~37 to ~34 Ma) (Amenábar et al., 2020; Douglas et al., 2014; Montes et al., 2013). · *d18Op*: Mean δ18Op values of silver phosphate crystals precipitated from shark tooth bioapatite. Specimens were run in triplicates, corrected, and standardized on the V-SMOW scale. · *sd*: Standard deviation of silver phosphate triplicate samples per specimen. · *Collection*: Institutional abbreviations of museum collections where shark tooth specimens are housed. Infrared spectra were obtained from a selected subset of tooth specimens in the care of the Swedish Natural History Museum (NRM-PZ; Stockholm, Sweden). **Grunenwald et al., 2014_CO3.csv** · *sample*: Sample code. · *material*: Material type of samples (i.e., standard material, bone, and enamel). · *v3CO3*: Area under the curves of Type-A and Type-B carbonate bands having maximum infrared absorbance at ~1410 (Type-B), ~1456 (Type-B), and ~1545 cm-1 (Type-A). · *v3PO4*: *AUC_v3PO4*: Area under the curve of the ν3PO4 and ν1PO4 bands where maximum absorbance is at ~1025 cm-1 and ~960 cm-1, respectively. · *v3CO3_v3PO4_ratio*: *v3CO3_v3PO4_ratio*: Ratio between area under the curves of carbonate and phosphate bands (i.e., *v3CO3* /*v3PO4*). · *CO3_wt*: Carbonate content measured via CO2 coulometry. Further details about the analytical measurements are found in Grunenwald et al. (2014). **4 “Bayes_FEST_Temperautre Estimates” directory** The folder includes the Bayesian approach used to estimate posterior seawater temperature, δ18Ow values from δ18Op of sharks bioapatite using a Bayesian approach modified after Griffiths et al. (2023). The original scripts used in Griffiths et al. (2023) are reposited here: [https://github.com/robintrayler/bayesian_phosphate](https://github.com/robintrayler/bayesian_phosphate). The directory includes: · The R project file “Bayes_FEST.Rproj”. This project file navigates through scripts for raw data analysis. · The “.Rproj.user” folder includes project-specific temporary files (e.g. auto-saved source documents, window-state, etc.) stored by the R project file “Bayes_FEST.Rproj”. The folder may be hidden depending on directory view options. · The “data” folder includes the spreadsheets for modeled seawater temperature and δ18Ow values (“CA_x3CO2.csv”) and δ18Op values of shark tooth bioapatite (“shark FEST d18Op.csv”) used as prior information for the Bayesian model. We refer to section 2.1 for the full description of spreadsheets. · The “R” folder includes customized functions for the Bayesian model stored in the “functions” directory and the script for data analysis (“01_model_sharks.R”). The script includes a comparison of paleothermometer equations after Kolodny et al. (1983), Lécuyer et al. (2013), Longinelli & Nuti (1973), and (Pucéat et al. (2010) using the bulk δ18Op shark tooth bioapatite, simulated seawater temperature and δ18Ow values as prior inputs. While all paleothermometers estimate similar posterior bulk δ18Op close to empirical values, temperature estimates using the Pucéat et al. (2010) method are often the highest, generating estimates ~8°C higher than other equations. We therefore used the Lécuyer et al. (2013) paleothermomether for temperature estimates using δ18Op of shark bioapatite grouped by taxa because it: 1\) Provides consistent posterior temperature estimates relative to other equations (Longinelli & Nuti, 1973, Kolodny et al., 1983). 2\) provides temperature values from fish tooth specimens consistent with estimates of co-existing bivalves or brachiopod carbonate shells. The script provides annotation through libraries, statistical analysis, figures, and tables. **4 Software** **4.1 R** R and R Studio (R Development Core Team, 2024; RStudio Team, 2024) are required to run scripts included in the "d18O data and maps", “FTIR data”, and “Bayes_FEST_Temperautre Estimates” directories, which were created using versions 4.4.1 and 2024.04.02, respectively. Install the following libraries before running scripts: “cowplot” (Wilke, 2024), “colorspace” (Zeileis et al., 2020), “DescTools” (Signorell, 2024), “lattice” (Sarkar, 2008), “flextable” (Gohel & Skintzos, 2024), “ggh4x” (van den Brand, 2024), “ggnewscale” (Campitelli, 2024), “ggpubr” (Kassambara, 2023a), “ggspatial” (Dunnington, 2023), “ggstance” (Henry et al., 2024), “ggstar” (Xu, 2022), “greekLetters” (Kévin Allan Sales Rodrigues, 2023), “gridExtra” (Auguie, 2017), “mapdata” (code by Richard A. Becker & version by Ray Brownrigg., 2022); “mapproj” (for R by Ray Brownrigg et al., 2023), “maps” (code by Richard A. Becker et al., 2023), “ncdf4” (Pierce, 2023), “oce” (Kelley & Richards, 2023), “rasterVis” (Oscar Perpiñán & Robert Hijmans, 2023), “RColorBrewer” (Neuwirth, 2022), “rnaturalearth” (Massicotte & South, 2023), “rnaturalearthhires” (South et al., 2024),”rstatix” (Kassambara, 2023b), “scales” (Wickham et al., 2023), “tidyverse” (Wickham et al., 2019), “viridisLite” (Garnier et al., 2023). **4.2 Python** Python scripts, including “d18O Analysis Script.ipynb” and “NetCDF Plotting.ipynb”, utilize the Jupyter Notebook interactive ‘platform and are executed using Python version 3.9.16. Install the following libraries before running scripts: “xarray” (Hoyer & Joseph, 2017), “matplotlib” (Hunter, 2007), “cartopy” (Met Office, 2015). **5 References** Amenábar, C. R., Montes, M., Nozal, F., & Santillana, S. (2020). Dinoflagellate cysts of the la Meseta Formation (middle to late Eocene), Antarctic Peninsula: Implications for biostratigraphy, palaeoceanography and palaeoenvironment. *Geological Magazine*, *157*(3), 351–366. [https://doi.org/10.1017/S0016756819000591](https://doi.org/10.1017/S0016756819000591) Auguie, B. (2017). gridExtra: Miscellaneous Functions for “Grid” Graphics. Retrieved from [https://cran.r-project.org/package=gridExtra](https://cran.r-project.org/package=gridExtra) van den Brand, T. (2024). ggh4x: Hacks for “ggplot2.” Retrieved from [https://cran.r-project.org/package=ggh4x](https://cran.r-project.org/package=ggh4x) Campitelli, E. (2024). ggnewscale: Multiple Fill and Colour Scales in “ggplot2.” Retrieved from [https://cran.r-project.org/package=ggnewscale](https://cran.r-project.org/package=ggnewscale) code by Richard A. Becker, O. S., & version by Ray Brownrigg., A. R. W. R. (2022). mapdata: Extra Map Databases. Retrieved from [https://cran.r-project.org/package=mapdata](https://cran.r-project.org/package=mapdata) code by Richard A. Becker, O. S., version by Ray Brownrigg. Enhancements by Thomas P Minka, A. R. W. R., & Deckmyn., A. (2023). maps: Draw Geographical Maps. Retrieved from [https://cran.r-project.org/package=maps](https://cran.r-project.org/package=maps) Douglas, P. M. J., Affek, H. P., Ivany, L. C., Houben, A. J. P., Sijp, W. P., Sluijs, A., et al. (2014). Pronounced zonal heterogeneity in Eocene southern high-latitude sea surface temperatures. *Proceedings of the National Academy of Sciences of the United States of America*, *111*(18), 6582–6587. [https://doi.org/10.1073/pnas.1321441111](https://doi.org/10.1073/pnas.1321441111) Dunnington, D. (2023). ggspatial: Spatial Data Framework for ggplot2. Retrieved from [https://cran.r-project.org/package=ggspatial](https://cran.r-project.org/package=ggspatial) Engelbrecht, A., Mörs, T., Reguero, M. A., & Kriwet, J. (2016a). A new sawshark, Pristiophorus laevis, from the Eocene of Antarctica with comments on Pristiophorus lanceolatus. *Historical Biology*, *29*(6), 841–853. [https://doi.org/10.1080/08912963.2016.1252761](https://doi.org/10.1080/08912963.2016.1252761) Engelbrecht, A., Mörs, T., Reguero, M. A., & Kriwet, J. (2016b). Revision of Eocene Antarctic carpet sharks (Elasmobranchii, Orectolobiformes) from Seymour Island, Antarctic Peninsula. *Journal of Systematic Palaeontology*, *15*(12), 969–990. [https://doi.org/10.1080/14772019.2016.1266048](https://doi.org/10.1080/14772019.2016.1266048) Engelbrecht, A., Mörs, T., Reguero, M. A., & Kriwet, J. (2017a). Eocene squalomorph sharks (Chondrichthyes, Elasmobranchii) from Antarctica. *Journal of South American Earth Sciences*, *78*, 175–189. [https://doi.org/10.1016/j.jsames.2017.07.006](https://doi.org/10.1016/j.jsames.2017.07.006) Engelbrecht, A., Mörs, T., Reguero, M. A., & Kriwet, J. (2017b). New carcharhiniform sharks (Chondrichthyes, Elasmobranchii) from the early to middle Eocene of Seymour Island, Antarctic Peninsula. *Journal of Vertebrate Paleontology*, *37*(6). [https://doi.org/10.1080/02724634.2017.1371724](https://doi.org/10.1080/02724634.2017.1371724) Engelbrecht, A., Mörs, T., Reguero, M. A., & Kriwet, J. (2019). Skates and rays (Elasmobranchii, Batomorphii) from the Eocene La Meseta and Submeseta formations, Seymour Island, Antarctica. *Historical Biology*, *31*(8), 1028–1044. [https://doi.org/10.1080/08912963.2017.1417403](https://doi.org/10.1080/08912963.2017.1417403) for R by Ray Brownrigg, D. M. P., Minka, T. P., & transition to Plan 9 codebase by Roger Bivand. (2023). mapproj: Map Projections. Retrieved from [https://cran.r-project.org/package=mapproj](https://cran.r-project.org/package=mapproj) Garnier, Simon, Ross, Noam, Rudis, Robert, et al. (2023). {viridis(Lite)} - Colorblind-Friendly Color Maps for R. [https://doi.org/10.5281/zenodo.4678327](https://doi.org/10.5281/zenodo.4678327) Gohel, D., & Skintzos, P. (2024). flextable: Functions for Tabular Reporting. Retrieved from [https://cran.r-project.org/package=flextable](https://cran.r-project.org/package=flextable) Griffiths, M. L., Eagle, R. A., Kim, S. L., Flores, R. J., Becker, M. A., IV, H. M. M., et al. (2023). Endothermic physiology of extinct megatooth sharks. *Proceedings of the National Academy of Sciences*, *120*(27), e2218153120. [https://doi.org/10.1073/PNAS.2218153120](https://doi.org/10.1073/PNAS.2218153120) Grunenwald, A., Keyser, C., Sautereau, A. M., Crubézy, E., Ludes, B., & Drouet, C. (2014). Revisiting carbonate quantification in apatite (bio)minerals: A validated FTIR methodology. *Journal of Archaeological Science*, *49*(1), 134–141. [https://doi.org/10.1016/j.jas.2014.05.004](https://doi.org/10.1016/j.jas.2014.05.004) Henry, L., Wickham, H., & Chang, W. (2024). ggstance: Horizontal “ggplot2” Components. Retrieved from [https://cran.r-project.org/package=ggstance](https://cran.r-project.org/package=ggstance) Hoyer, S., & Joseph, H. (2017). xarray: N-D labeled Arrays and Datasets in Python. *Journal of Open Research Software*, *5*(1), 17. [https://doi.org/10.5334/jors.148](https://doi.org/10.5334/jors.148) Hunter, J. D. (2007). Matplotlib: A 2D graphics environment. *Computing in Science & Engineering*, *9*(3), 90–95. [https://doi.org/10.1109/MCSE.2007.55](https://doi.org/10.1109/MCSE.2007.55) Ivany, L. C., Lohmann, K. C., Hasiuk, F., Blake, D. B., Glass, A., Aronson, R. B., & Moody, R. M. (2008). Eocene climate record of a high southern latitude continental shelf: Seymour Island, Antarctica. *Bulletin of the Geological Society of America*, *120*(5–6), 659–678. [https://doi.org/10.1130/B26269.1](https://doi.org/10.1130/B26269.1) Judd, E. J., Ivany, L. C., DeConto, R. M., Halberstadt, A. R. W., Miklus, N. M., Junium, C. K., & Uveges, B. T. (2019). Seasonally Resolved Proxy Data From the Antarctic Peninsula Support a Heterogeneous Middle Eocene Southern Ocean. *Paleoceanography and Paleoclimatology*, *34*(5), 787–799. [https://doi.org/10.1029/2019PA003581](https://doi.org/10.1029/2019PA003581) Kassambara, A. (2023a). ggpubr: “ggplot2” Based Publication Ready Plots. Retrieved from [https://cran.r-project.org/package=ggpubr](https://cran.r-project.org/package=ggpubr) Kassambara, A. (2023b). rstatix: Pipe-Friendly Framework for Basic Statistical Tests. Retrieved from [https://cran.r-project.org/package=rstatix](https://cran.r-project.org/package=rstatix) Kelley, D., & Richards, C. (2023). oce: Analysis of Oceanographic Data. Retrieved from [https://cran.r-project.org/package=oce](https://cran.r-project.org/package=oce) Kévin Allan Sales Rodrigues. (2023). greekLetters: routines for writing Greek letters and mathematical symbols on the RStudio and RGui. Retrieved from [https://cran.r-project.org/package=greekLetters](https://cran.r-project.org/package=greekLetters) Kolodny, Y., Luz, B., & Navon, O. (1983). Oxygen isotope variations in phosphate of biogenic apatites, I. Fish bone apatite-rechecking the rules of the game. *Earth and Planetary Science Letters*, *64*(3), 398–404. [https://doi.org/10.1016/0012-821X(83)90100-0](https://doi.org/10.1016/0012-821X\(83\)90100-0) Kriwet, J. (2005). Additions to the Eocene selachian fauna of Antarctica with comments on Antarctic selachian diversity. *Journal of Vertebrate Paleontology*, *25*(1), 1–7. [https://doi.org/10.1671/0272-4634(2005)025\[0001:ATTESF\]2.0.CO;2](https://doi.org/10.1671/0272-4634\(2005\)025[0001:ATTESF]2.0.CO;2) Kriwet, J., Engelbrecht, A., Mörs, T., Reguero, M., & Pfaff, C. (2016). Ultimate Eocene (Priabonian) chondrichthyans (Holocephali, Elasmobranchii) of Antarctica. *Journal of Vertebrate Paleontology*, *36*(4). [https://doi.org/10.1080/02724634.2016.1160911](https://doi.org/10.1080/02724634.2016.1160911) Larocca Conte, G., Lopes, L. E., Mine, A. H., Trayler, R. B., & Kim, S. L. (2024). SPORA, a new silver phosphate precipitation protocol for oxygen isotope analysis of small, organic-rich bioapatite samples. *Chemical Geology*, *651*, 122000. [https://doi.org/10.1016/J.CHEMGEO.2024.122000](https://doi.org/10.1016/J.CHEMGEO.2024.122000) Lécuyer, C., Amiot, R., Touzeau, A., & Trotter, J. (2013). Calibration of the phosphate δ18O thermometer with carbonate-water oxygen isotope fractionation equations. *Chemical Geology*, *347*, 217–226. [https://doi.org/10.1016/j.chemgeo.2013.03.008](https://doi.org/10.1016/j.chemgeo.2013.03.008) Long, D. J. (1992). Sharks from the La Meseta Formation (Eocene), Seymour Island, Antarctic Peninsula. *Journal of Vertebrate Paleontology*, *12*(1), 11–32. [https://doi.org/10.1080/02724634.1992.10011428](https://doi.org/10.1080/02724634.1992.10011428) Longinelli, A. (1965). Oxygen isotopic composition of orthophosphate from shells of living marine organisms. *Nature*, *207*(4998), 716–719. [https://doi.org/10.1038/207716a0](https://doi.org/10.1038/207716a0) Longinelli, A., & Nuti, S. (1973). Revised phosphate-water isotopic temperature scale. *Earth and Planetary Science Letters*, *19*(3), 373–376. [https://doi.org/10.1016/0012-821X(73)90088-5](https://doi.org/10.1016/0012-821X\(73\)90088-5) Marramá, G., Engelbrecht, A., Mörs, T., Reguero, M. A., & Kriwet, J. (2018). The southernmost occurrence of Brachycarcharias (Lamniformes, Odontaspididae) from the Eocene of Antarctica provides new information about the paleobiogeography and paleobiology of Paleogene sand tiger sharks. *Rivista Italiana Di Paleontologia e Stratigrafia*, *124*(2), 283–297. Massicotte, P., & South, A. (2023). rnaturalearth: World Map Data from Natural Earth. Retrieved from [https://cran.r-project.org/package=rnaturalearth](https://cran.r-project.org/package=rnaturalearth) Met Office. (2015). Cartopy: a cartographic python library with a Matplotlib interface. Exeter, Devon. Retrieved from [https://scitools.org.uk/cartopy](https://scitools.org.uk/cartopy) Mine, A. H., Waldeck, A., Olack, G., Hoerner, M. E., Alex, S., & Colman, A. S. (2017). Microprecipitation and δ18O analysis of phosphate for paleoclimate and biogeochemistry research. *Chemical Geology*, *460*(March), 1–14. [https://doi.org/10.1016/j.chemgeo.2017.03.032](https://doi.org/10.1016/j.chemgeo.2017.03.032) Montes, M., Nozal, F., Santillana, S., Marenssi, S., & Olivero, E. (2013). Mapa Geológico de Isla Marambio (Seymour), Antártida, escala 1:20,000. *Serie Cartográfica*. Neuwirth, E. (2022). RColorBrewer: ColorBrewer Palettes. Retrieved from [https://cran.r-project.org/package=RColorBrewer](https://cran.r-project.org/package=RColorBrewer) Oscar Perpiñán, & Robert Hijmans. (2023). rasterVis. Retrieved from [https://oscarperpinan.github.io/rastervis/](https://oscarperpinan.github.io/rastervis/) Pierce, D. (2023). ncdf4: Interface to Unidata netCDF (Version 4 or Earlier) Format Data Files. Retrieved from [https://cran.r-project.org/package=ncdf4](https://cran.r-project.org/package=ncdf4) Pucéat, E., Joachimski, M. M., Bouilloux, A., Monna, F., Bonin, A., Motreuil, S., et al. (2010). Revised phosphate-water fractionation equation reassessing paleotemperatures derived from biogenic apatite. *Earth and Planetary Science Letters*, *298*(1–2), 135–142. [https://doi.org/10.1016/j.epsl.2010.07.034](https://doi.org/10.1016/j.epsl.2010.07.034) R Development Core Team. (2024). A language and environment for statistical computing. R Foundation for Statistical Computing, Vienna, Austria. Vienna, Austria. RStudio Team. (2024). RStudio: Integrated Development for R. Boston, MA: RStudio, PBC. Retrieved from [http://www.rstudio.com/](http://www.rstudio.com/). Sarkar, D. (2008). *Lattice: Multivariate Data Visualization with R*. New York: Springer. Retrieved from [http://lmdvr.r-forge.r-project.org](http://lmdvr.r-forge.r-project.org) Shemesh, A. (1990). Crystallinity and diagenesis of sedimentary apatites. *Geochimica et Cosmochimica Acta*, *54*(9), 2433–2438. [https://doi.org/10.1016/0016-7037(90)90230-I](https://doi.org/10.1016/0016-7037\(90\)90230-I) Signorell, A. (2024). DescTools: Tools for Descriptive Statistics. Retrieved from [https://cran.r-project.org/package=DescTools](https://cran.r-project.org/package=DescTools) South, A., Michael, S., & Massicotte, P. (2024). rnaturalearthhires: High Resolution World Vector Map Data from Natural Earth used in rnaturalearth. Retrieved from [https://github.com/ropensci/rnaturalearthhires](https://github.com/ropensci/rnaturalearthhires) Trayler, R. B., Landa, P. V., & Kim, S. L. (2023). Evaluating the efficacy of collagen isolation using stable isotope analysis and infrared spectroscopy. *Journal of Archaeological Science*, *151*, 105727. [https://doi.org/10.1016/j.jas.2023.105727](https://doi.org/10.1016/j.jas.2023.105727) Wickham, H., Averick, M., Bryan, J., Chang, W., McGowan, L. D., François, R., et al. (2019). Welcome to the {tidyverse}. *Journal of Open Source Software*, *4*(43), 1686. [https://doi.org/10.21105/joss.01686](https://doi.org/10.21105/joss.01686) Wickham, H., Pedersen, T. L., & Seidel, D. (2023). scales: Scale Functions for Visualization. Retrieved from [https://cran.r-project.org/package=scales](https://cran.r-project.org/package=scales) Wilke, C. O. (2024). cowplot: Streamlined Plot Theme and Plot Annotations for “ggplot2.” Retrieved from [https://cran.r-project.org/package=cowplot](https://cran.r-project.org/package=cowplot) Xu, S. (2022). ggstar: Multiple Geometric Shape Point Layer for “ggplot2.” Retrieved from [https://cran.r-project.org/package=ggstar](https://cran.r-project.org/package=ggstar) Zeileis, A., Fisher, J. C., Hornik, K., Ihaka, R., McWhite, C. D., Murrell, P., et al. (2020). {colorspace}: A Toolbox for Manipulating and Assessing Colors and Palettes. *Journal of Statistical Software*, *96*(1), 1–49. [https://doi.org/10.18637/jss.v096.i01](https://doi.org/10.18637/jss.v096.i01) Zhu, J., Poulsen, C. J., Otto-Bliesner, B. L., Liu, Z., Brady, E. C., & Noone, D. C. (2020). Simulation of early Eocene water isotopes using an Earth system model and its implication for past climate reconstruction. *Earth and Planetary Science Letters*, *537*, 116164. [https://doi.org/10.1016/j.epsl.2020.116164](https://doi.org/10.1016/j.epsl.2020.116164) Eocene climate cooling, driven by the falling pCO2 and tectonic changes in the Southern Ocean, impacted marine ecosystems. Sharks in high-latitude oceans, sensitive to these changes, offer insights into both environmental shifts and biological responses, yet few paleoecological studies exist. The Middle-to-Late Eocene units on Seymour Island, Antarctica, provide a rich, diverse fossil record, including sharks. We analyzed the oxygen isotope composition of phosphate from shark tooth bioapatite (δ18Op) and compared our results to co-occurring bivalves and predictions from an isotope-enabled global climate model to investigate habitat use and environmental conditions. Bulk δ18Op values (mean 22.0 ± 1.3‰) show no significant changes through the Eocene. Furthermore, the variation in bulk δ18Op values often exceeds that in simulated seasonal and regional values. Pelagic and benthic sharks exhibit similar δ18Op values across units but are offset relative to bivalve and modeled values. Some taxa suggest movements into warmer or more brackish waters (e.g., Striatolamia, Carcharias) or deeper, colder waters (e.g., Pristiophorus). Taxa like Raja and Squalus display no shift, tracking local conditions in Seymour Island. The lack of difference in δ18Op values between pelagic and benthic sharks in the Late Eocene could suggest a poorly stratified water column, inconsistent with a fully opened Drake Passage. Our findings demonstrate that shark tooth bioapatite tracks the preferred habitat conditions for individual taxa rather than recording environmental conditions where they are found. A lack of secular variation in δ18Op values says more about species ecology than the absence of regional or global environmental changes. See methods in Larocca Conte, G., Aleksinski, A., Liao, A., Kriwet, J., Mörs, T., Trayler, R. B., Ivany, L. C., Huber, M., Kim, S. L. (2024). Eocene Shark Teeth From Peninsular Antarctica: Windows to Habitat Use and Paleoceanography. Paleoceanography and Paleoclimatology, 39, e2024PA004965.

Aim: We propose that forest trees create a vertical dimension for ecological niche variation that generates different regimes of climatic exposure, which in turn drives species elevation distributions. We test this hypothesis by statistically modelling the vertical and elevation distributions and microclimate exposure of rainforest ants. Location: Wet Tropics Bioregion, Australia Methods: We conducted 60 ground-to-canopy surveys to determine the vertical (tree) and elevation distributions, and microclimate exposure of ants (101 species) at 15 sites along four mountain ranges. We statistically modelled elevation range size as a function of ant species’ vertical niche breadth and exposure to temperature variance for 55 species found at two or more trees. Results: We found a positive association between vertical niche and elevation range of ant species: for every 3 m increase in vertical niche breadth our models predict a ~150% increase in mean elevation range size. Temperature variance increased with vertical height along the arboreal gradient and ant species exposure to temperature variance explained some of the variation in elevation range size. Main Conclusions: We demonstrate that arboreal ants have broader elevation ranges than ground-dwelling ants and are likely to have increased resilience to climatic variance. The capacity of species to expand their niche by climbing trees could influence their ability to persist over broader elevation ranges. We propose that wherever vertical layering exists - from oceans to forest ecosystems - vertical niche breadth is a potential mechanism driving macrogeographic distribution patterns and resilience to climate change. Data_collections.csv Main survey collections data in a site by species matrix showing all data for all sites surveyed. Tuna baited vials were placed every three metres from ground to canopy in trees at elevation sites at four subregion mountain ranges of the Australian Wet Tropics Bioregion. Note data file includes empty vials that lacked ants. Microclimate_AthertonTemp.csv This file contains Atherton Uplands temperature data from ibuttons deployed at one tree per elevation (200, 400, 600, 800, 1000) at every three metres in height in Dec-Jan 2017- 2018 set to record every half hour. See file Metadata for details of column names and data values.

Patterns of movement of marine species can reflect strategies of reproduction and dispersal, species’ interactions, trophodynamics, and susceptibility to change, and thus critically inform how we manage populations and ecosystems. On coral reefs, the density and diversity of metazoan taxa is greatest in dead coral and rubble, which is suggested to fuel food webs from the bottom-up. Yet, biomass and secondary productivity in rubble is predominantly available in some of the smallest individuals, limiting how accessible this energy is to higher trophic levels. We address the bioavailability of motile coral reef cryptofauna based on small-scale patterns of emigration in rubble. We deployed modified RUbble Biodiversity Samplers (RUBS) and emergence traps in a shallow rubble patch at Heron Island, Great Barrier Reef, to detect community-level differences in the directional influx of motile cryptofauna under five habitat accessibility regimes. The mean density (0.13–4.5 ind.cm-3) and biomass (0.14–5.2 mg.cm-3) of cryptofauna were high and varied depending on microhabitat accessibility. Emergent zooplankton represented a distinct community (dominated by the Appendicularia and Calanoida) with the lowest density and biomass, indicating constraints on nocturnal resource availability. Mean cryptofauna density and biomass were greatest when interstitial access within rubble was blocked, driven by the rapid proliferation of small harpacticoid copepods from the rubble surface, leading to trophic simplification. Individuals with high biomass (e.g., decapods, gobies, and echinoderms) were greatest when interstitial access within rubble was unrestricted. Treatments with a closed rubble surface did not differ from those completely open, suggesting that top-down predation does not diminish rubble-derived resources. Our results show that conspecific cues and species’ interactions (e.g., competition and predation) within rubble are most critical in shaping ecological outcomes within the cryptobiome. These findings have implications for prey accessibility through trophic and community size structuring in rubble, which may become increasingly relevant as benthic reef complexity shifts in the Anthropocene. We address the bioavailability of coral reef cryptofauna in rubble based on small-scale patterns of emigration. We adapted the accessibility of Rubble Biodiversity Samplers (RUBS), models used to standardise biodiversity sampling in rubble (Wolfe and Mumby 2020), to explore the local movement patterns of rubble-dwelling fauna, with inference to predation processes within and beyond the cryptobenthos. Five treatments were developed to detect community-level differences in the directional influx of motile cryptofauna under various habitat accessibility regimes. Four of these treatments were developed by modifying accessibility into RUBS (https://www.thingiverse.com/thing:4176644/files) to understand limitations on the directional influx and movement of cryptofauna within coral rubble patches using four treatments; (1) open (completely accessible), (2) interstitial access (top closed), (3) surficial access (sides and bottom closed), and (4) raised (above rubble substratum). The fifth treatment involved a series of emergence plankton traps, designed to target demersal cryptofauna that vertically migrate from within the rubble benthos at night, given emergent zooplankton biomass and diversity are greatest at night. Fieldwork was conducted over several weeks (11th September to 5th October 2021) in a shallow (~3–5 m depth) reef slope site on the southern margin of Heron Island (-23˚26.845’ S, 151˚54.732’ E), Great Barrier Reef, Australia (Fig. 1). All collections were conducted under the Great Barrier Reef Marine Park Authority permit G20/44613.1.

This file contains the AiNU data used for the article entitled by Physics-based material parameters extraction from perovskite experiments via Bayesian optimization (https://arxiv.org/abs/2402.11101).

Literature search and screening We searched for relevant literature with publication month and years Jan 2000- Nov 2022 in two databases: Web of Science (https://www.webofscience.com/; The Core Collection) and Scopus (https://www.scopus.com). We used the same nested Boolean (i.e., AND between different groups of search terms, OR within groups of similar search terms and NOT for excluding search terms) search string in the title, abstract and keywords fields for both Web of Science (TS) and Scopus (TITLE-ABS-KEY) (complete search strings in the supplementary material, Appendix S1). We targeted the relevant deer species for the boreal and temperate forests (i.e., Alces alces, Capreolus capreolus, Cervus spp., Dama dama, Odocoileus spp., Rangifer tarandus; for distribution maps, see Fig. S2), by using a combination of Latin and common names that we combined with geographical constraints based on names of biogeographical regions, countries, and states. We combined this search string with climate related variables (temperature, precipitation etc., Appendix S2). From here on, we refer to Cervus elaphus as red deer, and C. canadensis as wapiti. We refer to R. tarandus living in Europe and Asia as reindeer but as caribou when living in North America. We restrained the search by language (English) and document type (peer-reviewed papers). Our aim was to be as least exclusive as possible, but this led to some unexpected irrelevant documents. We therefore added exclusion terms to filter out non-targeted biogeographical regions and scientific fields. We did not exclude any topical part of our search because it would be impossible to make a coherent pre-emptive list of terms to exclude. The search hits from Web of Science and Scopus were merged and cleaned of duplicates, resulting in 8154 unique papers. Screening of papers was conducted using Rayyan (Ouzzani et al. 2016), a free web application for reviewing articles. Decisions on exclusion or inclusion were first made by reading the title and abstract of each article and determining their conformity to the criteria targeted by the search terms: right topic (i.e., in context of climate change), species (Cervidae excluding semi-domestic reindeer), geography (boreal and temperate zones), language (English) and type of study (new, or new synthesis of, empirical temporal data on deer response to climate). We included papers of migratory caribou residing in forest for larger parts of the year. Note that papers did not have to specify a climate change context to be included. It was sufficient that it contained temporal data on deer and weather variations. Given the controversies surrounding definitions of climate change, rather few papers proclaim having documented climate change and a stricter criterion would have excluded almost all papers. The robustness of the exclusion criteria and the individual screener divergence of the first screening were tested before the actual screening was done. Fifty randomly drawn papers were reviewed by all authors individually without conferring. The papers were randomly distributed among authors. The discrepancies were rather few (13 out of 49 papers (27%) had at least 1 person with a different opinion than the others). After discussing each of these cases in detail, the basis for coherent decision making was improved. To verify the improvement, another control procedure was applied for the remaining screening: 289 papers were each read by two to four authors. The result of this control screening showed 18 (6%) conflicting decisions. Screening of the remaining 7815 papers was done by the authors one by one and assigned equally among readers according to alphabetic order by the first author of the papers. The first screening finally generated 556 papers possibly relevant for the review. All papers with conflicting decisions in the test and control screenings were included among the 556. The possibly relevant papers were then equally divided between the authors. These papers were read completely and again scrutinized for conformation to criteria, resulting in a final list of 218 papers relevant for review. Data from these papers were then tabulated and systemized per demographics (species, location, season, etc.), deer responses and climate factor. Further details on this data collection are specified in Appendix S3. The table here in Dryad includes the detailed tabulations used to produce Table 1, Figure 1, Figure in the main article, and Table S3 in the Appendix.  Climate change causes far-reaching disruption in nature, where tolerance thresholds already have been exceeded for some plants and animals. In the short-term, deer may respond to climate through individual physiological and behavioral responses. Over time, individual responses can aggregate to the population level and ultimately lead to evolutionary adaptations. We systematically reviewed literature (published 2000-2022) to summarize the effect of temperature, rainfall, snow, combined measures (e.g., the North Atlantic Oscillation) and extreme events, on deer species inhabiting boreal and temperate forests in terms of their physiology, spatial use and population dynamics. We targeted deer species which inhabit relevant biomes in North America, Europe and Asia: moose, roe deer, elk, red deer, sika deer, fallow deer, white-tailed deer, mule deer, caribou and reindeer. Our review (218 papers) shows that many deer populations will likely benefit in-part from warmer winters, but hotter and drier summers may exceed their physiological tolerances. We found support for deer expressing both morphological, physiological, and behavioral plasticity in response to climate variability. For example, some deer species can limit the effects of harsh weather conditions by modifying habitat use and daily activity patterns, while the physiological responses of female deer can lead to long-lasting effects on population dynamics. We identified 20 patterns, among which some illustrate antagonistic pathways, suggesting that detrimental effects will cancel out some of the benefits of climate change. Our findings highlight the influence of local variables (eg. population density and predation) for how deer will respond to climatic conditions. We identified several knowledge gaps, such as studies regarding the potential impact on these animals of extreme weather events, snow type and wetter autumns. The patterns we have identified in this literature review should help managers understand how populations of deer may be affected by regionally projected futures regarding temperature, rainfall and snow. # Literature review protocol: Climate change and deer in boreal and temperate regions [https://doi.org/10.5061/dryad.jh9w0vtmd](https://doi.org/10.5061/dryad.jh9w0vtmd) ## Description of the data and file structure We systematically reviewed literature (published 2000-2022) to summarize the effect of temperature, rainfall, snow, combined measures (e.g., the North Atlantic Oscillation) and extreme events, on deer species inhabiting boreal and temperate forests in terms of their physiology, spatial use and population dynamics. We targeted deer species which inhabit relevant biomes in North America, Europe and Asia: moose, roe deer, elk, red deer, sika deer, fallow deer, white-tailed deer, mule deer, caribou and reindeer. After screening, 218 articles remained. The data made available here pertains to these articles. ### Files and variables #### File: Felton\_et\_al\_2024\_GCB\_Protocol\_literature\_review\_Dryad 30 aug no hidden columns.xlsx **Description:** protocol for tabulating relevant information from published literature. ##### Variables * Column B-G: Climatic variables that the studies assessed (temperature, rainfall, snow, combined measures, extreme climatic events) * Column H: animal species * Column I: extreme events * Column K-AF: registration whether information is presented that relate to the three larger topics of the review (Physiology, Spatial use, Population dynamics) and to any of the 20 Patterns Found, which are summarised in Table 2 in the main article. Abbreviations refer to details of such patterns, which are explained in the heading of Table 2 in the main article. * Blank cells = no relevant information exist. Data was derived from the following sources: * We searched for relevant literature with publication month and years Jan 2000- Nov 2022 in two databases: Web of Science ([https://www.webofscience.com/](https://www.webofscience.com/); The Core Collection) and Scopus ([https://www.scopus.com](https://www.scopus.com/)).