Soil property variation mapping through data mining of soil category maps

Spatial information on soil properties is an important input to hydrological models. In current hydrological modelling practices, soil property information is often derived from soil category maps by the linking method in which a representative soil property value is linked to each soil polygon. Limited by the area‐class nature of soil category maps, the derived soil property variation is discontinuous and less detailed than high resolution digital terrain or remote sensing data. This research proposed dmSoil, a data‐mining‐based approach to derive continuous and spatially detailed soil property information from soil category maps. First, the soil–environment relationships are extracted through data mining of a soil map. The similarity of the soil at each location to different soil types in the soil map is then estimated using the mined relationships. Prediction of soil property values at each location is made by combining the similarities of the soil at that location to different soil types and the representative soil property values of these soil types. The new approach was applied in the Raffelson Watershed and Pleasant Valley in the Driftless Area of Wisconsin, United States to map soil A horizon texture (in both areas) and depth to soil C horizon (in Pleasant Valley). The property maps from the dmSoil approach capture the spatial gradation and details of soil properties better than those from the linking method. The new approach also shows consistent accuracy improvement at validation points. In addition to the improved performances, the inputs for the dmSoil approach are easy to prepare, and the approach itself is simple to deploy. It provides an effective way to derive better soil property information from soil category maps for hydrological modelling. Copyright © 2014 John Wiley & Sons, Ltd.     

Discharge‐sediment processes of the Zhadang glacier on the Tibetan Plateau measured with a high frequency data acquisition system

In high elevation cold regions of the Tibetan Plateau, suspended sediment transfer from glacier meltwater erosion is one of the important hydrological components. The Zhadang glacier is a typical valley‐type glacier in the Nyainqentanglha Mountains on the Tibetan Plateau. To make frequent and long period records of meltwater runoff and sediment processes in the very high elevation and isolated regions, an automatic system was installed near the glacier snout (5400 m a.s.l) in August 2013, to measure the transient discharge and sediment processes at 5‐min interval, which is shorter than the time span for the water flow to traverse the catchment from the farthest end to the watershed outlet. Diurnal variations of discharge, and suspended sediment concentration (SSC) were recorded at high frequency for the Zhadang glacier, before suspended sediment load (SSL) was computed. Hourly SSC varied from the range of 0.2 kg/m³ to 0.5 kg/m³ (at 8:00–9:00) to the range of 2.0 kg/m³ to 4.0 kg/m³ (at 17:00–18:00). The daily SSL was 32.24 t during the intense ablation period. Hourly SSC was linearly correlated with discharge (r = 0.885**, n = 18, p < 0.01). A digit‐eight hysteresis loop was observed for the sediment transport in the glacier area. Air temperature fluctuations influence discharge, and then result in the sediment variations. The results of this study provide insight into the responses of suspended sediment delivery processes with a high frequency data in the high elevation cold regions. Copyright © 2016 John Wiley & Sons, Ltd.     

A nonstationary extreme value distribution for analysing the cessation of karst spring discharge

The effects of climate change and population growth in recent decades are leading us to consider their combined and potentially extreme consequences, particularly regarding hydrological processes, which can be modeled using a generalized extreme value (GEV) distribution. Most of the GEV models were based on a stationary assumption for hydrological processes, in contrast to the nonstationary reality due to climate change and human activities. In this paper, we present the nonstationary generalized extreme value (NSGEV) distribution and use it to investigate the risk of Niangziguan Springs discharge decreasing to zero. Rather than assuming the location, scale, and shape parameters to be constant as one might do for a stationary GEV distribution analysis, the NSGEV approach can reflect the dynamic processes by defining the GEV parameters as functions of time. Because most of the GEV model is designed to evaluate maxima (e.g. flooding, represented by positive numbers), and spring discharge cessation is a ?minima’, we deduced an NSGEV model for minima by applying opposite numbers, i.e. negative instead of positive numbers. The results of the model application to Niangziguan Springs showed that the probability of zero discharge at Niangziguan Springs will be 1/80 in 2025, and 1/10 in 2030. After 2025, the rate of decrease in spring discharge will accelerate, and the probability that Niangziguan Springs will cease flowing will dramatically increase. The NSGEV model is a robust method for analysing karst spring discharge. Copyright © 2013 John Wiley & Sons, Ltd.     

A new standardized Palmer drought index for hydro‐meteorological use

Accepting the concept of standardization introduced by the standardized precipitation index, similar methodologies have been developed to construct some other standardized drought indices such as the standardized precipitation evapotranspiration index (SPEI). In this study, the authors provided deep insight into the SPEI and recognized potential deficiencies/limitations in relating to the climatic water balance it used. By coupling another well‐known Palmer drought severity index (PDSI), we proposed a new standardized Palmer drought index (SPDI) through a moisture departure probabilistic approach, which allows multi‐scalar calculation for accurate temporal and spatial comparison of the hydro‐meteorological conditions of different locations. Using datasets of monthly precipitation, temperature and soil available water capacity, the moisture deficit/surplus was calculated at multiple temporal scales, and a couple of techniques were adopted to adjust corresponding time series to a generalized extreme value distribution out of several candidates. Results of the historical records (1900–2012) for diverse climates by multiple indices showed that the SPDI was highly consistent and correlated with the SPEI and self‐calibrated PDSI at most analysed time scales. Furthermore, a simple experiment of hypothetical temperature and/or precipitation change scenarios also verified the effectiveness of this newly derived SPDI in response to climate change impacts. Being more robust and preferable in spatial consistency and comparability as well as combining the simplicity of calculation with sufficient accounting of the physical nature of water supply and demand relating to droughts, the SPDI is promising to serve as a competent reference and an alternative for drought assessment and monitoring. Copyright © 2013 John Wiley & Sons, Ltd.     

Watershed delineation using hydrographic features and a DEM in plain river network region

Watershed delineation is a required step when conducting any spatially distributed hydrological modelling. Automated approaches are often proposed to delineate a watershed based on a river network extracted from the digital elevation model (DEM) using the deterministic eight‐neighbour (D8) method. However, a realistic river network cannot be derived from conventional DEM processing methods for a large flat area with a complex network of rivers, lakes, reservoirs, and polders, referred to as a plain river network region (PRNR). In this study, a new approach, which uses both hydrographic features and DEM, has been developed to address the problems of watershed delineation in PRNR. It extracts the river nodes and determines the flow directions of the river network based on a vector‐based hydrographic feature data model. The river network, lakes, reservoirs, and polders are then used to modify the flow directions of grid cells determined by D8 approach. The watershed is eventually delineated into four types of catchments including lakes, reservoirs, polders, and overland catchments based on the flow direction matrix and the location of river nodes. Multiple flow directions of grid cells are represented using a multi‐direction encoding method, and multiple outflows of catchments are also reflected in the topology of catchments. The proposed approach is applied to the western Taihu watershed in China. Comparisons between the results obtained from the D8 approach, the ‘stream burning’ approach, and those from the proposed approach clearly demonstrate an improvement of the new approach over the conventional approaches. This approach will benefit the development of distributed hydrological models in PRNR for the consideration of different types and multiple inlets and outlets of catchments. Copyright © 2015 John Wiley & Sons, Ltd.     

求解GM(1,1)模型新总体最小二乘算法

在GM(1,1)模型中系数矩阵和观测向量都是由原始序列组成的。系数矩阵中同样是有误差的,与观测向量中的误差一样,亦来源于原始序列,即它们误差同源。不同位置的相同元素应该有相同的改正数,采用传统总体最小二乘求解则不能达到此目的。针对这一缺陷,推导了一种新的总体最小二乘算法;并且通过算例验证了新方法的可行性和有效性。     

Si and O self-diffusion in hydrous forsterite and iron-bearing olivine from the perspective of defect chemistry

We discuss the experimental results of silicon and oxygen self-diffusion coefficients in forsterite and iron-bearing olivine from the perspective of defect chemistry. Silicon diffusion is dominated by V_O ^··-associated V_Si″″, whereas oxygen diffusion is dominated by hopping of V_O ^·· under anhydrous conditions, and by (OH)_O ^· under hydrous conditions. By considering the charge neutrality condition of [(OH)_O ^·] = 2[V_Me″] in hydrous forsterite and iron-bearing olivine, we get D _Si ∝ (\(C_{{{\text{H}}_{2} {\text{O}}}}\))^1/3 and D _O ∝ (\(C_{{{\text{H}}_{2} {\text{O}}}}\))⁰, which explains the experimental results of water effects on oxygen and silicon self-diffusion rates (Fei et al. in Nature 498:213–215, 2013; J Geophys Res 119:7598–7606, 2014). The \(C_{{{\text{H}}_{2} {\text{O}}}}\) dependence of creep rate in the Earth’s mantle should be close to that given by Si and O self-diffusion coefficients obtained under water unsaturated conditions.