Hydrological impact and potential flooding of convective rain cells in a semi-arid environment

Abstract

t Spatio-temporal storm properties have a large impact on catchment hydrological response. The sensitivity of simulated flash floods to convective rain-cell characteristics is examined for an extreme storm event over a 94 km² semi-arid catchment in southern Israel. High space–time resolution weather radar data were used to derive and model convective rain cells that then served as input into a hydrological model. Based on alterations of location, direction and speed of a major rain cell, identified as the flooding cell for this case, the impacts on catchment rainfall and generated flood were examined. Global sensitivity analysis was applied to identify the most important factors affecting the flash flood peak discharge at the catchment outlet. We found that the flood peak discharge could be increased three-fold by relatively small changes in rain-cell characteristics. We assessed that the maximum flash flood magnitude that this single rain cell can produce is 175 m³/s, and, taking into account the rest of the rain cells, the flash flood peak discharge can reach 260 m³/s.

Editor Z.W. Kundzewicz; Guest editor R.J. Moore

Citation Morin, E. and Yakir, H., 2013. Hydrological impact and potential flooding of convective rain cells in a semi-arid environment. Hydrological Sciences Journal, 59 (7), 1275–1284. http://dx.doi.org/10.1080/02626667.2013.841315     

Radar rainfall estimation for flash flood forecasting in small urban watersheds

Radar rainfall estimation for flash flood forecasting in small, urban catchments is examined through analyses of radar, rain gage and discharge observations from the 14.3 km² Dead Run drainage basin in Baltimore County, Maryland. The flash flood forecasting problem pushes the envelope of rainfall estimation to time and space scales that are commensurate with the scales at which the fundamental governing laws of land surface processes are derived. Analyses of radar rainfall estimates are based on volume scan WSR-88D reflectivity observations for 36 storms during the period 2003–2005. Gage-radar analyses show large spatial variability of storm total rainfall over the 14.3 km² basin for flash flood producing storms. The ability to capture the detailed spatial variation of rainfall for flash flood producing storms by WSR-88D rainfall estimates varies markedly from event to event. As spatial scale decreases from the 14.3 km² scale of the Dead Run watershed to 1 km² (and the characteristic time scale of flash flood producing rainfall decreases from 1 h to 15 min) the predictability of flash flood response from WSR-88D rainfall estimates decreases sharply. Storm to storm variability of multiplicative bias in storm total rainfall estimates is a dominant element of the error structure of radar rainfall estimates, and it varies systematically over the warm season and with flood magnitude. Analyses of the 7 July 2004 and 28 June 2005 storms illustrate microphysical and dynamical controls on radar estimation error for extreme flash flood producing storms.     

Development of a thresholds approach for real‐time flash flood prediction in complex geomorphological river basins

Flash floods represent one of the deadliest and costliest natural disasters worldwide. The hydrological analysis of a flash flood event contributes in the understanding of the runoff creation process. This study presents the analysis of some flash flood events that took place in a complex geomorphological Mediterranean River basin. The objective of the present work is to develop the thresholds for a real‐time flash flood forecasting model in a complex geomorphological watershed, based on high‐frequency data from strategically located hydrological and meteorological telemetric stations. These stations provide hourly real‐time data which were used to determine hydrological and meteorological parameters. The main characteristics of various hydrographs specified in this study were the runoff coefficients, lag time, time to peak, and the maximum potential retention. The estimation of these hydrometeorological parameters provides the necessary information in order to successfully manage flash floods events. Especially, the time to peak is the most significant hydrological parameter that affects the response time of an oncoming flash flood event. A study of the above parameters is essential for the specification of thresholds which are related to the geomorphological characteristics of the river basin, the rainfall accumulation of an event, the rainfall intensity, the threshold runoff through karstic area, the season during which the rainfall takes place and the time intervals between the rainstorms that affect the soil moisture conditions. All these factors are combined into a real‐time‐threshold flash flood prediction model. Historical flash flood events at the downstream are also used for the validation of the model. An application of the proposed model is presented for the Koiliaris River basin in Crete, Greece. Copyright © 2011 John Wiley & Sons, Ltd.     

Fully coupled free‐surface flow and sediment transport modelling of flash floods in a desert stream in the Mojave Desert,California

For most of the year, a dry‐bed desert wash is void of water flow. Intensive rain events, however, could trigger significant flash floods that bring about highly complicated hydrodynamics and morphodynamics processes within a desert stream. We present a fully coupled three‐phase flow model of air, water, and sediment to simulate numerically the propagation of a flash flood in a field‐scale fluvial desert stream, the so‐called Tex Wash located in the Mojave Desert, California, United States. The turbulent flow of the flash flood is computed using the three‐dimensional unsteady Reynolds‐averaged Navier–Stokes equations closed with the shear stress transport k ? ω model. The free surface of the flash flood at the interface of air and water phases is computed with the level‐set method, which enables instantaneous tracking of the water surface as the flash flood propagates over the dry bed of the desert stream. The evolution of the desert fluvial stream's morphology, due to the action of the propagating flash flood on the mobile bed, is calculated using a Eulerian morphodynamics model based on the curvilinear immersed boundary method. The capabilities of the proposed numerical framework are demonstrated by applying it to simulate a flash flood event in a 0.65‐km ‐long reach of the Tex Wash, the intricate channel morphology of which is obtained using light imaging detection and ranging technology. The simulated region of the stream includes a number of bridge foundations. The simulation results of the model for the flash flood event revealed the formation of a highly complex flow field and scour patterns within the stream. Moreover, our simulation results showed that most scour processes take place during the steady phase of the flash flood, that is, after the flash flood fills the stream. The transient phase of the flash flood is rather short and contributes to a very limited amount of erosion within the desert stream.     

Reconstructing extreme post‐wildfire floods: a comparison of convective and mesoscale events

In much of western United States destructive floods after wildfire are frequently caused by localized, short‐duration convective thunderstorms; however, little is known about post‐fire flooding from longer‐duration, low‐intensity mesoscale storms. In this study we estimate and compare peak flows from convective and mesoscale floods following the 2012 High Park Fire in the ungaged 15.5 km² Skin Gulch basin in the northcentral Colorado Front Range. The convective storm on 6 July 2012 came just days after the wildfire was contained. Radar data indicated that the total rainfall was 20–47 mm, and the maximum rainfall intensities (upwards of 50 mm h^?1) were concentrated over portions of the watershed that burned at high severity. The mesoscale storm on 9–15 September 2013 produced 220–240 mm of rain but had maximum 15‐min intensities of only 25–32 mm h^?1. Peak flows for each flood were estimated using three independent techniques. Our best estimate using a 2D hydraulic model was 28 m³ s^?1 km^?2 for the flood following the convective storm, placing it among the largest rainfall‐runoff floods per unit area in the United States. In contrast, the flood associated with the mesoscale flood was only 6 m³ s^?1 km^?2, but the long‐duration flood caused extensive channel incision and widening, indicating that this storm was much more geomorphically effective. The peak flow estimates for the 2013 flood had a higher relative uncertainty and this stemmed from whether we used pre‐ or post‐flood channel topography. The results document the extent to which a high and moderate severity forest fire can greatly increase peak flows and alter channel morphology, illustrate how indirect peak flow estimates have larger errors than is generally assumed, and indicate that the magnitude of post‐fire floods and geomorphic change can be affected by the timing, magnitude, duration, and sequence of rainstorms. Copyright © 2017 John Wiley & Sons, Ltd.     

The relationship between rainfall inputs and flood generation in south–east Spain

Rainfall and flood data are relatively sparse in semi‐arid areas; hence there have been relatively few investigations into the relationships between rainfall inputs and flood generation in these environments. Previous work has shown that flood properties are influenced by a combination of precipitation characteristics including amount, intensity, duration and spatial distribution. Therefore floods may be produced by high intensity, short duration storms, or longer duration, low intensity rainfall. Most of this research has been undertaken in small catchments in either hyper‐arid or relatively high rainfall Mediterranean climates. This paper presents results from a 6 year data record in south‐east Spain from research conducted in two basins, the Rambla Nogalte (171 km²) and the Rambla de Torrealvilla (200 km²). Data cover an area of approximately 500 km² and an annual average rainfall of 300 mm. At coarse temporal resolutions gauges spread over large areas record similar patterns of rainfall, although spells of rain show much more complexity; pulses of rain within storms can vary considerably in total rainfall, intensity and duration over the same area. The analysis for south‐east Spain shows that most storms occur over a period of less than 24 h, but that the number of rainfall events declines as the duration exceeds 8 h. This is at odds with data on floods for the study area suggesting that they are produced by storms lasting longer than 18 h. However, one flood event was produced by a very short (15 min) storm with high intensity rainfall. Most floods tended to occur in May/June or September, which coincides with wetter months of the year (September, October, December and May). Floods are also more highly related to the total rainfall occurring in a spell of rain, than to intensity. The complexity of storm rainfall increases with the storm total, which makes it difficult to generalize on the importance of rainfall intensity for flood generation. Copyright © 2007 John Wiley & Sons, Ltd.     

Correlations between deep convection and lightning activity on a global scale

Satellite observations of cloud top temperature and lightning flash distribution are used to examine the relationship between deep convection and lightning activity over the tropical regions of the northern and southern hemispheres. In agreement with previous work, the analysis of the results shows that, in the summer of both hemispheres, the lightning activity in continental deep convective storms is more intense than that in marine deep convective storms by a factor of between 7 and 10. Furthermore, it was observed that on average the daily lightning rate per 1°×1° grid cell for the southern hemisphere (SH) is about 20% greater than that of the northern hemisphere (NH), which can be attributed to a larger fractional cover by deep convective clouds in the SH. By using a set of independent indicators, it is shown that deep convection and lightning activity over land are well correlated (with correlation coefficients of 0.8 and 0.6 for NH and SH, respectively). This suggests the capacity for observations to act as a possible method of monitoring continental deep convective clouds, which play a key role in regulating the Earth’s climate. Since lightning can be monitored easily from ground networks and satellites, it could be a useful tool for validating the performance of model convective schemes and for monitoring changes in climate parameters.     

Testing Taylor’s hypothesis in Amazonian rainfall fields during the WETAMC/LBA experiment

Taylor’s hypothesis (TH) for rainfall fields states that the spatial correlation of rainfall intensity at two points at the same instant of time can be equated with the temporal correlation at two instants of time at some fixed location. The validity of TH is tested in a set of 12 storms developed in Rondonia, southwestern Amazonia, Brazil, during the January–February 1999 Wet Season Atmospheric Meso-scale Campaign. The time Eulerian and Lagrangian Autocorrelation Functions (ACF) are estimated, as well as the time-averaged space ACF, using radar rainfall rates of storms spanning between 3.2 and 23 h, measured at 7–10-min time resolution, over a circle of 100 km radius, at 2 km spatial resolution. TH does not hold in 9 out of the 12 studied storms, due to their erratic trajectories and very low values of zonal wind velocity at 700 hPa, independently from underlying atmospheric stability conditions. TH was shown to hold for 3 storms, up to a cutoff time scale of 10–15 min, which is closely related to observed features of the life cycle of convective cells in the region. Such cutoff time scale in Amazonian storms is much shorter than the 40 min identified in mid-latitude convective storms, due to much higher values of CAPE and smaller values of storm speed in Amazonian storms as compared to mid-latitude ones, which in turn contribute to a faster destruction of the rainfall field isotropy. Storms satisfying TH undergo smooth linear trajectories over space, and exhibit the highest negative values of maximum, mean and minimum zonal wind velocity at 700 hPa, within narrow ranges of atmospheric stability conditions. Non-dimensional parameters involving CAPE (maximum, mean and minimum) and CINE (mean) are identified during the storms life cycle, for which TH holds: CAPE mean/CINE mean = [30–35], CAPE max/CINE mean = [32–40], and CAPE min/CINE mean = [22–28]. These findings are independent upon the timing of storms within the diurnal cycle. Also, the estimated Eulerian time ACF’s decay faster than the time-averaged space and the Lagrangian time ACF’s, irrespectively of TH validity. The Eulerian ACF’s exhibit shorter e-folding times, reflecting smaller correlations over short time scales, but also shorter scale of fluctuation, reflecting less persistence in time than over space. No significant associations (linear, exponential or power law) were found between estimated e-folding times and scale of fluctuation, with all estimates of CAPE and CINE. Secondary correlation maxima appear between 50 and 70 min in the Lagrangian time ACF’s for storms satisfying TH. No differences were found in the behavior of each of the three ACF’s for storms developed during either the Easterly or Westerly zonal wind regimes which characterize the development of meso-scale convective systems over the region. These results have important implications for modelling and downscaling rainfall fields over tropical land areas.     

The June 2013 Alberta catastrophic flooding event – part 2: fine‐scale precipitation and associated features

Data obtained from a variety of sources including the Canadian Lightning Detection Network, weather radars, weather stations and operational numerical weather model analyses were used to address the evolution of precipitation during the June 2013 southern Alberta flood. The event was linked to a mid‐level closed low pressure system to the west of the region and a surface low pressure region initially to its south. This configuration brought warm, moist unstable air into the region that led to dramatic, organized convection with an abundance of lightning and some hail. Such conditions occurred in the southern parts of the region whereas the northern parts were devoid of lightning. Initially, precipitation rates were high (extreme 15‐min rainfall rates up to 102 mm h^?1 were measured) but decreased to lower values as the precipitation shifted to long‐lived stratiform conditions. Both the convective and stratiform precipitation components were affected by the topography. Similar flooding events, such as June 2002, have occurred over this region although the 2002 event was colder and precipitation was not associated with substantial convection over southwest Alberta. Copyright © 2016 Her Majesty the Queen in Right of Canada. Hydrological Processes. © John Wiley & Sons, Ltd.     

Hydrological control of flooding: Tuscany,October 1992

In the study of flash-flood occurrence in small catchments the lack of flow measurements is often one of the main limiting factors. Prior to estimating the forecasting potentialities and techniques for such events, an accurate reconstruction of past event flood dynamics is first required. This issue is here addressed by analyzing, with the use of a distributed hydrological model, the hydrometeorological conditions in which a severe flash-flood occurred, on October 1992, on a 48 square kilometers catchment in the Arno basin. Such an event was caused by the persistence of intense convective clusters on the background of widespread rain bands of frontal origin. The distributed hydrological model here adopted is devoted to simulate the evolution and the variability of the primary processes involved in the runoff cycle. Together with the hydrological model structure, other particular aspects of the event reconstruction procedure are discussed: the managing and processing of the information coming from different sensors, with different temporal and spatial resolutions; the identification of local precipitation dynamics (frontal or convective) within small areas of integrated radar and rain gauges data fields; the interpolation of rain gauge data on the basis of the radar-estimated spatial correlation. The results of the distributed modeling, concerning the estimate of the flood wave at various sites, are compared with analogous results obtained with simpler lumped models.     

Using Lightning Data to Better Understand and Predict Flash Floods in the Mediterranean

One of the costliest natural hazards around the globe is flash floods, resulting from localized intense convective precipitation over short periods of time. Since intense convective rainfall (especially over the continents) is well correlated with lightning activity in these storms, a European Union FP6 FLASH project was realized from 2006 to 2010, focusing on using lightning observations to better understand and predict convective storms that result in flash floods. As part of the project, 23 case studies of flash floods in the Mediterranean region were examined. For the analysis of these storms, lightning data were used together with rainfall estimates in order to understand the storms?? development and electrification processes. In addition, these case studies were simulated using mesoscale meteorological models to better understand the local and synoptic conditions leading to such intense and damaging storms. As part of this project, tools for short-term predictions (nowcasts) of intense convection across the Mediterranean and Europe, and long-term forecasts (a few days) of the likelihood of intense convection, were developed and employed. The project also focused on educational outreach through a special Web site http://flashproject.org supplying real-time lightning observations, real-time experimental nowcasts, medium-range weather forecasts and educational materials. While flash floods and intense thunderstorms cannot be prevented, long-range regional lightning networks can supply valuable data, in real time, for warning the public, end-users and stakeholders of imminent intense rainfall and possible flash floods.     

The significance of rainstorm variations to shallow translational hillslope failure

Landsliding in eastern Scotland results from high-magnitude rainstorms generated under either cyclonic or anticyclonic conditions, particularly during the summer. Data from Aviemore indicated that cyclonic storms produce higher rainfall totals than anticyclonic storms, as well as being of longer duration and lower intensity. The distribution of rain during individual storms also varies with the synoptic conditions under which the storms are produced. These different rainfall characteristics produce different geomorphic responses, which can be investigated in detail using physical based modelling. In this paper, a physically based coupled hydrology–stability model is used to assess the significance of these rainfall characteristics to soil moisture response and slope instability for mature podsols. The results provide evidence that rainstorms of different synoptic origin produce varying hydrological response, involving both the extent and the timing of moisture content change. This affects the depth and timing of slope failure, with anticyclonic storms promoting a large, rapid response in the factor of safety at shallow depths within the soil. Cyclonic storms produce a more gradual response, with the region of probable failure being deeper. Futhermore, each of these storm types is associated with different rainfall distributions, and this is also shown to have a significant effect on the timing and depth of slope instability.     

Characteristics of rain events at an inland locality in northeastern Borneo,Malaysia

Understanding the intensity and duration of tropical rain events is critical to modelling the rate and timing of wet‐canopy evaporation, the suppression of transpiration, the generation of infiltration‐excess overland flow and hence to erosion, and to river responsiveness. Despite this central role, few studies have addressed the characteristics of equatorial rainstorms. This study analyses rainfall data for a 5 km² region largely comprising of the 4 km² Sapat Kalisun Experimental Catchment in the interior of northeastern Borneo at sampling frequencies from 1 min^?1 to 1 day^?1. The work clearly shows that most rainfall within this inland, forested area is received during regular short‐duration events (<15 min) that have a relatively low intensity (i.e. less than two 0·2 mm rain‐gauge tips in almost all 5 min periods). The rainfall appears localized, with significant losses in intergauge correlations being observable in minutes in the case of the typical mid‐afternoon, convective events. This suggests that a dense rain‐gauge network, sampled at a high temporal frequency, is required for accurate distributed rainfall‐runoff modelling of such small catchments. Observed rain‐event intensity is much less than the measured infiltration capacities, and thus supports the tenet of the dominance of quick subsurface responses in controlling river behaviour in this small equatorial catchment. Copyright © 2006 John Wiley & Sons, Ltd.     

Critical rainfall statistics for predicting watershed flood responses: rethinking the design storm concept

Recent advances have been made to modernize estimates of probable precipitation scenarios; however, researchers and engineers often continue to assume that rainfall events can be described by a small set of event statistics, typically average intensity and event duration. Given the easy availability of precipitation data and advances in desk‐top computational tools, we suggest that it is time to rethink the ‘design storm’ concept. Design storms should include more holistic characteristics of flood‐inducing rain events, which, in addition to describing specific hydrologic responses, may also be watershed or regionally specific. We present a sensitivity analysis of nine precipitation event statistics from observed precipitation events within a 60‐year record for Tompkins County, NY, USA. We perform a two‐sample Kolmogorov–Smirnov (KS) test to objectively identify precipitation event statistics of importance for two related hydrologic responses: (1) peak outflow from the Six Mile Creek watershed and (2) peak depth within the reservoir behind the Six Mile Creek Dam. We identify the total precipitation depth, peak hourly intensity, average intensity, event duration, interevent duration, and several statistics defining the temporal distribution of precipitation events to be important rainfall statistics to consider for predicting the watershed flood responses. We found that the two hydrologic responses had different sets of statistically significant parameters. We demonstrate through a stochastic precipitation generation analysis the effects of starting from a constrained parameter set (intensity and duration) when predicting hydrologic responses as opposed to utilizing an expanded suite of rainfall statistics. In particular, we note that the reduced precipitation parameter set may underestimate the probability of high stream flows and therefore underestimate flood hazard. Copyright © 2016 John Wiley & Sons, Ltd.     

Flood simulation using the gauge‐adjusted radar rainfall and physics‐based distributed hydrologic model

This paper reports the results of an investigation into flood simulation by areal rainfall estimated from the combination of gauged and radar rainfalls and a rainfall–runoff model on the Anseong‐cheon basin in the southern part of Korea. The spatial and temporal characteristics and behaviour of rainfall are analysed using various approaches combining radar and rain gauges: (1) using kriging of the rain gauge alone; (2) using radar data alone; (3) using mean field bias (MFB) of both radar and rain gauges; and (4) using conditional merging technique (CM) of both radar and rain gauges. To evaluate these methods, statistics and hyetograph for rain gauges and radar rainfalls were compared using hourly radar rainfall data from the Imjin‐river, Gangwha, rainfall radar site, Korea. Then, in order to evaluate the performance of flood estimates using different rainfall estimation methods, rainfall–runoff simulation was conducted using the physics‐based distributed hydrologic model, Vflo^?. The flood runoff hydrograph was used to compare the calculated hydrographs with the observed one. Results show that the rainfall field estimated by CM methods improved flood estimates, because it optimally combines rainfall fields representing actual spatial and temporal characteristics of rainfall. Copyright © 2008 John Wiley & Sons, Ltd.     

Changes to flood peaks of a mountain river: implications for analysis of the 2013 flood in the Upper Bow River,Canada

The mountain headwater Bow River at Banff, Alberta, Canada, was subject to a large flood in June 2013, over which considerable debate has ensued regarding its probability of occurrence. It is therefore instructive to consider what information long‐term streamflow discharge records provide about environmental change in the Upper Bow River basin above Banff. Though protected as part of Banff National Park, since 1885, the basin has experienced considerable climate and land cover changes, each of which has the potential to impact observations, and hence the interpretations of flood probability. The Bow River at Banff hydrometric station is one of Canada's longest‐operating reference hydrological basin network stations and so has great value for assessing changes in flow regime over time. Furthermore, the station measures a river that provides an extremely important water supply for Calgary and irrigation district downstream and so is of great interest for assessing regional water security. These records were examined for changes in several flood attributes and to determine whether flow changes may have been related to landscape change within the basin as caused by forest fires, conversion from grasslands to forest with fire suppression, and regional climate variations and/or trends. Floods in the Upper Bow River are generated by both snowmelt and rain‐on‐snow (ROS) events, the latter type which include flood events generated by spatially and temporally large storms such as occurred in 2013. The two types of floods also have different frequency characteristics. Snowmelt and ROS flood attributes were not correlated significantly with any climate index or with burned area except that snowmelt event duration correlated negatively to the Pacific Decadal Oscillation. While there is a significant negative trend in all floods over the past 100 years, when separated based on generating process, neither snowmelt floods nor large ROS floods associated with mesoscale storms show any trends over time. Despite extensive changes to the landscape of the basin and in within the climate system, the flood regime remains unchanged, something identified at smaller scales in the region but never at larger scales. Copyright © 2016 John Wiley & Sons, Ltd.     

The Hat Yai 2000 flood: the worst flood in Thai history

Hat Yai, the largest commercial and tourist city in southern Thailand, is subjected to regular flood events, primarily during the northeast monsoon period. Flooding in this region is recognized as a serious disaster in terms of frequency, rate of risk, and affected areas. The monsoon of 21–25 November 2000 caused extremely heavy rain in the southern part of Thailand, resulting in a great flood occupying Hat Yai. This caused significant damage. Therefore, the use of both structural and non‐structural measures is mandatory to reduce the economic losses and the risk for society. This paper investigates two modelling approaches for flood prevention and mitigation of Hat Yai city. First, a hard computing approach by a physically distributed model was applied to study the flood behaviour in a two‐dimensional floodplain flow. Second, a soft computing approach using a neuro‐genetic algorithm was used to develop a flood‐forecasting tool. It was found that the great flood of 2000 can be simulated well by the FLO‐2D model. Computed discharges and flood level in the floodplain are close to the observed data. Countermeasures using diversion canals are guaranteed to accelerate the floodwater drainage to Songkla Lake, significantly reducing the flood impact to the people. In addition, the flood forecasting technique developed in this study can give satisfactory results. This would be very useful as a flood‐warning tool for the community Copyright © 2005 John Wiley & Sons, Ltd.     

The cold rain‐on‐snow event of June 2013 in the Canadian Rockies — characteristics and diagnosis

The June 2013 flood in the Canadian Rockies featured rain‐on‐snow (ROS) runoff generation at alpine elevations that contributed to the high streamflows observed during the event. Such a mid‐summer ROS event has not been diagnosed in detail, and a diagnosis may help to understand future high discharge‐producing hydrometeorological events in mountainous cold regions. The alpine hydrology of the flood was simulated using a physically based model created with the modular cold regions hydrological modelling platform. The event was distinctive in that, although at first, relatively warm rain fell onto existing snowdrifts inducing ROS melt; the rainfall turned to snowfall as the air mass cooled and so increased snowcover and snowpacks in alpine regions, which then melted rapidly from ground heat fluxes in the latter part of the event. Melt rates of existing snowpacks were substantially lower during the ROS than during the relatively sunny periods preceding and following the event as a result of low wind speeds, cloud cover and cool temperatures. However, at the basin scale, melt volumes increased during the event as a result of increased snowcover from the fresh snowfall and consequent large ground heat contributions to melt energy, causing snowmelt to enhance rainfall–runoff by one fifth. Flow pathways also shifted during the event from relatively slow sub‐surface flow prior to the flood to an even contribution from sub‐surface and fast overland flow during and immediately after the event. This early summer, high precipitation ROS event was distinctive for the impact of decreased solar irradiance in suppressing melt rates, the contribution of ground heat flux to basin scale snowmelt after precipitation turned to snowfall, the transition from slow sub‐surface to fast overland flow runoff as the sub‐surface storage saturated and streamflow volumes that exceeded precipitation. These distinctions show that summer, mountain ROS events should be considered quite distinct from winter ROS and can be important contributors to catastrophic events. Copyright © 2016 John Wiley & Sons, Ltd.     

Quantifying the role of the snowpack in generating water available for run‐off during rain‐on‐snow events from snow pillow records

Rain‐on‐snow events have generated major floods around the world, particularly in coastal, mountainous regions. Most previous studies focused on a limited number of major rain‐on‐snow events or were based primarily on model results, largely due to a lack of long‐term records from lysimeters or other instrumentation for quantifying event water balances. In this analysis, we used records from five automated snow pillow sites in south coastal British Columbia, Canada, to reconstruct event water balances for 286 rain‐on‐snow events over a 10‐year period. For large rain‐on‐snow events (event rainfall >40 mm), snowmelt enhanced the production of water available for run‐off (WAR) by approximately 25% over rainfall alone. For smaller events, a range of antecedent and meteorological factors influenced WAR generation, particularly the antecedent liquid water content of the snowpack. Most large events were associated with atmospheric rivers. Rainfall dominated WAR generation during autumn and winter events, whereas snowmelt dominated during spring and summer events. In the majority of events, the sensible heat of rain contributed less than 10% of the total energy consumed by snowmelt. This analysis illustrated the importance of understanding the amount of rainfall occurring at high elevations during rain‐on‐snow events in mountainous regions.     

Lightning activity and precipitation structure of hailstorms

By using the cloud-to-ground (CG) lightning location data from the lightning detection network of He- nan Province, surface Doppler radar data and standard orbit data of PR, TMI and LIS on TRMM satellite, the spatial and temporal characteristics of CG lightning flashes in 10 severe hailstorms are analyzed. The results show that the percentage of CG lightning in these hailstorms is high with an average value of 45.5%. There is a distinct increase in CG flash rate during the rapid development stage of hailstorms. The hailstone falling corresponds to an active positive flash period, and the increase of CG flash rate is generally accompanied with a decrease of –CG flash rate. The flash rate declines rapidly during the dissipating stage of hailstorms. The precipitation structure and lightning activity in two typical hail- storms are studied in detail. It is found that strong convective cells with reflectivity greater than 30dBZ mainly are situated in the front region of hailstorms, whereas the trailing stratiform region is in the rear part of the hailstorms. The maximum heights of echo top are higher than 14 km. Convective rain con- tributes much more rainfall to the total than stratiform rain, and the convective rain takes about 85% and 97% of the total in the two cases, respectively. Total lightning in the hailstorms is very active with the flash rate up to 183 fl/min and 55 fl/min, respectively. The results also indicate that most lightning flashes occurred in the echo region greater than 30 dBZ and its immediate periphery. The probability of lightning occurrence is 20 times higher in the convective region than in the stratiform region. The result suggests that the lightning information is helpful to the identification of convective rain region. The linear relationship between flash rate and ice water content is disclosed primarily.