Global Climate Datasets Research Articles

Combining global precipitation data and machine learning to predict flood peaks in ungauged areas with similar climate

Increasing flood risk due to urbanization and climate change poses a significant challenge to societies at global scale. Hydrologic information that is required for understanding flood processes and for developing effective warning procedures is currently lacking in most parts of the world. Procedures that can combine global climate dataset from satellite and reanalysis with fast and low computational cost prediction systems, are attractive solutions for addressing flood predictions in ungauged areas. This work develops and tests a prediction framework that relies on two fundamental components. First, meteorological data from global datasets (IMERG and ERA5-Land) provide key input variables and second, ML models trained in the data-rich contiguous US, are applied in climatically similar regions in other parts of the world. Catchments in Australia, Brazil, Chile, Switzerland, and Great Britain were used as pseudo-ungauged regions for testing. Results indicate acceptable performance for both IMERG and ERA5-Land forced models with relative difference in flood peak prediction within 30 % and similar overall performance to locally trained ML models. Specific climate regions for which ML models have revealed good performance include Mediterranean climates like the US West Coast, subtropical areas like the Southern Atlantic Gulf, and mild temperate regions like the Mid-Atlantic Basin. This work highlights the potential of combining global precipitation dataset with pre-trained ML models in data-rich areas, for flood prediction in ungauged areas with similar climate.

Physical quantity characteristics of severe aircraft turbulence near convective clouds over Australia

Using FY–2G satellite data, Aircraft Meteorological Data Relay (AMDAR) downlink data, and the fifth generation European Centre for Medium–Range Weather Forecasts (ECMWF) atmospheric reanalysis of the global climate (ERA5) dataset, we analyzed the circulation, thermal, and dynamic characteristics of the convective aircraft turbulence over Australia in the Southern Hemisphere. The results show that the near-convective clouds turbulence (NCCT) in the Southern Hemisphere mostly occurred in front of deep warm high–pressure ridges in mid–latitude regions and on the left side of the axis of the subtropical westerly jet stream. The isotherms in this area were relatively dense (i.e., large gradients), and the wind speed was high, with strong horizontal and/or vertical cyclonic wind shear. In addition, the NCCT usually occurred near the zero–divergence or zero–vorticity line and in areas with large vertical wind speed gradients. There were also strong vertical and horizontal wind shears in this area, which could easily trigger severe turbulence. Furthermore, the NCCT in the Southern Hemisphere mostly occurred at the intersection of cold and warm temperature advections (i.e., near the zero–temperature advection line), and the turbulence point was located near the high–altitude frontal zone where there was a strong gradient of cold and warm advections. There was temperature inversion with pseudo–equivalent potential temperatures in the middle and lower troposphere on the warmer side of the turbulence point. The unstable stratification of cold air at the top and warm air at the bottom was conducive to triggering convection from the ground, forming strong convective clouds, and causing severe turbulence.

ASSESSMENT OF SOLAR PV OUTPUT PERFORMANCE WITH VARYING TILT ANGLES AND WEATHER DATA FROM ERA5: CASE OF MUNTINLUPA CITY, PHILIPPINES

Abstract. Solar photovoltaic (PV) technology has been gaining popularity in the Philippines as an alternative source of sustainable energy. In such technology, module tilt angle and weather conditions are among the system parameters that have substantial impacts on PV system performance. Previous studies have considered either tilt angles or weather conditions, but not the combined impact of these two parameters on solar PV power output. The objective of this study is to examine the effects of weather variables and tilt angle on the output of solar photovoltaic (PV) systems in Muntinlupa City, the Philippines. Three 120W monocrystalline solar PV panels were used and set up to three different tilt angles (i.e., 5°, 10°, and 15°). The fifth generation of the ECMWF's global climate and weather reanalysis (ERA5) dataset was used to gather hourly weather information such as surface solar radiation, wind speed, wind direction, relative humidity, ambient temperature, and total precipitation. Three principal components (PC), which together account for 95% of the variability, were identified using principal components analysis (PCA), which was used to address multicollinearity among the weather parameters. To assess the effects of tilt angle, time, and PCs on solar PV production, analysis of covariance (ANCOVA) was carried out. Results show that all weather variables, except for wind speed and total precipitation, have a significant impact on solar PV production with configuration producing the best results. Moreover, a significant difference in mean solar PV production was observed among the three tilt angles. From 6:00 AM to 2:00 PM, solar PV output gradually increases and declines thereafter. Outputs this study can help in optimizing the design and configuration of solar PV systems in the Philippines by considering weather variables and module tilt angle. Lastly, this study provides useful information for system designers, installers, and policymakers in improving energy generation and utilization, encouraging the use of renewable energy sources, and advancing sustainable energy objectives.

Harnessing deep learning to forecast local microclimate using global climate data

Microclimate is a complex non-linear phenomenon influenced by both global and local processes. Its understanding holds a pivotal role in the management of natural resources and the optimization of agricultural procedures. This phenomenon can be effectively monitored in local areas by employing models that integrate physical laws and data-driven algorithms relying on climate data and terrain conformation. Climate data can be acquired from nearby meteorological stations when available, but in their absence, global climate datasets describing 10 km-scale areas are often utilized. The present research introduces an innovative microclimate model that combines physical laws and deep learning to reproduce temperature and relative humidity variations at the meter-scale within a study area located in the Lombardian foothills. The model is exploited to perform a comparative study investigating whether employing the global climate dataset ERA5 as input reduces model’s accuracy in reproducing the microclimate variations compared to using data collected by the Lombardy Regional Environment Protection Agency (ARPA) from a nearby meteorological station. The comparative analysis shows that using local meteorological data as inputs provides more accurate results for microclimate modeling. However, in situations where local data is not available, the use of global climate data remains a viable and reliable approach.

A Multi-Framework of Google Earth Engine and GEV for Spatial Analysis of Extremes in Non-Stationary Condition in Southeast Queensland, Australia

The frequency and severity of extremes, including extreme precipitation events, extreme evapotranspiration and extreme water storage deficit events, are changing. Thus, the necessity for developing a framework that estimates non-stationary conditions is urgent. The aim of this paper is to develop a framework using the geeSEBAL platform, Generalised Extreme Value (GEV) models and spatiotemporal analysis techniques that incorporate the physical system in terms of cause and effect. Firstly, the geeSEBAL platform has enabled the estimation of actual evapotranspiration (ETa) with an unprecedented level of spatial-temporal resolution. Following this, the Non-stationary Extreme Value Analysis (NEVA) approach employs the Bayesian method using a Differential Evolution Markov Chain technique to calculate the frequency and magnitude of extreme values across the parameter space. Station and global climate datasets have been used to analyse the spatial and temporal variation of rainfall, reference evapotranspiration (ETo), ETa and water storage (WS) variables in the Lockyer Valley located in Southeast Queensland (SEQ), Australia. Frequency analysis of rainfall, ETa, and water storage deficit for 14 stations were performed using a GEV distribution under stationary and non-stationary assumptions. Comparing the ETa, ETo and ERA5 rainfall with station data showed reasonable agreement as follows: Pearson correlation of 0.59–0.75 for ETa, RMSE of 45.23–58.56 mm for ETa, Pearson correlation of 0.96–0.97 for ETo, RMSE of 73.13–87.73 mm for ETo and Pearson correlation of 0.87–0.92 for rainfall and RMSE of 37.53–57.10 mm for rainfall. The lower and upper uncertainty bounds between stationary and non-stationary conditions for rainfall station data of Gatton varied from 550.98 mm (stationary) to 624.97 mm (non-stationary), and for ERA5 rainfall datasets, 441.30 mm (stationary) to 450.77 mm (non-stationary). The results demonstrate that global climate datasets underestimate the difference between stationary and non-stationary conditions by 9.47 mm compared to results of 73.99 mm derived from station data. Similarly, the results demonstrate less variation between stationary and non-stationary conditions in water storage, followed by a sharp variation in rainfall and moderate variation in evapotranspiration. The findings of this study indicate that neglecting the non-stationary condition in some hydrometeorological variables can lead to underestimating their amounts. This framework can be applied to any geographical area for estimating extreme conditions, providing valuable insights for infrastructure planning and design, risk assessment and disaster management.

CONUS404: The NCAR–USGS 4-km Long-Term Regional Hydroclimate Reanalysis over the CONUS

Abstract A unique, high-resolution, hydroclimate reanalysis, 40-plus-year (October 1979–September 2021), 4 km (named as CONUS404), has been created using the Weather Research and Forecasting Model by dynamically downscaling of the fifth-generation European Centre for Medium-Range Weather Forecasts (ECMWF) atmospheric reanalysis of the global climate dataset (ERA5) over the conterminous United States. The paper describes the approach for generating the dataset, provides an initial evaluation, including biases, and indicates how interested users can access the data. The motivation for creating this National Center for Atmospheric Research (NCAR)–U.S. Geological Survey (USGS) collaborative dataset is to provide research and end-user communities with a high-resolution, self-consistent, long-term, continental-scale hydroclimate dataset appropriate for forcing hydrological models and conducting hydroclimate scientific analyses over the conterminous United States. The data are archived and accessible on the USGS Black Pearl tape system and on the NCAR supercomputer Campaign storage system.

DOCU-CLIM: A global documentary climate dataset for climate reconstructions

Documentary climate data describe evidence of past climate arising from predominantly written historical documents such as diaries, chronicles, newspapers, or logbooks. Over the past decades, historians and climatologists have generated numerous document-based time series of local and regional climates. However, a global dataset of documentary climate time series has never been compiled, and documentary data are rarely used in large-scale climate reconstructions. Here, we present the first global multi-variable collection of documentary climate records. The dataset DOCU-CLIM comprises 621 time series (both published and hitherto unpublished) providing information on historical variations in temperature, precipitation, and wind regime. The series are evaluated by formulating proxy forward models (i.e., predicting the documentary observations from climate fields) in an overlapping period. Results show strong correlations, particularly for the temperature-sensitive series. Correlations are somewhat lower for precipitation-sensitive series. Overall, we ascribe considerable potential to documentary records as climate data, especially in regions and seasons not well represented by early instrumental data and palaeoclimate proxies.

CHELSA-W5E5: daily 1 km meteorological forcing data for climate impact studies

Abstract. Current changes in the world's climate increasingly impact a wide variety of sectors globally, from agriculture and ecosystems to water and energy supply or human health. Many impacts of climate on these sectors happen at high spatio-temporal resolutions that are not covered by current global climate datasets. Here we present CHELSA-W5E5 (https://doi.org/10.48364/ISIMIP.836809.3, Karger et al., 2022): a climate forcing dataset at daily temporal resolution and 30 arcsec spatial resolution for air temperatures, precipitation rates, and downwelling shortwave solar radiation. This dataset is a spatially downscaled version of the 0.5∘ W5E5 dataset using the CHELSA V2 topographic downscaling algorithm. We show that the downscaling generally increases the accuracy of climate data by decreasing the bias and increasing the correlation with measurements from meteorological stations. Bias reductions are largest in topographically complex terrain. Limitations arise for minimum near-surface air temperatures in regions that are prone to cold-air pooling or at the upper extreme end of surface downwelling shortwave radiation. We further show that our topographically downscaled climate data compare well with the results of dynamical downscaling using the Weather Research and Forecasting (WRF) regional climate model, as time series from both sources are similarly well correlated to station observations. This is remarkable given the lower computational cost of the CHELSA V2 algorithm compared to WRF and similar models. Overall, we conclude that the downscaling can provide higher-resolution climate data with increased accuracy. Hence, the dataset will be of value for a wide range of climate change impact studies both at global level and for applications that cover more than one region and benefit from using a consistent dataset across these regions.

Effect of Socioeconomic Variables in Predicting Global Fire Ignition Occurrence

Fires are a pervasive feature of the terrestrial biosphere and contribute large carbon emissions within the earth system. Humans are responsible for the majority of fire ignitions. Physical and empirical models are used to estimate the future effects of fires on vegetation dynamics and the Earth’s system. However, there is no consensus on how human-caused fire ignitions should be represented in such models. This study aimed to identify which globally available predictors of human activity explain global fire ignitions as observed by satellites. We applied a random forest machine learning framework to state-of-the-art global climate, vegetation, and land cover datasets to establish a baseline against which influences of socioeconomic data (cropland fraction, gross domestic product (GDP), road density, livestock density, grazed lands) on fire ignition occurrence were evaluated. Our results showed that a baseline random forest without human predictors captured the spatial patterns of fire ignitions globally, with hotspots over Sub-Saharan Africa and South East Asia. Adding single human predictors to the baseline model revealed that human variables vary in their effects on fire ignitions and that of the variables considered GDP is the most vital driver of fire ignitions. A combined model with all human predictors showed that the human variables improve the ignition predictions in most regions of the world, with some regions exhibiting worse predictions than the baseline model. We concluded that an ensemble of human predictors can add value to physical and empirical models. There are complex relationships between the variables, as evidenced by the improvement in bias in the combined model compared to the individual models. Furthermore, the variables tested have complex relationships that random forests may struggle to disentangle. Further work is required to detangle the complex regional relationships between these variables. These variables, e.g., population density, are well documented to have substantial effects on fire at local and regional scales; we determined that these variables may provide more insight at more continental scales.

The global historical climate database HCLIM

There is a growing need for past weather and climate data to support science and decision-making. This paper describes the compilation and construction of a global multivariable (air temperature, pressure, precipitation sum, number of precipitation days) monthly instrumental climate database that encompasses a substantial body of the known early instrumental time series. The dataset contains series compiled from existing databases that start before 1890 (though continuing to the present) as well as a large amount of newly rescued data. All series underwent a quality control procedure and subdaily series were processed to monthly mean values. An inventory was compiled, and the collection was deduplicated based on coordinates and mutual correlations. The data are provided in a common format accompanied by the inventory. The collection totals 12452 meteorological records in 118 countries. The data can be used for climate reconstructions and analyses. It is the most comprehensive global monthly climate dataset for the preindustrial period so far.

Statistical Characteristics of Blocking High in the Ural Mountains during Winters and Relationship with Changes in Sea Surface Temperature and Sea Ice

A blocking high in the Ural Mountains, which is recognized as the third major blocking high area in the northern hemisphere, describes a deep warm high-pressure system superimposed on the westerly belt. Based on the ERA-5 daily reanalysis data (the fifth-generation European Centre for Medium-Range Weather Forecasts atmospheric reanalysis global climate dataset) and using the Tibaldi and Molteni (TM) method, we selected 43 blocking high events in the Ural Mountains during the extended winters of 1979–2020 and analyzed their atmospheric circulation characteristics and influencing factors. Our findings revealed a downward trend in the frequency of occurrence of blocking highs in the Ural Mountains in winter, most of them were short-lived; furthermore, the frequency and duration of these occurrences generally followed a 3–4 years oscillating cycle. The synthetic results of the geopotential height (HGT) anomaly field and the surface air temperature (SAT) anomaly field of these 43 extended wintertime blocking high events in the Ural Mountains region showed that during the development of a blocking high, the central intensity of the positive anomalies in the Ural Mountains region first increased and then weakened, while the central intensity and meridional span of the negative anomalies in the Eurasian mid-latitudes of the SAT anomaly field increased continuously. In addition, abnormally high sea surface temperature (SST) in the North Atlantic sea area and abnormal reduction of sea ice (SI) in the Barents-Kara Sea and the Chukchi Sea in autumn had a significant impact on the wintertime formation of Ural Mountains blocking highs. In contrast, in autumn, the abnormal reduction of SI in the Barents-Kara and Chukchi Seas might also have led to the westward positioning of Ural Mountains blocking highs.

Gridded 20-year climate parameterization of Africa and South America for a stochastic weather generator (CLIGEN)

ABSTRACT CLIGEN is a stochastic weather generator that creates statistically representative timeseries of daily and sub-daily point-scale weather variables from observed monthly statistics and other parameters. CLIGEN precipitation timeseries are used as climate input for various risk-assessment modelling applications as an alternative to observe long-term, high temporal resolution records. Here, we queried gridded global climate datasets (TerraClimate, ERA5, GPM-IMERG, and GLDAS) to estimate various 20-year climate statistics and obtain complete CLIGEN input parameter sets with coverage of the African and South American continents at 0.25 arc degree resolution. The estimation of CLIGEN precipitation parameters was informed by a ground-based dataset of >10,000 locations worldwide. The ground observations provided target values to fit regression models that downscale CLIGEN precipitation input parameters. Aside from precipitation parameters, CLIGEN’s parameters for temperature, solar radiation, etc. were in most cases directly calculated according to the original global datasets. Cross-validation for estimated precipitation parameters quantified errors that resulted from applying the estimation approach in a predictive fashion. Based on all training data, the RMSE was 2.23 mm for the estimated monthly average single-event accumulation and 4.70 mm/hr for monthly maximum 30-min intensity. This dataset facilitates exploration of hydrological and soil erosional hypotheses across Africa and South America.

Density Altitude: Climatology of Daily Maximum Values and Evaluation of Approximations for General Aviation

Abstract Density altitude (DA) is an aviation parameter that helps determine specific aircraft performance characteristics for the expected atmospheric conditions. However, there are currently no detailed graphical tools for general aviation (GA) pilot education demonstrating the spatial and temporal variation of DA to help improve situational awareness. In this study, the fifth-generation European Centre for Medium-Range Weather Forecasts atmospheric reanalysis of the global climate (ERA5) dataset is used to construct a 30-yr monthly climatology of DA for the conterminous United States. Several DA characteristics are also investigated, including the effect of humidity on DA, the determination of reasonable worst-case conditions, and the applicability of two DA rules of thumb (ROTs). Maximum values of DA (worst aircraft performance) occur during July, reaching 3600 m over areas with high surface elevations. Humidity, while tertiary to the effects of temperature and pressure, causes the DA to increase from their dry values by more than 140 m as far north as the U.S.-Canada border. The dry DA ROT performs well for all conditions outside of strong tropical cyclones, where GA flights would not be expected. The ROT to correct for the effects of humidity performs well except in high elevations or when the dewpoint temperatures fall outside the applicable range of ≥5°C. When applied outside this range, in some situations, DA errors can be greater than if no humidity correction were applied. Therefore, a new ROT to correct for humidity is introduced here that extends the applicable dewpoint temperature range to ≥−28°C and reduces errors in estimated DA. Significance Statement The impacts of density altitude on aircraft performance have led to numerous general aviation (GA) accidents. This study helps GA pilots better understand the spatial and temporal variability in density altitude, thereby increasing their situational awareness during flight planning. This study also evaluates commonly used approximations to estimate density altitude, so pilots can understand the situations where these approximations are (in)applicable. Results suggest the need for a humidity correction approximation when dewpoint temperatures are <5°C, which is introduced in this study.

The fast‐changing climate reality of Arequipa, Peru

AbstractDespite the availability of global and continental climate datasets and climate change information, locally relevant quantification of historic trends in climate variables is still lacking in developing countries, especially at local scales. This is particularly true in the Department of Arequipa, Peru. An arid region with a booming population, substantial mining activities, and large irrigated agriculture, which is highly susceptible to climate change. This study aims to evaluate climate trends from 1988 to 2017 in the Arequipa Department and provide information that can facilitate stakeholders' adaptation to the rapid‐changing climate. The daily precipitation (Prec), and maximum (Tmax) and minimum (Tmin) daily air temperature data used in this study came from the Servicio Nacional de Meteorología e Hidrología del Perú (SENAMHI) and the National Ocean and Atmospheric Administration's (NOAA) Global Summary of the Day (GSOD). Data passed through a quality checking process for removal of implausible data, data gap filling, and inhomogeneity detection. The Mann–Kendall test, at a significance level of 0.10, was used to determine trends and the Theil–Sen slope (Sen's slope) was used to estimate the magnitude of the change. Sen's slope was also calculated spatially, using the gridded Arequipa Climate Maps (ACM) dataset. Results indicate that precipitation seasonality has been increasing, as the observed increase in annual precipitation is happening mostly in the rainy season (December–March) and the start and end of the rainy season are delayed. Positive temperature trends were dominant in the whole region. Tmin is increasing more than Tmax, especially at higher altitudes. Exceptions to increasing temperatures were found in areas influenced by irrigation projects that underwent great expansion. The effect of increasing temperature on glaciers was evaluated by mapping the change in the area with average annual temperature below 0°C between the decades of 1988–1997 and 2008–2017, which reduced by 73.2%, with small areas disappearing and larger contiguous areas shrinking.

Segmentation of global climate dataset into contiguous spatial units having quantitatively homogeneous climates

AbstractDiscretization of global climate simplifies the spatial variability of climates. The Köppen–Geiger classification (KGC), and more recently, clustering of global climate grids, discretize the landmass into several zones corresponding to qualitatively distinct climate types. However, quantitatively, in terms of local climatograms, such zones have significant climatic variability. In this paper, we propose a segmentation of the global climatic grid as a means to delineate contiguous segments—spatial units characterized by quantitatively homogeneous climates. The level of segments' climatic homogeneity is controlled by a single parameter δ whose value determines a maximum allowable level of segments' inhomogeneity. Segmentation is implemented using a graph‐based segmentation algorithm customized to the climatic data. Application of segmentation with three different values of δ is demonstrated on the WorldClim 1.4 2.5 arc minutes resolution data. Major uses for climate segmentation are: (a) precise delineation of highly homogeneous, local climatic units, (b) checking to what degree delineation of vegetation zones can be based on climate alone, and (c) data compression of large climatic grids. The first use is demonstrated using the island of Great Britain as an example. The second use is demonstrated by the juxtaposition of KGC zones with segments. We discuss the third use in the context of using segmentation as a pre‐processing to climate clustering.

Assessment on spatial extent of arid and semi-arid climatic zones of India using GIS


 In this study, a quantitative assessment of spatial extent of arid and semi-arid climatic zones of India was performed for the period from 1988 to 2018 using potential evapo-transpiration (PET) calculated by Modified Penman Method, estimated from global climate data sets. Climatic water balances computed for 625 stations across the country are used for classifying to bio-climate types based on moisture index and areas falling under arid climatic zones in India are delineated using ArcGIS 10.5. It was noticed a considerable changes in the country’s arid and semi-arid climatic zones between the two periods; 1992 and 2018. Overall, there has been a net percent change in hyper arid, typic arid and semi arid (dry) areas is 5.62, 1.62 and 7.17 percent, respectively. Dryness and wetness are increasing in different parts of the country. There is also change in rainfall, PET and moisture index over a period of time which is vital to determine aridity pattern of any region. Thus, results are of great significance for studying the assessment of temporal and spatial dry climatic water balance of India, which can help immensely in the management of water resources and sustainability of crop production under changing climatic conditions.

Land suitability determination for citrus cultivation using a GIS-based multi-criteria analysis in Mersin, Turkey

Determination of suitable locations for agricultural products is an important subject to study, as it is essential to identify optimal areas for crop productivity and sustainable farming activities. This study proposes a site selection assessment model for citrus cultivation using Geographic Information System (GIS) Based Multi-Criteria Decision Analysis (MCDA) in the Mersin province of Turkey. The purpose of the study is to create a Multi-Criteria Assessment-Based Spatial Decision Support System which takes into account the opinions of local citrus farmers as well as those of experts. 14 preference factors and 10 exclusion factors were determined for the analysis of land suitability for citrus farming. The weighting of preference factors for citrus cultivation was determined using the Analytic Hierarchy Process, and a land suitability map was generated using the weighted overlay method in the GIS environment. Based on the opinions of farmers and experts, higher weights were assigned to the average temperature in the flowering period, the average temperature in the winter season, and elevation, as compared to all other preference factors. The results also demonstrate that WorldClim, a global climate dataset, is an essential source of data in such studies, as climate conditions play a crucial role in potential citrus expansion. The final land suitability map shows that the most suitable and moderately suitable areas for citrus cultivation have surface areas of 85 824 ha and 85 430 ha, respectively. The Receiver Operating Characteristic (ROC) curve was used to evaluate validation of the land suitability for citrus and the Area Under the Curve (AUC) value was calculated at 0.87, indicating substantial consistency. This study proposes a methodology for the determination of suitable sites for citrus cultivation which relies on realistic data. The land suitability map produced in the study will be essential for the creation of sustainable agriculture strategies in the eastern Mediterranean region.

Higher plant photosynthetic capability in autumn responding to low atmospheric vapor pressure deficit

It has been long established that the terrestrial vegetation in spring has stronger photosynthetic capability than in autumn. However, this study challenges this consensus by comparing photosynthetic capability of terrestrial vegetation between the spring and autumn seasons based on measurements of 100 in situ eddy covariance towers over global extratropical ecosystems. At the majority of these sites, photosynthetic capability, indicated by light use efficiency (LUE) and apparent quantum efficiency, is significantly higher in autumn than in spring, due to lower atmosphere vapor pressure deficit (VPD) at the same air temperature. Seasonal VPD differences also substantially explain the interannual variability of the differences in photosynthetic capability between spring and autumn. We further reveal that VPD in autumn is significantly lower than in spring over 74.14% of extratropical areas, based on a global climate dataset. In contrast, LUE derived from a data-driven vegetation production dataset is significantly higher in autumn in over 61.02% of extratropical vegetated areas. Six Earth system models consistently projected continuous larger VPD values in spring compared with autumn, which implies that the impacts on vegetation growth will long exist and should be adequately considered when assessing the seasonal responses of terrestrial ecosystems to future climate conditions.

Biophysical permafrost map indicates ecosystem processes dominate permafrost stability in the Northern Hemisphere

The stability of permafrost is of fundamental importance to socio-economic well-being and ecological services, involving broad impacts to hydrological cycling, global budgets of greenhouse gases and infrastructure safety. This study presents a biophysical permafrost zonation map that uses a rule-based geographic information system (GIS) model integrating global climate and ecological datasets to classify and map permafrost regions (totaling 19.76 × 106 km2, excluding glaciers and lakes) in the Northern Hemisphere into five types: climate-driven (CD) (19% of area), CD/ecosystem-modified (41%), CD/ecosystem protected (3%), ecosystem-driven (29%), and ecosystem-protected (8%). Overall, 81% of the permafrost regions in the Northern Hemisphere are modified, driven, or protected by ecosystems, indicating the dominant role of ecosystems in permafrost stability in the Northern Hemisphere. Permafrost driven solely by climate occupies 19% of permafrost regions, mainly in High Arctic and high mountains areas, such as the Qinghai–Tibet Plateau. This highlights the importance of reducing ecosystem disturbances (natural and human activity) to help slow permafrost degradation and lower the related risks from a warming climate.

A statistics-based reconstruction of high-resolution global terrestrial climate for the last 800,000 years

Curated global climate data have been generated from climate model outputs for the last 120,000 years, whereas reconstructions going back even further have been lacking due to the high computational cost of climate simulations. Here, we present a statistically-derived global terrestrial climate dataset for every 1,000 years of the last 800,000 years. It is based on a set of linear regressions between 72 existing HadCM3 climate simulations of the last 120,000 years and external forcings consisting of CO2, orbital parameters, and land type. The estimated climatologies were interpolated to 0.5° resolution and bias-corrected using present-day climate. The data compare well with the original HadCM3 simulations and with long-term proxy records. Our dataset includes monthly temperature, precipitation, cloud cover, and 17 bioclimatic variables. In addition, we derived net primary productivity and global biome distributions using the BIOME4 vegetation model. The data are a relevant source for different research areas, such as archaeology or ecology, to study the long-term effect of glacial-interglacial climate cycles for periods beyond the last 120,000 years.