西北太平洋浪流相互作用对有效波高的影响研究

西北太平洋强流区会对海浪的特征和分布产生显著的影响,尤其是研究台风过程中海流与海浪的相互作用具有重要的研究意义。本文以ROMS海洋模式和SWAN海浪模式为基础,构建了浪流耦合模式系统,对2013年10月6-17日间的台风“丹娜丝”、“百合”、“韦帕”过程中西北太平洋浪流相互作用中海流对有效波高的影响进行了研究。通过对比模式模拟有效波高与浮标观测资料,发现耦合后的有效波高比非耦合结果更接近观测值,耦合模式中海流的存在对有效波高的分布有明显的影响。研究表明,特别是在有效波高峰值处,海流引起的有效波高增大最大可达1 m。海浪浪向及流向的空间分布以及中国近海浮标处浪向与流向的时间序列表明,流向与浪向反向时,海流的影响造成有效波高增大;二者同向时,有效波高减小。海流对有效波高的调整会沿着海浪传播的方向传播相当一段距离。在西北太平洋的海浪场计算中,引入海流的耦合模式计算结果对改善强流区海浪预报具有重要意义,并且海流的模拟精度对于高精度的海浪预报非常重要。     

围填海对海洋水动力与生态环境的影响

近10年来,中国海岸带围填海活动呈现出规模大、速度快的发展态势。围填海能带来显著的经济效益,但对海洋环境与生态的负面影响也不可忽视。针对围填海对海洋环境和生态的影响及作用机制,分别从水动力和生态系统两个方面进行了概述。围填海改变了海洋的自然几何属性(原始岸线、地形地貌、海湾面积),引起水动力环境的变化(潮汐系统和海湾水交换能力),进而影响了海湾的环境容量;围填海破坏了生物栖息地、导致生物多样性的丧失,影响到生态系统结构与功能的稳定性;水动力与生物多样性的变化可显著影响到生物地球化学过程,加速富营养化进程,恶化水质,增加生态灾害风险。目前,围填海后的生态修复策略主要有增加生物量、建设自然保护区、退陆还海3种方式;而生态补偿策略则多基于“生态系统服务功能与生境面积的大小为线性关系”,通过对其经济价值的量化后进行生态补偿与实施相关政策。国际上,生态系统服务功能的量化参数逐步纳入实际管理,并在线性关系研究的基础上,逐步纳入一些非线性的理念,使生态补偿机制更为合理化;而我国对于围填海生态效应的定量化研究及科学理论在管理政策中的实际应用仍亟待提高。整体而言,全面、准确地评估围填海对海洋环境与生态的影响离不开自然科学与社会科学的交叉与融合。     

Paratacamite from Gypsum Veins in Sandstones of the Kuqa Basin,Xinjiang, China: Implications for a New Epigenetic Cu Enrichment Mechanism

The Kuqa Basin filled with Paleogene evaporite series is located in the northeast of the Tarim Basin, Xinjiang, China. It is famous for sandstone‐hosted Cu deposits formed by synsedimentary processes. However, our recent studies reveal that there has been another Cu mineralization mechanism in this basin. Field investigations show that there is a close relationship among faults, salt domes, and brine. Cu deposits are mainly located in two east–west‐trending anticlinal belts in the basin, adjacent to salt domes in the belts. Cu minerals in gypsum veins of the Jidike and Kangcun formations have been investigated by SEM, EDS, and X‐ray diffraction methods. The occurrence of paratacamite in gypsum veins has been reported to coexist with glauberite and halite in the joint planes of sandstones. In addition, it occurs accompanying residual crystal salt encrustation in limestone fractures, or in sandstones in dendritic form. These features indicate that the surface‐Cu enrichment in the Kuqa Basin might have originated from Cu‐bearing brine in the underlying evaporite units, which migrated upward along fractures. In addition, the presence of H₂S in the east–west fault belt in the Kuqa Basin, and the discovery of surface sulfur, calcium carbonate, and covellite, suggest thermochemical sulfate reduction near salt domes in the deeper parts of the rock units. This process resulted in the generation of reduced brine and provided a favorable environment for Cu enrichment. Therefore, the surface‐Cu mineralization near salt domes is interpreted to be the result of Cu‐bearing brine migrating upward to the surface along faults (or joints) following the intrusion of deep salt domes. The geological evidence indicates the presence of reducing brine and Cu‐bearing brine near the salt dome in the deeper rocks of the Kuqa Basin, thus making the intrusive contact zone of the salt dome a favorable site for the epigenetic enrichment of Cu. Our study demonstrates that Cu enrichment in the Kuqa Basin resulted not only from synsedimentary deposition but also through epigenetic enrichment associated with salt dome intrusion and brine‐rich fluids.     

Assessment of small-diameter shallow wells for managed aquifer recharge at a site in southern Styria,Austria

An approach to establish the recharge component of managed aquifer recharge (MAR) has recently been proposed that uses small-diameter shallow wells installed using relatively inexpensive drilling methods such as direct push. As part of further development of that approach, a generalized procedure is presented for a technical and economic assessment of the approach’s potential in comparison to other systems. Following this procedure, the use of small-diameter wells was evaluated both experimentally and numerically for a site located in southern Styria, Austria. MAR is currently done at the site using a horizontal pipe infiltration system, and system expansion has been proposed with a target rate of 12 l/s using small-diameter wells as one possible option. A short-duration single-well field recharge experiment (recharge rate 1.3–3.5 l/s) was performed (recharge by gravity only). Numerical modeling of the injection test was used to estimate hydraulic conductivity (K). Quasi-steady-state, single-well recharge simulations for different locations, as well as a long-term transient simulation, were performed using the K value calibrated from the field injection test. Results indicate that a recharge capacity of 4.1 l/s was achievable with a maximum head rise of 0.2 m at the injection well. Finally, simulations were performed for three different well fields (4, 6 and 8 wells, respectively) designed to infiltrate a target rate of 12 l/s. The experimental and numerical assessments, supported by a cost analysis of the small-diameter wells, indicate that the small-diameter wells are a viable, cost-effective recharge approach at this and other similar sites.     

A multiobjective steepest descent method with applications to optimal well control

Multiobjective optimization deals with mathematical optimization problems where two or more objective functions (cost functions) are to be optimized (maximized or minimized) simultaneously. In most cases of interest, the objective functions are in conflict, i.e., there does not exist a decision (design) vector (vector of optimization variables) at which every objective function takes on its optimal value. The solution of a multiobjective problem is commonly defined as a Pareto front, and any decision vector which maps to a point on the Pareto front is said to be Pareto optimal. We present an original derivation of an analytical expression for the steepest descent direction for multiobjective optimization for the case of two objectives. This leads to an algorithm which can be applied to obtain Pareto optimal points or, equivalently, points on the Pareto front when the problem is the minimization of two conflicting objectives. The method is in effect a generalization of the steepest descent algorithm for minimizing a single objective function. The steepest-descent multiobjective optimization algorithm is applied to obtain optimal well controls for two example problems where the two conflicting objectives are the maximization of the life-cycle (long-term) net-present-value (NPV) and the maximization of the short-term NPV. The results strongly suggest the multiobjective steepest-descent (MOSD) algorithm is more efficient than competing multiobjective optimization algorithms.     

Stability analysis of an ancient landslide considering shear strength reduction behavior of slip zone soil

In this paper, the stability of an ancient landslide during the first impounding of a nearby reservoir is investigated through the analyses of the shear strength reduction behavior of slip zone soil. In view of the experimental observations, an empirical strain-dependent soil model is established and is then incorporated in finite element analyses. The numerical analysis results show that the failing sections progressively develop due to the soil strength declines from peak toward residual, and the shear zone propagates within the front slope. It is demonstrated in the numerical results that the toe weighting measure has a significant effect on restraining the shear displacements of the soils and preventing the progressive failure of the landslide. The field observations further confirmed the stability condition of the reinforced landslide.     

Two-dimensional landslide dynamic simulation based on a velocity-weakening friction law

Surprisingly, hypermobility (high velocity and long run-out) is a remarkable feature of large landslides and is still poorly understood. In this paper, a velocity-weakening friction law is incorporated into a depth-averaged landslide model for explaining the higher mobility mechanism of landslides. In order to improve the precision of the calculation, a coupled numerical method based on the finite volume method is proposed to solve the model equations. Finally, several numerical tests are performed to verify the stability of the algorithm and reliability of the model. The comparison between numerical results and experimental data indicates that the presented model can predict the movement of landslide accurately. Considering the effect of velocity-weakening friction law, the presented model can better reflect the hypermobility of landslide than the conventional Mohr–Coulomb friction model. This work shows that the application of a universal velocity-weakening friction law is effective in describing the hypermobility of landslide and predicting the extent of landslides.     

Characterization,dissolution and solubility of the hydroxypyromorphite–hydroxyapatite solid solution [(Pb<Subscript>x</Subscript>Ca<Subscript>1−x</Subscript>)<Subscript>5</Subscript>(PO<Subscript>4</Subscript>)<Subscript>3</Subscript>OH] at 25 °C and pH 2–9

Background

The interaction between Ca-HAP and Pb²⁺ solution can result in the formation of a hydroxyapatite–hydroxypyromorphite solid solution [(Pb_xCa_1?x)₅(PO₄)₃(OH)], which can greatly affect the transport and distribution of toxic Pb in water, rock and soil. Therefore, it’s necessary to know the physicochemical properties of (Pb_xCa_1?x)₅(PO₄)₃(OH), predominantly its thermodynamic solubility and stability in aqueous solution. Nevertheless, no experiment on the dissolution and related thermodynamic data has been reported.

Results

Dissolution of the hydroxypyromorphite–hydroxyapatite solid solution [(Pb_xCa_1?x)₅(PO₄)₃(OH)] in aqueous solution at 25 °C was experimentally studied. The aqueous concentrations were greatly affected by the Pb/(Pb + Ca) molar ratios (X_Pb) of the solids. For the solids with high X_Pb [(Pb_0.89Ca_0.11)₅(PO₄)₃OH], the aqueous Pb²⁺ concentrations increased rapidly with time and reached a peak value after 240–720 h dissolution, and then decreased gradually and reached a stable state after 5040 h dissolution. For the solids with low X_Pb (0.00–0.80), the aqueous Pb²⁺ concentrations increased quickly with time and reached a peak value after 1–12 h dissolution, and then decreased gradually and attained a stable state after 720–2160 h dissolution.

Conclusions

The dissolution process of the solids with high X_Pb (0.89–1.00) was different from that of the solids with low X_Pb (0.00–0.80). The average K _sp values were estimated to be 10^?80.77±0.20 (10^?80.57–10^?80.96) for hydroxypyromorphite [Pb₅(PO₄)₃OH] and 10^?58.38±0.07 (10^?58.31–10^?58.46) for calcium hydroxyapatite [Ca₅(PO₄)₃OH]. The Gibbs free energies of formation (ΔG _f ^o ) were determined to be ?3796.71 and ?6314.63 kJ/mol, respectively. The solubility decreased with the increasing Pb/(Pb + Ca) molar ratios (X_Pb) of (Pb_xCa_1?x)₅(PO₄)₃(OH). For the dissolution at 25 °C with an initial pH of 2.00, the experimental data plotted on the Lippmann diagram showed that the solid solution (Pb_xCa_1?x)₅(PO₄)₃(OH) dissolved stoichiometrically at the early stage of dissolution and moved gradually up to the Lippmann solutus curve and the saturation curve for Pb₅(PO₄)₃OH, and then the data points moved along the Lippmann solutus curve from right to left. The Pb-rich (Pb_xCa_1?x)₅(PO₄)₃(OH) was in equilibrium with the Ca-rich aqueous solution.

Open image in new window

Graphical abstractLippmann diagrams for dissolution of the hydroxypyromorphite–hydroxyapatite solid solution [(Pb_xCa_1?x)₅(PO₄)₃OH] at 25??C and an initial pH of 2.00.