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.     

基于投影变换的线阵旋转扫描相机标定方法

针对线阵旋转扫描相机标定数据难以获取及模型解算不稳定的问题,提出了一种基于投影变换的线阵旋转扫描相机标定方法。首先构建了一种新的线阵旋转扫描相机成像模型,该模型可以根据相机旋转平台参数、线阵旋转扫描相机成像平面与其切平面的位置关系将旋转扫描图像投影为框幅式图像;然后根据构建出的成像模型,采用直接线性变换法和非线性优化方法进行相机标定和参数优化。试验结果表明,此方法简单易行,具有较高的精度。     

红透山矿床变质岩系褶皱构造特征浅析

红透山式铜-锌矿床形成于太古代克拉通裂谷槽盆中海底火山喷发-沉积环境。因此，矿床的空间展布严格受特定的含矿岩系所控制。槽盆构造中的岩层褶皱构造最发育，矿床赋存在具多期褶皱构造叠加的大型同倾向形褶皱内，褶皱的层间虚脱是控矿有利构造空间，可以借此预测矿体侧伏方向及赋存部位。     

模拟地震的弹簧滑块模型的混沌运动

我们构造了一种自由度为２的双弹簧滑块系统，讨论模型参数的不均匀性和滑块间的相互作用与复杂现象的关系。假定滑块在所受力小于其静摩擦力时保持静止，因此，只有发生了滑动，两滑块所受力之差才改变。如果我们以两滑块受力之差为变量，就可以构造一种一维的映射。这种映射直观地显示了系统的演化，反映了方程终态解的形态。结果表明，用混沌运动来描述地震现象可能是较为合理的。     

Analysis of snowpack dynamics during the spring melt season for a sub‐alpine site using point measurements and numerical modeling

Snowpack dynamics through October 2014–June 2017 were described for a forested, sub‐alpine field site in southeastern Wyoming. Point measurements of wetness and density were combined with numerical modeling and continuous time series of snow depth, snow temperature, and snowpack outflow to identify 5 major classes of distinct snowpack conditions. Class (i) is characterized by no snowpack outflow and variable average snowpack temperature and density. Class (ii) is characterized by short durations of liquid water in the upper snowpack, snowpack outflow values of 0.0008–0.005 cm hr^?1, an increase in snowpack temperature, and average snow density between 0.25–0.35 g cm^?3. Class (iii) is characterized by a partially saturated wetness profile, snowpack outflow values of 0.005–0.25 cm hr^?1, snowpack temperature near 0 °C, and average snow density between 0.25–0.40 g cm^?3. Class (iv) is characterized by strong diurnal snowpack outflow pattern with values as high as 0.75 cm hr^?1, stable snowpack temperature near 0 °C, and stable average snow density between 0.35–0.45 g cm^?3. Class (v) occurs intermittently between Classes (ii)–(iv) and displays low snowpack outflow values between 0.0008–0.04 cm hr^?1, a slight decrease in temperature relative to the preceding class, and similar densities to the preceding class. Numerical modeling of snowpack properties with SNOWPACK using both the Storage Threshold scheme and Richards' equation was used to quantify the effect of snowpack capillarity on predictions of snowpack outflow and other snowpack properties. Results indicate that both simulations are able to predict snow depth, snow temperature, and snow density reasonably well with little difference between the 2 water transport schemes. Richards' equation more accurately simulates the timing of snowpack outflow over the Storage Threshold scheme, especially early in the melt season and at diurnal timescales.     

多尺度非静力通用模式框架的设计策略

针对目前国际上数值预报模式的最新发展趋势,综述了新形势下多尺度非静力通用数值预报模式框架的设计策略。目前世界上开发通用模式主要有两个途经,一是建立一个离散化方案和源程序代码共享的、全球和有限区通用的模式,二是建立一个单一的全球可变分辨率的模式。还从通用模式方程组的选取策略、模式网格属性的构造,时间积分方案、空间离散方案,垂直坐标的选取等方面进行了分析。     

Under‐ice salinity and stable isotope distribution of Saroma‐ko Lagoon,Hokkaido, northern Japan

In winter, lakes and lagoons at high altitudes or high latitudes have interesting hydrological cycles that differ from those in other seasons or in other regions, because water surfaces are covered with ice. Hydrological balances of lakes and lagoons are complex dynamic systems, and to elucidate them, isotopic tracers of water have been used as effective tools along with observations of precipitation, evaporation, inflows, and outflows. Here, to understand hydrological processes during freezing periods in the brackish Saroma‐ko Lagoon, Hokkaido, northern Japan, we examined horizontal and vertical distributions of salinity and isotope compositions of lagoon water and ice in 2005 and 2006. Horizontal and vertical gradients of salinity and isotope compositions were observed from the river mouth to the sea channel, and factors determining these distributions were considered. The mixing of freshwater and seawater and a freezing effect were presumed to be factors in relationships between salinity and isotopes and in relationships between surface waters and ice just above the water. A simple box model for water balance was constructed based on these putative factors to reproduce the distributions of salinity and isotope compositions of surface waters and ice. An evaluation of the model revealed that this hydrological system is controlled primarily by horizontal advection of the epilimnion, freshwater influx, and the ice growth rate. Copyright © 2009 John Wiley & Sons, Ltd.     

Lateral tidal mixing in the Antarctic marginal seas

Tidal mixing plays an important role in the modification of dense water masses around the Antarctic continent. In addition to the vertical (diapycnal) mixing in the near-bottom layers, lateral mixing can also be of relevance in some areas. A numerical tide simulation shows that lateral tidal mixing is not uniformly distributed along the shelf break. In particular, strong mixing occurs all along the Ross Sea and Southern Weddell Sea shelf breaks, while other regions (e.g., the western Weddell Sea) are relatively quiet. The latter regions correspond surprisingly well to areas where indications for cross-shelf exchange of dense water masses have been found. The results suggest that lateral tidal mixing may account for the relatively small contribution of Ross Sea dense water masses to Antarctic Bottom Water.     

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.