Rain on snow induced urban floods in cold maritime climate: Risk,indicators and trends

Urban floods pose a societal and economical risk. This study evaluated the risk and hydro-meteorological conditions that cause pluvial flooding in coastal cities in a cold climate. Twenty years of insurance claims data and up to 97 years of meteorological data were analysed for Reykjavík, Iceland (64.15°N; <100 m above sea level). One third of the city's wastewater collection system is combined, and pipe grades vary from 0.5% to 10%. Results highlight semi-intensive rain (<7 mm/h; ≤3 year return period) in conjunction with snow and frozen ground as the main cause for urban flood risk in a climate which undergoes frequent snow and frost cycles (avg. 13 and 19 per season, respectively). Floods in winter were more common, more severe and affected a greater number of neighbourhoods than during summer. High runoff volumes together with debris remobilized with high winds challenged the capacity of wastewater systems regardless of their age or type (combined vs. separate). The two key determinants for the number of insurance claims were antecedent frost depth and total precipitation volume per event. Two pluvial regimes were particularly problematic: long duration (13–25 h), late peaking rain on snow (RoS), where snowmelt enhanced the runoff intensity, elongated and connected independent rainfall into a singular, more voluminous (20–76 mm) event; shorter duration (7–9 h), more intensive precipitation that evolved from snow to rain. Closely timed RoS and cooling were believed to trigger frost formation. A positive trend was detected in the average seasonal snow depth and volume of rain and snowmelt during RoS events. More emphasis, therefore, needs to be placed on designing and operating urban drainage infrastructure with regard to RoS co-acting with frozen ground. Furthermore, more detailed, routine monitoring of snow and soil conditions is important to predict RoS flood events.     

Application of the numerical manifold method for stress wave propagation across rock masses

In this paper, the numerical manifold method (NMM) is extended to study wave propagation across rock masses. First, improvements to the system equations, contact treatment, and boundary conditions of the NMM are performed, where new system equations are derived based on the Newmark assumption of the space–time relationship, the edge‐to‐edge contact treatment is further developed for the NMM to handle stress wave propagation across discontinuities, and the viscous non‐reflection boundary condition is derived based on the energy minimisation principle. After the modification, numerical comparisons between the original and improved NMM are presented. The results show that the original system equations result in artificial numerical damping, which can be overcome by the Newmark system equations. Meanwhile, the original contact scheme suffers some calculation problems when modelling stress wave propagation across a discontinuity, which can be solved by the proposed edge‐to‐edge contact scheme. Subsequently, the influence of the mesh size and time step on the improved NMM for stress wave propagation is studied. Finally, 2D wave propagation is modelled, and the model's results are in good agreement with the analytical solution. Copyright © 2013 John Wiley & Sons, Ltd.     

Cracking risk of partially saturated porous media—Part I: Microporoelasticity model

Drying of deformable porous media results in their shrinkage, and it may cause cracking provided that shrinkage deformations are hindered by kinematic constraints. This is the motivation to develop a thermodynamics‐based microporoelasticity model for the assessment of cracking risk in partially saturated porous geomaterials. The study refers to 3D representative volume elements of porous media, including a two‐scale double‐porosity material with a pore network comprising (at the mesoscale) 3D mesocracks in the form of oblate spheroids, and (at the microscale) spherical micropores of different sizes. Surface tensions prevailing in all interfaces between solid, liquid, and gaseous matters are taken into account. To establish a thermodynamics‐based crack propagation criterion for a two‐scale double‐porosity material, the potential energy of the solid is derived, accounting—in particular—for mesocrack geometry changes (main original contribution) and for effective micropore pressures, which depend (due to surface tensions) on the pore radius. Differentiating the potential energy with respect to crack density parameter yields the thermodynamical driving force for crack propagation, which is shown to be governed by an effective macrostrain. It is found that drying‐related stresses in partially saturated mesocracks reduce the cracking risk. The drying‐related effective underpressures in spherical micropores, in turn, result in a tensile eigenstress of the matrix in which the mesocracks are embedded. This way, micropores increase the mesocracking risk. Model application to the assessment of cracking risk during drying of argillite is the topic of the companion paper (Part II). Copyright © 2009 John Wiley & Sons, Ltd.     

Monitoring over‐winter soil water dynamics in a freezing and snow‐covered environment using a thermally insulated tensiometer

Water potential below a frozen soil layer was continuously monitored over an entire winter period (using thermally insulated tensiometers sheltered in a heated chamber) along with other soil, snow and atmospheric variables. In early winter, the freezing front advanced under a thin snow cover, inducing upward soil water flow in the underlying unfrozen soil. The freezing front started to retreat when the snow cover became thick enough to insulate the soil, resulting in the reversal of the flow direction in the unfrozen zone. These data provide a clear illustration of soil water dynamics, which have rarely been monitored with a tensiometer. Copyright © 2005 John Wiley & Sons, Ltd.     

Investigation of directional hydraulic fracturing based on true tri-axial experiment and finite element modeling

The initiation and propagation of directional hydraulic fracturing (DHF) was investigated based on true tri-axial experiment and finite element modeling. The influences of notch angle, notch length and injection rate on the DHF were investigated. The initiation and propagation of DHF was modeled by a 3D nonlinear finite element method. A comparison between experimental investigation and numerical modeling results indicates that there is a good correlation between unbalanced force (UF) and fracturing. UF can be used to predict the hydraulic fracture initiation and propagation.     

Anomalous Western Pacific Subtropical High during Late Summer in Weak La Nia Years: Contrast between 1981 and 2013

Both 1981 and 2013 were weak La Niña years with a similar sea surface temperature (SST) anomaly in the tropical Pacific, yet the western Pacific subtropical high (WPSH) during August exhibited an opposite anomaly in the two years. A comparison indicates that, in the absence of a strong SST anomaly in the tropics, the cold advection from Eurasian high latitudes and the convection of the western Pacific warm pool play important roles in influencing the strength and position of the WPSH in August. In August 1981, the spatial pattern of 500 hPa geopotential height was characterized by a meridional circulation with a strong ridge in the Ural Mountains and a deep trough in Siberia, which provided favorable conditions for cold air invading into the lower latitudes. Accordingly, the geopotential height to the north of the WPSH was reduced by the cold advection anomaly from high latitudes, resulting in an eastward retreat of the WPSH. Moreover, an anomalous cyclonic circulation in the subtropical western Pacific, excited by enhanced warm pool convection, also contributed to the eastward retreat of the WPSH. By contrast, the influence from high latitudes was relatively weak in August 2013 due to a zonal circulation pattern over Eurasia, and the anomalous anticyclonic circulation induced by suppressed warm pool convection also facilitated the westward extension of the WPSH. Therefore, the combined effects of the high latitude and tropical circulations may contribute a persistent anomaly of the WPSH in late summer, despite the tropical SST anomaly being weak.     

5—6月东北冷涡持续活动及强度特征

利用1960—2012年5—6月NCEP/NCAR逐日再分析资料,基于冷涡经典定义,采取客观识别方法检索东北冷涡活动过程,根据东北冷涡活动时空变化特征给出东北冷涡持续活动过程标准, 通过冷涡强度指数进行定量化分析,该指数对冷涡持续活动过程具有较好表征意义。冷涡活动强对应5月乌拉尔山阻塞高压、贝加尔湖阻塞高压和6月鄂霍次克海阻塞高压活动频繁。通过强弱指数年合成,得到6月强指数年冷涡系统较深厚,集中于对流层中高层,冷心结构明显,具有一定大气斜压特征; 高层存在冷中心,低层有冷空气活动,中高层西风带呈明显的上游分流和下游汇合特征,分汇流之间呈东北高、西南低的偶极子阻塞形势;弱指数年冷涡系统较浅薄,主要集中在对流层中低层,冷心结构不明显,不存在阻塞形势。     

长江中下游地区暖区暴雨特征分析

利用2007到2013年5-9月间常规和非常规资料以及6 h一次的NCEP 1°×1°再分析资料,将长江中下游地区暖区暴雨按天气形势划分为冷锋前暖区暴雨、暖切变暖区暴雨以及副热带高压边缘暖区暴雨三种类型。统计表明暖区暴雨一般发生在距离切变线(锋线)100~300 km的暖区内。主要结论包括:(1)冷锋型降水强度偏弱且分布均匀,集中在5、6月;暖切变型发生次数最多且强度大,主要发生在6、7月长江中下游地区的偏南部;副热带高压边缘型发生次数最少但强度较大,发生在7、8月。暖区暴雨的发生次数及强度在大别山、皖南山区较为集中。(2)暖区暴雨中短时强降水贡献大。(3)冷锋背景下的暖区暴雨一般产生在锋前低压槽中,暴雨落区与高低空急流耦合有紧密联系;暖切变型以低层暖切变线为主要天气背景,地面常有弱静止锋,暖区对流活动与中尺度急流结构、地形强迫等因素存在较高的相关性;副热带高压边缘暖区暴雨与局地的水汽积累和对流不稳定条件的发展有密切关系。据此建立三类暖区暴雨的概念模型。     

中亚低涡背景下南疆西部两次强冰雹环境场对比分析

利用常规资料、14:00加密探空、NCEP1°×1°再分析资料,对比分析2013年6月18日和2014年6月23日南疆西部两次强冰雹环境场及物理机制,表明两次冰雹环境场有相同之处:在环流经向度较大的中亚低涡背景下,前期有明显降水,傍晚前后由中亚低涡后部西北气流下产生,高空干冷平流、低层暖湿环境、较大的垂直温度递减率及低层辐合线或切变线为冰雹天气提供了强不稳定层结和动力触发条件,为冷平流强迫类型。但两次冰雹中亚低涡位置、强度、物理机制变化等有所不同:"6·18"中亚低涡压至南疆西部,低涡强,西北风大,而"6·23"中亚低涡在巴尔喀什湖附近,位置偏北,低涡较弱,西北风较小;"6·18"偏西北气流、风(垂直)切变、强回波伸展高度、层结不稳定、水汽及动力等变化均较"6·23"明显偏强;"6·18"的0℃层和-20℃层的高度均比"6·23"的低300~400 m;"6·18"低层南疆盆地偏东风日变化明显,即傍晚到夜间增强,山前东西风辐合促使上升运动发展,并有明显的中尺度垂直环流圈,而"6·23"低层为西北风风速辐合及与偏东风的切变。     

京津冀地区中α尺度飑线过程中大风特征分析及成因初探

应用常规天气观测资料、地面加密自动气象站资料、大风灾情报资料、京津冀地区7部多普勒天气雷达组网观测资料及VDRAS资料,从多个角度对2013年8月4日京津冀地区一次飑线过程产生的大范围大风天气过程进行了分析,结果显示:此次过程是在高空冷空气南下、低层暖湿气流北上、系统前倾及位势不稳定的有利层结条件下,由多单体风暴演变为中α尺度的强飑线所致。飑线形成于低层垂直切变加强、冷池合并之后;大风主要发生在飑线主体回波中,其次是主体回波前和中前,主体回波后很少发生。大风发生的位置取决于飑线结构中气流的性质,气流的性质与冷池前进的程度和对流的强度关系密切。大风大部分由下沉冷气流产生,少数为近地面上升暖气流导致。大风发生的范围和强度与低层风垂直切变的强度呈正比,大范围低层风垂直切变的加强增强了飑线入流和出流的强度,是大范围大风、局部强风形成的重要原因。大风发生站次与冷池的强度和范围密切相关,冷池的加强和范围的扩大加强了后侧冷入流和前侧暖入流的强度和范围,也是大范围大风形成的重要因素。