J. I. López-Moreno

A New Global 0.5° Gridded Dataset (1901–2006) of a Multiscalar Drought Index: Comparison with Current Drought Index Datasets Based on the Palmer Drought Severity Index

Abstract A monthly global dataset of a multiscalar drought index is presented and compared in terms of spatial and temporal variability with the existing continental and global drought datasets based on the Palmer drought severity index (PDSI). The presented dataset is based on the standardized precipitation evapotranspiration index (SPEI). The index was obtained using the Climatic Research Unit (CRU) TS3.0 dataset at a spatial resolution of 0.5°. The advantages of the new dataset are that (i) it improves the spatial resolution of the unique global drought dataset at a global scale; (ii) it is spatially and temporally comparable to other datasets, given the probabilistic nature of the SPEI; and, in particular, (iii) it enables the identification of various drought types, given the multiscalar character of the SPEI. The dataset is freely available on the Web page of the Spanish National Research Council (CSIC) in three different formats [network Common Data Form (netCDF), binary raster, and plain text].

Impact of climate and land use change on water availability and reservoir management: Scenarios in the Upper Aragón River, Spanish Pyrenees

Streamflows in a Mediterranean mountain basin in the central Spanish Pyrenees were projected under various climate and land use change scenarios. Streamflow series projected for 2021-2050 were used to simulate the management of the Yesa reservoir, which is critical to the downstream supply of irrigation and domestic water. Streamflows were simulated using the Regional Hydro-Ecologic Simulation System (RHESSys). The results show that increased forest cover in the basin could decrease annual streamflow by 16%, mainly in early spring, summer and autumn. Regional climate models (RCMs) project a trend of warming and drying in the basin for the period 2021-2050, which will cause a 13.8% decrease in annual streamflow, mainly in late spring and summer. The combined effects of forest regeneration and climate change are expected to reduce annual streamflows by 29.6%, with marked decreases affecting all months with the exception of January and February, when the decline will be moderate. Under these streamflow reduction scenarios it is expected that it will be difficult for the Yesa reservoir to meet the current water demand, based on its current storage capacity (476 hm(3)). If the current project to enlarge the reservoir to a capacity of 1059 hm(3) is completed, the potential to apply multi-annual streamflow management, which will increase the feasibility of maintaining the current water supply. However, under future climate and land cover scenarios, reservoir storage will rarely exceed half of the expected capacity, and the river flows downstream of the reservoir is projected to be dramatically reduced.

Mapping the annual evolution of snow depth in a small catchment in the Pyrenees using the long-range terrestrial laser scanning

This paper presents the methodology used to develop snow depth distribution maps for a small catchment in the Central Spanish Pyrenees covering 55 ha in a 1:10,000 scale. The Main Map was obtained using LiDAR (light detection and ranging) technology from a long-range Terrestrial Laser Scanner (TLS) in six field surveys undertaken during the 2012 winter–spring period. This technique enabled the acquisition of information at a very high resolution concerning the spatial variability of snow cover, providing snow depth information for remote areas where data acquisition is complex and hazardous. We describe the methodological steps and the quality assessment applied in developing the maps. Comparison with manual measurements confirmed the reliability of the snow depth maps, including areas located at large distances from the scanner (800 m). This method provides a promising tool for future investigations of snow dynamics in mountainous environments.

Spatial Predictions of Extreme Wind Speeds over Switzerland Using Generalized Additive Models

The purpose of this work is to present a methodology aimed at predicting extreme wind speeds over Switzerland. Generalized additive models are used to regionalize wind statistics for Swiss weather stations using a number of variables that describe the main physiographical features of the country. This procedure enables one to present the results for Switzerland in the form of a map that provides the 98th percentiles of daily maximum wind speeds (W98) at a 10-m anemometer height for cells with a 50-m grid interval. This investigation comprises three major steps. First, meteorological data recorded by the weather stations was gathered to build local wind statistics at each station. Then, data describing the topographic and landscape characteristics of the country were prepared using geographic information systems (GIS). Third, appropriate regression models were selected to make spatially explicit predictions of extreme wind speeds in Switzerland. The predictions undertaken in this study provide realistic values of the W98. The effects of topography on the results are particularly conspicuous. Wind speeds increase with altitude and are greatest on mountain peaks in the Alps, as would be intuitively expected. Relative errors between observations and model results calculated for the meteorological stations do not exceed 30%, and only 12 out of 70 stations exhibit errors that exceed 20%. The combination of GIS techniques and statistical models used to predict a highly uncertain variable, such as extreme wind speed, yields interesting results that can be extended to other fields, such as the assessment of storm damage on infrastructures.

Canopy influence on snow depth distribution in a pine stand determined from terrestrial laser data

In this study, we analyzed the effects of the forest canopy and trunks of a pine stand in the central Spanish Pyrenees on the snow depth (SD) distribution. Using LiDAR technology with a terrestrial laser scanner (TLS), high-resolution data on the SD distribution were acquired during the 2011-2012 and 2012-2013 snow seasons, which were 2 years having very contrasting climatic and snow accumulation conditions. Average SD evolution in open and canopy areas was characterized. Principal component analysis was applied to identify days having similar spatial patterns of SD distribution. There was a clear contrast in the temporal variability of the snowpack in different areas of the forest stand, corresponding generally to beneath the canopy, and in open sites. The canopy and openings showed markedly different accumulation and melting, with higher snow accumulation found in openings. Differences ranged from 14 to 80% reduction (average 49%) in the SD beneath the canopy relative to open sites. The difference in SD between open and canopy areas increased throughout the snow season. The surveyed days were classified in terms of SD distribution, and included days associated with: high SD, low SD, intense melting conditions and periods when the SD distribution was driven by wind conditions. The SD increased with distance from the trunks to a distance of 3.5-4.5 m, coinciding with the average size of the crown of individual trees.

Hydrological and depositional processes associated with recent glacier recession in Yanamarey catchment, Cordillera Blanca (Peru)

In this study, we investigate changes in the glaciated surface and the formation of lakes in the headwater of the Querococha watershed in Cordillera Blanca (Peru) using 24 Landsat images from 1975 to 2014. Information of glacier retreat was integrated with available climate data, the first survey of recent depositional dynamics in proglacial Yanamarey Lake (4600m a.s.l.), and a relatively short hydrological record (2002-2014) at the outlet of Yanamarey Lake. A statistically significant temperature warming (0.21°C decade-1 for mean annual temperature) has been detected in the region, and it caused a reduction of the glacierized area since 1975 from 3.5 to 1.4km-2. New small lakes formed in the deglaciated areas, increasing the flooded area from1.8ha in 1976 to 2.8ha in 2014. A positive correlation between annual rates of glacier recession and runoff was found. Sediment cores revealed a high sedimentation rate (>1cmyr-1) and two contrasted facies, suggesting a shift toward a reduction of meltwater inputs and higher hydrological variability likely due to an increasing role of precipitation on runoff during the last decades. Despite the age control uncertainties, the main transition likely occurred around 1998-2000, correlating with the end of the phase with maximum warming rates and glacier retreat during the 1980s and 1990s, and the slowing down of expansion of surface lake-covered surface. With this hydrological - paleolimnological approach we have documented the association between recent climate variability and glacier recession and the rapid transfer of hydroclimate signal to depositional and geochemical processes in high elevation Andean environments. This, study also alerts about water quality risks as proglacial lakes act as secondary reservoirs that trap trace and minor elements in high altitude basins.

Spatio-temporal snowmelt variability across the headwaters of the Southern Rocky Mountains

Understanding the rate of snowmelt helps inform how water stored as snow will transform into streamflow. Data from 87 snow telemetry (SNOTEL) stations across the Southern Rocky Mountains were used to estimate spatio-temporal melt factors. Decreases in snow water equivalent were correlated to temperature at these monitoring stations for eight half-month periods from early March through late June. Time explained 70% of the variance in the computed snow melt factors. A residual linear correlation model was used to explain subsequent spatial variability. Longitude, slope, and land cover type explained further variance. For evergreen trees, canopy density was relevant to find enhanced melt rates; while for all other land cover types, denoted as non-evergreen, lower melt rates were found at high elevation, high latitude and north facing slopes, denoting that in cold environments melting is less effective than in milder sites. A change in the temperature sensor about mid-way through the time series (1990 to 2013) created a discontinuity in the temperature dataset. An adjustment to the time series yield larger computed melt factors.

Average monthly and annual climate maps for Bolivia

This study presents monthly and annual climate maps for relevant hydroclimatic variables in Bolivia. We used the most complete network of precipitation and temperature stations available in Bolivia, which passed a careful quality control and temporal homogenization procedure. Monthly average maps at the spatial resolution of 1 km were modeled by means of a regression-based approach using topographic and geographic variables as predictors. The monthly average maximum and minimum temperatures, precipitation and potential exoatmospheric solar radiation under clear sky conditions are used to estimate the monthly average atmospheric evaporative demand by means of the Hargreaves model. Finally, the average water balance is estimated on a monthly and annual scale for each 1 km cell by means of the difference between precipitation and atmospheric evaporative demand. The digital layers used to create the maps are available in the digital repository of the Spanish National Research Council.

Application of terrestrial laser scanner techniques for monitoring dynamic geomorphological processes: snow accumulation and ice masses in mountain areas

This paper presents the application of a long range terrestrial laser scanner for monitoring snow accumulation and evolution of glaciers in the Pyrenees. With this study we are showing the great potential of the methodology presented for studies in mountain environments, particularly in research related to the cryosphere monitoring in relation to climate variability and the complex topography typical of mountain environments. We present a comprehensive