Trends in indices of daily precipitation extremes in Catalonia (NE Spain), 1951-2003

The aim of this paper is to quantitatively characterize the climatology of daily precipitation indices in Catalonia (northeastern Iberian Peninsula) from 1951 to 2003. This work has been performed analyzing a subset of the ETCCDI (Expert Team on Climate Change Detection and Indices) precipitation indices calculated from a new interpolated dataset of daily precipitation, namely SPAIN02, regular at 0.2° horizontal resolution (around 20 km) and from two high-quality stations: the Ebro and Fabra observatories. Using a jack-knife technique, we have found that the sampling error of the SPAIN02 regional averaged is relatively low. The trend analysis has been implemented using a Circular Block Bootstrap procedure applicable to non-normal distributions and autocorrelated series. A running trend analysis has been applied to analyze the trend persistence. No general trends at a regional scale are observed, considering the annual or the seasonal regional averaged series of all the indices for all the time windows considered. Only the consecutive dry days index (CDD) at annual scale shows a locally coherent spatial trend pattern; around 30% of the Catalonia area has experienced an increase of around 2¿3 days decade¿1. The Ebro and Fabra observatories show a similar CDD trend, mainly due to the summer contribution. Besides this, a significant decrease in total precipitation (around ¿10 mm decade¿1) and in the index "highest precipitation amount in five-day period" (RX5DAY, around ¿5 mm decade¿1), have been found in summer for the Ebro observatory.

Enhanced radar precipitation estimates using a combined clutter and beam blockage correction technique

Weather radar observations are currently the most reliable method for remote sensing of precipitation. However, a number of factors affect the quality of radar observations and may limit seriously automated quantitative applications of radar precipitation estimates such as those required in Numerical Weather Prediction (NWP) data assimilation or in hydrological models. In this paper, a technique to correct two different problems typically present in radar data is presented and evaluated. The aspects dealt with are non-precipitating echoes - caused either by permanent ground clutter or by anomalous propagation of the radar beam (anaprop echoes) - and also topographical beam blockage. The correction technique is based in the computation of realistic beam propagation trajectories based upon recent radiosonde observations instead of assuming standard radio propagation conditions. The correction consists of three different steps: 1) calculation of a Dynamic Elevation Map which provides the minimum clutter-free antenna elevation for each pixel within the radar coverage; 2) correction for residual anaprop, checking the vertical reflectivity gradients within the radar volume; and 3) topographical beam blockage estimation and correction using a geometric optics approach. The technique is evaluated with four case studies in the region of the Po Valley (N Italy) using a C-band Doppler radar and a network of raingauges providing hourly precipitation measurements. The case studies cover different seasons, different radio propagation conditions and also stratiform and convective precipitation type events. After applying the proposed correction, a comparison of the radar precipitation estimates with raingauges indicates a general reduction in both the root mean squared error and the fractional error variance indicating the efficiency and robustness of the procedure. Moreover, the technique presented is not computationally expensive so it seems well suited to be implemented in an operational environment.

The effect of scale in daily precipitation hazard assessment

Daily precipitation is recorded as the total amount of water collected by a rain-gauge in 24 h. Events are modelled as a Poisson process and the 24 h precipitation by a Generalised Pareto Distribution (GPD) of excesses. Hazard assessment is complete when estimates of the Poisson rate and the distribution parameters, together with a measure of their uncertainty, are obtained. The shape parameter of the GPD determines the support of the variable: Weibull domain of attraction (DA) corresponds to finite support variables as should be for natural phenomena. However, Fréchet DA has been reported for daily precipitation, which implies an infinite support and a heavy-tailed distribution. Bayesian techniques are used to estimate the parameters. The approach is illustrated with precipitation data from the Eastern coast of the Iberian Peninsula affected by severe convective precipitation. The estimated GPD is mainly in the Fréchet DA, something incompatible with the common sense assumption of that precipitation is a bounded phenomenon. The bounded character of precipitation is then taken as a priori hypothesis. Consistency of this hypothesis with the data is checked in two cases: using the raw-data (in mm) and using log-transformed data. As expected, a Bayesian model checking clearly rejects the model in the raw-data case. However, log-transformed data seem to be consistent with the model. This fact may be due to the adequacy of the log-scale to represent positive measurements for which differences are better relative than absolute

Mesoscale numerical analysis of the historical November 1982 heavy precipitation event over Andorra (Eastern Pyrenees)

From 6 to 8 November 1982 one of the most catastrophic flash-flood events was recorded in the Eastern Pyrenees affecting Andorra and also France and Spain with rainfall accumulations exceeding 400 mm in 24 h, 44 fatalities and widespread damage. This paper aims to exhaustively document this heavy precipitation event and examines mesoscale simulations performed by the French Meso-NH non-hydrostatic atmospheric model. Large-scale simulations show the slow-evolving synoptic environment favourable for the development of a deep Atlantic cyclone which induced a strong southerly flow over the Eastern Pyrenees. From the evolution of the synoptic pattern four distinct phases have been identified during the event. The mesoscale analysis presents the second and the third phase as the most intense in terms of rainfall accumulations and highlights the interaction of the moist and conditionally unstable flows with the mountains. The presence of a SW low level jet (30 m s-1) around 1500 m also had a crucial role on focusing the precipitation over the exposed south slopes of the Eastern Pyrenees. Backward trajectories based on Eulerian on-line passive tracers indicate that the orographic uplift was the main forcing mechanism which triggered and maintained the precipitating systems more than 30 h over the Pyrenees. The moisture of the feeding flow mainly came from the Atlantic Ocean (7-9 g kg-1) and the role of the Mediterranean as a local moisture source was very limited (2-3 g kg-1) due to the high initial water vapour content of the parcels and the rapid passage over the basin along the Spanish Mediterranean coast (less than 12 h).

Spatiotemporal variability in radial growth of trees in the Central Pyrenees: climatic and site influences.

To understand how tree growth has responded to recent climate warming, an understanding of the tree-climate-site complex is necessary. To achieve this, radial growth variability among 204 trees established before 1850 was studied in relation to both climatic and site factors. Seventeen forest stands were sampled in the Spanish Central Pyrenees. Three species were studied: Pinus uncinata, Abies alba, and Pinus sylvestris. For each tree, a ring-width residual chronology was built. All trees cross-dated well, indicating a common influence of the regional climate. For the 1952-1993 period, the radial growth of all species, especially P. uncinata, was positively correlated with warm Novembers during the year before ring formation and warm Mays of the year the annual ring formed. Differences in species-stand elevation modulated the growth-climate associations. Radial growth in P. uncinata at high elevation sites was reduced when May temperatures were colder and May precipitation more abundant. In the 20th century, two contrasting periods in radial growth were observed: one (1900-1949) with low frequency of narrow and wide rings, low mean annual sensitivity, and low common growth variation; and another (1950-1994) with the reverse characteristics. The increased variability in radial growth since the 1950s was observed for all species and sites, which suggests a climatic cause. The low shared variance among tree chronologies during the first half of the 20th century may result from a"relaxation" of the elevation gradient, allowing local site conditions to dominate macroclimatic influence. These temporal trends may be related to the recently reported increase of climatic variability and warmer conditions. This study emphasizes the need to carefully assess the relationships between radial growth and site conditions along ecological gradients to improve dendroclimatic reconstructions.

Daily precipitation concentration across Europe 1971-2010

Daily Precipitation Concentration Index (CI) was used in this paper to investigate the statistical structure of daily precipitation across Europe based on 530 daily rainfall series for the period 1971-2010. Annual CI shows a North- West to South-East gradient (excluding Turkey and Greece). The same gradient is also observed in winter, spring and autumn, while in summer the gradient is North-South. Highest annual and seasonal daily concentrations of rainfall were detected in the western Mediterranean basin, mainly along Spanish and French coastlands. Latitude and distance from the sea seems to play a major role on spatial CI distribution; at subregional scale also relief plays an important role. The Mann-Kendall test did not identify uniform significant pattern in temporal trend across Europe for 1971-2010 period. The only broad areas with increasing annual and seasonal CI values are located in northern and south-western France and northern coastlands of the Iberian Peninsula. This findings suggest that daily precipitation distribution has not significantly changed during the 1971-2010 over Europe.

Mesoscale numerical analysis of the historical November 1982 heavy precipitation event over Andorra (Eastern Pyrenees)

From 6 to 8 November 1982 one of the most catastrophic flash-flood events was recorded in the Eastern Pyrenees affecting Andorra and also France and Spain with rainfall accumulations exceeding 400 mm in 24 h, 44 fatalities and widespread damage. This paper aims to exhaustively document this heavy precipitation event and examines mesoscale simulations performed by the French Meso-NH non-hydrostatic atmospheric model. Large-scale simulations show the slow-evolving synoptic environment favourable for the development of a deep Atlantic cyclone which induced a strong southerly flow over the Eastern Pyrenees. From the evolution of the synoptic pattern four distinct phases have been identified during the event. The mesoscale analysis presents the second and the third phase as the most intense in terms of rainfall accumulations and highlights the interaction of the moist and conditionally unstable flows with the mountains. The presence of a SW low level jet (30 m s-1) around 1500 m also had a crucial role on focusing the precipitation over the exposed south slopes of the Eastern Pyrenees. Backward trajectories based on Eulerian on-line passive tracers indicate that the orographic uplift was the main forcing mechanism which triggered and maintained the precipitating systems more than 30 h over the Pyrenees. The moisture of the feeding flow mainly came from the Atlantic Ocean (7-9 g kg-1) and the role of the Mediterranean as a local moisture source was very limited (2-3 g kg-1) due to the high initial water vapour content of the parcels and the rapid passage over the basin along the Spanish Mediterranean coast (less than 12 h).

Hydrological evidence for a North Atlantic oscillation during the Little Ice Age outside its range observed since 1850.

An annual-resolved precipitation reconstruction for the last 800 yr in Southern Spain has been performed using stable carbon isotope (δ13C) of Pinus nigra tree rings. The reconstruction exhibits high- to low-frequency variability and distinguishes a Little Ice Age (LIA, AD 13501850) characterized by lower averaged rainfall than both in the transition from the Medieval Climate Anomaly to the LIA and in the 20th century. The driest conditions are recorded during the Maunder solar Minimum (mid 17thearly 18th centuries), in good agreement with the Spanish documentary archive. Similar linkage between solar activity (maximum/minimum) and precipitation (increase/decrease) is observed throughout the entire LIA. Additionally, the relationship between the hydrological pattern in the Iberian Peninsula and Morocco during the LIA suggests different spatial distribution of precipitation in the south-eastern sector of the North Atlantic region such as it is known currently. Whereas in the instrumental record the precipitation evolves similarly in both regions and opposite to the North Atlantic oscillation (NAO) index, the coldest periods of the LIA shows a contrasting pattern with drier conditions in the South of Spain and wetter in Northern Africa. We suggest an extreme negative NAO conditions, accompanied by a southward excursion of the winter rainfall band beyond that observed in the last century, can explain this contrast. The sustained NAO conditions could have been triggered by solar minima and higher volcanic activity during the LIA.