This dataset provides estimate of the potential increase in soil organic carbon within the top 30 cm of soil in croplands after 20 years, following implementation of better land managment practices under a high sequestration scenario. The per pixel values here take in to consideration the percent of each pixel which is classified as cropland (from the GLC-Share/GLC-02 dataset), and values have been converted to total tonnes of carbon (x 100) per pixel.<br/><br>See: <a href="https://doi.org/10.1038/s41598-017-15794-8">Zomer, R.J., Bossio, D.A., Sommer, R., Verchot, L.V., 2017. Global Sequestration Potential of Increased Organic Carbon in Cropland Soils. Scientific Reports 7, 15554</a>.<br/>For descriptions of sequrestion scenarions see: <a href="https://doi.org/10.1016/j.jenvman.2014.05.017">Sommer, R., Bossio, D., 2014. Dynamics and climate change mitigation potential of soil organic carbon sequestration. Journal of Environmental Management 144, 83–87</a>.<br/>

This data layer combines estimates of pollution coming from commercial shipping and from ports. As such, it is a combination of the shipping and port volume data layers, with the port volume data plumed to estimate pollution from commercial ports (with exponential decline in intensity from the port). Ocean-based pollution is assumed to derive from commercial and recreational ship activity. No data on global recreational ship activity currently exist, and therefore we modelled this driver to oceans using a combination of the commercial shipping traffic data and port data. The shipping data provide an estimate of the occurrence of ships at a particular location, and therefore an estimate of the amount of pollution they produce (via fuel leaks, oil discharge, waste disposal, etc.) that is unique from their contribution to ship strikes, etc. described above. We recognize that ocean currents can disperse this pollution into untraveled regions, but small-scale oceanography is known for only a few select locations around the world, and pollutants are likely to be most concentrated in high traffic areas. The dispersal of port-derived pollution was modelled as a diffusive plume with a maximum distance of 100 km. These plumes were not clipped to shallow regions as was done for Invasive Species.Raw stressor data from "Benjamin Halpern, Melanie Frazier, John Potapenko, Kenneth Casey, Kellee Koenig, et al. 2015. Cumulative human impacts: raw stressor data (2008 and 2013). Knowledge Network for Biocomplexity. doi:10.5063/F1S180FS."

Non-point source inorganic pollution was modelled with global 1 km2 impervious surface area data (http://www.ngdc.noaa.gov/dmsp/) under the assumption that most of this pollution comes from urban runoff. These data will not capture point-sources of pollution or nonpoint sources where paved roads do not exist (e.g., select places in developing countries). These values were then aggregated to the watershed and distributed to the pour point (i.e., stream and river mouths) for the watershed with raster statistics (i.e., aggregation by watershed). Finally, spread of the driver values into coastal waters at each pour point was modelled with a cost-path surface on the basis of a decay function that assigns a fixed amount of the driver (0.5% of the value in the previous cell) in the initial cell and then evenly distributes the remaining amount of driver in all adjacent and ‘unvisited’ cells, repeated until a minimum threshold (0.05% of global maximum) is reached. This approach to modelling river plumes is diffusive and so allows drivers to wrap around headlands and islands. Raw stressor data from "Benjamin Halpern, Melanie Frazier, John Potapenko, Kenneth Casey, Kellee Koenig, et al. 2015. Cumulative human impacts: raw stressor data (2008 and 2013). Knowledge Network for Biocomplexity. doi:10.5063/F1S180FS."

Input data are from 2007-2010 FAO statistics, the most recent available. FAO data on annual country-level fertilizer were used were available, averaged over the time periods. Missing values were filled using a linear regression model of fertilizers as a function of pesticides (gaps: N=4; regression: R2 = 0.72) when pesticide data were available or agricultural GDP (gaps: N=22; regression: R2 = 0.62) when not. These country-level average fertilizer values were then dasymetrically distributed over a country’s landscape using global land cover data from 2009, derived from the Moderate Resolution imaging Spectroradiometer (MODIS) instrument at ~500m resolution. Finally, spread of the driver values into coastal waters at each pour point was modelled with a cost-path surface on the basis of a decay function that assigns a fixed amount of the driver (0.5% of the value in the previous cell) in the initial cell and then evenly distributes the remaining amount of driver in all adjacent and ‘unvisited’ cells, repeated until a minimum threshold (0.05% of global maximum) is reached. This approach to modelling river plumes is diffusive and so allows drivers to wrap around headlands and islands. Raw stressor data from "Benjamin Halpern, Melanie Frazier, John Potapenko, Kenneth Casey, Kellee Koenig, et al. 2015. Cumulative human impacts: raw stressor data (2008 and 2013). Knowledge Network for Biocomplexity. doi:10.5063/F1S180FS."

Anomaly for the period 2041-2060 compared to climatological data (1979-2013) on precipitation and temperature data based on two different scenarios (RCP4.5 and RCP8.5). The layer is calculated at UNEP/GRID-Geneva from the layers on annual mean temperature and annual precipitations provided in the products CHELSA V1.2 and CHELSA-[CMIP5]. CHELSA-[CMIP5] is a delta change climatological dataset for the years 2041-2060 and 2061- 2080 for mean monthly maximum temperatures, mean monthly minimum temperatures, monthly precipitation amounts, and several derived parameters. We use the delta change method by B-spline interpolation of anomalies (deltas) of the respective CMIP5 GCM dataset. Anomalies were interpolated between all CMIP5 grid cells and are then added (for temperature variables) or multiplied (in case of precipitation) to high resolution climate data from CHELSA V1.2. This method has the assumption that climate only varies on the scale of the coarser (CMIP5) dataset, and the spatial pattern (from CHELSA) is consistent over time. CHELSA- [CMIP5] does not take changing wind patterns, or temperature lapse rates into account, but rather expects them to be constant over time, and similar to the long term averages. CHELSA V1.2 (http://chelsa-climate.org/) is a high resolution (30 arc sec, ~1 km) climate data set for the earth land surface areas. It includes monthly and annual mean temperature and precipitation patterns for the time period 1979-2013. Methods are described in http://chelsa-climate.org/wp-admin/download-page/CHELSA_tech_specification.pdf. CHELSA Version 1.2 is licensed under a Creative Commons Attribution 4.0 International License. Specifications: High resolution (30 arcsec, ~1 km) Precipitation & Temperature Climatologies for the years 1979 – 2013 Incorporation of topoclimate (e.g. orographic rainfall & wind fields). All products of CHELSA are in a geographic coordinate system referenced to the WGS 84 horizontal datum, with the horizontal coordinates expressed in decimal degrees. The CHELSA layer extents (minimum and maximum latitude and longitude) are a result of the coordinate system inherited from the 1-arc-second GMTED2010 data which itself inherited the grid extent from the 1-arc-second SRTM data. Note that because of the pixel center referencing of the input GMTED2010 data the full extent of each CHELSA grid as defined by the outside edges of the pixels differs from an integer value of latitude or longitude by 0.000138888888 degree (or 1/2 arc-second). Users of products based on the legacy GTOPO30 product should note that the coordinate referencing of CHELSA (and GMTED2010) and GTOPO30 are not the same. In GTOPO30, the integer lines of latitude and longitude fall directly on the edges of a 30-arc-second pixel. Thus, when overlaying CHELSA with products based on GTOPO30 a slight shift of 1/2 arc-second will be observed between the edges of corresponding 30-arc-second pixels. To redistribute the data, please cite the following peer reviewed articles: <a href="https://www.nature.com/articles/sdata2017122"target=_blank>Karger, D.N., Conrad, O., Böhner, J., Kawohl, T., Kreft, H., Soria-Auza, R.W., Zimmermann, N.E., Linder, H.P. & Kessler, M. (2017) Climatologies at high resolution for the earth’s land surface areas. Scientific Data 4, 170122.</a> <a href="https://doi.org/10.5061/dryad.kd1d4"target=_blank>Karger, D.N., Conrad, O., Böhner, J., Kawohl, T., Kreft, H., Soria-Auza, R.W., Zimmermann, N.E., Linder, H.P., Kessler, M. (2017) Data from: Climatologies at high resolution for the earth’s land surface areas. Dryad Digital Repository. </a> CHELSA – Climatologies at high resolution for the Earth land surface areas. Version 1.2

CHELSA V1.2 (http://chelsa-climate.org/) is a high resolution (30 arc sec, ~1 km) climate data set for the earth land surface areas. It includes monthly and annual mean temperature and precipitation patterns for the time period 1979-2013. Methods are described in http://chelsa-climate.org/wp-admin/download-page/CHELSA_tech_specification.pdf. CHELSA Version 1.2 is licensed under a Creative Commons Attribution 4.0 International License. Specifications: High resolution (30 arcsec, ~1 km) Precipitation & Temperature Climatologies for the years 1979 – 2013 Incorporation of topoclimate (e.g. orographic rainfall & wind fields). All products of CHELSA are in a geographic coordinate system referenced to the WGS 84 horizontal datum, with the horizontal coordinates expressed in decimal degrees. The CHELSA layer extents (minimum and maximum latitude and longitude) are a result of the coordinate system inherited from the 1-arc-second GMTED2010 data which itself inherited the grid extent from the 1-arc-second SRTM data. Note that because of the pixel center referencing of the input GMTED2010 data the full extent of each CHELSA grid as defined by the outside edges of the pixels differs from an integer value of latitude or longitude by 0.000138888888 degree (or 1/2 arc-second). Users of products based on the legacy GTOPO30 product should note that the coordinate referencing of CHELSA (and GMTED2010) and GTOPO30 are not the same. In GTOPO30, the integer lines of latitude and longitude fall directly on the edges of a 30-arc-second pixel. Thus, when overlaying CHELSA with products based on GTOPO30 a slight shift of 1/2 arc-second will be observed between the edges of corresponding 30-arc-second pixels. To redistribute the data, please cite the following peer reviewed articles: <a href="https://www.nature.com/articles/sdata2017122"target=_blank>Karger, D.N., Conrad, O., Böhner, J., Kawohl, T., Kreft, H., Soria-Auza, R.W., Zimmermann, N.E., Linder, H.P. & Kessler, M. (2017) Climatologies at high resolution for the earth’s land surface areas. Scientific Data 4, 170122.</a> <a href="https://doi.org/10.5061/dryad.kd1d4"target=_blank>Karger, D.N., Conrad, O., Böhner, J., Kawohl, T., Kreft, H., Soria-Auza, R.W., Zimmermann, N.E., Linder, H.P., Kessler, M. (2017) Data from: Climatologies at high resolution for the earth’s land surface areas. Dryad Digital Repository. </a>

Vessel identity and location information was obtained using two approaches. (1) Over the past 20 years, 10-20% of the vessel fleet has voluntarily participated in collecting meteorological data for the open ocean, which includes location at the time of measurement, as part of the Volunteer Observing System (VOS). (2) In order to improve maritime safety, in 2002 the International Maritime Organization SOLAS agreement required all vessels over 300 gross tonnage (GT) and vessels carrying passengers to equip Automatic Identification System (AIS) transceivers, which use the Global Positioning System (GPS) to precisely locate vessels. Eight broad classes of vessels were taken into account separately: authority, cargo, fishing, high-speed, passenger, pleasure, support, tanker and an ‘other’ class. The vessel classes which move globally (cargo, tanker, and passenger) are required to carry AIS transceivers, and in these three classes 60-70% of the total vessel fleet was observed using AIS. The resulting data layer is primarily composed of these vessel classes in both the AIS and VOS data sources, and is almost exclusively these ship types in the open ocean. We used a simple linear average of the two data sources, producing a final model resolved for the whole ocean at a resolution of 0.1 decimal degrees (~11km). Data have limited observation frequency, leading to gaps that when directly interpolated with geodesic paths, create invalid routes which cross land masses. Routing model was used to create a visibility graph of the oceans, creating valid potential movement paths. These movement paths are based on the assumption that mariners will prefer great circle distances when possible. Raw stressor data from "Benjamin Halpern, Melanie Frazier, John Potapenko, Kenneth Casey, Kellee Koenig, et al. 2015. Cumulative human impacts: raw stressor data (2008 and 2013). Knowledge Network for Biocomplexity. doi:10.5063/F1S180FS."

Vessel identity and location information was obtained using two approaches. (1) Over the past 20 years, 10-20% of the vessel fleet has voluntarily participated in collecting meteorological data for the open ocean, which includes location at the time of measurement, as part of the Volunteer Observing System (VOS). (2) In order to improve maritime safety, in 2002 the International Maritime Organization SOLAS agreement required all vessels over 300 gross tonnage (GT) and vessels carrying passengers to equip Automatic Identification System (AIS) transceivers, which use the Global Positioning System (GPS) to precisely locate vessels. A single year sample of the VOS data was used for analysis. These data ignores vessel type, and included observations from only 12% of the vessel fleet. The ships included are a spatially- and statistically-biased sample of the population, making the modelled results somewhat misleading. Data have limited observation frequency, leading to gaps that when directly interpolated with geodesic paths, create invalid routes which cross land masses. Routing model was used to create a visibility graph of the oceans, creating valid potential movement paths. These movement paths are based on the assumption that mariners will prefer great circle distances when possible.Raw stressor data from "Benjamin Halpern, Melanie Frazier, John Potapenko, Kenneth Casey, Kellee Koenig, et al. 2015. Cumulative human impacts: raw stressor data (2008 and 2013). Knowledge Network for Biocomplexity. doi:10.5063/F1S180FS."

Input data are from FAO statistics published in the time period 2003-2006. FAO data on annual country-level fertilizer were used were available, averaged over the time periods. Missing values were filled using a linear regression model of fertilizers as a function of pesticides (gaps: N=4; regression: R2 = 0.72) when pesticide data were available or agricultural GDP (gaps: N=22; regression: R2 = 0.62) when not. These country-level average fertilizer values were then dasymetrically distributed over a country’s landscape using global land cover data from 2005. Finally, spread of the driver values into coastal waters at each pour point was modelled with a cost-path surface on the basis of a decay function that assigns a fixed amount of the driver (0.5% of the value in the previous cell) in the initial cell and then evenly distributes the remaining amount of driver in all adjacent and ‘unvisited’ cells, repeated until a minimum threshold (0.05% of global maximum) is reached. This approach to modelling river plumes is diffusive and so allows drivers to wrap around headlands and islands. Raw stressor data from "Benjamin Halpern, Melanie Frazier, John Potapenko, Kenneth Casey, Kellee Koenig, et al. 2015. Cumulative human impacts: raw stressor data (2008 and 2013). Knowledge Network for Biocomplexity. doi:10.5063/F1S180FS."

This data layer combines estimates of pollution coming from commercial shipping and from ports. As such, it is a combination of the shipping and port volume data layers, with the port volume data plumed to estimate pollution from commercial ports (with exponential decline in intensity from the port). Ocean-based pollution is assumed to derive from commercial and recreational ship activity. No data on global recreational ship activity currently exist, and therefore we modelled this driver to oceans using a combination of the commercial shipping traffic data and port data. The shipping data provide an estimate of the occurrence of ships at a particular location, and therefore an estimate of the amount of pollution they produce (via fuel leaks, oil discharge, waste disposal, etc.) that is unique from their contribution to ship strikes, etc. described above. We recognize that ocean currents can disperse this pollution into untraveled regions, but small-scale oceanography is known for only a few select locations around the world, and pollutants are likely to be most concentrated in high traffic areas. The dispersal of port-derived pollution was modelled as a diffusive plume with a maximum distance of 100 km. These plumes were not clipped to shallow regions as was done for Invasive Species.Raw stressor data from "Benjamin Halpern, Melanie Frazier, John Potapenko, Kenneth Casey, Kellee Koenig, et al. 2015. Cumulative human impacts: raw stressor data (2008 and 2013). Knowledge Network for Biocomplexity. doi:10.5063/F1S180FS."