Search

On November 15, 2017, a high intensity convective storm, reaching 300 mm in 13 h in the core zone of the event, hit the western part of the region of Attica in Greece, causing a catastrophic flash flood in the town of Mandra and a tragic loss of 24 people, making it the most deadly flood in the country, in a period of 40 years. The research team surveyed the area during and after the flood using a combination of systematic ground and aerial observations with the aid of an unmanned aerial vehicle (UAV), aiming to reconstruct the basic physical and hydrological characteristics of the flood and its impacts. The analysis produced detailed flood extent and depth maps that provided a comprehensive description of the physical characteristics of floodwaters across the inundated area. Peak discharge was estimated, using a UAV-derived digital surface model, at two locations, corresponding to the two main tributaries and indicated an impressive hydrological response, between 9 and 10 m3/s/km2. Impact analysis on the basis of these observations showed an extensive diversity, including effects in geomorphology, vegetation, buildings, infrastructure and human population. Analysis of meteorological, botanical and geomorphological evidence lead to the conclusion that this flash flood was a very rare event. Results demonstrate that the combination of aerial and ground observations allow an enhanced and holistic reconstruction of a flash flood and its impacts with high accuracy, leading to the conclusion that the approach used has a significant potential in many aspects of flood disaster investigations.

Extreme atmospheric and surface conditions over the vast Saharan desert result in the development of one of the deepest atmospheric boundary layers (Saharan Atmospheric Boundary Layer, SABL) on the planet. The land cover in the Sahara mainly consists of wide interchanging areas of stony and sandy desert intercepted by narrower bare rock formations. Significant changes in surface albedo, surface temperature, turbulence fluxes and wind speed have been observed with remote measurement platforms like aircraft and satellites. This was a motivation to simulate the convective boundary layer (CBL) that develops over the heterogeneous Saharan surface with a two-way coupled system of a large eddy simulation code (LES) and a land surface model (LSM). Initial and boundary conditions are provided by airborne observations and mesoscale atmospheric simulations. In order to investigate the effect of surface heterogeneity on the vertical structure of the SABL, a large scale (larger than 10 km) idealized warm surface anomaly is simulated and analyzed. A land strip with a low surface albedo produces almost doubled fluxes in density and stronger convergence near the ground than over the surrounding brighter surfaces. First and second order turbulence statistics reveal that strong thermals penetrate the inversion layer above the CBL producing a vigorous exchange between unstable and overlying stable air. Furthermore, a convective internal boundary layer (CIBL) is formed over the warm strip. The downwind development of the CIBL is studied and CBL depth equations from literature are revised. The average turbulent structure of the SABL reveals flow changes locally and downwind of the surface heterogeneity (strip), which can be more significant than the dispersive fluxes due to surface heterogeneity, and may have an important effect on processes such as dust uplift and its long range transport that influence weather and climate globally.

The performance of the Tropical Rainfall Measuring Mission (TRMM) Precipitation Radar (PR) rainfall estimation algorithm is assessed, locally, in Crete island, south Greece, using data from a 2D-video disdrometer and a ground-based, X-band, polarimetric radar. A three-parameter, normalized Gamma drop size distribution is fitted to the disdrometer rain spectra; the latter are classified in stratiform and convective rain types characterized by different relations between distribution parameters. The method of moments estimates more accurately the distribution parameters than the best fit technique, which exhibits better agreement with and is more biased by the observed droplet distribution at large diameter values. Power laws between the radar reflectivity factor ( Z ) and the rainfall rate ( R ) are derived from the disdrometer data. A significant diversity of the prefactor and the exponent of the estimated power laws is observed, depending on the scattering model and the regression technique. The Z – R relationships derived from the disdrometer data are compared to those obtained from TRMM-PR data. Generally, the power laws estimated from the two datasets are different. Specifically, the greater prefactor found for the disdrometer data suggests an overestimation of rainfall rate by the TRMM-PR algorithm for light and moderate stratiform rain, which was the main rain type in the disdrometer dataset. Finally, contemporary data from the TRMM-PR and a ground-based, X-band, polarimetric radar are analyzed. Comparison of the corresponding surface rain rates for a rain event with convective characteristics indicates a large variability of R in a single TRMM-PR footprint, which typically comprises several hundreds of radar pixels. Thus, the coarse spatial resolution of TRMM-PR may lead to miss of significant high local peaks of convective rain. Also, it was found that the high temporal variability of convective rain may introduce significant errors in the estimation of bias of the satellite rainfall estimates with respect to data from ground-based radars.

In this work the capability of reliable rainfall measurements with small weather radars in complex terrain for flood forecasting purposes is examined. Rain measurements were carried out during winter–spring 2007 with a mobile X-band dual-polarization radar in the northwestern mountainous part of the island of Crete in Greece. In this area a 2D-video disdrometer and a network of raingauges was installed for radar calibration and evaluation of rainfall measurements, respectively. The complex terrain of the experimental site may significantly reduce the performance of rain measurements due to ground clutter and partial or total beam blockage. A beam blockage algorithm using high resolution terrain data was applied in order to find areas with significant terrain effects and estimate correction of the radar measurements. Rain attenuation correction was based on modern sophisticated algorithms using differential phase measurements. The accuracy of rainfall estimation from standard or polarimetric algorithms at plan position indicator (PPI) scans was examined for high-temporal resolution (1 min) rainfall rates and accumulated rainfall values for winter and spring time rain events. For high elevation measurements, which were required in order to avoid terrain effects in large areas of interest, a correction for the vertical-profile-of-reflectivity (VPR) was also applied to the radar data. An average VPR model used in the corresponding correction of reflectivity was constructed based on range–height indicator (RHI) scans. It was concluded that quantitative high resolution X-band radar observations of rainfall in complex terrain is possible after careful application of all the above processing steps.