Analysis of runoff in ungauged mountain watersheds in Sichuan, China using kinematic-wave-based GIUH model

Floods are one of the most common natural hazards occurring all around the world. However, the knowledge of the origins of a food and its possible magnitude in a given region remains unclear yet. This lack of understanding is particularly acute in mountainous regions with large degrees in Sichuan Province, China, where runoff is seldom measured. The nature of streamflow in a region is related to the time and spatial distribution of rainfall quantity and watershed geomorphology. The geomorphologic characteristics are the channel network and surrounding landscape which transform the rainfall input into an output hydrograph at the outlet of the watershed. With the given geomorphologic properties of the watershed, theoretically the hydrological response function can be determined hydraulically without using any recorded data of past rainfall or runoff events. In this study, a kinematic-wave-based geomorphologic instantaneous unit hydrograph (KW-GIUH) model was adopted and verified to estimate runoff in ungauged areas. Two mountain watersheds, the Yingjing River watershed and Tianquan River watershed in Sichuan were selected as study sites. The geomorphologic factors of the two watersheds were obtained by using a digital elevation model (DEM) based on the topographic database obtained from the Shuttle Radar Topography Mission of US’s NASA. The tests of the model on the two watersheds were performed both at gauged and ungauged sites. Comparison between the simulated and observed hydrographs for a number of rainstorms at the gauged sites indicated the potential of the KW-GIUH model as a useful tool for runoff analysis in these regions. Moreover, to simulate possible concentrated rainstorms that could result in serious flooding in these areas, synthetic rainfall hyetographs were adopted as input to the KW-GIUH model to obtain the flow hydrographs at two ungauged sites for different return period conditions. Hydroeconomic analysis can be performed in the future to select the optimum design return period for determining the flood control work.     

Watershed classification by remote sensing indices: A fuzzy c-means clustering approach

Determining the relatively similar hydrological properties of the watersheds is very crucial in order to readily classify them for management practices such as flood and soil erosion control. This study aimed to identify homogeneous hydrological watersheds using remote sensing data in western Iran. To achieve this goal, remote sensing indices including SAVI, LAI, NDMI, NDVI and snow cover, were extracted from MODIS data over the period 2000 to 2015. Then, a fuzzy method was used to clustering the watersheds based on the extracted indices. A fuzzy c-mean (FCM) algorithm enabled to classify 38 watersheds in three homogeneous groups. The optimal number of clusters was determined through evaluation of partition coefficient, partition entropy function and trial and error. The results indicated three homogeneous regions identified by the fuzzy c-mean clustering and remote sensing product which are consistent with the variations of topography and climate of the study area. Inherently, the grouped watersheds have similar hydrological properties and are likely to need similar management considerations and measures.     

甘东南及邻近地区地震目录的聚类分析

利用2013-01-01~2021-11-01期间甘肃甘东南及邻近地区（32°~36°N,102°~106°E）地震目录中的11 659次地震震中位置数据,使用硬聚类方法中的K-means和软聚类方法高斯混合模型GMM聚类方法对地震原始目录以及地震精定位目录的空间位置进行聚类分析。为确定最佳聚类数,使用AIC和BIC模型选择法,最终将原始目录和精定位目录分别聚类成6个和14个地震群。结果显示,精定位地震目录和GMM聚类方法结合可以更好地从地震大数据中找到具有不同空间分布特征的地震群。     

遥感及再分析降水产品在缺资料干旱内陆盆地的适用性评估

山区降水较集中,但降水测站多位于山谷或人口密集区,代表性差。遥感和再分析降水产品能提供时空分布连续的数据,不受地形条件限制。柴达木盆地中心属干旱荒漠区,水是制约该区开发的首要条件,其四周属高寒山区,降水相对较多,但降水监测十分薄弱。为获取该区相对精确的降水时空分布信息,本文评估了4套高分辨率降水产品（CMADS、TRMM、GPM和MSWEP）的适用性。首先基于地面站点数据评估它们在不同时空尺度上的精度,并分析它们在柴达木盆地的空间分布和年内分配特征。然后,以盆地东南隅的无测站山区香日德河流域为研究区,利用降水产品驱动SWAT模型来评估它们的分布式水文模拟适用性。结果表明：① MSWEP在年、月尺度上与站点降水的吻合程度最高（R ≥ 0.79,PBIAS = 0.5%）,其次是GPM和TRMM,CMADS精度最低（R ≥ 0.64,PBIAS = 5.8%）;② 从降水精度与站点高程的关系来看,降水产品在相对低海拔区容易高估站点降水,而在相对高海拔区常低估实际降水;③ 在香日德河流域,MSWEP（NSE = 0.64）在基准期（2009—2012年）的径流模拟表现明显好于其它降水产品（NSE = 0.36~0.59）,变化期（2013—2016年）表现最好的是CMADS（NSE = 0.75,其余产品NSE = 0.53~0.68）。本研究可为缺资料干旱山区获取精确的降水时空信息和后续水资源的科学管理与规划提供重要支撑。     

Hydrological Simulation Using TRMM and CHIRPS Precipitation Estimates in the Lower Lancang-Mekong River Basin

Satellite-based products with high spatial and temporal resolution provide useful precipitation information for data-sparse or ungauged large-scale watersheds. In the Lower Lancang-Mekong River Basin, rainfall stations are sparse and unevenly distributed, and the transboundary characteristic makes the collection of precipitation data more difficult, which has restricted hydrological processes simulation. In this study, daily precipitation data from four datasets(gauge observations, inverse distance weighted(IDW) data, Tropical Rainfall Measuring Mission(TRMM) estimates, and Climate Hazards Group InfraRed Precipitation with Stations(CHIRPS) estimates), were applied to drive the Soil and Water Assessment Tool(SWAT) model, and then their capability for hydrological simulation in the Lower Lancang-Mekong River Basin were examined. TRMM and CHIRPS data showed good performances on precipitation estimation in the Lower Lancang-Mekong River Basin, with the better performance for TRMM product. The Nash-Sutcliffe efficiency(NSE) values of gauge, IDW, TRMM, and CHIRPS simulations during the calibration period were 0.87, 0.86, 0.95, and 0.93 for monthly flow, respectively, and those for daily flow were 0.75, 0.77, 0.86, and 0.84, respectively. TRMM and CHIRPS data were superior to rain gauge and IDW data for driving the hydrological model, and TRMM data produced the best simulation performance. Satellite-based precipitation estimates could be suitable data sources when simulating hydrological processes for large data-poor or ungauged watersheds, especially in international river basins for which precipitation observations are difficult to collect. CHIRPS data provide long precipitation time series from 1981 to near present and thus could be used as an alternative precipitation input for hydrological simulation, especially for the period without TRMM data. For satellite-based precipitation products, the differences in the occurrence frequencies and amounts of precipitation with different intensities would affect simulation results of water balance components, which should be comprehensively considered in water resources estimation and planning.     

一种面向栅格的空间-属性双重约束聚类方法

针对栅格数据,传统聚类方法大都基于专题属性进行聚类,分裂了栅格对象的空间特性与专题属性,而兼顾空间与专题属性的现有空间聚类方法又存在算法复杂、参数设置多等问题,因此本文提出了一种面向栅格的空间-属性双重约束聚类算法（A Raster-oriented Clustering Method with Space-Attribute Constraints, RoCMSAC）。RoCMSAC利用栅格数据空间邻域和空间连通特性,重新定义栅格簇的相似性度量准则,通过属性均质簇生成,空间相邻栅格簇合并和空间邻近栅格簇合并3个步骤对栅格数据进行空间-属性双重约束聚类。利用太平洋海域海表温度栅格数据对算法的可行性以及有效性进行验证,并与现有算法进行对比分析。通过实例验证与对比发现：① RoCMSAC方法能够保证栅格簇空间域的邻近性和属性域的均质性;② RoCMSAC方法可发现复杂形状的栅格簇,且算法时间复杂度低,需输入参数较少。     

Revealing storage-area relationship of open water in ungauged subalpine wetland — Napahai in Northwest Yunnan, China

Mountains and plateaus in Southwest China contain many subalpine and alpine wetlands, with significant hydroecological functions. But ungauged or poorly gauged conditions limit the study and understanding of hydrological regimes of these wetland types. This study selects an ungauged subalpine wetland-Napahai in Northwest Yunnan, China-as a case for developing a practical approach to revealing its storage-area relationship of open water. A Trimble R8 GNSS (Global Navigation Satellites Systems) RTK (Real-time Kinematic system) and sonar fathometer were used to survey fine-resolution elevation data and generate a digital elevation model of the Napahai Wetland. Forty-four Landsat images from 1987 to 2011 were collected, and the Normalized Difference Water Index was used to classify open water features in the area. The area of open water in Napahai was calculated for each phase. With these data and a developed conceptual model, the storage of open water for each phase was estimated using ArcGIS tools. Both storage and area of open water showed significant intra-annual and inter-annual variations. In the rainy season, the monthly change of average storage of open water in Napahai showed about 1-2 months lag behind mean monthly rainfall. The storage-area relationship of open water was well fit by a power function equation (R 2 ≈0.91, n=44). This study indicates that if detailedelevations are available for similarly ungauged subalpine wetlands in Southwest China, researchers can use this practical approach to estimate multi- temporal areas and storages and reveal the storage-area relationship of open water in the wetlands. The study provided valuable information of this case wetland for optimizing its hydro-ecological managements and a new method to wetland researchers and managers for the hydrological study of similarly ungauged wetland complex.     

浮动车轨迹数据聚类的有向密度方法

为了充分挖掘浮动车轨迹数据的潜在特性,本文在OPTICS空间密度聚类算法基础上,提出了一种有向密度的快速聚类方法(D-OPTICS)。该方法通过扇形空间邻域计算其有向密度信息,并基于方向信息约束其密度可连通性,通过有向可达距离曲线生成数据基本簇,最后,通过空间网格及类簇聚合等优化方法,实现其大规模浮动车轨迹数据的快速聚类处理。通过有向时空数据的聚类分析,发现浮动车轨迹的时空分布特性,以提取复杂路网的结构信息。本文以福州市大规模浮动车轨迹数据,对D-OPTICS进行了系统实验,分析表明,该算法可实现浮动车轨迹数据的快速有向密度聚类分析,有助于挖掘发现时空轨迹数据的分布规律,且基于聚类结果提取了福州市区复杂路网的有向拓扑结构图。同时,与DBSCAN及OPTICS等传统的密度聚类算法进行性能对比,实验表明,D-OPTICS算法能更好地支持大规模浮动车轨迹数据的处理要求。     

Characteristics of clustering debris flows in Wenchuan earthquake zone

Clustering debris-flow events, namely many debris flows simultaneously triggered by a regional rainstorm in a large-scale mountainous area, occurred in four regions of Wenchuan earthquake stricken areas in 2008 and 2010. The characteristics of the clustering debris flows are examined with regard to triggering rainfall, formation process, and relationship with the earthquake by field survey and remote sensing interpretation. It is found that the clustering events occurred nearly at the same time with the local peak rainstorms, and the rainfall intensity-duration bottom limit line for clustering debris flows is higher than the worldwide line. It means that more rainfall is needed for the occurrence of the clustering debris flows. Four kinds of major formation processes for these debris flows are summarized： tributary-dominated, mainstream- dominated, transformation from slope failures, and mobilization or liquefaction of landslide. The four regions has a spatial correlation with the strong- quake-influenced zone with the peak ground acceleration = 0.2 g and the seismic intensity 〉 X.     

A multiscale spatio-temporal framework to regionalize annual precipitation using k-means and self-organizing map technique

Determination of homogenous precipitation-based regions is a very important task in effective management of water resources. The present study tried to propose an effective precipitation-based regionalization methodology by conjugating both temporal pre-processing and spatial clustering approaches in a way to take advantage of multiscale properties of precipitation time series. Annual precipitation data of 51 years (1960-2010) for 31 rain gauges (RGs) were collected and used in proposed clustering approaches. Discreet wavelet transform (DWT) was used to capture the time-frequency attributes of the time series and multiscale regionalization was performed by using k-means and Self Organizing Maps (SOM) clustering techniques. Daubechies function (db) was selected as mother wavelet to decompose the precipitation time series. Also, proper boundary extensions and decomposition level were applied. Different combinations of the approximation (A) and detail (D) coefficients were used to determine the input dataset as a basis of spatial clustering. The proposed model’s efficiency in spatial clustering stage was verified using three different indexes namely, Silhouette Coefficient (SC), Dunn index and Davis Bouldin index (DB). Results approved superior performance of k-means technique in comparison to SOM. It was also deduced that DWT-based regionalization methodology showed improvements in comparison to historical-based models. Cross mutual information was used to investigate the RGs of cluster 3’s homogeneousness in DWT-k-means approach. Results of non-linear correlation approach verified homogeneity of cluster 3. Verifications based on mean annual precipitation values of rain gauges in each cluster also approved the capability of multiscale approach in precipitation regionalization.     

集合卡尔曼滤波在流域水文模型流量预报中的应用

针对流域降雨入渗过程,引入集合卡尔曼滤波(EnKF)理论,视整个边坡流域为一个随机动态系统,将边坡流域流量观测值作为系统的输出,用集合卡尔曼滤波模型来描述系统的状态;结合流域流量计算方法,实现水文模型参数的随机动态估计,在有效获得待估参数的同时还给出估计值的不确定性.通过数值算例表明,集合卡尔曼滤波可以有效地对含噪声的量测数据进行处理,能够跟踪水文模型的动态变化.相对于常用最优化算法,集合卡尔曼滤波同时给出反演结果和先验知识的后验分布,显示出更好的实时性和可靠性.     

EFFECTS OF SEASONALITY ON STREAMFLOW AND WATER QUALITY OF THE PINANG RIVER IN PENANG ISLAND, MALAYSIA

For the Pinang River, originating in the western highlands of Penang Island, the nature, sources and extent of pollution were studied. The river water samples collected at five selected sites were analyzed for various physical and chemical parameters, namely temperature, DO, BOD, COD, SS, pH, ammoniac nitrogen (AN), and conductance. Long-term data of rainfall and temperature were analyzed to determine the seasonal variations of the streamflow.The streamflow during the dry season is extremely low compared to the wet season, thus concentrations of contaminants derived from point pollution source increase due to lack of rainfall and runoff events. On the contrary, in the predominantly urban and agricultural catchments, non-point pollution source increases during rainy season through seepage and runoff. Effects of seasonal variations consequently deterrnine the quantity and quality of the water parameters.The Jelutong River, the Dondang River and the Air Itam River carry the seepage from widely urban and residential areas to the main Pinang River systems. Water quality of the Pinang River at different points assessed by the water quality indices was compared. According to the quality indices during the study period, water quality in the upper reaches of the river is medium to good. It dwindled in the plains, due to the seepage from urban areas and discharges from the industrial and agricultural lands.     

Comparison of the rise of water level in the typical catchments,Three Gorges Reservoir area

Water level is an important index for studying hydrologic processes. Water level rise processes were studied in three catchments (catchment I, II, III in Chen Jiagou watershed in the Three Gorge Reservoir Area) with different areas to provide useful information to inform data extension from a gauged-catchment to an ungauged catchment. The results showed that there are seasonal changes in the dominant driving mode of the rise of the water level. The rise of the water level in March is likely mainly driven by the mode of stored-full runoff, and in September or October, it is mainly driven by Horton-flow. The correlation coefficients of all indexes were significant among the three catchments, suggesting that these catchments have similarities and that water level data extension is likely to be completed successfully between the large catchment (III-Catchment) and the small catchment (I-Catchment). It was confirmed that there is good similarity between the 0.6 km2 and 6 km2 catchments, and the data correlation is good between the catchments with the area differences in the Three Gorges Reservoir Area. In addition, the rise processes of the water level in the catchments were not only different under the same rain conditions, but this difference could also change with the rain condition.     

Probabilistic Properties of a Curve Number:A Case Study for Small Polish and Slovak Carpathian Basins

The proper determination of the curve number(CN) in the SCS-CN method reduces errors in predicting runoff volume. In this paper the variability of CN was studied for 5 Slovak and 5 Polish Carpathian catchments. Empirical curve numbers were applied to the distribution fitting. Next, theoretical characteristics of CN were estimated. For 100-CN the Generalized Extreme Value(GEV) distribution was identified as the best fit in most of the catchments. An assessment of the differences between the characteristics estimated from theoretical distributions and the tabulated values of CN was performed. The comparison between the antecedent runoff conditions(ARC) of Hawkins and Hjelmfelt was also completed. The analysis was done for various magnitudes of rainfall. Confidence intervals(CI) were helpful in this evaluation. The studies revealed discordances between the tabulated and estimated CNs. The tabulated CNs were usually lower than estimated values; therefore, an application of the median value and the probabilistic ARC of Hjelmfelt for wet runoff conditions is advisable. For dry conditions the ARC of Hjelmfelt usually better estimated CN than ARC of Hawkins did, but in several catchments neither the ARC of Hawkins nor Hjelmfelt sufficiently depicted the variability in CN.     

一种基于空间-拓扑结构相似性的复杂轨迹聚类算法

复杂的面状空间实体如海洋涡旋、环流和降雨过程在运动过程中会产生更复杂的轨迹,即具有分支结构的复杂轨迹。为了挖掘这类复杂轨迹的运动模式特征,本文从复杂轨迹的拓扑结构和空间特征出发,创新性地提出复杂轨迹的空间-拓扑结构相似性度量算法（Spatial-Topological Similarity Measurement, STSM）,该算法是基于图同构算法VF2改进的。首先STSM算法将复杂轨迹用带有节点和边的图结构表达,并将空间信息融入图结构的节点属性中,通过匹配复杂轨迹之间所有最大公共子结构,找到匹配结构中节点之间一一对应的关系,利用加权的欧式距离计算复杂轨迹匹配结构中点对之间的空间距离。然后,基于STSM相似性算法进行层次聚类分析,旨在发现复杂轨迹之间相似的拓扑结构在空间上的聚集模式。最后,利用1993-2016年长时间序列的中国南海冷涡复杂轨迹验证方法的有效性,并对比分析复杂轨迹拓扑结构相似性算法CSM。结果表明：单纯用拓扑结构相似性算法CSM进行聚类分析,不能充分挖掘空间的聚集模式,因为不同空间位置也存在拓扑结构相似的轨迹。而本文提出的STSM算法将南海冷涡复杂轨迹分为5类,第一类分布在南海北部、第二类分布在南海中部、其他三类交错在南海南部。这种聚集模式在一定程度上反映了冷涡的生成和演化过程在南海北部、中部、南部的差异性,同时也表明了冷涡移动在南海南部存在更为复杂的异质性。因此,本文提出的方法可以有效地从复杂轨迹数据中发现其演化过程的潜在聚集模式,为认识这类复杂动态现象的时空演化特征提供了一种新的方法。     

一种逐级合并OD流向时空联合聚类算法

现有OD流向聚类多将O点和D点相分离或者将OD流向看作4维空间的数据点进行聚类处理,忽视了流向长度、方向、时间对流向聚类的影响。本文以流向作为研究对象,提出一种基于流向间相似性度的逐级合并OD流向时空联合聚类算法。首先在充分研究OD流向的空间信息和时间信息的基础上,构建合理的OD流向间时空相似性度量方法,对OD流向间的时空相似性进行量化;然后提出逐级合并OD流向聚类策略,优化类簇合并的顺序,以减少层次聚类的时间开销,实现OD流向的时空联合聚类。以成都市的滴滴出行OD数据和纽约市出租车数据为例对本文方法进行了验证,结果表明：① 本算法聚类获得的流向类簇不仅带有空间特征还具备时间特征;② 在不同参数下本方法可以得到不同时空尺度的聚类结果;③ 与现有较高水平的流向聚类算法相对比,本文方法的聚类效果更好。这体现在流向类簇内部的流向之间有着充分的相似性,以及本文方法不仅可以提取出显著的流向类簇,还可以提取出非热点区域之间的流向类簇。本算法顾及空间因素和时间因素,可以通过调整时空相似性度量方法中的时间参数和空间参数以实现不同时空尺度的流向聚类,这使得从不同时空角度研究城市居民出行模式成为可能。本文提出的OD流向时空联合聚类算法从联合时间信息和空间信息的角度获得对运动数据的新见解,有助于合理全面地研究居民的移动模式、区域之间的空间联系、已知出行结构的确定以及出行目的的探索,是后续一系列分析工作的基础。     

Full 2D Hydrodynamic Modelling of Rainfall-induced Flash Floods

Mountain catchments are prone to flash flooding due to heavy rainfall. Enhanced understanding of the generation and evolution processes of flash floods is essential for effective flood risk management. However, traditional distributed hydrological models based on kinematic and diffusion wave approximations ignore certain physical mechanisms of flash floods and thus bear excessive uncertainty. Here a hydrodynamic model is presented for flash floods based on the full two-dimensional shallow water equations incorporating rainfall and infiltration. Laboratory experiments of overland flows were modelled to illustrate the capability of the model. Then the model was applied to resolve two observed flash floods of distinct magnitudes in the Lengkou catchment in Shanxi Province, China. The present model is shown to be able to reproduce the flood flows fairly well compared to the observed data. The spatial distribution of rainfall is shown to be crucial for the modelling of flash floods. Sensitivity analyses of the model parameters reveal that the stage and discharge hydrographs are more sensitive to the Manning roughness and initial water content in the catchment than to the Green-Ampt head. Most notably, as the flash flood augments due to heavier rainfall, the modelling results agree with observed data better, which clearly characterizes the paramount role of rainfall in dictating the floods. From practical perspectives, the proposed model is more appropriate for modelling large flash floods.     

Impact of Meteorological Drought on Streamflow Drought in Jinghe River Basin of China

Under global climate change, drought has become one of the most serious natural hazards, affecting the ecological environ- ment and human life. Drought can be categorized as meteorological, agricultural, hydrological or socio-economic drought. Among the different categories of drought, hydrological drought, especially streamflow drought, has been given more attention by local govern- ments, researchers and the public in recent years. Identifying the occurrence of streamflow drought and issuing early warning can pro- vide timely information for effective water resources management. In this study, streamflow drought is detected by using the Standard- ized Runoff Index, whereas meteorological drought is detected by the Standardized Precipitation Index. Comparative analyses of fre- quency, magnitude, onset and duration are conducted to identify the impact of meteorological drought on streamflow drought. This study focuses on the Jinghe River Basin in Northwest China, mainly providing the following findings. 1） Eleven meteorological droughts and six streamflow droughts were indicated during 1970 and 1990 after pooling using Inter-event time and volume Criterion method. 2） Streamflow drought in the Jinghe River Basin lagged meteorological drought for about 127 days. 3） The frequency of streamflow drought in Jinghe River Basin was less than meteorological drought. However, the average duration of streamflow drought is longer. 4） The magnitude of streamflow drought is greater than meteorological drought. These results not only play an important theo- retical role in understanding relationships between different drought categories, but also have practical implications for streamflow drought mitigation and regional water resources management.     

面向山洪预警的水雨情监测站网布设方法研究

山洪灾害是中国高频发、高死亡率的自然灾害之一。水雨情站网的合理布设及优化,有利于捕获区域暴雨、洪水情势变化的时空异质性,可显著提高中小流域山洪预警的精度,增强山洪灾害防御能力。本文以山洪灾害高发的福建省顺昌县为例,提出了面向山洪预警的水雨情站网布设方法。县内现状雨量和水位站网监测密度分别为37 km²/站和76 km²/站,主要分布在平原主干河流地区,山洪灾害重点防治区内站网布设不足,小流域暴雨山洪监测和预警能力较弱。针对上述问题,综合分析流域降雨时空特征、历史山洪灾害与山洪灾害预警预报需求,对研究区水雨情站网进行了合理性分析和布设研究,建议增设雨量站3座、水位站3座,其中一座水位站同时监测降雨过程,调整后县内雨量站和水位站的监测密度达到34 km²/站和68 km²/站。本文研究对山洪灾害高发区的水雨情站网布设具有参考和指导意义。     

加拿大Hudson海湾地区的水文研究

Hudson海湾的河川径流量在整个加拿大河川径流量中占有很大的比例，所有河流排泄量约30％流人了Hudson和James海湾。研究了该地区河川径流与当地气候变化之间的关系；把河川向James海湾的排泄与曼尼托巴省Churchill地区海水位的变化联系起来，讨论了其对Hudson海湾再循环的影响；根据曼尼托巴省Churchill地区海水位的观测资料重点讨论了该地区海水位在历史时期的变化及未来的预测变化等三个重要水文问题。最后指出了今后在Hudson海湾地区进行水文研究的方向。