GIS Spatial-Temporal Modeling of Water Systems in Greater Toronto Area, Canada

Modeling landscape with high-resolution digital elevation model (DEM) in a geographic information system can provide essential morphological and structural information for modeling surface processes such as geomorphologic process and water systems. This paper introduces several DEM-based spatial analysis processes applied to characterize spatial distribution and their interactions of ground and surface water systems in the Greater Toronto Area (GTA), Canada. The stream networks and drainage basin systems were derived from the DEM with 30 m resolution and the regularities of the surface stream and drainage patterns were modeled from a statistical/multifractal point of view. Together with the elevation and slope of topography, other attributes defmed from modeling the stream system, and drainage networks were used to associate geological, hydrological and topographical features to water flow in river systems and the spatial locations of artesian aquifers in the study area. Stream flow data derived from daily flow measurements recorded at river gauging stations for multi-year period were decomposed into “drainage-area dependent“ and “drainage-area independent“ flow components by two-step “frequency“ and “spatial“ analysis processes. The latter component was further demonstrated to relate most likely to the ground water discharge. An independent analysis was conducted to model the distribution of aquifers with information derived from the records of water wells. The focus was given on quantification of the likelihood of ground water discharge to river and ponds through flowing wells, springs and seepages. It has been shown that the Oak Ridges Moraine (ORM) is a unique glacial deposit that serves as a recharge layer and that the aquifers in the ORM underlain by Hilton Tills and later deposits exposed near the steep slope zone of the ridges of ORM provide significant discharge to the surface water systems (river flow and ponds) through flowing wells, springs and seepages. Various statistics (cross- and auto-correlation coefficients, fractal R/S exponent) were used in conjunction with GIS to demonstrate the influence of land types, topography and geometry of drainage basins on short- and long-term persistence of river flows as well as responding time to precipitation events. The current study has provided not only insight in understanding the interaction of water systems in the GTA, but also a base for further establishment of an on-line GIS system for predicting spatial-temporal changes of river flow and groundwater level in the GTA.     

基于物方影像匹配和概率松弛的断面自动提取

提出了基于物方影像匹配的断面提取,为了提高匹配的可靠性,引入了概率松弛法,最终获得了高精度的断面.     

基于通量均衡基准的地壳垂直形变场分析模型

根据全球通量均衡假说建立了地壳垂直形变场的分析模型，在球面近似下求得了该模型的球谐函数级数解并拟合了两个不同区域的地壳垂直形变场，编绘了等值线图，继而在密切平面坐标系下导出了该模型的Fourier级数解。结果表明，该方法不仅具有几何意义，而且具有一定的地球动力学意义。     

Basal sediment evacuation by subglacial meltwater: suspended sediment transport from Haut Glacier d'Arolla,Switzerland

Proglacial suspended sediment transport was monitored at Haut Glacier d'Arolla, Switzerland, during the 1998 melt season to investigate the mechanisms of basal sediment evacuation by subglacial meltwater. Sub‐seasonal changes in relationships between suspended sediment transport and discharge demonstrate that the structure and hydraulics of the subglacial drainage system critically influenced how basal sediment was accessed and entrained. Under hydraulically inefficient subglacial drainage at the start of the melt season, sediment availability was generally high but sediment transport increased relatively slowly with discharge. Later in the melt season, sediment transport increased more rapidly with discharge as subglacial meltwater became confined to a spatially limited network of channels following removal of the seasonal snowpack from the ablation area. Flow capacity is inferred to have increased more rapidly with discharge within subglacial channels because rapid changes in discharge during highly peaked diurnal runoff cycles are likely to have been accommodated largely by changes in flow velocity. Basal sediment availability declined during channelization but increased throughout the remainder of the monitored period, resulting in very efficient basal sediment evacuation over the peak of the melt season. Increased basal sediment availability during the summer appears to have been linked to high diurnal water pressure variation within subglacial channels inferred from the strong increase in flow velocity with discharge. Basal sediment availability therefore appears likely to have been increased by (1) enhanced local ice‐bed separation leading to extra‐channel flow excursions and[sol ]or (2) the deformation of basal sediment towards low‐pressure channels due to a strong diurnally reversing hydraulic gradient between channels and areas of hydraulically less‐efficient drainage. Copyright © 2005 John Wiley & Sons, Ltd.     

华蓥山隧道岩溶发育强度垂直分带

论述了华蓥山隧道的岩溶发育特征,根据岩溶发育强度的划分依据和划分标准,结合岩溶发育的历史,采用类比判别法、主因素判别法和综合判别法对岩溶发育强度作了垂直分带,据此预测的溶洞规模与实际情况吻合。     

Semi‐analytical solution to one‐dimensional consolidation for unsaturated soils with semi‐permeable drainage boundary under time‐dependent loading

This paper presents semi‐analytical solutions to Fredlund and Hasan's one‐dimensional consolidation of unsaturated soils with semi‐permeable drainage boundary under time‐dependent loadings. Two variables are introduced to transform two coupled governing equations of pore‐water and pore‐air pressures into an equivalent set of partial differential equations, which are easily solved by the Laplace transform. The pore‐water pressure, pore‐air pressure and settlement are obtained in the Laplace domain. Crump's method is adopted to perform the inverse Laplace transform in order to obtain semi‐analytical solutions in time domain. It is shown that the present solutions are more general and have a good agreement with the existing solutions from literatures. Furthermore, the current solutions can also be degenerated into conventional solutions to one‐dimensional consolidation of unsaturated soils with homogeneous boundaries. Finally, several numerical examples are provided to illustrate consolidation behavior of unsaturated soils under four types of time‐dependent loadings, including instantaneous loading, ramp loading, exponential loading and sinusoidal loading. Parametric studies are illustrated by variations of pore‐air pressure, pore‐water pressure and settlement at different values of the ratio of air–water permeability coefficient, depth and loading parameters. Copyright © 2017 John Wiley & Sons, Ltd.     

Age of Gimli beach of Lake Agassiz based on new OSL dating

Gimli beach in Manitoba is one of the lowest elevation beaches in the southern Lake Agassiz basin, and is a distinct ridge composed of bedded sand and gravel that rises above the lake plain and extends for more than 40 km. Ten new optically stimulated luminescence (OSL) ages from Gimli beach yield ages mostly ranging from 9.7 ± 0.7 to 10.5 ± 0.8 ka (average 10.3 ± 0.5 ka), which is older by 0.6 to >1.0 ka than age estimates of previous researchers. Two of our new OSL ages are notably older than the others, dating to ~11.3 ± 0.8 and 13.9 ± 1.0 ka, which we attribute to poorly bleached sands. We ascribe an age of about 10 ka to Gimli beach, which is several centuries before overflow from Lake Agassiz and its vast drainage basin shifted from the western Great Lakes to glacial Lake Ojibway and the St. Lawrence Valley.     

Vertical vibration of a rigid strip footing on viscoelastic half-space

This paper presents an analytical method for modeling the dynamic response of a rigid strip footing subjected to vertical-only loads. The footing is assumed to rest on the surface of a viscoelastic half-space; therefore, effects of hysteretic soil damping on the impedance of the foundation and the generated ground vibrations are considered in the solution. To solve the mixed boundary value problem, we use the Fourier transform to cast a pair of dual integral equations providing contact stresses, which are solved by means of Jacobi orthogonal polynomials. The resulting soil and footing displacements and stresses are obtained by means of the Fourier inverse transform. The solution provides more realistic estimates of footing impedance, compared to existing solutions for elastic soil, as well as of the attenuation of ground vibrations with distance of the footing. The latter is important for the estimation of machine vibration effects on nearby structures and installations.     

联合GPS和GRACE数据分析陕甘宁地区地表垂直形变

采用甘肃省CORS网和中国大陆构造环境监测网络中共48个台站的GPS观测数据，解算得到观测台站的垂直位移，并与GRACE时变重力场Mascon模型解CSR RL05M数据计算得到的垂直形变进行比较，分析区域地表垂直形变特征。结果表明，研究区内台站垂直形变存在局部特征，甘肃庆阳和平凉地区垂直形变与其他地区存在明显差异，相关系数、均方根减少量和周年信号减少量均高于其他地区；扣除趋势项后，观测台站GPS垂直位移与GRACE垂直形变时间序列相关系数均值为0.72，GPS和GRACE周年信号振幅均值分别为6.00 mm和3.70 mm，周年信号减少量和均方根误差减少量均值分别为0.51和0.29；研究区内GPS垂直位移和GRACE垂直形变时间序列一致性较强，GRACE垂直形变能有效解释50%以上的GPS垂直位移周年信号，GPS垂直位移时间序列包含的非构造形变中平均约29%来源于环境负载变化所引起的负荷形变。