The data‐driven approach as an operational real‐time flood forecasting model

Accurate water level forecasts are essential for flood warning. This study adopts a data‐driven approach based on the adaptive network–based fuzzy inference system (ANFIS) to forecast the daily water levels of the Lower Mekong River at Pakse, Lao People's Democratic Republic. ANFIS is a hybrid system combining fuzzy inference system and artificial neural networks. Five ANFIS models were developed to provide water level forecasts from 1 to 5 days ahead, respectively. The results show that although ANFIS forecasts of water levels up to three lead days satisfied the benchmark, four‐ and five‐lead‐day forecasts were only slightly better in performance compared with the currently adopted operational model. This limitation is imposed by the auto‐ and cross‐correlations of the water level time series. Output updating procedures based on the autoregressive (AR) and recursive AR (RAR) models were used to enhance ANFIS model outputs. The RAR model performed better than the AR model. In addition, a partial recursive procedure that reduced the number of recursive steps when applying the AR or the RAR model for multi‐step‐ahead error prediction was superior to the fully recursive procedure. The RAR‐based partial recursive updating procedure significantly improved three‐, four‐ and five‐lead‐day forecasts. Our study further shows that for long lead times, ANFIS model errors are dominated by lag time errors. Although the ANFIS model with the RAR‐based partial recursive updating procedure provided the best results, this method was able to reduce the lag time errors significantly for the falling limbs only. Improvements for the rising limbs were modest. Copyright © 2011 John Wiley & Sons, Ltd.     

The Early Cretaceous Jacuí Group,a newly discovered volcaniclastic–epiclastic accumulation at the top of the Paraná Basin,southern Brazil

The presence of volcaniclastic rocks related to the silicic magmatism within the Serra Geral Formation has been a matter of long-standing debate. In this paper, we present extensive documentation that supports the presence and abundance of these rocks in the Jacuí Group, a newly discovered volcaniclastic and epiclastic accumulation in southern Brazil. The Jacuí Group is composed of two interfingered stratigraphic units, the Volta Alegre and Tupanciretã formations, and it represents the uppermost stratigraphic unit of the Paraná Basin. The Volta Alegre Formation is primarily composed of resedimented volcaniclastic tuffites, the pyroclasts which were sourced from the Santa Maria subgroup of the Palmas-type of the Serra Geral Formation. The Tupanciretã Formation is composed of fluvial and aeolian deposits transported towards the north–northwest. Deposition of the Jacuí Group began in the Early Cretaceous (∼132 Ma) and was coeval with the acidic volcanism of the Santa Maria subgroup. This group was deposited in a probable interior sag basin that represents either the beginning of the extension in the inner part of the continent that subsequently migrated to the east or the far-field impact of extensional processes that preceded the break-up of Gondwana and the opening of the South Atlantic Ocean.     

江苏一次锢囚状MCS和相关中涡旋MCV的观测分析

利用常规地面和高空气象观测资料,结合气象卫星云图和雷达回波,分析了2009年6月14日15—23时(北京时,下同),造成江苏强对流天气的一个中尺度对流系统(MCS)的锢囚状特征的形成过程及其垂直结构。地面中尺度分析表明,雷暴高压东侧在飑前倒槽北端发展的闭合低压环流的东南气流将暖湿空气输送到冷性雷暴高压的北侧形成东南一西北向的暖舌,从而形成锢囚状的结构。长三角探空网资料的垂直结构分析表明,在对流层下部地面到850 hPa为冷性的雷暴高压,在对流层中部700 hPa为冷性的α中尺度涡旋(MCV),而500 hPa已转变为暖性的MCV。静力学关系可以说明MCV仅仅存在于700~500 hPa的原因和MCS下冷上暖的热力结构密切相关。     

Geochemistry of Paraná-Etendeka basalts from Misiones,Argentina: Some new insights into the petrogenesis of high-Ti continental flood basalts

The Early Cretaceous (∼135–131 Ma) Paraná-Etendeka continental flood basalts, preserved in bulk in the Paraná basin of southern Brazil and vicinity, have been divided into low-Ti and high-Ti types that govern the southern and northern halves of the basin, respectively. We have examined a new sample set from the southern margin of the northern high-Ti segment of Paraná basalts in Misiones, northeastern Argentina. These basalts are strongly to moderately enriched in TiO₂ (2–4 wt.%), have relatively high Ti/Y (300–500), low MgO (3.5–6.5 wt.%), and high Fe (FeO(tot) 12–14 wt.%) and belong to the Pitanga and Paranapanema magma types of Peate et al. (1992). Nd and Sr isotope compositions are quite unvarying with ε_Nd (at 133 Ma) values of −4.6 to −3.6 and initial ⁸⁷Sr/⁸⁶Sr of 0.7054–0.7059 and show no variation with fractionation. Compared to high-Ti lavas in the central and northern parts of the Paraná high-Ti basalt segment, the lavas from Misiones are similar to those in the northeastern magin of the basin but less radiogenic in initial Nd isotope composition than those in the central part. This variation probably reflects mixed EM1-EM2 source components in the sublithospheric mantle. A polybaric melt model of a sublithospheric mantle source at the garnet lherzolite-spinel lherzolite transition is compatible with the observed Ti budget of the Pitanga and Paranapanema lavas, regardless of the Nd isotope composition of their purported source.     

Rainfall flood hazard at nuclear power plants in India

Flood hazard evaluation is an important input for Nuclear Power Plants external events safety studies. In the present study, flood hazard at various nuclear sites in India due to rainfall has been evaluated. Hazard estimation is a statistical procedure by which rainfall intensity versus occurrence frequency is estimated from historical records of rainfall data and extrapolated with asymptotic extreme value distribution. Rainfall data needed for flood hazard assessment are daily annual maximum rainfall (24?h data). The observed data points have been fitted using Gumbel, power law and exponential distribution, and return period has been estimated. To study the stationarity of rainfall data, a moving window estimate of the parameters has been performed. The rainfall pattern is stationary in both coastal and inland regions over the period of observation. The coastal regions show intense rainfall and higher variability than inland regions. Based on the plant layout, catchment area and drainage capacity, the prototype fast breeder reactor (PFBR) site is unlikely to be flooded.     

融雪洪水跳跃变异对其参数估计不确定性影响分析

利用Pettitt非参数检验法和Mann-Kendall非参数趋势检验法,分析年最大洪峰流量序列的非一致性,确定序列的变异形式,采用“分解-合成”理论对其进行一致性修正,得到过去、现状两种条件下年最大洪峰流量序列,根据贝叶斯理论对序列一致性修正前后参数不确定性进行估计,并对其预报区间优良性进行评价。研究结果表明：年最大洪峰流量序列变异点发生在1993年,序列整体上升趋势不显著,在1957-1993年子序列呈显著下降趋势,而1994-2006年子序列变化趋势不显著,跳跃变异为序列主要变异形式;给出了实测、还原及还现序列参数后验分布估计值及95%置信区间,将其结合优化适线法进行P-Ⅲ型频率分析,得到修正前后设计频率年最大洪峰流量预报区间估计值;还原、还现序列与实测序列相比,预报区间覆盖率均提高24%,平均带宽分别减少39.59%、23.17%,平均偏移幅度分别减少28.45%、11.39%。通过对非一致性年最大洪峰流量序列还原/还现计算,可减小参数估计不确定性对其计算产生的影响,从而提高预报区间的可靠性。     

金沙江“11·3”白格堰塞湖溃决洪水事件在奔子栏—石鼓段的地貌作用和沉积特征

古堰塞湖溃决洪水事件的重建是当前地学研究的热点问题之一,寻找足够的可参考的现代溃决洪水事件案例是顺利开展这项工作的基础。2018年11月13日发生在金沙江干流的白格堰塞湖超万年一遇的溃决洪水事件(学术界称之为“11·3”白格堰塞湖溃决洪水事件)就是一个难得的样本。这次溃决事件发生在枯水期,洪峰完全由溃决洪水产生,没有叠加其他来源,对评估流域地貌和沉积体系对堰塞溃决事件的响应有很好的参考价值。本文以溃决洪水事件受灾最为严重的奔子栏—石鼓段为研究区,通过详细的野外调查和初步的水力学估算发现“11·3”白格堰塞湖溃决洪水事件在奔子栏—石鼓段的地貌作用主要表现为洪水淹没区的岸坡塌岸和沉积物堆积,未发现明显的基岩侵蚀。沉积物主要由分选良好的具水平纹层的砂组成。受金沙江较低的河床比降影响,洪水产生的基底剪切应力较弱在27~142 N/m²,不能悬浮和搬运直径5 cm以上的砾石,也不能产生明显的磨蚀和冲(撞)击作用。在发生塌岸的部分段落,洪水沉积物中有砾石坠入,甚至会出现类似浊流沉积的层序。这些现象的发现对深入理解堰塞湖溃决洪水的复杂地貌过程和沉积特征有重要参考意义。     

2018年一次罕见早春飑线大风过程演变和机理分析

2018年3月4—5日,华南、江南等地发生了一次大范围强对流过程,发生时间早,落区范围广,多地伴有雷暴大风、冰雹、短时强降水等剧烈对流天气,尤其飑线在江西境内造成了严重大风灾害。基于大气环流和雷达回波发展演变特征,将该次过程分为初始、发展和减弱三个阶段：初始阶段西风槽前西南急流造成的低压倒槽为强对流提供大尺度触发条件;发展阶段对流活动位于槽前暖区中,飑线在江西造成极端大风;入夜后,冷锋南下,对流进入减弱阶段。环境场及对流参数诊断表明江西中北部低层高温高湿,中层干冷,温度垂直递减率大,有利于产生雷暴大风。南昌探空长时间序列分析表明温湿要素气候态异常,与历史同期比,低层明显偏暖偏湿,中层偏干,有利于极端对流天气发生。综合多源观测资料和雷达资料分析中小尺度特征,本次江西飑线过程特点及成因包括：（1）受引导气流和前向传播共同作用,飑线移动速度快。（2）自动站分析显示飑锋后雷暴高压强,与锋前暖低压作用造成强密度流,有利于产生大范围直线型大风;（3）通过对比飑线弓状回波南北段回波结构差异表明,飑线后侧中层干后向入流促使降水粒子相变,剧烈降温形成的强下沉运动(下击暴流)是导致极端大风的主要原因,后部层云区下沉气流增强雷暴高压加之动量下传作用对雷暴大风有增幅作用。     

内蒙古东南部暴雨洪涝灾害风险评估与区划——以通辽市为例

文章选取通辽市各气象站点气象数据、基础地理信息数据、历史灾情数据、第二次土地调查数据及相关社会经济等统计数据,基于自然灾害风险原理,利用数理统计、层次分析法及空间叠加分析等方法,提出了通辽市暴雨洪涝灾害风险评价指标体系,得到通辽市暴雨洪涝灾害风险评估空间分布图。结果表明:高风险区主要分布在通辽市北部山区、中东部地区及偏南部地区,包括扎鲁特旗鲁北镇、科尔沁区、科左中旗东南部及库伦旗东南部;而科左后旗大部地区易形成由短时强降水引发的洪涝,为次高风险区;低风险区及次低风险区主要分布在通辽市中部、西南部地区,包括开鲁县、奈曼旗北部、科左中旗西部。     

FLOOD MODELLING WITH GIS-DERIVED DISTRIBUTED UNIT HYDROGRAPHS

The concept of a spatially distributed unit hydrograph is based on the fact that the unit hydrograph can be derived from the time–area curve of a watershed by the S-curve method. The time–area diagram is a graph of cumulative drainage area contributing to discharge at the watershed outlet within a specified time of travel. Accurate determination of the time–area diagram is made possible by using a GIS. The GIS is used to describe the connectivity of the links in the watershed flow network and to calculate distances and travel times to the watershed outlet for various points within the watershed. Overland flow travel times are calculated by the kinematic wave equation for time to equilibrium; channel flow times are based on the Manning and continuity equations. To account for channel storage, travel times for channel reaches are increased by a percentage depending on the channel reach length and geometry. With GIS capability for rainfall mapping, the assumption of a uniform spatial rainfall distribution is no longer necessary; hence the term, spatially distributed unit hydrograph. An example of the application for the Waiparous Creek in the Alberta Foothills is given. IDRISI is used to develop a simple digital elevation model of the 229 km² watershed, using 1 km × 1 km grid cells. A grid of flow directions is developed and used to create an equivalent channel network. Excess rainfall for each 1 km × 1 km cell is individually computed by the Soil Conservation Service (SCS) runoff curve method and routed through the equivalent channel network to obtain the time–area curve. The derived unit hydrograph gave excellent results in simulating an observed flood hydrograph. The distributed unit hydrograph is no longer a lumped model, since it accounts for internal distribution of rainfall and runoff. It is derived for a watershed without the need for observed rainfall and discharge data, because it is essentially a geomorphoclimatic approach. As such, it allows the derivation of watershed responses (hydrographs) to inputs of various magnitudes, thus eliminating the assumption of proportionality of input and output if needed. The superposition of outputs is retained in simulating flood hydrographs by convolution, since it has been shown that some non-linear systems satisfy the principle of superposition. The distributed unit hydrograph appears to be a very promising rainfall runoff model based on GIS technology.