A methodology for the classification of convective structures using meteorological radar: Application to heavy rainfall events on the Mediterranean coast of the Iberian Peninsula

During the period 1996-2000, forty-three heavy rainfall events have been detected in the Internal Basins of Catalonia (Northeastern of Spain). Most of these events caused floods and serious damage. This high number leads to the need for a methodology to classify them, on the basis of their surface rainfall distribution, their internal organization and their physical features. The aim of this paper is to show a methodology to analyze systematically the convective structures responsible of those heavy rainfall events on the basis of the information supplied by the meteorological radar. The proposed methodology is as follows. Firstly, the rainfall intensity and the surface rainfall pattern are analyzed on the basis of the raingauge data. Secondly, the convective structures at the lowest level are identified and characterized by using a 2-D algorithm, and the convective cells are identified by using a 3-D procedure that looks for the reflectivity cores in every radar volume. Thirdly, the convective cells (3-D) are associated with the 2-D structures (convective rainfall areas). This methodology has been applied to the 43 heavy rainfall events using the meteorological radar located near Barcelona and the SAIH automatic raingauge network.

Anomalous rainfall and associated atmospheric circulation in the northeast Spanish Mediterranean area and its relationship to sediment fluidization events in a lake

Changes in the dynamics of sediment transport in a Mediterranean lake (sediment fluidization events) are linked to atmospheric circulations patterns (trough monthly precipitation). In the basins of Lake Banyoles, located in the northeast of Spain, water enters mainly through subterranean springs, and associated fluctuations in the vertical migration of sediment distribution (fluidization events) present episodic behavior as a result of episodic rainfall in the area. The initiation of the fluidization events takes place when the monthly rainfall is ∼2.7 times greater than the mean monthly rainfall of the rainiest months in the area, especially in spring (April and May), October, and December. The duration of these events is found to be well correlated with the accumulated rainfall of the preceding 10 months before the process initiation. The rainfall, in turn, is mainly associated with six atmospheric circulation patterns among the 19 fundamental circulations that emerged in an earlier study focused on significant rainfall days in Mediterranean Spain. Among them, accentuated surface lows over the northeast of Spain, general northeasterly winds by low pressure centered to the east of Balearic Islands and short baroclinic waves over the Iberian Peninsula, with easterly flows over the northeastern coast of Spain, are found the most relevant atmospheric circulations that drive heavy rainfall events

Predictive ability of severe rainfall events over Catalonia for the year 2008

This paper analyses the predictive ability of quantitative precipitation forecasts (QPF) and the so-called "poor-man" rainfall probabilistic forecasts (RPF). With this aim, the full set of warnings issued by the Meteorological Service of Catalonia (SMC) for potentially-dangerous events due to severe precipitation has been analysed for the year 2008. For each of the 37 warnings, the QPFs obtained from the limited-area model MM5 have been verified against hourly precipitation data provided by the rain gauge network covering Catalonia (NE of Spain), managed by SMC. For a group of five selected case studies, a QPF comparison has been undertaken between the MM5 and COSMO-I7 limited-area models. Although MM5's predictive ability has been examined for these five cases by making use of satellite data, this paper only shows in detail the heavy precipitation event on the 9¿10 May 2008. Finally, the "poor-man" rainfall probabilistic forecasts (RPF) issued by SMC at regional scale have also been tested against hourly precipitation observations. Verification results show that for long events (>24 h) MM5 tends to overestimate total precipitation, whereas for short events (¿24 h) the model tends instead to underestimate precipitation. The analysis of the five case studies concludes that most of MM5's QPF errors are mainly triggered by very poor representation of some of its cloud microphysical species, particularly the cloud liquid water and, to a lesser degree, the water vapor. The models' performance comparison demonstrates that MM5 and COSMO-I7 are on the same level of QPF skill, at least for the intense-rainfall events dealt with in the five case studies, whilst the warnings based on RPF issued by SMC have proven fairly correct when tested against hourly observed precipitation for 6-h intervals and at a small region scale. Throughout this study, we have only dealt with (SMC-issued) warning episodes in order to analyse deterministic (MM5 and COSMO-I7) and probabilistic (SMC) rainfall forecasts; therefore we have not taken into account those episodes that might (or might not) have been missed by the official SMC warnings. Therefore, whenever we talk about "misses", it is always in relation to the deterministic LAMs' QPFs.

Effect of radar rainfall time resolution on the predictive capability of a distributed hydrologic model

The performance of a hydrologic model depends on the rainfall input data, both spatially and temporally. As the spatial distribution of rainfall exerts a great influence on both runoff volumes and peak flows, the use of a distributed hydrologic model can improve the results in the case of convective rainfall in a basin where the storm area is smaller than the basin area. The aim of this study was to perform a sensitivity analysis of the rainfall time resolution on the results of a distributed hydrologic model in a flash-flood prone basin. Within such a catchment, floods are produced by heavy rainfall events with a large convective component. A second objective of the current paper is the proposal of a methodology that improves the radar rainfall estimation at a higher spatial and temporal resolution. Composite radar data from a network of three C-band radars with 6-min temporal and 2 × 2 km2 spatial resolution were used to feed the RIBS distributed hydrological model. A modification of the Window Probability Matching Method (gauge-adjustment method) was applied to four cases of heavy rainfall to improve the observed rainfall sub-estimation by computing new Z/R relationships for both convective and stratiform reflectivities. An advection correction technique based on the cross-correlation between two consecutive images was introduced to obtain several time resolutions from 1 min to 30 min. The RIBS hydrologic model was calibrated using a probabilistic approach based on a multiobjective methodology for each time resolution. A sensitivity analysis of rainfall time resolution was conducted to find the resolution that best represents the hydrological basin behaviour.

Characteristics of 2D convective structures in Catalonia (NE Spain): an analysis using radar data and GIS

Flood simulation studies use spatial-temporal rainfall data input into distributed hydrological models. A correct description of rainfall in space and in time contributes to improvements on hydrological modelling and design. This work is focused on the analysis of 2-D convective structures (rain cells), whose contribution is especially significant in most flood events. The objective of this paper is to provide statistical descriptors and distribution functions for convective structure characteristics of precipitation systems producing floods in Catalonia (NE Spain). To achieve this purpose heavy rainfall events recorded between 1996 and 2000 have been analysed. By means of weather radar, and applying 2-D radar algorithms a distinction between convective and stratiform precipitation is made. These data are introduced and analyzed with a GIS. In a first step different groups of connected pixels with convective precipitation are identified. Only convective structures with an area greater than 32 km2 are selected. Then, geometric characteristics (area, perimeter, orientation and dimensions of the ellipse), and rainfall statistics (maximum, mean, minimum, range, standard deviation, and sum) of these structures are obtained and stored in a database. Finally, descriptive statistics for selected characteristics are calculated and statistical distributions are fitted to the observed frequency distributions. Statistical analyses reveal that the Generalized Pareto distribution for the area and the Generalized Extreme Value distribution for the perimeter, dimensions, orientation and mean areal precipitation are the statistical distributions that best fit the observed ones of these parameters. The statistical descriptors and the probability distribution functions obtained are of direct use as an input in spatial rainfall generators.

The role of different factors related to social impact of heavy events: considerations about the intensity thresholds in densely populated areas

In the assessment of social impact caused by meteorological events, factors of different natures need to be considered. Not only does hazard itself determine the impact that a severe weather event has on society, but also other features related to vulnerability and exposure. The requests of data related to insurance claims received in meteorological services proved to be a good indicator of the social impact that a weather event causes, according to studies carried out by the Social Impact Research Group, created within the framework of the MEDEX project. Taking these requests as proxy data, diverse aspects connected to the impact of heavy rain events have been studied. The rainfall intensity, in conjunction with the population density, has established itself as one of the key factors in social impact studies. One of the conclusions we obtained is that various thresholds of rainfall should be applied for areas of varying populations. In this study, the role of rainfall intensity has been analysed for a highly populated urban area like Barcelona. A period without significant population changes has been selected for the study to minimise the effects linked to vulnerability and exposure modifications. First, correlations between rainfall recorded in different time intervals and requests were carried out. Afterwards, a method to include the intensity factor in the social impact index was suggested based on return periods given by intensity duration frequency (IDF) curves.

Proposal of three thermodynamic variables to discriminate between storms associated with hail and storms with intense rainfall in Catalonia (Proposta de tres variables termodinàmiques per discriminar entre tempestes amb pedregades i tempestes amb pluges intenses)

Between late spring and early fall, the development of storms is common in Catalonia. Despite the fact that they usually produce heavy showers of short duration, they can also involve severe weather with ice pellets or hail. While the latter usually affect inland regions, and there are numerous publications on these cases; the analysis of events affecting the coast and causing damage to public and private properties is not so well developed. The aim of this study is to provide additional thermodynamic indicators that help differentiate storms with hail from storms without hail, considering cases that have affected various regions of Catalonia, mainly coastal areas. The aim is to give more information to improve prognosis and the ability to detail information in these situations. The procedure developed involved the study of several episodes of heavy rainfall and hail that hit Catalonia during the 2003-2009 period, mainly in the province of Girona, and validated the proposal during the campaign of late summer and fall of 2009, as well as 2012. For each case, several variables related to temperature, humidity and wind were analyzed at different levels of the atmosphere, while the information provided by the radio sounding in Barcelona was also taken into account. From this study, it can be concluded that the temperature difference between 500 hPa and 850 hPa, the humidity in the lower layers of the atmosphere and the LI index are good indicators for the detection of storms with associated hail.

Towards a database on societal impact of Mediterranean floods in the framework of the HYMEX project

The NW Mediterranean region experiences every year heavy rainfall and flash floods that occasionally produce catastrophic damages. Less frequent are floods that affect large regions. Although a large number of databases devoted exclusively to floods or considering all kind of natural hazards do exist, usually they only record catastrophic flood events. This paper deals with the new flood database that is being developed within the framework of HYMEX project. Results are focused on four regions representative of the NW sector of Mediterranean Europe: Catalonia, Spain; the Balearic Islands, Spain; Calabria, Italy; and Languedoc-Roussillon, Midi-Pyrenées and PACA, France. The common available 30-yr period starts in 1981 and ends in 2010. The paper shows the database structure and criteria, the comparison with other flood databases, some statistics on spatial and temporal distribution, and an identification of the most important events. The paper also provides a table that includes the date and affected region of all the catastrophic events identified in the regions of study, in order to make this information available for all audiences.

Torrential rainfall in northeast of the Iberian Peninsula: synoptic patterns and WeMO influence

This study presents a catalogue of synoptic patterns of torrential rainfall in northeast of the Iberian Peninsula (IP). These circulation patterns were obtained by applying a T-mode Principal Component Analysis (PCA) to a daily data grid (NCEP/NCAR reanalysis) at sea level pressure (SLP). The analysis made use of 304 days which recorded >100 mm in one or more stations in provinces of Barcelona, Girona and Tarragona (coastland area of Catalonia) throughout the 1950-2005 period. The catalogue comprises 7 circulation patterns showing a great variety of atmospheric conditions and seasonal or monthly distribution. Likewise, we computed the mean index value of the Western Mediterranean Oscillation index (WeMOi) for the synoptic patterns obtained by averaging all days grouped in each pattern. The results showed a clear association between the negative values of this teleconnection index and torrential rainfall in northeast of the IP. We therefore put forward the WeMO as an essential tool for forecasting heavy rainfall in northeast of Spain

Predicting the potential for lightning activity in mediterranean storms based on the WRF model dynamic and microphysical fields

A new parameter is introduced: the lightning potential index (LPI), which is a measure of the potential for charge generation and separation that leads to lightning flashes in convective thunderstorms. The LPI is calculated within the charge separation region of clouds between 0 C and 20 C, where the noninductive mechanism involving collisions of ice and graupel particles in the presence of supercooled water is most effective. As shown in several case studies using the Weather Research and Forecasting (WRF) model with explicit microphysics, the LPI is highly correlated with observed lightning. It is suggested that the LPI may be a useful parameter for predicting lightning as well as a tool for improving weather forecasting of convective storms and heavy rainfall.

Analysis of warm season thunderstorms using an object-oriented tracking method based on radar and total lightning data

Monitoring thunderstorms activity is an essential part of operational weather surveillance given their potential hazards, including lightning, hail, heavy rainfall, strong winds or even tornadoes. This study has two main objectives: firstly, the description of a methodology, based on radar and total lightning data to characterise thunderstorms in real-time; secondly, the application of this methodology to 66 thunderstorms that affected Catalonia (NE Spain) in the summer of 2006. An object-oriented tracking procedure is employed, where different observation data types generate four different types of objects (radar 1-km CAPPI reflectivity composites, radar reflectivity volumetric data, cloud-to-ground lightning data and intra-cloud lightning data). In the framework proposed, these objects are the building blocks of a higher level object, the thunderstorm. The methodology is demonstrated with a dataset of thunderstorms whose main characteristics, along the complete life cycle of the convective structures (development, maturity and dissipation), are described statistically. The development and dissipation stages present similar durations in most cases examined. On the contrary, the duration of the maturity phase is much more variable and related to the thunderstorm intensity, defined here in terms of lightning flash rate. Most of the activity of IC and CG flashes is registered in the maturity stage. In the development stage little CG flashes are observed (2% to 5%), while for the dissipation phase is possible to observe a few more CG flashes (10% to 15%). Additionally, a selection of thunderstorms is used to examine general life cycle patterns, obtained from the analysis of normalized (with respect to thunderstorm total duration and maximum value of variables considered) thunderstorm parameters. Among other findings, the study indicates that the normalized duration of the three stages of thunderstorm life cycle is similar in most thunderstorms, with the longest duration corresponding to the maturity stage (approximately 80% of the total time).

GIS-based debris flow source and runout susceptibility assessment from DEM data - a case study in NW Nicaragua

In October 1998, Hurricane Mitch triggered numerous landslides (mainly debris flows) in Honduras and Nicaragua, resulting in a high death toll and in considerable damage to property. The potential application of relatively simple and affordable spatial prediction models for landslide hazard mapping in developing countries was studied. Our attention was focused on a region in NW Nicaragua, one of the most severely hit places during the Mitch event. A landslide map was obtained at 1:10 000 scale in a Geographic Information System (GIS) environment from the interpretation of aerial photographs and detailed field work. In this map the terrain failure zones were distinguished from the areas within the reach of the mobilized materials. A Digital Elevation Model (DEM) with 20 m×20 m of pixel size was also employed in the study area. A comparative analysis of the terrain failures caused by Hurricane Mitch and a selection of 4 terrain factors extracted from the DEM which, contributed to the terrain instability, was carried out. Land propensity to failure was determined with the aid of a bivariate analysis and GIS tools in a terrain failure susceptibility map. In order to estimate the areas that could be affected by the path or deposition of the mobilized materials, we considered the fact that under intense rainfall events debris flows tend to travel long distances following the maximum slope and merging with the drainage network. Using the TauDEM extension for ArcGIS software we generated automatically flow lines following the maximum slope in the DEM starting from the areas prone to failure in the terrain failure susceptibility map. The areas crossed by the flow lines from each terrain failure susceptibility class correspond to the runout susceptibility classes represented in a runout susceptibility map. The study of terrain failure and runout susceptibility enabled us to obtain a spatial prediction for landslides, which could contribute to landslide risk mitigation.

Scale effect in hazard assessment - application to daily rainfall

Daily precipitation is recorded as the total amount of water collected by a rain-gauge in 24h. Events are modelled as a Poisson process and the 24h precipitation by a Generalized 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. We use the fact that a log-scale is better suited to the type of variable analyzed to overcome this inconsistency, thus showing that using the appropriate natural scale can be extremely important for proper hazard assessment. The approach is illustrated with precipitation data from the Eastern coast of the Iberian Peninsula affected by severe convective precipitation. The estimation is carried out by using Bayesian techniques

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).

Preface: Forecast and projection in climate scenario of Mediterranean intense events: uncertainties and propagation on environment (the MEDUP project)

The Mediterranean basin is a particularly vulnerable region to climate change, partly due to its quite unique character that results both from physiographic conditions and societal development. The region features indeed a near-closed sea surrounded by very urbanised littorals and mountains from which numerous rivers originate. This results in a lot of interactions and feedbacks between oceanic-atmospheric-hydrological processes that play a predominant role on climate and extreme events that frequently cause heavy dam- ages and human losses in the Mediterranean ...