Enhanced identification of a hydrologic model using streamflow and satellite water storage data: A multicriteria sensitivity analysis and optimization approach

Hydrologic model development and calibration have continued in most cases to focus only on accurately reproducing streamflows. However, complex models, for example, the so‐called physically based models, possess large degrees of freedom that, if not constrained properly, may lead to poor model performance when used for prediction. We argue that constraining a model to represent streamflow, which is an integrated resultant of many factors across the watershed, is necessary but by no means sufficient to develop a high‐fidelity model. To address this problem, we develop a framework to utilize the Gravity Recovery and Climate Experiment's (GRACE) total water storage anomaly data as a supplement to streamflows for model calibration, in a multiobjective setting. The VARS method (Variogram Analysis of Response Surfaces) for global sensitivity analysis is used to understand the model behaviour with respect to streamflow and GRACE data, and the BORG multiobjective optimization method is applied for model calibration. Two subbasins of the Saskatchewan River Basin in Western Canada are used as a case study. Results show that the developed framework is superior to the conventional approach of calibration only to streamflows, even when multiple streamflow‐based error functions are simultaneously minimized. It is shown that a range of (possibly false) system trajectories in state variable space can lead to similar (acceptable) model responses. This observation has significant implications for land‐surface and hydrologic model development and, if not addressed properly, may undermine the credibility of the model in prediction. The framework effectively constrains the model behaviour (by constraining posterior parameter space) and results in more credible representation of hydrology across the watershed.     

Geochemical evolution of groundwater in hard-rock aquifers of South India using statistical and modelling techniques

ABSTRACT

Multivariate statistical analysis and inverse geochemical modelling techniques were employed to deduce the mechanism of groundwater evolution in the hard-rock terrain of Telangana, South India. Q-mode hierarchical cluster analysis (HCA) and principal component analysis (PCA) were used to extract the hydrogeochemical characteristics and classify the groundwater samples into three principal groups. Use of thermodynamic stability diagrams and inverse geochemical modelling in PHREEQC identified the chemical reactions controlling hydrogeochemistry of each of the groups obtained from statistical analysis. The model output showed that a few phases are governing the water chemistry in this area and the geochemical reactions responsible for evolution of groundwater chemistry along the flow path are (i) dissolution of evaporite minerals (dolomite, halite); (ii) dissolution of primary silicate minerals (albite, anorthite, K-feldspar, biotite); (iii) precipitation of secondary silicate minerals (kaolinite, quartz, gibbsite, Ca-montmorillonite) along with anhydrite and calcite; and (iv) reverse ion exchange processes.     

中国宏观经济的国际评价

本文以世界平均水平为参照系，对中国的宏观经济特征进行概括和评价，中国在经济规模上是发展很快的大国，在经济水平上是不断提高的究国，在经济结构上具有相对完整的国民经济体系和工业经济体系，中国宏观经济的主要问题是第三产业和城镇化严重滞后，经济效益相对低下，行为开支比例太高，而科技教育投入明显不足。     

P–T–t–d evolution of orogenic middle crust of the Roc de Frausa Massif (Eastern Pyrenees): a result of horizontal crustal flow and Carboniferous doming?

Structural, petrological and textural studies are combined with phase equilibria modelling of metapelites from different structural levels of the Roc de Frausa Massif in the Eastern Pyrenees. The pre‐Variscan lithological succession is divided into the Upper, Intermediate and Lower series by two orthogneiss sheets and intruded by Variscan igneous rocks. Structural analysis reveals two phases of Variscan deformation. D1 is marked by tight to isoclinal small‐scale folds and an associated flat‐lying foliation (S1) that affects the whole crustal section. D2 structures are characterized by tight upright folds facing to the NW with steep NE–SW axial planes. D2 heterogeneously reworks the D1 fabrics, leading to an almost complete transposition into a sub‐vertical foliation (S2) in the high‐grade metamorphic domain. All structures are affected by late open to tight, steeply inclined south‐verging NW–SE folds (F3) compatible with steep greenschist facies dextral shear zones of probable Alpine age. In the micaschists of the Upper series, andalusite and sillimanite grew during the formation of the S1 foliation indicating heating from 580 to 640 °C associated with an increase in pressure. Subsequent static growth of cordierite points to post‐D1 decompression. In the Intermediate series, a sillimanite–biotite–muscovite‐bearing assemblage that is parallel to the S1 fabric is statically overgrown by cordierite and K‐feldspar. This sequence points to ~1 kbar of post‐D1 decompression at 630–650 °C. The Intermediate series is intruded by a gabbro–diorite stock that has an aureole marked by widespread migmatization. In the aureole, the migmatitic S1 foliation is defined by the assemblage biotite–sillimanite–K‐feldspar–garnet. The microstructural relationships and garnet zoning are compatible with the D1 pressure peak at ~7.5 kbar and ~750 °C. Late‐ to post‐S2 cordierite growth implies that F2 folds and the associated S2 axial planar leucosomes developed during nearly isothermal decompression to <5 kbar. The Lower series migmatites form a composite S1–S2 fabric; the garnet‐bearing assemblage suggests peak P–T conditions of >5 kbar at suprasolidus conditions. Almost complete consumption of garnet and late cordierite growth points to post‐D2 equilibration at <4 kbar and <750 °C. The early metamorphic history associated with the S1 fabric is interpreted as a result of horizontal middle crustal flow associated with progressive heating and possible burial. The upright F2 folding and S2 foliation are associated with a pressure decrease coeval with the intrusion of mafic magma at mid‐crustal levels. The D2 tectono‐metamorphic evolution may be explained by a crustal‐scale doming associated with emplacement of mafic magmas into the core of the dome.     

Water balance components estimation under scenarios of land cover change in the Vea catchment,West Africa

ABSTRACT

The need for a detailed investigation of the Vea catchment water balance components cannot be overemphasized due to its accelerated land-cover dynamics and the associated impacts on the hydrological processes. This study assessed the possible consequences of land-use change scenarios (i.e. business as usual, BAU, and afforestation for the year 2025) compared to the 2016 baseline on the Vea catchment’s water balance components using the Soil and Water Assessment Tool (SWAT) model. The data used include daily climate and discharge, soil and land use/land cover maps. The results indicate that the mean annual water yield may increase by 9.1% under the BAU scenario but decrease by 2.7% under the afforestation scenario; actual evapotranspiration would decrease under BAU but increase under afforestation; and groundwater recharge may increase under both scenarios but would be more pronounced under the afforestation scenario. These outcomes highlight the significance of land-cover dynamics in water resource management and planning at the catchment.     

Is there a link between deep‐sea biodiversity and ecosystem function?

Studying the diversity‐ecosystem function relationship in the deep sea is of primary importance in the face of biodiversity loss and for our understanding of how the deep sea functions. Results from the first study of diversity‐ecosystem function relationships in the deep sea (Danovaro et al. 2008; Current Biology, 18, 1–8) are unexpected and show an exponential relationship between deep‐sea nematode diversity and ecosystem function and efficiency, although this relationship appears largely restricted to relatively low diversities [ES(51) <25]. Here, we investigate the relationship between nematode diversity and several independent measures/proxies of ecosystem function (sediment community oxygen consumption, bacterial biomass, bacterial extracellular enzyme activity) and efficiency (ratio of bacterial/nematode carbon to organic C content of the sediment) on the New Zealand continental slope. Nematode diversity at our study sites was relatively high [ES(51) = 30–42], and there was no relationship between species/functional diversity and ecosystem function/efficiency after accounting for the effects of water depth and food availability. Our results are consistent with a breakdown of the exponential diversity‐function relationship at high levels of diversity, which may be due to increased competition or greater functional redundancy. Future studies need to take into account as many environmental factors and as wide a range of diversities as possible to provide further insights into the diversity‐ecosystem function relationship in the largest ecosystem on Earth.     

信息量与相似度约束下的网络地图服务缩略图自动生成算法

缩略图可显著提高图片和视频等图形图像资源的展现效率,有效改善信息检索过程中的用户体验。地图服务是一种耦合空间和尺度信息的图像资源,其设计创作、检索筛选和存储管理均需要缩略图的支撑。设计精巧的缩略图带给人良好、鲜活的第一印象,帮助用户实现高效的交互与探索;粗糙凌乱的缩略图则让用户产生抵触情绪,打消其继续探索地图服务的主动性。本文借鉴视频关键帧的思路并结合地图表达特点,提出了地图服务关键位置和关键尺度的概念,并设计了相应的量化指标和提取算法,重点解决了缩略图视觉效果不佳、地图信息量不足和自动化程度不高等问题。算法基于信息量识别地图中内容丰富的关键位置,然后利用跨比例尺相似度判断发生显著变化的关键尺度,最后自动筛选出若干张代表地图服务内容的缩略图。试验表明本方法提取的缩略图数量适中、内容丰富、代表性强,可高度还原地图服务的内容构成和外观样式。     

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.     

Spatial Breakdown Point of Variogram Estimators

In the context of robust statistics, the breakdown point of an estimator is an important feature of reliability. It measures the highest fraction of contamination in the data that an estimator can support before being destroyed. In geostatistics, variogram estimators are based on measurements taken at various spatial locations. The classical notion of breakdown point needs to be extended to a spatial one, depending on the construction of most unfavorable configurations of perturbation. Explicit upper and lower bounds are available for the spatial breakdown point in the regular unidimensional case. The difficulties arising in the multidimensional case are presented on an easy example in IR² , as well as some simulations on irregular grids. In order to study the global effects of perturbations on variogram estimators, further simulations are carried out on data located on a regular or irregular bidimensional grid. Results show that if variogram estimation is performed with a 50% classical breakdown point scale estimator, the number of initial data likely to be contaminated before destruction of the estimator is roughly 30% on average. Theoretical results confirm the previous statement on data in IR^d , d 1.     

A rooftop tuned mass damper frame

This article documents the analytical study and feasibility of placing a tuned mass damper in the form of a limber rooftop moment frame atop relatively stiff structures to reduce seismic acceleration response. Six existing structures were analytically studied using a suite of time history and response spectra records. The analyses indicate that adding mass in conjunction with a limber frame results in an increase in the fundamental period of each structure. The fundamental period increase generally results in a decrease in seismic acceleration response for the same time history and response spectra records. Owing to the limber nature of the rooftop frames, non‐linear analysis methods were required to evaluate the stability of the rooftop tuned mass damper frame. The results indicate the addition of a rooftop tuned mass damper frame reduces the seismic acceleration response for most cases although acceleration response can increase if the rooftop frame is not tuned to accommodate the specific structure's dynamic behaviour and localized soil conditions. Appropriate design of the rooftop tuned mass damper frame can result in decreased seismic acceleration response. This translates to safer structures if used as a retrofit measure or a more economical design if used for new construction. Copyright © 2003 John Wiley & Sons, Ltd.