Abstract
This study presents a new global gridded dataset of bioclimatic indicators at 0.5° by 0.5° resolution for historical and future conditions. The dataset, called CMCC-BioClimInd, provides a set of 35 bioclimatic indices, expressed as mean values over each time interval, derived from post-processing both climate reanalysis for historical period (1960–1999) and an ensemble of 11 bias corrected CMIP5 simulations under two greenhouse gas concentration scenarios for future climate projections along two periods (2040–2079 and 2060–2099). This new dataset complements the availability of spatialized bioclimatic information, crucial aspect in many ecological and environmental wide scale applications and for several disciplines, including forestry, biodiversity conservation, plant and landscape ecology. The data of individual indicators are publicly available for download in the commonly used Network Common Data Form 4 (NetCDF4) format.
Measurement(s) | bioclimatic indicators |
Technology Type(s) | Geographic Information System • computational modeling technique |
Factor Type(s) | temporal interval • geographic location |
Sample Characteristic - Environment | climate system |
Sample Characteristic - Location | global |
Machine-accessible metadata file describing the reported data: https://doi.org/10.6084/m9.figshare.12927443
Similar content being viewed by others
Background & Summary
Climate change impacts, affecting primarily ecosystems’ functions and consequently human sectors, have become a crucial topic within the scientific community (https://www.ipcc.ch/working-group/wg2/) and most recently across the whole society and sustainable development efforts (https://sustainabledevelopment.un.org/sdg13). Average climate patterns, mainly represented by intra-annual (monthly to seasonal) temperature and precipitation cycle, directly influence the distribution, abundance and interactions of biological species, with consolidated evidence of how observed and expected variations in climate conditions can undermine the ecosystems’ ecological equilibrium1,2,3.
Under faster and faster environmental modifications over lands4, climate datasets and efficient processing chains applied on them allow answering many urgent questions of biogeographical sciences about climate change impacts on living organisms (e.g. animals, plants, bacteria), even through the interactions they have with the surrounding natural resources, like water and soil. For example, to reproduce and model, respectively, current and future habitats’ ranges by mean of Species Distribution Modelling (SDM), large amount of high quality and up-to-date environmental (especially climate) data is required, to be operated into further indicators proxy of typical climate settings across the domains of interest5,6,7,8,9.
During the long history of scientific research on the relationships between climate and Earth’s communities, numerous meteorological variables and/or derived indices have been formulated, calculated and applied to explain the geographic distribution of natural populations along climate gradients, characterized by intra-annual patterns of temperature and precipitation10. Such variables and indices are also known as bioclimatic indicators (hereafter BioClimInd). They mainly result from primary - observed or modelled - climate fields (e.g. minimum, maximum and mean temperature, precipitation amount) and contribute to delineate the bioclimatic “envelope” for species in terms of favourable environmental conditions11,12, also referred to as “suitability”7,9.
Considering the increasing interest for spatially-explicit assessments, the computation of BioClimInd for the historical period can benefit from gridded climate products building on meteorological stations’ data spatialized through either regression-based interpolation or climate reanalysis13. Instead, the primary data source to derive BioClimInd for future time horizons are the outputs of climate model simulations that, despite the substantial progress occurred in the last few decades, are still affected by both systematic and random errors preventing their direct use in climate impact studies without affecting their reliability13,14,15. To reduce these errors and related effects, different approaches can be used14: improving the physics in the models; quantifying the uncertainty by using a multi-model ensemble; or removing the model biases through climate data post-processing16. The first option requires an extensive investment in research efforts while the second one needs a huge computing power. Although lacking a sound physical basis, the last option, also know as “bias-correction”, has become increasingly common in the climate impact community. Thanks to its high implementation simplicity and low computational demand along with growing spatialized databases on global and regional climate observations to be used as benchmark17.
Under the raising availability of climate information, several global datasets of bioclimatic indicators have become popular as reference point for the scientific community. Worthy of mention are undoubtedly WorldClim18,19, CHELSA20,21, CliMond22,23, ecoClimate24, ENVIREM25 and MERRAclim26. WorldClim, the most cited dataset and recently (January 2020) updated to the 2.1 version, provides 19 BioClimInd at very high spatial resolution (up to 30 arc-seconds, ca. 0.9 km at the equator) for recent climate (1970–2000) and also for the future (2021–2040, 2041–2060, 2061–2080, 2081–2100) using 9 General Circulation and Earth System Models (GCMs and ESMs) as part of the Coupled Model Intercomparison Project Phase 6 (CMIP6). Climate time series have been downscaled, bias-corrected and covering the four Shared Socio-economic Pathways (SSPs) that drive simulations in support of the 6th Assessment Report from the Intergovernmental Panel on Climate Change (IPCC AR6) (https://pcmdi.llnl.gov/CMIP6/). CHELSA dataset (in the last version 1.2) has the same spatial resolution and indicators as WorldClim for the climatological period 1979–2013 and future periods 2041–2060 and 2061–2080 for a large set of future simulations under CMIP5 (Coupled Model Intercomparison Project Phase 5); CHELSA substantially differs from WorldClim over the mountain regions especially for rainfall-based indicators since the basic algorithm incorporates further orographic predictors20. Recently, still within the CHELSA project, a high resolution dataset was released with information on precipitation and temperatures from 2006 to 210021. CliMond (in the latest 1.2 version), is a dataset of 40 BioClimInd for current and short to long-term future time frames from 2 GCMs with 10’ and 0.5° spatial resolution (ca. 9.2 and 55 km at the equator, respectively). While CliMond uses also radiation and pan evaporation as input variables, it is based on “old” greenhouse gas concentration scenarios (A1B and A2) supporting the IPCC 4th Assessment Report. ENVIREM provides both 16 BioClimInd at multiple spatial resolutions (up to 30 arc-seconds) from three global climate models for several periods in the past and the code to generate them for the future. EcoClimate provides data starting from Pliocene up to 2080–2100 for the same 19 BioClimInd in WorldClim but at 0.5° spatial resolution and for 9 among GCMs/ESMs from CMIP5 and PMIP3 (Paleoclimate Modelling Intercomparison Project Phase 3). MerraClim was recently added to these datasets and provides the 19 WorldClim indicators at three spatial resolutions (from 2.5 to 10 arc-minutes) for three recent decades (1980s, 1990s and 2000s) starting from hourly data of 2-meter air temperature and specific humidity (instead of precipitation). Even when the above datasets are provided interpolated to very high spatial resolution (few kilometers), the represented climate physical processes remain consistent to native map units (tens of thousands to millions of hectares), so that the datasets are more suitable to serve studies encompassing very broad domains (large regions/continents up to the entire globe) and hosting, in terms of biological organization, generalists species and ecosystems’ mosaics27. In case of investigations over sub-regional to local (e.g. mountain area) scale, some studies looking at species with more restricted habitats or at smaller patch ecosystems started to exploit simulations from regional climate models, like those making part of the Coordinated Regional Downscaling Experiment (CORDEX), their further downscaling and/or similar ensemble initiatives28,29,30,31. Moreover, aiming at CORDEX-based applications at global level would raise the issue of heterogeneity due to different combinations of global and regional models, greenhouse gas concentration scenarios and spatial resolutions (0.11°, 0.22° and/or 0.44°) across the CORDEX domains32. To give an example, only CORDEX domains of North America, Mediterranean and Europe are currently covered at 0.11° (ca. 12 km) resolution. Finally, CORDEX data require bias-correction as well, and only in one case (Europe domain) adjusted data have been made available to the community.
In this work, a global dataset of 35 BioClimInd is presented, named CMCC-BioClimInd, with a spatial resolution of 0.5°, both for the historical period (1960–1999) and for two future time horizons (2040–2079 and 2060–2099) from, respectively, post-processing of climate reanalysis and an ensemble of 11 CMIP5 climate simulations33,34. Simulations were bias-corrected benefiting of the trend-preserving approach adopted within the InterSectoral Impact-Model Intercomparison Project (ISIMIP)35 (see Methods). The CMCC-BioClimInd dataset contributes to widening the availability of spatial information useful to the community by (1) providing an ensemble of bioclimatic indicators for the historical and future time frames (e.g. Figs 1 and 2) (2) adopting models and/or other analysis methods for robust (i.e. taking into account uncertainty) climate change impacts’ assessments, at broad scale and in a wide range of research fields such as wildlife ecology, natural resources’ conservation and management, climate impacts’ mitigation. On the one hand, the exploitation of indicators instead of just raw climate variables enables easier inferring of relationships between the studied topic (species occurrence, resources availability etc.) and the climate regime to support decision for complex systems36; on the other hand, using the ensemble allows considering the variability across simulations due to the different models’ physics and the uncertain future development pathways37.
Methods
Input data
The CMCC-BioClimInd were elaborated, for the historical period, starting from the daily time series of temperature and precipitation available from the Water and Global Change (WATCH) forcing dataset (WFD) (http://www.eu-watch.org/data_availability). The WFD is a twentieth century meteorological dataset based on the European Centre for Medium-range Weather Forecasts (ECMWF) Re-Analysis (ERA-40)38 interpolated to a 0.5° by 0.5° grid, with successive elevation-based correction of surface meteorological variables plus monthly bias-correction from the Climatic Research Unit (CRU) gridded observational dataset39,40. In this dataset, the years from 1960 to 1999 were considered as reference period.
For the future time frame, daily time series of precipitation and temperature are available, thanks to the coordinated climate model simulation experiments in the frame of CMIP533,34. These simulations are forced by multiple Representative Concentration Pathways (RCPs) formulated to support the IPCC AR541,42. It is known, as stated above, how raw climate model simulations are however affected by biases that, to serve impact studies, need to be reduced14. To achieve this, the ISIMIP Fast-Track (FT) initiative (https://www.isimip.org/about/#simulation-rounds) provides bias-corrected daily time series of temperature and precipitation globally downscaled (through bilinear interpolation) at 0.5° grid resolution for five Earth System Models (ESMs) participating in CMIP5: GFDL-ESM2M, HadGEM2-ES, IPSL-CM5A-LR, MIROC-ESM-CHEM, NorESM1-M. Bias-correction is based on a trend-preserving method35 (https://github.com/ISIMIP/BC), developed to adjust the probability distribution over the reference period 1960–1999, while preserving the long-term trend in the data.
For the purpose of this work, this method was applied also to the ESM developed and used by the Foundation Euro-Mediterranean Center on Climate Change (CMCC-CESM43), and whose raw data are available on the CMIP5 data archive (https://esgf-node.llnl.gov/search/cmip5/). For the five different ESMs already bias-corrected within the ISIMIP FT initiative, the RCP4.544 and the RCP8.545 were selected, while projections under the RCP8.5 are available from CMCC-CESM runs.
Bioclimatic indicators
Here, a description of BioClimInd (numbered from 1 to 35 with the prefix “Bio” and with a short name) is provided building on the definitions, in some case slightly revised, of the U.S. Geological Survey (USGS)46, the one in the ANUCLIM User Guide (https://fennerschool.anu.edu.au/files/anuclim61.pdf), and those in47 and/or48,49. Based on the indicator considered, source variables can be the bias-corrected mean temperature (TG), maximum temperature (TX), minimum temperature (TN) or precipitation (Pr), directly from daily series as described in the Section “Input data”. Before computing any indicator, bias-corrected daily series of maximum, minimum and mean temperature were converted from degrees Kelvin to Celsius by subtracting 273.15 and indicated as Tx, Tn and Tg, respectively, while precipitation was converted from kg/m2/sec to mm/day by multiplying by 86400 (seconds in a day) and the new variable is then indicated with P. The BioClimInd short name, input variables and definition sources are reported in Table 1, formulas are described in detail in Supplementary Table 1. In the extended description below, the long-term or climatological period (or simply “period”) refers to the 40-year time frames used for the analysis, i.e. 1960–1999 for the historical period (to match the reference period for the bias-correction) and 2040–2079 and 2060–2099 for the future periods, in order to enable comparison between future and historical time frames with the same length. When calendar months are mentioned, it means that a month is considered at climatological level (i.e. its 40-year average). The approach for the definition of wettest/driest/coldest/warmest months or quarters identifies these periods year by year allowing to consider inter-annual variability of intra-annual cycles, as also interestingly investigated for different “climatological” windows by50.
- Bio1 Annual mean temperature
This indicator, expressed in °C, indicates the total amount of energy inputs for the ecosystems in a year. The annual average of daily mean temperature (Tg) is first computed for each year in the considered long-term period, and then it is further averaged among all years in the period.
- Bio2 Mean diurnal range
The daily fluctuations of temperatures have a strong influence on ecosystems51. This indicator is expressed in °C and is calculated by averaging, within the considered period, the daily differences between the maximum (Tx) and minimum (Tn) temperature.
- Bio3 Isothermality
This indicator quantifies how large the day-to-night temperatures oscillate relatively to the annual oscillations among extreme (warmest and coldest) months. Isothermality is the ratio, expressed in %, of Bio2/Bio7 (see below for Bio7).
- Bio4 Temperature seasonality
This indicator, expressed in °C, measures the temperature change throughout the year. Based on the formulation developed by Hijmans in WorldClim18, the average of daily mean temperature (Tg) is calculated for each calendar month in the selected period, and then the Standard Deviation is computed among the 12 monthly values obtained. The larger the value of the Standard Deviation, the greater the variability of temperature within the year.
- Bio5 Maximum temperature of warmest month
For this indicator, expressed in °C, first the monthly mean of daily maximum temperature (Tx) is calculated, then for each year in the period the warmest month, i.e. the one with the highest mean of daily Tx, is selected and the climatological mean of these values is finally computed. Note that the warmest month could be different from year to year.
- Bio6 Minimum temperature of coldest month
For this indicator, expressed in °C, first the monthly mean of daily minimum temperature (Tn) is calculated, then for each year in the period the coldest month, i.e. the one with the lowest mean of daily Tn, is selected and the climatological mean of these values is finally computed. Note that the coldest month could be different from year to year.
- Bio7 Temperature annual range
Similar to Bio2, also the temperature fluctuations between the warmest vs. coldest month within a year play a crucial role on ecosystems. This indicator measures the range of temperature between extreme (warmest and coldest) months. Practically, Bio7 is the difference between Bio5 and Bio6, expressed in °C, and is also used to calculate Bio3.
- Bio8 Mean temperature of wettest quarter
To calculate this indicator, expressed in °C, first the wettest quarter (i.e. that with the highest precipitation amount calculated by summing daily P) of each year in the whole period is identified (considering a quarter belonging to the year of its central month), and then the mean temperature among all wettest quarters is calculated by averaging daily Tg. Note that the wettest quarter could be different from year to year.
- Bio9 Mean temperature of driest quarter
To calculate this indicator, expressed in °C, first the driest quarter (i.e. that with the lowest precipitation amount calculated by summing daily P) of each year in the whole period is identified (considering a quarter belonging to the year of its central month), and then the mean temperature among all driest quarters is calculated by averaging daily Tg. Note that the driest quarter could be different from year to year.
- Bio10 Mean temperature of warmest quarter
To calculate this indicator, expressed in °C, first the warmest quarter (i.e. that with the highest average of daily mean temperature Tg) of each year in the whole period is identified (considering a quarter belonging to the year of its central month), and then the mean temperature among all warmest quarters is calculated by averaging daily Tg. Note that the warmest quarter could be different from year to year.
- Bio11 Mean temperature of coldest quarter
To calculate this indicator, expressed in °C, first the coldest quarter (i.e. that with the lowest average of daily mean temperature Tg) of each year in the whole period is identified (considering a quarter belonging to the year of its central month), and then the mean temperature among all coldest quarters is calculated by averaging daily Tg. Note that the coldest quarter could be different from year to year.
- Bio12 Annual precipitation
This is a widely used indicator, representing the total amount of water inputs to the ecosystems and to their water cycle. This indicator is expressed in millimeter (mm) per year and it is derived by averaging, along the whole period, the annual sum of daily precipitation amounts P.
- Bio13 Precipitation of wettest month
For this indicator, expressed in mm, first the monthly sum of daily precipitation P is calculated, then for each year in the period the wettest month, i.e. the one with the highest precipitation amount, is identified and then the precipitation amount among all wettest months is averaged. Note that the wettest month could be different from year to year.
- Bio14 Precipitation of driest month
For this indicator, expressed in mm, first the monthly sum of daily precipitation P is calculated, then for each year in the period the driest month, i.e. the one with the lowest precipitation amount, is identified and then the precipitation amount among all driest months is averaged. Note that the driest month could be different from year to year.
- Bio15 Precipitation seasonality
Species distribution can be heavily influenced by precipitation variability during the year46. The precipitation seasonality, expressed in %, is the ratio between the standard deviation and the mean of 12 values representing the monthly average precipitation as calculated at calendar month level within the considered period. To avoid division by 0, the denominator is increased by 1.
- Bio16 Precipitation of wettest quarter
To calculate this indicator, expressed in mm, first the wettest quarter (i.e. that with the highest precipitation amount calculated by summing daily P) of each year in the whole period is identified (considering a quarter belonging to the year of its central month) (see Bio8), and then the precipitation amount is averaged among all wettest quarters. Note that the wettest quarter could be different from year to year.
- Bio17 Precipitation of driest quarter
To calculate this indicator, expressed in mm, first the driest quarter (i.e. that with the lowest precipitation amount calculated by summing daily P) of each year in the whole period is identified (considering a quarter belonging to the year of its central month) (see Bio9), and then the precipitation amount is averaged among all driest quarters. Note that the driest quarter could be different from year to year.
- Bio18 Precipitation of warmest quarter
To calculate this indicator, expressed in mm, first the warmest quarter (i.e. that with the highest average of daily mean temperature Tg) of each year in the whole period is identified (considering a quarter belonging to the year of its central month) (see Bio10), and then the precipitation amount is averaged among all warmest quarters. Note that the warmest quarter could be different from year to year.
- Bio19 Precipitation of coldest quarter
To calculate this indicator, expressed in mm, first the coldest quarter (i.e. that with the lowest average of daily mean temperature Tg) of each year in the whole period is identified (considering a quarter belonging to the year of its central month) (see Bio11), and then the precipitation amount is averaged among all coldest quarters. Note that the coldest quarter could be different from year to year.
- Bio20 Ellenberg quotient
This indicator is firstly described by Ellenberg in 196352 and it relates the temperature of the warmest month of the year to the annual precipitation. It is calculated in °C/mm as the ratio Bio28/Bio12 multiplied by 1000. This is a simple measure of humidity in respect to continentality.
- Bio21 Yearly positive temperature
Following the Rivas-Martinez bioclimatic approach47, this indicator, expressed in °C, consists in the sum of the monthly average temperature for those months whose long-term average of daily temperature Tg is higher than 0 °C.
- Bio22 Sum of annual temperature
This indicator is similar to Bio21 from which it differs because it is the sum of the monthly average temperature for all months expressed in °C. In large portion of the globe this data results smaller than Bio21. It is also approximately equivalent to Bio1 multiplied by 12, but has been included in this dataset for a direct comparison with Bio21, Bio25 and Bio26.
- Bio23 Ombrothermic index
This indicator, calculated in mm/°C, is ten times the ratio between the average precipitation amount (in mm) for months with positive average of daily Tg (Bio24) and the sum of average temperature for the same months (Bio21)47.
- Bio24 Yearly positive precipitation
Still following47 this indicator, expressed in mm, represents the sum of the long-term average of total precipitation (cumulated of daily P) for those months whose long-term average of daily temperature Tg is higher than 0 °C.
- Bio25 Modified Kira coldness index
Described by Kira48,49, and partially modified, this indicator describes the amount of energy inputs in the coldest portion of the year. It consists of the sum of the monthly average temperature for those months whose long-term average of daily temperature Tg is lower than 5 °C.
- Bio26 Modified Kira warmth index
Described by Kira48,49, and partially modified, this indicator describes the amount of energy inputs in the warmest portion of the year. It consists of the sum of the monthly average temperature for those months whose long-term average of daily temperature Tg is higher than 5 °C.
- Bio27 Simplified continentality index
Similar to Bio7, but based on Tg instead of Tx and Tn, this is an another measure of temperature extremes. It is the difference, in °C, between the average Tg of warmest (Bio28) and of the coldest (Bio29) months from the years of the considered period. Partially modified from Driscoll index53
- Bio28 Mean temperature of warmest month
For this indicator, expressed in °C, first the warmest months as in Bio5 are extracted for each year in the considered period, and then the average among daily mean temperature (Tg) is calculated among these months.
- Bio29 Mean temperature of coldest month
For this indicator, expressed in °C, first the coldest months as in Bio6 are extracted for each year in the considered period, and then the average among daily mean temperature (Tg) is calculated among these months.
- Bio30 Mean temperature of driest month
For this indicator, expressed in °C, first the driest months as in Bio14 are extracted for each year in the considered period, and then the average among daily mean temperature (Tg) is calculated among these months.
- Bio31 Mean temperature of wettest month
For this indicator, expressed in °C, first the wettest months as in Bio13 are extracted for each year in the considered period, and then the average among daily mean temperature (Tg) is calculated among these months.
- Bio32 Modified Thermicity index
Revised from what was described by Rivas-Martinez47, this indicator, expressed in tents of °C, is the sum of annual mean temperature (Bio1), the average of maximum temperatures of warmest month (Bio5) and the average of minimum temperatures of coldest month (Bio6).
- Bio33 Ombrothermic index of summer and the previous month
Modified from Rivas-Martinez47, this indicator, expressed in mm/°C, is suited for the Northern Hemisphere, dividing the sum of P and the average of Tg for the late spring/summer period (May, June, July and August)
- Bio34 Potential Evapotranspiration Hargreaves
This indicator calculates the potential evapotranspiration (PET, mm) through the formula of Hargreaves and Samani54, that requires as input the time series of monthly mean of Tg and of daily temperature range (this calculated from daily input of Tx and Tn). First the PET is cumulated for each year in the considered period and then averaged among years.
- Bio35 Potential Evapotranspiration Thornthwaite
This indicator calculates the potential evapotranspiration (PET, mm) through the formula of Thornthwaite55, that requires as input the time series of monthly mean of Tg. First the PET is cumulated for each year in the considered period and then averaged among years.
Data Records
The complete CMCC-BioClimInd dataset56 is available through PANGAEA (https://doi.org/10.1594/PANGAEA.904278). It consists in 805 files in NetCDF4 format with a 0.5° by 0.5° grid resolution and global coverage (except Antarctica). Files represent 35 bioclimatic indicators calculated for a 40-years historical interval (1960–1999) under climate reanalysis and for two future 40-years time intervals (2040–2079 and 2060–2099) under 6 ESMs’ projections and 2 RCPs (except RCP4.5 for CMCC-CESM) (Table 2). The size of the array for each map is 720 (longitudes) x 360 (latitudes) grid points for a total of 259200 points, whose 67415 differ from NULL as representing land areas. The individual file names follow this structure: “BIOx institute-model rcp yyyy zz.nc”
where x identifies the code of the bioclimatic indicator, institute-model (see Table 3, first column) identifies the source of data (short name for historical observation or for the producing Center and/or ESMs), rcp is the Representative Concentration Pathway (RCP 4.5 or RCP 8.5), yyyy is the starting year of time interval and zz the last two digits of the ending year of the time interval.
Technical Validation
Comparison of CMCC-BioClimInd 1.0 and WorldClim 2.0 datasets for the historical period
To check the quality of our dataset the first 19 BioClimInd were compared with the same indicators from WorldClim for the historical period, ensuring independence of datasets as relying on different climate data sources (reanalysis-based for CMCC-BioClimInd, interpolation from station data for WorldClim). The WorldClim Bioclimatic variables were downloaded from http://www.worldclim.org/ version 2 at a 10’ resolution, then resampled to a 0.5° resolution in ESRI ArcGIS Desktop 10.1 using the bilinear resampling technique (Fig. 3). Subsequently, for each indicators only the 64515 valid (different from no data) grid cells from both CMCC-BioClimInd and WorldClim datasets were considered in the comparison. It is necessary to remark that historical CMCC-BioClimInd refer to 1960–1999 period whereas WorldClim to 1970–2000.
Although the two datasets overlap in time for 30 out of 40 years (1970–1999), Fig. 4 shows a good correlation between them. For 11 BioClimInd the R2 is above 0.9, for 7 between 0.75 and 0.9, only the comparison of Bio19 between the two datasets has shown a lower level of correlation. In particular, for the basic indicators Annual Mean Temperature (Bio1) and Annual precipitation (Bio12) the level of correlation is excellent (i.e. R2 = 0.991 - 0.917).
In Fig. 5 are reported maps of spatial differences for Bio1 and Bio12; main divergences in annual temperature can be found in the North-Eastern portion of Siberia, Arabian peninsula, western Africa where BioClimInd results colder, in contrast with a large portion of Greenland and Himalayan mountains where BioClimInd results sensibly hotter. However, in most part of the globe, temperature differences are within ±2.5°C. For precipitation, the differences, in terms of lower BioCLimInd values are recognized in almost all desert areas (i.e. example Sahara, Arabian peninsula, Karakum, Gobi, Kalahari etc.) although it is important to underline that even with high percentage values the differences in absolute terms (mm) are minimal but appear high as compared to low reference values; WorlClim results sensibly drier in many areas of central Africa, Southern America and Himalayan mountains.
It is worth noticing that the above mentioned differences are mainly located in areas where the variables estimate is less accurate due to the paucity of ground observations57 and some artifacts may arise from the interpolation function used when creating the spatial gridded dataset. Therefore, one of the possible explanations for the differences between WorldClim and CMCC-BioClimInd, could be the different weight given to observations with respect to climate model data when creating the datasets. In Table 4 the main descriptive statistics for the 19 common indicators and the results of the Paired Samples T-Test are reported.
Further remarks
After the above necessary comparison, the authors want to emphasize that the entire dataset presented in this work is based on elaborations of climate reanalysis for the historical period, whose validation effort are available in literature40,58,59. Hence, it can be assumed that outputs are at least as robust as the datasets from which they originate and therefore do not need further validation. For future projections, the dataset uses bias-corrected outputs of climate models and it is recommended to exploit this dataset under an ensemble perspective37. In conclusion it should be emphasized that combining several parameters for the creation of indices (see Table 1) inevitably increases the level of uncertainty16 due to error propagation.
Usage Notes
When using CMCC-BioclimInd dataset or part of it, please cite this manuscript. Besides data sources for the dataset presented in this manuscript, data from further climate simulations can be integrated, e.g. considering also RCPs 2.6 and 6.0 and/or additional models as available from other ISIMIP rounds. Data are in NetCDF4 format, but other common interoperable formats (i.e. ESRI grid, GEOTIFF) can be provided, by contacting the authors, for a sub-selection of the dataset, whose size will be defined based on the specific requests from interested users and on the processing time needed. For any questions, suggestions or request of collaboration regarding CMCC-BioClimInd please contact the corresponding author. Data are freely available under the Creative Commons License: CC BY
Code availability
All the BioClimInd calculations were conducted exploiting NetCDF data manipulation and analysis tools, namely CDO (Climate Data Operators) and NCO (NetCDF Operators) combined through DOS Batch (.bat) commands and scripts available through Github at https://github.com/CMCC-Foundation/BioClimInd.
References
Walther, G. et al. Ecological responses to recent climate change. Nature 416, 389–395 (2011).
Ehrlén, J. & Morris, F. M. Predicting changes in the distribution and abundance of species under environmental change. Ecology Letters 18, 303–314 (2015).
Valiente-Banuet, A. et al. Beyond species loss: the extinction of ecological interactions in a changing world. Functional Ecology 29, 299–307 (2015).
The Intergovernamental Panel on Climate Change - IPCC Special Report on Climate Change and Land (2019).
Ancillotto, L., Santini, L., Ranc, N., Maiorano, L. & Russo, D. Extraordinary range expansion in a common bat: the potential roles of climate change and urbanisation. The Science of Nature 103, 15 (2016).
Deblauwe, V. et al. Remotely sensed temperature and precipitation data improve species distribution modelling in the tropics. Global Ecology and Biogeography 25, 443–454 (2016).
Noce, S., Collalti, A. & Santini, M. Likelihood of changes in forest species suitability, distribution, and diversity under future climate: The case of Southern Europe. Ecology and Evolution 7, 9358–9375 (2017).
Marchi, M., Sinjur, I., Bozzano, M. & Westergren, M. Evaluating WorldClim Version 1 (1961–1990) as the Baseline for Sustainable Use of Forest and Environmental Resources in a Changing Climate. Sustainability 11, 3043 (2019).
Noce, S., Caporaso, L., & Santini, M. Climate change and geographic ranges: the implications for Russian forests. Frontiers in Ecology and Evolution 19, (2019)
Yim, Y. & Kira, T. Distribution of forest vegetation and climate in the Korean peninsula: I. distribution of some indices of thermal climate. Japanese Journal of Ecology 25, 77–88 (1975).
Pearson, R. G. & Dawson, T. P. Bioclimate envelope models: what they detect and what they hide – response to Hampe. Global Ecology and Biogeography 13, 471–473 (2004).
Araujo, M. & Luoto, M. The importance of biotic interactions for modelling species distributions under climate change. Global Ecology and Biogeography 16, 743–753 (2007).
Hu, Z., Hu, Q., Zhang, C., Chen, X. & Li, Q. Evaluation of reanalysis, spatially interpolated and satellite remotely sensed precipitation data sets in central Asia. Journal of Geophysical Research: Atmospheres 121, 5648–5663 (2016).
Ehret, U., Zehe, E., Warrach-Sagi, K. & Liebert, J. HESS Opinions “Should we apply bias correction to global and regional climate model data?”. Hydrology and Earth System Sciences 9, 3391–3404 (2012).
Galmarini, S. et al. Adjusting climate model bias for agricultural impact assessment: How to cut the mustard Climate Services 13, 65–69 (2019).
Maraun, D. et al. Towards process-informed bias correction of climate change simulations. Nature Climate Change 7, 764–773 (2017).
Maraun, D. Bias Correcting Climate Change Simulations-a Critical Review. Current Climate Change Reports 2, 211–220 (2016).
Hijmans, R. J., Cameron, S. E., Parra, J. L., Jones, P. G. & Jarvis, A. Very high resolution interpolated climate surfaces for global land areas. International Journal of Climatology 25, 1965–1978 (2012).
Fick, S. & Hijmans, R. WorldClim 2: new 1-km spatial resolution climate surfaces for global land areas. International Journal of Climatology 37, 4302–4315 (2017).
Karger, D.N. et al. Climatologies at high resolution for the earth’s land surface areas. Scientific Data 4, (2017).
Karger, D. N., Schmatz, D. R., Dettling, G. & Zimmermann, N. E. High-resolution monthly precipitation and temperature time series from 2006 to 2100. Scientific Data 7, 1–10 (2020).
Kriticos, D. et al. CliMond: Global high-resolution historical and future scenario climate surfaces for bioclimatic modelling. Methods in Ecology and Evolution 3, 53–64 (2012).
Kriticos, D., Jarosik, B. & Ota, N. Extending the suite of bioclim variables: a proposed registry system and case study using principal components analysis. Methods in Ecology and Evolution 5, 956–960 (2012).
Lima-Ribeiro, M. S. et al. EcoClimate: a database of climate data from multiple models for past, present, and future for macroecologists and biogeographers. Biodiversity Informatics 10, 1–21 (2015).
Title, P. & Bemmels, J. ENVIREM: an expanded set of bioclimatic and topographic variables increases flexibility and improves performance of ecological niche modeling. Ecography 41, 291–307 (2018).
Vega, G.C., Pertierra, L.R., Ollala-Tarraga, M.A. Data from: MERRAclim, a high-resolution global dataset of remotely sensed bioclimatic variables for ecological modelling. Scientific Data 4 (2017).
Poiani, K. A., Ritcher, B. D., Anderson, M. G. & Ritcher, H. E. Biodiversity Conservation at Multiple Scales: Functional Sites, Landscapes, and Networks. BioScience 50, 133–146 (2000).
Tshwene-Maychaza, B. & Aguirre-Gutiérrez, J. Climatic Drivers of Plant Species Distributions Across Spatial Grains in Southern Africa Tropical Forests. Frontiers in Forests and Global Change 2, 69 (2019).
Trivedi, M. R., Berry, P. B., Morecroft, M. D. & Dawson, T. P. Spatial scale affects bioclimate model projections of climate change impacts on mountain plants. Global Change Biology 14, 1089–1103 (2008).
Priya, P. & Bofu, Y. Interannual variations in rainfall of different intensities in South West of Western Australia. International Journal of Climatology 40, 3052–3071 (2019).
Pantano, V. C., Penalba, O. C., Spescha, L. B. & Murphy, G. M. Assessing how accumulated precipitation and long dry sequences impact the soil water storage. International Journal of Climatology 37, 4316–4326 (2017).
Spinoni, J. et al. Future Global Meteorological Drought Hot Spots: A Study Based on CORDEX Data. Journal of Climate 33, 3635–3661 (2020).
van Vuuren, D. P. et al. The representative concentration pathways: an overview. Climate Change 109, 5–31 (2011).
Taylor, K. E., Stouffer, R. J. & Meehl, G. A. An overview of CMIP5 and the experiment design. Bulletin of the American Meteorological Society 93, 485–498 (2012).
Hempel, S., Frieler, K., Warszawski, L., Schewe, J. & Piontek, F. A trend-preserving bias correction – the ISIMIP approach. Earth System Dynamics 4, 219–236 (2013).
Kenney, M.A., Janetos, A.C., & Gerst, M.D. A framework for national climate indicators Climatic Change (2018).
Collins, M. Ensembles and probabilities: a new era in the prediction of climate change. Philosophical Transactions of The Royal Society 365, 1957–1970 (2007).
Uppala, S. M. et al. The ERA-40 re-analysis. Quarterly Journal of the Royal Meteorological Society: A journal of the atmospheric sciences, applied meteorology and physical oceanography 131, 2961–3012 (2005).
Piani, C. et al. Statistical bias correction of global simulated daily precipitation and temperature for the application of hydrological models. Journal of Hydrology 395, 199–215 (2010).
Weedon, G. et al. Creation of the WATCH Forcing Data and Its Use to Assess Global and Regional Reference Crop Evaporation over Land during the Twentieth Century. Journal of Hydrometeorology 12, 823–848 (2011).
IPCC Working Group I Contribution to the IPCC Fifth Assessment Report, Climate Change 2013: The Physical Science Basis. IPCC,AR5 (2013).
IPCC Contribution of Working Group II to the Fifth Assessment Report of the Intergovernmental Panel on Climate Change. IPCC,AR5 (2014).
Vichi, M. et al. Global and regional ocean carbon uptake and climate change: sensitivity to a substantial mitigation scenario. Climate Dynamics 37, 1929–1947 (2011).
Thomson, A. et al. RCP4.5: A pathway for stabilization of radiative forcing by 2100. Climatic Change 109, 77–94 (2011).
Riahi, K. et al. RCP 8.5- A scenario of comparatively high greenhouse gas emissions. Climatic Change 109, 33–57 (2011).
O’Donnell, M.A, Ignizio, D.A. Bioclimatic predictors for supporting ecological applications in the conterminous United States. U.S. Geological Survey Data Series 691 (2012).
Rivas-Martínez, S., Rivas Sáenz, S. & Penas, A. Worldwide bioclimatic classification system Global Geobotany 1, 1–634 (2011).
Kira, T. A new classification of climate in eastern Asia as the basis for agricultural geography. Horticultural Institute: Kyoto University (1945).
Kira, T. Forest ecosystems of east and southeast Asia in a global perspective. Ecological Research 6, 185–200 (1991).
Monahan, W.B., Fisichelli, N.A. Climate Exposure of US National Parks in a New Era of Change Plos One 9, (2014).
Ehbrecht, M., Schall, P., Ammer, C., Fisher, M. & Seidel, D. Effects of structural heterogeneity on the diurnal temperature range in temperate forest ecosystems. Forest Ecology and Management 432, 860–867 (2019).
Ellenberg, H. Vegetation Mitteleuropas mit den Alpen: in kausaler, dynamischer und historischer Sicht. Ulmer (1963).
Driscoll, D. M. & Yee Fong, J. M. Continentality: a basic climatic parameter re-examined. International Journal of Climatology 12, 185–192 (1992).
Hargreaves, G. H. & Samani, Z. Estimating Potential Evapotranspiration Journal of the Irrigation and Drainage Division 108, 225–230 (1982).
Thornthwaite, C. An Approach toward a Rational Classification of Climate. Geographical Review 38, 55–94 (1948).
Noce, S., Caporaso, L., Santini, M. CMCC-BioClimInd. A new global dataset of bioclimatic indicators. PANGAEA https://doi.org/10.1594/PANGAEA.904278 (2019)
Harris, I., Osborn, T. J., Jones, P., Lister, D. Version 4 of the CRU TS monthly high-resolution gridded multivariate climate dataset Scientific Data 109 (2020).
Li, L., Ngondongondo, C. S., Xu, C. & Gong, L. Comparison of the global TRMM and WFD precipitation datasets in driving a large-scale hydrological model in southern Africa. Hydrology Research 44, 770–788 (2013).
Li, L., Xu, C., Zhang, Z. & Jain, S. K. Validation of a new meteorological forcing data in analysis of spatial and temporal variability of precipitation in India. Stochastic environmental research and risk assessment 28, 239–252 (2014).
Author information
Authors and Affiliations
Contributions
S.N. designed the research, created the data, wrote the paper; L.C. provided the input climatic data, wrote the paper; M.S. designed the research, created the data, wrote the paper.
Corresponding author
Ethics declarations
Competing interests
The authors declare no competing interests.
Additional information
Publisher’s note Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Supplementary information
Rights and permissions
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.
The Creative Commons Public Domain Dedication waiver http://creativecommons.org/publicdomain/zero/1.0/ applies to the metadata files associated with this article.
About this article
Cite this article
Noce, S., Caporaso, L. & Santini, M. A new global dataset of bioclimatic indicators. Sci Data 7, 398 (2020). https://doi.org/10.1038/s41597-020-00726-5
Received:
Accepted:
Published:
DOI: https://doi.org/10.1038/s41597-020-00726-5
This article is cited by
-
Projections of future bioclimatic indicators using bias-corrected CMIP6 models: a case study in a tropical monsoon region
Environmental Science and Pollution Research (2024)
-
Modeling the species occurrence probability and response of climate change on Himalayan Somalata plant under different Shared Socioeconomic Pathways
Environmental Monitoring and Assessment (2024)
-
Elevation transition of aquatic insects closely matches a thermal feature in the Yungas of Northwestern Argentina
Journal of Mountain Science (2024)
-
Habitat suitability of Opuntia ficus-indica (L.) MILL. (CACTACEAE): a comparative temporal evaluation using diverse bio-climatic earth system models and ensemble machine learning approach
Environmental Monitoring and Assessment (2024)
-
Modeling soil acidity (pH) dynamics under extreme agroclimatic conditions in Horro Guduru Wallaga Zone, northwestern Ethiopia
Environmental Geochemistry and Health (2024)