BUSCADOR RECOLECTA

Dataset are structured following well-established data formats Darwing Core. Three files are provided. The first file (Don_liz-gec_event_20221201) contains the information of the project, the institution and the description each event (time of occurrence, geographical coordinates, habitat type, etc…); the second file (Don_liz-gec_occ_20221201) contains the information of the occurrences of species recorded in each transect, taxonomic classification, geographoical coordinates of its observation, etc…; and the third file (Don_liz-gec_mof_20221201) provide information of the description of other variables measured during counts., The monitoring of lizards and geckos’ community in Doñana was initiated in 2005 as part of the monitoring program of natural resources and processes. One of the aims of this project was to obtain a temporal and continuous series of data of the presence and abundance of these species to detect changes and trends in their wild populations within the protected area. Lizard and gecko counts are collected annually by members of the monitoring team three times per year in the periods when reptile activity is high (two samplings in March-June and one sampling in September-October), with good environmental conditions (temperature between 17 and 25 ºC and absence of rain or strong wind conditions). The method used to record the presence and abundance of lizard and gecko species (kilometric index: number of individuals/km) are the transect censuses. Individuals are visually searched along seven transects (linear or circular) that are carried out by one trained person on foot. Each transect have an established length, but the length surveyed is different (averaged 1858 ± 51.86 meters). Each transect is located in dunes or mediterranean vegetation habitats, representative of Doñana ecosystems. Two transects runs through wooden footpaths (within the Natural Park and five transects are placed on sand-trails (4 within the National Park and one in Natural Park which contains a small part of wooden footpath). Linear transects have been surveyed in the outward track direction and after waiting 15 minutes it has been again performed on the way back. This information is valid to account for species presence in the area. However, we suggest to choose the first survey performed for statistical analyses that require independence of samples. Eight species can be potentially observed during the samplings: Mediterranean house gecko (Hemidactylus turcicus), Common Wall Gecko (Tarentola mauritanica), Fringe-fingered Lizard (Acanthodactylus erythrurus), Algerian Psammodromus (Psammodromus algirus), Western Psammodromus (Psammodromus occidentalis), Carbonell's wall lizard (Podarcis carbonelli), Andalusian wall lizard (Podarcis vaucheri), Ocellated lizard (Timon lepidus). Other reptile species present in Doñana have not been included in this study as the detection with this method is very low. For instance, to detect species with a fossorial behaviour (the Mediterranean Worm Lizard Blanus cinereus) or those cryptic as adders (Vipera latastei), it is required larger investment of survey which consist in looking under vegetation, burrows or logs (there are not stones in Doñana). Data recorded during the surveys include weather description (cloud cover, temperature, rain, or wind speed), species identification, number of individuals, sex and life stage of the reptiles when possible, as well as time and georeferenced data of the observation. Between 2005-2007 data was registered in Excel file and since 2008 data is recorded with the app CyberTracker (see protocol). The protocol used has been supervised by herpetological researchers and the data have been validated by the members who performed the transects., We acknowledge financial support from National Parks Autonomous Agency (OAPN) between 2002-2007; Singular Scientific and Technical Infrastructures from the Spanish Science and Innovation Ministry (ICTS-MICINN); Ministry of Agriculture, Livestock, Fisheries and Sustainable Development from the Regional Government of Andalusia (CAGPDES-JA) since 2007; and Doñana Biological Station from the Spanish National Research Council (EBD-CSIC) since all the study period (2005)., 1. Don_liz-gec_event_20221201: eventID, institutionCode, institutionID, datasetName, continent, country, countryCode, Province, Location, Locality, LocalityID, eventDate, eventTime, decimalLatitude, decimalLongitude, decimalLatitudeEnd, decimalLongitudeEnd, verbatimCoordinate_func, verbatimCoordinate, habitat, sampleSizeValue, sampleSizeUnit, samplingEffort, recordedBy, samplingProtocol.-- 2. Don_liz-gec_occ_20221201: RecordedBy, eventID, OccurrenceID, OcurrenceTime, decimalLatitude, decimalLongitude, basisOfRecord, individualCount, lifeStage, sex, OccurrenceRemarks, behavior, kingdom, Class, Family, scientificName, genus, specificEpithet, scientificNameAuthorship, taxonRank.-- 3. Don_liz-gec_mof_20221201: eventID, measurementID, measurementType, measurementValue, measurementUnit, measurementMethod., Peer reviewed

Dataset are structured following well-established data formats. Three files are provided. The first file (Meta-data) contains the information of each event (time of occurrence, geographical coordinates, Ecosystem, Sampling mehtod, etc…); the second file (Fish) contains the information of the occurrences of fish species recorded in each station, taxonomic classification, etc…; and the third file (Macroinvertebrates) provide information of the occurrences of macroinvertebrates recorded in each station, taxonomic classification, abundance clases, etc…, The monitoring of the macroinvertebrates and fish community in Doñana wetlands was initiated in 2004 as part of the Monitoring Program of Natural Resources and Processes. The aim was to obtain a temporal and continuous series of data in the abundance and distribution of macroinvertebrates and fish species to analyze the evolution of their numbers and estimates biodiversity values. Data were recorded annually between 2004-2019 by more than 2 members of the monitoring team which performed samplings in different locations twice per year in winter-spring and summer seasons when the study sites are flooded. The macroinvertebrates and fishes were sampled at the 140 stations classified according to their location (on either aeolian sands or marshland). Funnel traps were used as a sampling method. Between 5-9 funnel traps were randomly distributed (until 50 cm of depth) in each location, depending of the flooded area and depth. The traps were left for 24 hours and emptied the content into white sorting pans. Individuals were counted and identified until the maximun taxonomic level in the field and realease. During samplings, it was identified 66 and 16 families, of macroinvertebrates and fishes respectively. The most abundances were Notonectidae and Corixidae in macroinvertebrates, and Poecilidae and Cyprinidae in fishes. Data recorded during the surveys included species identification, number of individuals, sex and life stage (pupa, larvae, inmature, adult) of the organisms when possible, as well as the time and georreferenced data of the observation. Between 2004-2007 data was registered in Excel file and since 2008 data was recorded in CyberTracker sequence). The protocol used has been supervised by researchers and the data have been validated by the members who performed the sampling., We acknowledge financial support from National Parks Autonomous Agency (OAPN) between 2002-2007; Singular Scientific and Technical Infrastructures from the Spanish Science and Innovation Ministry (ICTS-MICINN); Ministry of Agriculture, Livestock, Fisheries and Sustainable Development from the Regional Government of Andalusia (CAGPDES-JA) since 2007; and Doñana Biological Station from the Spanish National Research Council (EBD-CSIC) since all the study period (2005)., 1. Metadata: Taxa group, Site ID, Site name, Country, y coordinate, x coordinate, Ecosystem River/lake name, Sampling method, Starting year, Ending year, 1st Name, 1st Mail, 2nd Name, 2nd Mail, 3rd Name, 3rd Mail.-- 2. Fish: Site ID, Sample ID, Sampling date, Taxon name, Taxon ID, Definition of abundance class, Abundance class.-- 3. Macroinvertebrates: Site ID, Sample ID, Sampling date, Taxon name, Taxon ID, 0+, 1+, Adult, All., Peer reviewed

Dataset are structured following well-established data formats. Two files are provided. The first file (icts-rbd-dungBe_event_20221107) contains the information of each event (time of occurrence, geographical coordinates, habitat, sampling effort, etc…); the second file (icts-rbd-dungBe_occ_20221107) contains the information of the occurrences of dung-beetles species recorded in each site, numbers of individual recorded and taxonomic classification., The monitoring of the roller dung-beetles (Scarabaeinidae) in Doñana, southwestern Spain, was initiated in 2004 as part of the Monitoring Program of Natural Resources and Processes. The aim was to obtain a temporal and continuous series of data in the abundance and distribution of two species (Scarabeus sacer and S. cicatricosus) present in the area. Data were recorded annually from 2004 to 2012 by members of the monitoring team which performed one sampling (between May and August) in different habitats (sand dunes, mediterranean schrublands, flooplain meadows, and marshlands). Dung-baited pitfalls traps were used as a method to obtain samples of individuals of the two species. These traps were plastic cilinder of 30 cm diameter x 20 cm high buried on the ground. A baited grill of 2cm x 2 cm mesh rested on top of the trap. Bait was fresh horse or cow feaces (250 g) collected around the area early in the morning the day before of trapping. Five pitfall traps were established at each site separated 15 m each other during 24 hours. Two checking were conducted every 12 hours after baiting to avoid the mortality of individuals. Individual of each species were counted and release after it. Data recorded during the surveys included species identification and number of individuals. Between 2004-2008 data was registered in Excel file and since 2008 data was recorded in CyberTracker sequence. The protocol used has been supervised by researchers and the data have been validated by the members who performed the sampling., We acknowledge financial support from National Parks Autonomous Agency (OAPN) between 2004-2007; Singular Scientific and Technical Infrastructures from the Spanish Science and Innovation Ministry (ICTS-MICINN); Ministry of Agriculture, Livestock, Fisheries and Sustainable Development from the Regional Government of Andalusia (CAGPDES-JA) since 2007; and Doñana Biological Station from the Spanish National Research Council (EBD-CSIC) since all the study period (2004-2012)., 1.icts-rbd-dungBeetles_event_20221107: intitutionCode, institutionID, datasetName, eventID, eventDate, eventTime, continent, country, countryCode, stateProvince, locality, decimalLatitude, decimalLongitude, habitat, eventRemarks, sampleSizeValue, sampleSizeUnit, sampleEffort, dynamicPropertiesEvents, recordyBy and scientificName 2.icts-rbd-dungBeetles_occ_20221107: eventID, OccurrenceID, basisOfRecords, individualCount, kingdom, class, family, scientificName, genus, specificEpithet and scientificNameAuthorship., Peer reviewed

Dataset are structured following well-established data formats. Three files are provided. The first file (Don_macroinv_ev_20221222) contains the information of each event (eventID, event date, geographical coordinates, sample effort, etc…); the second file (Don_macroinv_occ_20221222) contains the information of the occurrences of individuals recorded in each station and its taxonomic classification; and the third file (Don_macroinv_mof_20221222) provide information of biometric variable (weithg) of macroinvertebrates samples recorded in each occurrence., The monitoring of the macroinvertebrates community in Doñana wetlands was initiated in 2004 as part of the Monitoring Program of Natural Resources and Processes. The aim was to obtain a temporal and continuous series of data in the abundance and distribution of macroinvertebrates species to analyze the evolution of their numbers and estimates biodiversity values. Data were recorded annually between 2004-2019 by more than 2 members of the monitoring team which performed samplings in different locations twice per year in winter-spring and summer seasons when the study sites are flooded. The macroinvertebrates were sampled at the 139 stations classified according to their location (on either aeolian sands or marshland). Funnel traps were used as a sampling method. Between 5-9 funnel traps were randomly distributed (until 50 cm of depth) in each location, depending of the flooded area and depth. The traps were left for 24 hours and emptied the content into white sorting pans. Individuals were counted and identified until the maximun taxonomic level in the field and realease. During samplings, it was identified 65 families. The most abundances were Notonectidae and Corixidae. Data recorded during the surveys included species identification, number of individuals, sex and life stage (pupa, larvae, juvenile, adult) of the organisms when possible, as well as the time and georreferenced data of the observation. Between 2004-2007 data was registered in Excel file and since 2008 data was recorded in CyberTracker sequence). The protocol used has been supervised by researchers and the data have been validated by the members who performed the sampling., We acknowledge financial support from National Parks Autonomous Agency (OAPN) between 2002-2007; Singular Scientific and Technical Infrastructures from the Spanish Science and Innovation Ministry (ICTS-MICINN); Ministry of Agriculture, Livestock, Fisheries and Sustainable Development from the Regional Government of Andalusia (CAGPDES-JA) since 2007; and Doñana Biological Station from the Spanish National Research Council (EBD-CSIC) since all the study period (2005)., 1. Don_macroinv_ev_20221222: eventID, intitutionCode, datasetName, eventDate, year, month, day, continent, country, stateProvince, location, localityID, locality, sampleSizeUnit, sampleSizeEffort, DynamicPropiertiesEvent, eventRemarks, recordedBy.-- 2. Don_macroinv_occ_20221222: eventID, occurrenceID, basisOfRecord, individualCount, sex, lifeStage, kingdom, phylum, order, family, genus, specificEpithet, scientificName.-- 3. Don_macroinv_mof_20221222: ocurrenceID, measurementID, measurementValue, measurementUnit, measurementType, measurementAccuracy, measurementMethod., Peer reviewed

Dataset are structured following well-established data formats. Three files are provided. The first file (Don_fish_ev_20221222) contains the information of each event (eventID, event date, geographical coordinates, sample effort, etc…); the second file (Don_fish_occ_20221222) contains the information of the occurrences of fish species recorded in each station, taxonomic classification; and the third file (Don_fish_mof_20221222) provide information of the biometric variable (weight) of fish sample in each occurrence., The monitoring of the fish community in Doñana wetlands was initiated in 2004 as part of the Monitoring Program of Natural Resources and Processes. The aim was to obtain a temporal and continuous series of data in the abundance and distribution of fish species to analyze the evolution of their numbers and estimates biodiversity values. Data were recorded annually between 2004-2019 by more than 2 members of the monitoring team which performed samplings in different locations twice per year in winter-spring and summer seasons when the study sites are flooded. The fishes were sampled at the 139 stations classified according to their location (on either aeolian sands or marshland). Funnel traps were used as a sampling method. Between 5-9 funnel traps were randomly distributed (until 50 cm of depth) in each location, depending of the flooded area and depth. The traps were left for 24 hours and emptied the content into white sorting pans. Individuals were counted and identified until the maximun taxonomic level in the field and realease. During samplings, it was identified 15 families. The most abundances were Poecilidae and Cyprinidae. Data recorded during the surveys included species identification, number of individuals, sex and life stage (pupa, larvae, inmature, mature) of the organisms when possible, as well as the time and georreferenced data of the observation. Between 2004-2007 data was registered in Excel file and since 2008 data was recorded in CyberTracker sequence). The protocol used has been supervised by researchers and the data have been validated by the members who performed the sampling., We acknowledge financial support from National Parks Autonomous Agency (OAPN) between 2002-2007; Singular Scientific and Technical Infrastructures from the Spanish Science and Innovation Ministry (ICTS-MICINN); Ministry of Agriculture, Livestock, Fisheries and Sustainable Development from the Regional Government of Andalusia (CAGPDES-JA) since 2007; and Doñana Biological Station from the Spanish National Research Council (EBD-CSIC) since all the study period (2005)., 1. Don_fish_ev_20221222: eventID, intitutionCode, institutionID, datasetName, eventDate, year, month, day, country, stateProvince, location, localityID, locality, decimalLatitude, decimalLongitude, habitat, sampleSizeUnit, sampleSizeEffort, DynamicPropiertiesEvent, eventRemarks, recordeBy.-- 2. Don_fish_occ_20221222: eventID, occurrenceID, individualCount, sex, lifeStage, kingdom, phylum, order, family, genus, specificEpithet, scientificName.-- 3. Don_fish_mof_20221222: OccurrenceID, measurementID, measurementType, measurementValue, measurementUnit, measurementMethod., Peer reviewed

Time Series water turbidity derived from Landsat TM, ETM+ & OLI in the Path 202 Row 34 (Doñana). Also, the product and its metadata are freely available to consult or downloaded in the LAST-EBD Cartography Server: http://mercurio.ebd.csic.es/imgs/ Teh methodology is described in the paper: Empirical models to estimate water turbidity from reflectance data from TM or ETM+ Landsat sensors in shallow wetlands such as Doñana marshes. See the reference: Bustamante, J. et al. 2009. Predictive models of turbidity and water depth in the Doñana marshes using Landsat TM and ETM+ images. Journal of Environmental Management. 90:2219-2225.https://doi.org/10.1016/j.jenvman.2007.08.021, European Commission: ECOPOTENTIAL - ECOPOTENTIAL: IMPROVING FUTURE ECOSYSTEM BENEFITS THROUGH EARTH OBSERVATIONS (641762), Peer reviewed

Time Series of NDVI derived from Landsat TM, ETM+ & OLI in the Path 202 Row 34 (Doñana). Also, the product and its metadata are freely available to consult or downloaded in the LAST-EBD Cartography Server: http://mercurio.ebd.csic.es/imgs/, European Commission: ECOPOTENTIAL - ECOPOTENTIAL: IMPROVING FUTURE ECOSYSTEM BENEFITS THROUGH EARTH OBSERVATIONS (641762), Peer reviewed

Time Series of annual Hydroperiods derived from Landsat MSS, TM, ETM+ & OLI in the Path 202 Row 34 (Doñana) cover the period 1974-2018., European Commission: ECOPOTENTIAL - ECOPOTENTIAL: IMPROVING FUTURE ECOSYSTEM BENEFITS THROUGH EARTH OBSERVATIONS (641762), Peer reviewed

Time Series of flooded areas derived from Landsat TM, ETM+ & OLI in the Path 202 Row 34 (Doñana). Also, these products and its metadata are freely available to consult or downloaded in the LAST-EBD Cartography Server: http://mercurio.ebd.csic.es/imgs/ Methodology is described in this paper: Remote Sensing 8(9):775 · September 2016. DOI: 10.3390/rs8090775, European Commission: ECOPOTENTIAL - ECOPOTENTIAL: IMPROVING FUTURE ECOSYSTEM BENEFITS THROUGH EARTH OBSERVATIONS (641762), Peer reviewed

[Methods] Experiment locations and climate Trans-Mediterranean translocation of Posidonia oceanica fragments took place between Catalunya (Spain), Mallorca (Spain) and Cyprus in July 2018 and were monitored until July 2019 (Fig. 1). Sea surface temperature data for each transplant site were based on daily SST maps with a spatial resolution of 1/4°, obtained from the National Center for Environmental Information (NCEI, https://www.ncdc.noaa.gov/oisst ) (Reynolds et al. 2007). These maps have been generated through the optimal interpolation of Advanced Very High Resolution Radiometer (AVHRR) data for the period 1981-2019. Underwater temperature loggers (ONSET Hobo pro v2 Data logger) were deployed at the transplant sites in Catalunya, Mallorca and Cyprus and recorded hourly temperatures throughout the duration of the experiment (one year). In order to obtain an extended time series of temperature at transplant sites, a calibration procedure was performed comparing logger data with sea surface temperature from the nearest point on SST maps. In particular, SST data were linearly fitted to logger data for the common period. Then, the calibration coefficients were applied to the whole SST time series to obtain corrected-SST data and reconstruct daily habitat temperatures from 1981-2019. Local climate data was also compared to the global thermal distribution of P. oceanica to assess how representative experimental sites were of the thermal distribution of the species (Supplementary materials). Collectively, seawater temperatures from the three locations span the 16th - 99th percentile of temperatures observed across the global thermal distribution of P. oceanica. As such Catalunya, Mallorca and Cyprus are herein considered to represent the cool-edge, centre and warm-edge of P. oceanica distribution, respectively. Transplantation took place toward warmer climates and procedural controls were conducted within each source location, resulting in six source-to-recipient combinations (i.e. treatments, Fig. 1). Initial collection of P. oceanica, handling and transplantation was carried out simultaneously by coordinated teams in July 2018 (Table S1). Each recipient location was subsequently resampled four times over the course of the experiment, in August/September 2018 (T1), October 2018 (T2), April 2019 (T3) and May/June 2019 (T4, Table S1). Between 60-100 fragments were collected for each treatment. A fragment was defined as a section of P. oceanica containing one apical shoot connected with approximately five vertical shoots by approximately 10-15 cm of rhizome with intact roots. Collection occurred at two sites within each location, separated by approximately 1 km. Within sites, collections were conducted between 4 – 5 m depth and were spaced across the meadow to minimise the dominance of a single clone and damage to the meadow. Upon collection, fragments were transported for up to one hour back to the nearest laboratory in shaded seawater. Handling methods In the laboratory, fragments were placed into holding tanks with aerated seawater, at ambient temperature and a 14:10 light-dark cycle. All shoots were clipped to 25 cm length (from meristem to the tip of the longest leaves), to standardise initial conditions and reduce biomass for transportation. For transport by plane or ferry between locations, fragments were packed in layers within cool-boxes. Each layer was separated by frozen cool-packs wrapped in wet tea towels (rinsed in sea water). All fragments spent 12 hrs inside a cool-box irrespective of their recipient destination, including procedural controls (i.e. cool-cool, centre-centre and warm-warm) to simulate the transit times of the plants travelling furthest from their source location (Fig. 1a). On arrival at the destination, fragments were placed in holding tanks with aerated seawater at ambient temperature as described above in their recipient location for 48 hrs, prior to field transplantation. Measurement methods One day prior to transplantation, fragments were tagged with a unique number and attached to U-shaped peg with cable-ties. Morphological traits for each fragment were measured and included: 1) length of the longest apical leaf, width and number of leaves 2) total number of bite marks on leaves of three vertical shoots per fragment, 3) number of vertical shoots, 4) leaf count of three vertical shoots per fragment and 5) overall horizontal rhizome length. A subset (n=10) of fragments per treatment were marked prior transplantation to measure shoot growth. To do this, all shoots within a single fragment were pierced using a hypodermic needle. Two holes were pierced side-by-side at the base of the leaf/top of the meristem. Transplant methods All transplant sites were located in 4 – 5 m depth in area of open dead-matte, surrounded by P. oceanica meadow. In Mallorca and Cyprus, fragments were distributed between two sites, separated by approximately 1 km. In Catalunya, a lack of suitable dead matte habitat, meant that all fragments were placed in one site. Fragments were planted along parallel transects at 50 cm intervals and with a 50 cm gap between parallel transects (Fig. S1). Different treatments were mixed and deployed haphazardly along each transect. Resampling methods and herbivory On day 10 of the experiment, a severe herbivory event was recorded at both warm-edge translocation sites. Scaled photos of all fragments were taken at this time to record the effects of herbivory on transplants. At the end of each main sampling period (T0 – T1, T1-T2 and T3 – T4), all pierced fragments were collected and taken back to the laboratory to measure shoot growth. At T1, T2 and T3, additional sets of fragments (n = 10 per treatment) were marked using the piercing method to record growth in the subsequent time period. In addition, at T1 and T3, n = 20 shoots within the natural meadow at each site were marked to compare growth rates between the native meadow and transplants. Underwater shoot counts and a scaled photo was taken to record fragment survivorship, shoot mortality, bite marks, and shoot length among all remaining fragments within each site and sampling time. In the laboratory, morphological measurements (described above) were repeated on the collected fragments and growth of transplant and natural meadow shoots was measured. Growth (shoot elongation, cm d-1) of the marked shoots was obtained by measuring the length from the base of meristem to marked holes of each leaf (new growth) of the shoot and dividing the leaf elongation per shoot by the marking period (in days). For each shoot, total leaf length (cm shoot-1) and the number of new leaves was also recorded. The rate of new leaf production (new leaves shoot-1 d-1) was estimated dividing the number of new leaves produced per shoot and the marking period. New growth was dried at 60 ºC for 48 hrs to determine carbon and nitrogen content of the leaves, and carbon to nitrogen (C:N) ratios. Carbon and nitrogen concentrations in the new growth leaf tissue was measured at the beginning of the experiment and each subsequent time point for each treatment. Nutrient analyses were conducted at Unidade de Técnicas Instrumentais de Análise (University of Coruña, Spain) with an elemental analyser FlashEA112 (ThermoFinnigan). Underwater photos of shoots were analysed using ImageJ software (https://imagej.net). Maximum leaf length on each shoot in warm-edge transplant sites (cool-warm, centre-warm and warm-warm) were recorded for the initial (day 10) herbivore impact, T1, T2 and T3 time-points and related to transplant nutrient concentrations. Herbivore impact was estimated as the proportional change in length of the longest leaf relative to initial length at T0. Thermal stress Long term maximum temperatures were recorded as the average of annual maximum daily temperatures in each transplant site, averaged between years from 1981-2019. Maximum thermal anomalies were calculated as the difference between daily temperatures in a recipient site over the course of the experiment and the long-term maximum temperature in the source site for each corresponding population. ‘Heat stress’ and ‘recovery’ growth periods of the experiment were defined as T0 -T2 (July-October) and T2-T4 (November-June), respectively, corresponding to periods of positive and negative maximum thermal anomalies. Thermal anomalies experienced by the different transplant treatments were plotted using the ‘geom_flame’, function in the ‘HeatwavesR’ package (Schlegel & Smit 2018) of R (version 3.6.1, 2019) ., 1. The prevalence of local adaptation and phenotypic plasticity among populations is critical to accurately predicting when and where climate change impacts will occur. Currently, comparisons of thermal performance between populations are untested for most marine species or overlooked by models predicting the thermal sensitivity of species to extirpation. 2. Here we compared the ecological response and recovery of seagrass populations (Posidonia oceanica) to thermal stress throughout a year-long translocation experiment across a 2800 km gradient in ocean climate. Transplants in central and warm-edge locations experienced temperatures >29 ºC, representing thermal anomalies >5ºC above long-term maxima for cool-edge populations, 1.5ºC for central and <1ºC for warm-edge populations. 3. Cool, central and warm-edge populations differed in thermal performance when grown under common conditions, but patterns contrasted with expectations based on thermal geography. Cool-edge populations did not differ from warm-edge populations under common conditions and performed significantly better than central populations in growth and survival. 4. Our findings reveal that thermal performance does not necessarily reflect the thermal geography of a species. We demonstrate that warm-edge populations can be less sensitive to thermal stress than cooler, central populations suggesting that Mediterranean seagrasses have greater resilience to warming than current paradigms suggest., Australian Research Council, Award: DE200100900. Horizon 2020 Framework Programme, Award: 659246. Fundación BBVA., Peer reviewed