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.     

Aragonite and calcite preservation in sediments from Lake Iznik related to bottom lake oxygenation and water column depth

This study examines the forcing mechanisms driving long‐term carbonate accumulation and preservation in lacustrine sediments in Lake Iznik (north‐western Turkey) since the last glacial. Currently, carbonates precipitate during summer from the alkaline water column, and the sediments preserve aragonite and calcite. Based on X‐ray diffraction data, carbonate accumulation has changed significantly and striking reversals in the abundance of the two carbonate polymorphs have occurred on a decadal time scale, during the last 31 ka cal bp . Different lines of evidence, such as grain size, organic matter and redox sensitive elements, indicate that reversals in carbonate polymorph abundance arise due to physical changes in the lacustrine setting, for example, water column depth and lake mixing. The aragonite concentrations are remarkably sensitive to climate, and exhibit millennial‐scale oscillations. Extending observations from modern lakes, the Iznik record shows that the aerobic decomposition of organic matter and sulphate reduction are also substantial factors in carbonate preservation over long time periods. Lower lake levels favour aragonite precipitation from supersaturated waters. Prolonged periods of stratification and, consequently, enhanced sulphate reduction favour aragonite preservation. In contrast, prolonged or repeated exposure of the sediment–water interface to oxygen results in in situ aerobic organic matter decomposition, eventually leading to carbonate dissolution. Notably, the Iznik sediment profile raises the hypothesis that different states of lacustrine mixing lead to selective preservation of different carbonate polymorphs. Thus, a change in the entire lake water chemistry is not strictly necessary to favour the preservation of one polymorph over another. Therefore, this investigation is a novel contribution to the carbon cycle in lacustrine systems.     

一次湖北暴雪天气的诊断与模拟

利用NCEP GFS资料分析了2007年1月15—16日鄂东南地区降雪过程,对造成暴雪过程的天气系统发生、发展背景场进行分析。并利用中尺度数值天气模式WRF模拟了这次暴雪过程,探讨了其发生发展的机制。天气系统的背景分析表明,这次暴雪过程主要是受700 hPa西南急流和地面冷空气的共同影响而产生的,降水过程与西南急流的变化密切联系。WRF模式较好地再现了此次暴雪的过程。模拟结果表明西南急流的减弱和移出,对应着降雪的开始和停止;在西南急流的左侧,由于低层涡度的增加,使低空辐合、高空辐散,在连续性原理和动力机制约束下导致上升运动的加强是该次暴雪的形成机制。模式结果说明,产生暴雪的上升运动要远小于产生暴雨的上升运动,且在暴雪过程中,中层为上升运动,近地层和高层伴随着下沉运动。     

A study on feasibility of SI for identification of earthquake damages in Taiwan

Using the actual damage data of the strong earthquakes in Taiwan, the effectiveness of the earthquake risk indices, namely, peak ground acceleration (PGA) and spectral intensity (SI), is verified. PGA and Housner's [Housner GW. Spectrum intensity of strong-motion earthquakes. In: Proceedings of symposium on earthquakes and blast effects on structures. EERI, UCLA; 1952] definition of SI are directly compared. A three-parameter spectral intensity system with spectral intensities SI_a, SI_v, and SI_d in the acceleration, velocity, and displacement regions, respectively, is discussed. Here, the effectiveness of SI, SI_a, SI_v, and SI_d has been compared by using the available earthquake-damage data in Taiwan. Three period ranges, namely, 0.1–0.6, 0.6–1.6, and 1.6–3.0 s were used for structures of 1–6 stories, 7–20 stories, and 21 and more stories, respectively. The results indicate that the three-parameter system is a good risk index of the damage potential of earthquakes.     

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.     

Solution for a plane strain rough‐walled hydraulic fracture driven by turbulent fluid through impermeable rock

The impact of turbulent flow on plane strain fluid‐driven crack propagation is an important but still poorly understood consideration in hydraulic fracture modeling. The changes that hydraulic fracturing has experienced over the past decade, especially in the area of fracturing fluids, have played a major role in the transition of the typical fluid regime from laminar to turbulent flow. Motivated by the increasing preponderance of high‐rate, water‐driven hydraulic fractures with high Reynolds number, we present a semianalytical solution for the propagation of a plane strain hydraulic fracture driven by a turbulent fluid in an impermeable formation. The formulation uses a power law relationship between the Darcy‐Weisbach friction factor and the scale of the fracture roughness, where one specific manifestation of this generalized friction factor is the classical Gauckler‐Manning‐Strickler approximation for turbulent flow in a rough‐walled channel. Conservation of mass, elasticity, and crack propagation are also solved simultaneously. We obtain a semianalytical solution using an orthogonal polynomial series. An approximate closed‐form solution is enabled by a choice of orthogonal polynomials embedding the near‐tip asymptotic behavior and thus giving very rapid convergence; a precise solution is obtained with 2 terms of the series. By comparison with numerical simulations, we show that the transition region between the laminar and turbulent regimes can be relatively small so that full solutions can often be well approximated by either a fully laminar or fully turbulent solution.     

祁连山中东段地区中小地震震源机制和区域应力场

利用哥伦比亚大学 GCMT 目录给出的祁连山中东段地区中强地震震源机制资料,研究较大区域（34°-41°N,100°-106°E）的应力场;利用该地区布设的中法微震数字监测台网多年监测资料和甘肃数字监测台网资料,使用 P 波和 S 波初动及振幅比联合反演方法,反演中小地震震源机制解和发震应力场。结果表明,地区构造应力大致为北东40°-45°水平向压应力;景泰地区主压应力方向约北东45°,绝大多数地震为走滑型。天祝-古浪地区有相当部分的逆断层地震分布,主压应力方向约60°,P 轴仰角在10°左右优势分布,大致为水平应力场。这与大区域构造应力场和断层实际分布基本一致。     

动态城市交通网络优化模型研究及算法设计

在城市交通网络中,行车路线所耗时间是用户关心的关键因素之一.本文构建了基于时间因素的城市交通网络模型,为反映交通网络随时间动态变化的特性,网络模型中弧的权重是和时间统计变量相关.本文设计了交通网络模型算法,并对算法进行了详细的论述.针对实际交通网络的行驶线路选择人群中乐观型、保守型和折中型,给出了路线选择的简化模型.文...     

双向近场地震作用下HDR隔震梁桥地震响应研究

以HDR隔震梁桥多自由度（MDOF）模型和等效双线性单自由度（SDOF）模型为研究对象，以典型近场地震动作为输入，研究HDR支座双向耦合效应对HDR隔震梁桥地震响应的影响。研究结果表明:不考虑双向耦合效应的HDR支座滞回曲线呈典型双线性；考虑双向耦合效应的HDR支座滞回曲线面积小于不考虑双向耦合效应的HDR支座滞回曲线面积。不考虑双向耦合效应的顺桥向HDR支座位移峰值d_b大于考虑双向耦合效应时，但横桥向的结果相反。近场地震作用下，对梁桥进行HDR支座隔震设计时，忽略双向耦合效应计算得到的墩底剪力峰值和弯矩峰值均偏于保守。可忽略HDR支座双向耦合效应对HDR隔震梁桥近场地震能量的影响。     

Degradation and Detoxification of Chlorophenols in Continuous‐Flow Fixed‐Bed Aerobic Reactors

An indigenous bacterial strain of Delftia sp. capable of degrading 2,4‐dicholorophenol and an indigenous bacterial community that degrades 2,4,6‐trichlorophenol (TCP) were employed to inoculate continuous down‐flow fixed‐bed reactors. Continuous‐reactors were constructed from PVC employing hollow PVC cylinders as support material. Synthetic wastewater was prepared by dissolving the corresponding chlorophenol in non‐sterile groundwater. Biodegradation was evaluated by spectrophotometry, chloride release, GC, and microbial growth. Detoxification was evaluated by using Daphnia magna as test organism. Delftia sp. was able to remove an average of 95.6% of DCP. Efficiency in terms of chemical oxygen demand (COD) was of 88.9%. The indigenous bacterial community that degrades TCP reached an average efficiency of 96.5 and 91.6% in terms of compound and COD removal, respectively. In both cases stoichiometric removal of chloride and detoxification was achieved. When synthetic wastewater feed was cut off for 7 days, both reactors showed a fast recovery after inflow restarting, reaching average outlet concentration values within 36 h. The promising behavior of the microorganisms and the low cost of the reactors tested allow us to suggest their possible application to remediation processes.