Blessing and curse of bioclimatic variables: a comparison of different calculation schemes and datasets for species distribution modeling within the extended Mediterranean area

Bioclimatic variables (BCVs) are the most widely used predictors within the field of species distribution modeling, but recent studies imply that BCVs alone are not sufficient to describe these limits. Unfortunately, the most popular database, WorldClim, offers only a limited selection of bioclimatological predictors; thus, other climatological datasets should be considered, and, for data consistency, the BCVs should also be derived from the respective datasets. Here, we investigate how well the BCVs are represented by different datasets for the extended Mediterranean area within the period 1970–2020, how different calculation schemes affect the representation of BCVs, and how deviations among the datasets differ regionally. We consider different calculation schemes for quarters/months, the annual mean temperature (BCV-1), and the maximum temperature of the warmest month (BCV-5). Additionally, we analyzed the effect of different temporal resolutions for BCV-1 and BCV-5. Differences resulting from different calculation schemes are presented for ERA5-Land. Selected BCVs are analyzed to show differences between WorldClim, ERA5-Land, E-OBS, and CRU. Our results show that (a) differences between the two calculation schemes for BCV-1 diminish as the temporal resolution decreases, while the differences for BCV-5 increase; (b) with respect to the definition of the respective month/quarter, intra-annual shifts induced by the calculation schemes can have substantially different effects on the BCVs; (c) all datasets represent the different BCVs similarly, but with partly large differences in some subregions; and (d) the largest differences occur when specific month/quarters are defined by precipitation. In summary, (a) since the definition of BCVs matches different calculation schemes, transparent communication of the BCVs calculation schemes is required; (b) the calculation, integration, or elimination of BCVs has to be examined carefully for each dataset, region, period, or species; and (c) the evaluated datasets provide, except in some areas, a consistent representation of BCVs within the extended Mediterranean region

Environmental drivers of circum-Antarctic glacier and ice shelf front retreat over the last two decades

The safety band of Antarctica, consisting of floating glacier tongues and ice shelves, buttresses ice discharge of the Antarctic Ice Sheet. Recent disintegration events of ice shelves along with glacier retreat indicate a weakening of this important safety band. Predicting calving front retreat is a real challenge due to complex ice dynamics in a data-scarce environment that are unique for each ice shelf and glacier. We explore the extent to which easy-to-access remote sensing and modeling data can help to define environmental conditions leading to calving front retreat. For the first time, we present a circum-Antarctic record of glacier and ice shelf front change over the last two decades in combination with environmental variables such as air temperature, sea ice days, snowmelt, sea surface temperature, and wind direction. We find that the Antarctic Ice Sheet area decreased by −29 618 ± 1193 km2 in extent between 1997–2008 and gained an area of 7108 ± 1029 km2 between 2009 and 2018. Retreat concentrated along the Antarctic Peninsula and West Antarctica including the biggest ice shelves (Ross and Ronne). In several cases, glacier and ice shelf retreat occurred in conjunction with one or several changes in environmental variables. Decreasing sea ice days, intense snowmelt, weakening easterlies, and relative changes in sea surface temperature were identified as enabling factors for retreat. In contrast, relative increases in mean air temperature did not correlate with calving front retreat. For future studies a more appropriate measure for atmospheric forcing should be considered, including above-zero-degree days and temperature extreme events. To better understand drivers of glacier and ice shelf retreat, it is critical to analyze the magnitude of basal melt through the intrusion of warm Circumpolar Deep Water that is driven by strengthening westerlies and to further assess surface hydrology processes such as meltwater ponding, runoff, and lake drainage

Late quaternary climate, precipitation δ18O, and Indian monsoon variations over the Tibetan Plateau

The Himalaya–Tibet orogen contains one of the largest modern topographic and climate gradients on Earth. Proxy data from the region provide a basis for understanding Tibetan Plateau paleo climate and paleo elevation reconstructions. Paleo climate model comparisons to proxy data compliment sparsely located data and can improve climate reconstructions. This study investigates temporal changes in precipitation, temperature and precipitation δ18O(δ18Op) over the Himalaya–Tibet from the Last Glacial Maximum (LGM) to present. We conduct a series of atmospheric General Circulation Model (GCM, ECHAM5-wiso) experiments at discrete time slices including a Pre-industrial (PI, Pre-1850 AD), Mid Holocene (MH, 6 ka BP) and LGM (21 ka BP) simulations. Model predictions are compared with existing proxy records. Model results show muted climate changes across the plateau during the MH and larger changes occurring during the LGM. During the LGM surface temperatures are ∼2.0–4.0◦C lower across the Himalaya and Tibet, and >5.0◦C lower at the northwest and northeast edge of the Tibetan Plateau. LGM mean annual precipitation is 200–600 mm/yr lower over on the Tibetan Plateau. Model and proxy data comparison shows a good agreement for the LGM, but large differences for the MH. Large differences are also present between MH proxy studies near each other. The precipitation weighted annual mean δ18Op lapse rate at the Himalaya is about 0.4h/km larger during the MH and 0.2h/km smaller during the LGM than during the PI. Finally, rainfall associated with the continental Indian monsoon (between 70◦E–110◦E and 10◦N–30◦N) is about 44% less in the LGM than during PI times. The LGM monsoon period is about one month shorter than in PI times. Taken together, these results document significant spatial and temporal changes in temperature, precipitation, and δ18Op over the last ∼21 ka. These changes are large enough to impact interpretations of proxy data and the intensity of the Indian monsoon

Precipitation δ18O over the Himalaya-Tibet orogen from ECHAM5-wiso simulations: Statistical analysis of temperature, topography and precipitation

Variations in oxygen isotope compositions (δ18O) provide insight into modern climate and past changes in climate and topography. In addition, in regions such as Tibet, geologic archives of isotope ratios record climate change driven by plateau uplift and therefore also provide information about the surface uplift history. A good understanding of modern-day controls on δ18O is crucial for interpreting geologic δ18O in this context. We use the ECHAM5-wiso global atmospheric general circulation model to calculate δ18O in precipitation (δ18Op) for the present-day climate. In the region of the Tibetan Plateau, spatial variations of monthly means of δ18Op are statistically related to spatial variations of 2 m air temperature and precipitation rate, as well as to topography. The size and location of investigated regions are varied in our study to capture regional differences in these relationships and the processes governing the modern δ18Op. In addition to correlation analyses, a cross-validated stepwise multiple regression is carried out using δ18Op as the predictand, and topography and atmospheric variables (temperature and precipitation amount) as predictors. The 2 m air temperature and topography yield the highest spatial correlation coefficients of >0.9 and 90% of the δ18Op spatial variance in the same regions. The 2 m air temperature is the dominant predictor and contributes 93.6% to the total explained spatial variance on average. The results demonstrate that most of the δ18Op pattern on and around the Tibetan Plateau can be explained by variation in 2 m air temperature and altitude. Correlation of the dependent predictors indicate that in low-altitude regions where topography does not determine temperature variability, the high correlation of temperature and δ18Op may partially be explained by variations in precipitation rates

Modern precipitation δ18O and trajectory analysis over the Himalaya-Tibet Orogen from ECHAM5-wiso simulations

Variations in oxygen isotope ratios (δ18O) measured from modern precipitation and geologic archives provide a promising tool for understanding modern and past climate dynamics and tracking elevation changes over geologic time. In areas of extreme topography, such as the Tibetan Plateau, the interpretation of δ18O has proven challenging. This study investigates the climate controls on temporal (daily and 6 h intervals) and spatial variations in present-day precipitation δ18O (δ18Op) across the Tibetan Plateau using a 30 year record produced from the European Centre/Hamburg ECHAM5-wiso global atmospheric general circulation model (GCM). Results indicate spatial and temporal agreement between model-predicted δ18Op and observations. Large daily δ18Op variations of 25 to +5‰ occur over the Tibetan Plateau throughout the 30 simulation years, along with interannual δ18Op variations of ~2‰. Analysis of extreme daily δ18Op indicates that extreme low values coincide with extreme highs in precipitation amount. During the summer, monsoon vapor transport from the north and southwest of the plateau generally corresponds with high δ18Op, whereas vapor transport from the Indian Ocean corresponds with average to low δ18Op. Thus, vapor source variations are one important cause of the spatial-temporal differences in δ18Op. Comparison of GCM and Rayleigh Distillation Model (RDM)-predicted δ18Op indicates a modest agreement for the Himalaya region (averaged over 86°–94°E), confirming application of the simpler RDM approach for estimating δ18Op lapse rates across Himalaya

Future climate change and its impact on runoff generation from the debris-covered Inylchek glaciers, central Tian Shan, Kyrgyzstan

The heavily debris-covered Inylchek glaciers in the central Tian Shan are the largest glacier system in the Tarim catchment. It is assumed that almost 50% of the discharge of Tarim River are provided by glaciers. For this reason, climatic changes, and thus changes in glacier mass balance and glacier discharge are of high impact for the whole region. In this study, a conceptual hydrological model able to incorporate discharge from debris-covered glacier areas is presented. To simulate glacier melt and subsequent runoff in the past (1970/1971-1999/2000) and future (2070/2071-2099/2100), meteorological input data were generated based on ECHAM5/MPI-OM1 global climate model projections. The hydrological model HBV-LMU was calibrated by an automatic calibration algorithm using runoff and snow cover information as objective functions. Manual fine-tuning was performed to avoid unrealistic results for glacier mass balance. The simulations show that annual runoff sums will increase significantly under future climate conditions. A sensitivity analysis revealed that total runoff does not decrease until the glacier area is reduced by 43%. Ice melt is the major runoff source in the recent past, and its contribution will even increase in the coming decades. Seasonal changes reveal a trend towards enhanced melt in spring, but a change from a glacial-nival to a nival-pluvial runoff regime will not be reached until the end of this century.Publisher PDFPeer reviewe

Modelling wind speed across Zambia: Implications for wind energy

Wind energy is a key option in global dialogues about climate change mitigation. Here, we combined observations from surface wind stations, reanalysis datasets, and state‐of‐the‐art regional climate models from the Coordinated Regional Climate Downscaling Experiment (CORDEX Africa) to study the current and future wind energy potential in Zambia. We found that winds are dominated by southeasterlies and are rarely strong with an average speed of 2.8 m·s−1. When we converted the observed surface wind speed to a turbine hub height of 100 m, we found a ~38% increase in mean wind speed for the period 1981–2000. Further, both simulated and observed wind speed data show statistically significant increments across much of the country. The only areas that divert from this upward trend of wind speeds are the low land terrains of the Eastern Province bordering Malawi. Examining projections of wind power density (WPD), we found that although wind speed is increasing, it is still generally too weak to support large‐scale wind power generation. We found a meagre projected annual average WPD of 46.6 W·m−2. The highest WPDs of ~80 W·m−2 are projected in the northern and central parts of the country while the lowest are to be expected along the Luangwa valley in agreement with wind speed simulations. On average, Zambia is expected to experience minor WPD increments of 0.004 W·m−2 per year from 2031 to 2050. We conclude that small‐scale wind turbines that accommodate cut‐in wind speeds of 3.8 m·s−1 are the most suitable for power generation in Zambia. Further, given the limitations of small wind turbines, they are best suited for rural and suburban areas of the country where obstructions are few, thus making them ideal for complementing the government of the Republic of Zambia's rural electrification efforts.Wind speed is very slow in Zambia. It is increasing but remains unlikely to support large commercial wind farms especially not at the 100‐m hub height. Any efforts for wind power generation in Zambia should be towards ultra‐tall wind turbines fitted with larger rotors.Alexander von Humboldt‐Stiftun