• Energy Research
• US close
• DE close
• English close
• ZENODO close

This project systematically reviewed the literature for measurements of aboveground carbon stocks in monoculture plantation forests. The data compiled here are for monoculture (single-species) plantation forests, which are a subset of a broader review to identify empirical measurements of carbon stocks across all forest types. The database is structured similarly to that of the ForC (https://forc-db.github.io/) and GROA databases (https://github.com/forc-db/GROA). When using these data, please cite: Bukoski, J.J., Cook-Patton, S.C., Melikov, C., Ban, H., Liu, J.C., Harris, N., Goldman, E., and Potts, M.D. 2022. Rates and drivers of aboveground carbon accumulation in global monoculture plantation forests. Nature Communications 13(4206). doi: 10.1038/s41467-022-31380-7 The code for all analyses in Bukoski et al., 2022 (paper associated with this dataset) is available at https://github.com/jbukoski/GPFC (doi: 10.5281/zenodo.6588710).

This climate change impact data (future scenarios on temperature-induced GDP losses) and climate change mitigation cost data (REMIND model scenarios) is published under doi: 10.5281/zenodo.3541809 and used in this paper: Ueckerdt F, Frieler K, Lange S, Wenz L, Luderer G, Levermann A (2018) The economically optimal warming limit of the planet. Earth System Dynamics. https://doi.org/10.5194/esd-10-741-2019 Below the individual file contents are explained. For further questions feel free to write to Falko Ueckerdt (ueckerdt@pik-potsdam.de). Climate change impact data File 1: Data_rel-GDPpercapita-changes_withCC_per-country_all-RCP_all-SSP_4GCM.csv Content: Data of relative change in absolute GDP/CAP levels (compared to the baseline path of the respective SSP in the SSP database) for each country, RCP (and a zero-emissions scenario), SSP and 4 GCMs (spanning a broad range of climate sensitivity). Negative (positive) values indicate losses (gains) due to climate change. For figure 1a of the paper, this data was aggregated for all countries. File 2: Data_rel-GDPpercapita-changes_withCC_per-country_all-SSP_4GCM_interpolated-for-REMIND-scenarios.csv Content: Data of relative change in absolute GDP/CAP levels (compared to the baseline path of the respective SSP in the SSP database) for each country, SSP and 4 GCMs (spanning a broad range of climate sensitivity). The RCP (and a zero-emissions scenario) are interpolated to the temperature pathways of the ten REMIND model scenarios used for climate change mitigation costs. Hereby the set of scenarios for climate impacts and climate change mitigation are consistent and can be combined to total costs of climate change (for a broad range of mitigation action). File 3: Data_rel-GDPpercapita-changes_withCC_per-country_SSP2_12GCM_interpolated-for-REMIND-scenarios.csv Content: Same as file 2, but only for the SSP2 (chosen default scenario for the study) and for all 12 GCMs. Data of relative change in absolute GDP/CAP levels (compared to the baseline path of the respective SSP in the SSP database) for each country, SSP-2 and 12 GCMs (spanning a broad range of climate sensitivity). The RCP (and a zero-emissions scenario) are interpolated to the temperature pathways of the ten REMIND model scenarios used for climate change mitigation costs. Hereby the set of scenarios for climate impacts and climate change mitigation are consistent and can be combined to total costs of climate change (for a broad range of mitigation action). In addition, reference GDP and population data (without climate change) for each country until 2100 was downloaded from the SSP database, release Version 1.0 (March 2013, https://tntcat.iiasa.ac.at/SspDb/, last accessed 15Nov 2019). Climate change mitigation cost data The scenario design and runs used in this paper have first been conducted in [1] and later also used in [2]. File 4: REMIND_scenario_results_economic_data.csv File 5: REMIND_scenarios_climate_data.csv Content: A broad range of climate change mitigation scenarios of the REMIND model. File 4 contains the economic data of e.g. GDP and macro-economic consumption for each of the countries and world regions, as well as GHG emissions from various economic sectors. File 5 contains the global climate-related data, e.g. forcing, concentration, temperature. In the scenario description “FFrunxxx” (column 2), the code “xxx” specifies the scenario as follows. See [1] for a detailed discussion of the scenarios. The first dimension specifies the climate policy regime (delayed action, baseline scenarios): 1xx: climate action from 2010 5xx: climate action from 2015 2xx climate action from 2020 (used in this study) 3xx climate action from 2030 4x1 weak policy baseline (before Paris agreement) The second dimension specifies the technology portfolio and assumptions: x1x Full technology portfolio (used in this study) x2x noCCS: unavailability of CCS x3x lowEI: lower energy intensity, with final energy demand per economic output decreasing faster than historically observed x4x NucPO: phase out of investments into nuclear energy x5x Limited SW: penetration of solar and wind power limited x6x Limited Bio: reduced bioenergy potential p.a. (100 EJ compared to 300 EJ in all other cases) x6x noBECCS: unavailability of CCS in combination with bioenergy The third dimension specifies the climate change mitigation ambition level, i.e. the height of a global CO2 tax in 2020 (which increases with 5% p.a.). xx1 0$/tCO2 (baseline) xx2 10$/tCO2 xx3 30$/tCO2 xx4 50$/tCO2 xx5 100$/tCO2 xx6 200$/tCO2 xx7 500$/tCO2 xx8 40$/tCO2 xx9 20$/tCO2 xx0 5$/tCO2 For figure 1b of the paper, this data was aggregated for all countries and regions. Relative changes of GDP are calculated relative to the baseline (4x1 with zero carbon price). [1] Luderer, G., Pietzcker, R. C., Bertram, C., Kriegler, E., Meinshausen, M. and Edenhofer, O.: Economic mitigation challenges: how further delay closes the door for achieving climate targets, Environmental Research Letters, 8(3), 034033, doi:10.1088/1748-9326/8/3/034033, 2013a. [2] Rogelj, J., Luderer, G., Pietzcker, R. C., Kriegler, E., Schaeffer, M., Krey, V. and Riahi, K.: Energy system transformations for limiting end-of-century warming to below 1.5 °C, Nature Climate Change, 5(6), 519–527, doi:10.1038/nclimate2572, 2015.

Climate and vegetation change in a coastal marsh: two snapshots of groundwater dynamics and tidal flooding at Piermont Marsh, NY spanning 20 years We include water levels measured along a transect of groundwater wells in 1999 and 2019, statistical analyses of ground water data, tidal efficiency estimates, vegetation data from 1997, 2005, 2014, and 2018, measures of tide gauge data and sea level rise from the Battery, New York Harbor. We attach the following three groups of files: (1) Files related to data from Piermont Marsh, which includes water levels in wells, tide gauge data collected from the tidal channel, and vegetation data; (2) Files related to analysis of water levels at Piermont Marsh; (3) Files related to analysis of Battery tide gauge data, Battery tide predictions, and precipitation data ## Description of the data and file structure **(1) Files related to data from Piermont Marsh, which includes water levels in wells, tide gauge data collected from the tidal channel, and vegetation data** 1999PiermontWaterlevels.csv 2019PiermontWaterLevels.csv channel_1999.xls channel_2019.xls water_level_elevations.csv Vegetation.xls 1999PiermontWaterlevels.csv and 2019PiermontWaterLevels.csv - Water levels collected at Piermont marsh in groundwater wells, at 0-m, 6-m, 12-m, 18-m, 24-m, 36-m, and 48-m from a tidal channel. The files contain three fields: daytime, well, and elevation. The daytime is the date and time the water level was collected, hours in Eastern Daylight Time -4GMT. The well number refers to its location relative to the tidal channel, with #1 referring to 0-m, #2 referring to 6-m, #3 referring to 12-m, #4 referring to 18-m, #5 referring to 24-m, #6 referring to 36-m, and #7 referring to 48-m. The elevation field refers to the water level in meters relative to the NAVD88 datum. In 1999 water levels were collected 14 April 2019 - 26 May 2019. In 2019, water levels were collected 5 May 2019 - 30 June 2019. channel_1999.xls - This file shows the elevation of water level in the channel. There is a field for date and time, in GMT -4, and water level in meters relative to NGVD29. channel_2019.xls - This file shows the elevation of water level in the channel. There is a field for Date, Time, in GMT -4, absolute pressure in in mbar, temperature in degrees C, and water level in meters relative to NAVD88. water_level_elevations.csv - This csv file includes five fields. The first is "year" or the year collected (1999 or 2019). The second is "well" numbered 1-7. Well 1 is closest to the channel while 7 is the furthest from the channel. #1 referrs to 0-m from the channel, #2 referring to 6-m from the channel, #3 referring to 12-m from the channel, #4 referring to 18-m from the channel, #5 referring to 24-m from the channel, #6 referring to 36-m from the channel, and #7 referring to 48-m from the channel. The datetime field refers to the day and time the measure was made in day/month/year HH:MM AM/PM format. The next field is lunarcyle which refers to whether the measure was made during "spring" or "neap" tidal cycles. Spring was assigned to the tides the week of full or new moons, Neap was assigned to tides the week of the first and last quarter. The last is "elevation" and is the measure of water levels in meters relative to the NAVD88 datum. Vegetation.xls - This Excel file includes four sheets that each refer to a year of vegetation date - 1997, 2005, 2014, and 2017. The first field is "well" which has a number 1 through 7. The well number refers to its location relative to the tidal channel, with #1 referring to 0-m, #2 referring to 6-m, #3 referring to 12-m, #4 referring to 18-m, #5 referring to 24-m, #6 referring to 36-m, and #7 referring to 48-m. There is a field for latitude (lat) and longitude (long), which refers to the location of the shape in UTM, in meters, in the 18N. Cover refers to the plant cover type, area is the area of the polygon in square meters. **(2) Files related to analysis of water levels at Piermont Marsh** Distancefromsurface.R MinNeap_MarshSurface.csv MaxNeap_MarshSurface.csv MinSpring_MarshSurface.csv MaxSpring_MarshSurface.csv PiermontEfficiencyRggplot.csv Tidalefficiency.R The R file Distancefromsurface.R includes calculations of mean and variance of water levels, and as well as production of relevant figures. MinNeap_MarshSurface.csv file has low tide minimum water levels during neap tides (weeks centered on the moons first and third quarter). It includes the following fields: distance, year, water_elevation, marsh_elevation, and distance_surface. The field distance, is distance from the tidal channel, in meters. The field year, refers to is the year collected (1999 or 2019). The field water_elevation, is the elevation of the water level at low tide, in meters relative to the NGVD88 datum. The field marsh_elevation refers to the height of the marsh at that location, in meters relative to the NGVD88 datum. The field distance_surface is the difference between the marsh elevation and the water elevation. Positive values are values below the marsh surface, while negative values are values above the marsh surface. MaxNeap_MarshSurface.csv file has high tide maximum water levels during neap tides (weeks centered on the moons first and third quarter). It includes the following fields: distance, year, water_elevation, marsh_elevation, and distance_surface. The field distance, is distance from the tidal channel, in meters. The field year, refers to is the year collected (1999 or 2019). The field water_elevation, is the elevation of the water level at high tide, in meters relative to the NGVD88 datum. The field marsh_elevation refers to the height of the marsh at that location, in meters relative to the NGVD88 datum. The field distance_surface is the difference between the marsh elevation and the water elevation. Positive values are values below the marsh surface, while negative values are values above the marsh surface. MinSpring_MarshSurface.csv file has low tide minimum water levels during spring tides (weeks centered on the new and full moon). It includes the following fields: distance, year, water_elevation, marsh_elevation, and distance_surface. The field distance, is distance from the tidal channel, in meters. The field year, refers to is the year collected (1999 or 2019). The field water_elevation, is the elevation of the water level at low tide, in meters relative to the NGVD88 datum. The field marsh_elevation refers to the height of the marsh at that location, in meters relative to the NGVD88 datum. The field distance_surface is the difference between the marsh elevation and the water elevation. Positive values are values below the marsh surface, while negative values are values above the marsh surface. MaxSpring_MarshSurface.csv file has high tide maximum water levels during spring tides (weeks centered on the new and full moon). It includes the following fields: distance, year, water_elevation, marsh_elevation, and distance_surface. The field distance, is distance from the tidal channel, in meters. The field year, refers to is the year collected (1999 or 2019). The field water_elevation, is the elevation of the water level at high tide, in meters relative to the NGVD88 datum. The field marsh_elevation refers to the height of the marsh at that location, in meters relative to the NGVD88 datum. The field distance_surface is the difference between the marsh elevation and the water elevation. Positive values are values below the marsh surface, while negative values are values above the marsh surface. PiermontEfficiencyRggplot.csv - file lists the well number (1-7), distance (a number 1-14, which gives a unique identifier to each combination of well and year), year, which was the year the data was collected. The last field is efficiency. This field refers to the ratio between the change in water level over the course of a tidal cycle in the well to the change in the water level over the course of the tidal cycle at the Battery tide gauge, NYC. Tidalefficiency.R - file that plots and calculates tidal efficiency during 1999 and 2019 at each well. **(3) Files related to analysis of Battery tide gauge data, Battery tide predictions, and precipitation data** MSL_time.R 3348871.csv 3348873.csv Battery.csv Bat_wls.csv monthly.csv sin2.csv predictions.csv tide_l.csv wls.csv MSL_time.R - This R code uses several data files to conduct analysis of change over time in water levels and monthly anomalies in precipitation and water levels. All necessary packages are described. 3348871.csv and 3348873.csv - are weather data from Westchester County airport, station USW00094745 from 1997 to 2001 (3348873.csv) 2017 to 2022 (3348871.csv). The field station lists the station. The field Name is the name of the station, Westchester County Airport. The date is the day data was collected. AWND refers to Average daily wind speed in miles per hour. PGTM refers to peak gust time (hours and minutes, i.e., HHMM). PRCP refers to precipitation in inches, TMAX refers to the maximum daily temperature, in degrees Fahrenheit. TMIN refers to the minimum daily temperature, in degrees Fahrenheit. WDF2 is the direction of fastest 2-minute wind in degrees. WDF5 is the direction of fastest 5-second wind in degrees. WSF2 is the fastest 2-minute wind speed in miles per hour. WSF5 is the fastest 5-second wind speed in miles per hour. Missing data is replaced with -999. Battery.csv - all high tide levels for 1997 through 2022. The two fields are level, referring to high tide water levels in meters relative to the NAVD88 datum. The second field is year. Bat_wls.csv is monthly tide levels from the Battery tide gauge, NY. The year field refers to year including fraction. Mean high water (MHW) refers to monthly mean high water relative to the NAVD88 datum in meters. Mean sea level (MSL) refers to monthly mean sea level relative to the NAVD88 datum in meters. Mean tide level (MTL) refers to monthly mean tide level relative to the NAVD88 datum in meters.. Mean Low Water (MLW) refers to monthly mean low water relative to the NAVD88 datum in meters. monthly.csv - is mean high water and mean sea level from 1980-2022, by month. The field month refers to the month (January =1). MHW is monthly mean high water for all months, relative to the NAVD88 datum, and MSL is monthly mean sea level relative to the NAVD88 datum. sin2.csv is the monthly mean sea level at the Battery tide gauge (1980-2022), with a 1 year rolling window median smooth added. There are three fields, month, MSL, and year. Month is the number of months elapsed since January 1961. MSL is the monthly mean sea level in meters, relative to the NAVD88 datum, with a one year smoothing function applied. Year refers to the observation month, expressed in years and the fraction of years so January 1980 would be 1980, while February 1980 is depicted as 1980.083. predictions.csv - tide predictions for the Battery tide gauge, New York City. Fields are y, which stands for year, represented by year, including fractions representing months. High_p is the highest predicted tide of the month, in meters relative to the NAVD88 datum. MHW_p is the predicted mean high tide for the month relative to the NAVD88 datum. MLW_p is the predicted mean low tide for the month relative to the NAVD88 datum. MTL_p is the predicted mean tide level for the month relative to the NAVD88 datum. High_1 is the highest actual tide of the month, in meters relative to the NAVD88 datum. MHW_a is the actual mean high tide for the month relative to the NAVD88 datum. MLW_a is the actual mean low tide for the month relative to the NAVD88 datum. MTL_a is the actual mean tide level for the month relative to the NAVD88 datum. tide_l.csv is a file with the monthly mean high water (MHW_l), monthly mean tide level (MTL_l), and mean low water (MLW_l) for 1960 -2021. wls.csv is a file that has monthly water levels from 1999 to 2019, listing year (as a fraction, not just an integer for month), Highest, as the highest tide of the month in meters relative to the NAVD88 datum. MHW refers to the mean high water during the month in meters relative to the NAVD88 datum. MTL refers to the mean tidal level during the month in meters relative to the NAVD88 datum. MLW refers to the mean low water during the month in meters relative to the NAVD88 datum. ## Sharing/Access information Data was derived from the following external sources: * Vegetation shapefiles for the Hudson River NERR for 1997, 2005, and 2014, were obtained through personal request to Sarah Fernald, *Reserve Manager and Research Coordinator.* Files should be available through the Reserve website, although the link is not functional at this time: * The 2018 vegetation shapefiles were obtained from under the heading, [Hudson River Estuary tidal wetlands](https://data.gis.ny.gov/datasets/ee2723393f894e929dbd6dbdc84770de_0/explore?location=41.308770%2C-73.842410%2C9.10). * We acknowledge the NYS DEC Hudson River Estuary Program, NYS DEC Hudson River National Estuarine Research Reserve, and Cornell Institute for Resource Information Sciences for collection and curation of the Hudson River NERR vegetation data. * Tide gauge data and tide predictions for the Battery, NY were obtained from NOAA tides and currents website: * Precipitation data was obtained from the National Centers for Environmental Information, NOAA: . The station for which data was obtained was the Westchester County airport, station USW00094745. ## Code/Software We provide three R files, which we ran using R version 4.3.1 (2023-06-16), in R Studio 2022.02.1, Build 461. All required packages are described in the .R files. Distancefromsurface.R - This R code utilizes four data files that include low tides during spring tides, low tides during neap tides, high tides during spring tides, and high tides during neap files to compare average and variance in low and high tide water levels during 1999 and 2019 relative to the marsh surface and relative to the NAVD88 datum. Code is also included to produce plots. Tidalefficiency.R - file that plots and calculates tidal efficiency during 1999 and 2019 at each well. MSL_time.R - This R code uses several data files to conduct analysis of change over time in water levels and monthly anomalies in precipitation and water levels. Hydrological measurements were collected during the spring and summer of 1999 and 2019 in Piermont Marsh (coordinates 41.0361°, -73.9105°). These measurements covered a transect that was laid out perpendicular to a tidal channel. The objective of this study was to compare the current tidal flooding and groundwater table levels with the data from 1999. The goal was to assess the differences in tidal hydrology between these two distinct time periods, which also differed in terms of marsh and water level elevations. To determine groundwater levels and tidal flooding across the marsh, we installed seven water level loggers along a gradient, ranging from the tidal channel to the upland area. We constructed wells by suspending pressure transducers within 7.5 cm diameter perforated PVC pipes lined with screening to prevent sediment from entering the well. These wells were positioned one meter below the marsh surface, 0.6 meters above the soil surface, vented to the atmosphere, and only the section below the soil surface was perforated. Additionally, we installed concrete collars at the marsh surface around the wells to prevent preferential water flow down the well sides. These seven wells were placed along the original transect, perpendicular to the creek, with increasing distances (0 meters, 6 meters, 12 meters, 18 meters, 24 meters, 36 meters, and 48 meters). We installed and monitored the wells from May 5 to June 30, 2019, and from April 6 to May 26, 1999. In 2019, we measured the absolute elevation of the top of each well using RTK-enabled static GPS measurements from Leica GNSS GS14 rover units and static measurements with an AX1202 GG base station unit to reference water levels to the NAVD88 vertical datum. We measured reference water levels each time data was collected, which involved determining the distance from the top of the well to the water surface and converting it to elevation relative to the NAVD88 datum. To relate marsh elevation to water elevations, GPS surveys were conducted along the transect using a Leica GNSS GS14 rover unit. In 1999, elevation control for the wells and water levels was similarly measured using survey-grade GPS. We compared changes in the marsh water table with significant potential hydrological and vegetation changes that have occurred over the past 20 years. We calculated the rates of change in monthly water levels at Battery, NY for the period from 1999 to 2019 using two different methods. We modeled changes over time in monthly highest water levels, mean high water (MHW), mean tide level (MTL), and mean low water (MLW) using an ordinary least squares regression model with ARIMA errors to account for the autoregressive structure of tide data. We removed the annual cycle first using a curve with a 1-year periodicity. The ARIMA errors model was fitted using the "auto.arima" function from the "forecast" package. We calculated the squared correlation of fitted values to actual values to produce a pseudo-r2. For comparison, we calculated trends using ordinary least squares regression for the 1999-2019 period, although it's important to note that the temporal autocorrelation likely results in underestimated uncertainty. We obtained vegetation maps from the HRNERR for 1997, 2005, 2014, and 2018 to help assess changes in the coverage of plant species over time, as these changes could impact evapotranspiration and water table patterns. A 20-meter buffer zone was created around each well location, and the composition of vegetation within this buffer zone was quantified using QGIS version 3.30.2. While four time-points may not be sufficient for statistically identifying trends, we analyzed the changes observed. To put the measurement time periods in context and ensure that our selected seasons were not anomalous, we compared water levels in spring 1999 and 2019 relative to the astronomical cycles driving interannual sea level variability using data from the Battery tide gauge. We also compared spring high tide levels in 1999 and 2019 with surrounding years. The main astronomical cycles thought to influence tides include the 18.6-year lunar nodal cycle and the 4.4-year subharmonic of the 8.85-year lunar perigee cycle. As our 1999 and 2019 measurements were collected during slightly different time periods (April/May 1999 vs. May/June 2019), we also examined mean monthly water levels (1980-2022) from the NOAA Battery tidal gauge to identify potential artifacts. We obtained rainfall data from spring 1999 and 2019 from the nearest precipitation monitoring station (Westchester airport) to determine whether the measurements were made during an unusually wet or dry period. The sampling periods were 20 years apart, so they occurred at approximately the same point in the 18.6-year lunar nodal cycle. Pressure transducer data was processed using HOBOware Pro (Version 3.7.16, Onset Computer Corporation, Bourne, MA) with reference water levels collected in the field. The data were corrected for atmospheric pressure using the HOBOware barometric compensation assistant, using data from the Hudson River National Estuarine Research Reserve. Raw water elevation data from 1999 was analyzed in conjunction with the 2019 data. Water level data from 1999 were converted from the NVGD29 to NAVD 88 datum using NOAA VDatum v4.0.1 prior to analysis. Well seven's transducer experienced three brief malfunctions from May 30 to June 3, 2019, resulting in inaccurate elevation measurements for a total of 19.5 hours. These data were excluded from the analysis. In 1999, well seven also experienced malfunctions, which were corrected by Montalto into smoothed six-hour increments using average water elevation measurements and calculated error, calibrated using regression. No other well transducers appeared to have malfunctioned. Groundwater hydrology plays an important role in coastal marsh biogeochemical function, in part because groundwater dynamics drive the zonation of macrophyte community distribution. Changes that occur over time, such as sea level rise and shifts in habitat structure are likely altering groundwater dynamics and eco-hydrological zonation. We examined tidal flooding and marsh water table dynamics in 1999 and 2019 and mapped shifts in plant distributions over time, at Piermont Marsh, a brackish tidal marsh located along the Hudson River Estuary near New York City. We found evidence that the marsh surface was flooded more frequently in 2019 than in 1999, and that tides were propagating further into the marsh in 2019, although marsh surface elevation gains were largely matching that of sea level rise. The changes in groundwater hydrology that we observed are likely due to the high tide rising at a rate that is greater than that of mean sea level. In addition, we reported on changes in plant cover by P. australis, which has displaced native marsh vegetation at Piermont Marsh. Although P. australis has increased in cover, wrack deposition and plant die off associated Superstorm Sandy allowed for native vegetation to rebound in part of our focus area. These results suggest that climate change and plant community composition may interact to shape ecohydrologic zonation. Considering these results, we recommend that habitat models consider tidal range expansion and groundwater hydrology as metrics when predicting the impact of sea level rise on marsh resilience.

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.

This dataset contains the underlying data for the following publication: Doukas, H., Spiliotis, E., Jafari, M. A., Giarola, S. & Nikas, A. (2021). Low-cost emissions cuts in container shipping: Thinking inside the box. Transportation Research Part D: Transport and Environment, 94, 102815, https://doi.org/10.1016/j.trd.2021.102815.

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.

Leaching dataset of dissolved organic carbon (DOC) and nitrogen (DON), nitrate (NO3+) and ammonium (NH4+) were collected from 6 cropping treatments (corn, switchgrass, miscanthus, native grass mix, restored prairie and poplar) established in the Bioenergy Cropping System Experiment (BCSE) which is a part of Great Lakes Bioenergy Research Center (www.glbrc.org) and Long Termn Ecological Research (LTER) program (www.lter.kbs.msu.edu). The site is located at the W.K. Kellogg Biological Station (42.3956° N, 85.3749° W and 288 m above sea level), 25 km from Kalamazoo in southwestern Michigan, USA. Prenart soil water samplers made of Teflon and silica (http://www.prenart.dk/soil-water-samplers/) were installed in blocks 1 and 2 of the BCSE (Fig. S1), and Eijkelkamp soil water samplers made of ceramic (http://www.eijkelkamp.com) were installed in blocks 3 and 4 (there were no soil water samplers in block 5). All samplers were installed at 1.2 m depth at a 45° angle from the soil surface, approximately 20 cm into the unconsolidated sand of the 2Bt2 and 2E/Bt horizons. Beginning in 2009, soil water was sampled at weekly to biweekly intervals during non-frozen periods (April to November) by applying 50 kPa of vacuum for 24 hours, during which water was collected in glass bottles. During the 2009 and 2010 sampling periods we obtained fewer soil water samples from blocks 1 and 2 where Prenart lysimeters were installed. We observed no consistent differences between the two sampler types in concentrations of the analytes reported here. Depending on the volume of leachate collected, water samples were filtered using either 0.45 µm pore size, 33-mm-dia. cellulose acetate membrane filters when volumes were <50 ml, or 0.45 µm, 47-mm-dia. Supor 450 membrane filters for larger volumes. Samples were analyzed for NO3-, NH4+, total dissolved nitrogen (TDN), and DOC. The NO3- concentration was determined using a Dionex ICS1000 ion chromatograph system with membrane suppression and conductivity detection; the detection limit of the system was 0.006 mg NO3--N L-1. The NH4+ concentration in the samples was determined using a Thermo Scientific (formerly Dionex) ICS1100 ion chromatograph system with membrane suppression and conductivity detection; the detection limit of the system was similar. The DOC and TDN concentrations were determined using a Shimadzu TOC-Vcph carbon analyzer with a total nitrogen module (TNM-1); the detection limit of the system was ~0.08 mg C L-1 and ~0.04 mg N L-1. DON concentrations were estimated as the difference between TDN and dissolved inorganic N (NO3- + NH4+) concentrations. The NH4+ concentrations were only measured in the 2013-2015 crop-years, but they were always small relative to NO3- and thus their inclusion or lack of it was inconsequential to the DON estimation. Leaching rates were estimated on a crop-year basis, defined as the period from planting or emergence of the crop in the year indicated through the ensuing year until the next year’s planting or emergence. For each sampling point, the concentration was linearly interpolated between sampling dates during non-freezing periods (April through November). The concentrations in the unsampled winter period (December through March) were also linearly interpolated based on the preceding November and subsequent April samples. Solute leaching (kg ha-1) was calculated by multiplying the daily solute concentration in pore-water (mg L -1) by the modeled daily drainage rates (m3 ha-1) from the overlying soil. The drainage rates were obtained using the SALUS (Systems Approach for Land Use Sustainability) model (Basso and Ritchie, 2015). SALUS simulates yield and environmental outcomes in response to weather, soil, management (planting dates, plant population, irrigation, nitrogen fertilizer application, tillage), and crop genetics. The SALUS water balance sub-model simulates surface run-off, saturated and unsaturated water flow, drainage, root water uptake, and evapotranspiration during growing and non-growing seasons (Basso and Ritchie, 2015). Drainage amounts and rates simulated by SALUS have been validated with measurements using large monolith lysimeters at a nearby site at KBS (Basso and Ritchie, 2005). On days when SALUS predicted no drainage, the leaching was assumed to be zero. The volume-weighted mean concentration for an entire crop-year was calculated as the sum of daily leaching (kg ha-1) divided by the sum of daily drainage rates (m3 ha-1). Weather data for the model were collected at the nearby KBS LTER meteorological station (lter.kbs.msu.edu). Leaching losses of dissolved organic carbon (DOC) and nitrogen (DON) from agricultural systems are important to water quality and carbon and nutrient balances but are rarely reported; the few available studies suggest linkages to litter production (DOC) and nitrogen fertilization (DON). In this study we examine the leaching of DOC, DON, NO3-, and NH4+ from no-till corn (maize) and perennial bioenergy crops (switchgrass, miscanthus, native grasses, restored prairie, and poplar) grown between 2009 and 2016 in a replicated field experiment in the upper Midwest U.S. Leaching was estimated from concentrations in soil water and modeled drainage (percolation) rates. DOC leaching rates (kg ha-1 yr-1) and volume-weighted mean concentrations (mg L-1) among cropping systems averaged 15.4 and 4.6, respectively; N fertilization had no effect and poplar lost the most DOC (21.8 and 6.9, respectively). DON leaching rates (kg ha-1 yr-1) and volume-weighted mean concentrations (mg L-1) under corn (the most heavily N-fertilized crop) averaged 4.5 and 1.0, respectively, which was higher than perennial grasses (mean: 1.5 and 0.5, respectively) and poplar (1.6 and 0.5, respectively). NO3- comprised the majority of total N leaching in all systems (59-92%). Average NO3- leaching (kg N ha-1 yr-1) under corn (35.3) was higher than perennial grasses (5.9) and poplar (7.2). NH4+ concentrations in soil water from all cropping systems were relatively low (<0.07 mg N L-1). Perennial crops leached more NO3- in the first few years after planting, and markedly less after. Among the fertilized crops, the leached N represented 14-38% of the added N over the study period; poplar lost the greatest proportion (38%) and corn was intermediate (23%). Requiring only one third or less of the N fertilization compared to corn, perennial bioenergy crops can substantially reduce N leaching and consequent movement into aquifers and surface waters. readme files are given that describe the data table