Beaver dam analogues drive heterogeneous groundwater–surface water interactions

Beaver dam analogues (BDAs) are a cost-effective stream restoration approach that leverages the recognized environmental benefits of natural beaver dams on channel stability and local hydrology. Although natural beaver dams are known to exert considerable influence on the hydrologic conditions of a stream system by mediating geomorphic processes, nutrient cycling, and groundwater–surface water interactions, the impacts of beaver-derived restoration methods on groundwater–surface water exchange are poorly characterized. To address this deficit, we monitored hyporheic exchange fluxes and streambed porewater biogeochemistry across a sequence of BDAs installed along a central Wyoming stream during the summer of 2019. Streambed fluxes were quantified by heat tracing methods and vertical hydraulic gradients. Biogeochemical activity was evaluated using major ion porewater chemistry and principal component analysis. Vertical fluxes of approximately 1.0 m/day were observed around the BDAs, as was the development of spatially heterogeneous zones of nitrate production, groundwater upwelling, and anaerobic reduction. Strong contrasts in hyporheic zone processes were observed across BDAs of differing sizes. This suggests that structures may function with size-dependent behaviour, only altering groundwater–surface water interactions after a threshold hydraulic step height is exceeded. Patterns of hyporheic exchange and biogeochemical cycling around the studied BDAs resemble those around natural beaver dams, suggesting that BDAs may provide comparable benefits to channel complexity and near-stream function over a 1-year period.     

有限元法与伪谱法混合求解弹性波动方程

在地震波场数值模拟中,有限差分法、有限元法和伪谱法都是常用的基本方法,但它们各有不同的适应性和优缺点,如有限差分法、有限元法都存在减弱网格频散和提高计算效率的矛盾,而伪谱法的网格频散小且计算效率高.有限差分法和伪谱法在处理地表结构复杂或地表剧烈起伏以及地下结构复杂的情况时存在较大的难度,而有限元法可较为理想地拟合起伏地表和任意弯曲界面,且可方便地处理自由边界条件和界面边界条件.尝试将有限元法和伪谱法相结合,形成地震波场数值模拟的一种混合方法,利用二者的优点,克服二者的缺点,达到既减弱网格频散又提高计算精度和效率的目的.并采用所谓的‘过度区域‘技术解决两种不同算法的衔接问题.模拟实例表明,给出的混合模拟方法不失为弹性波场数值模拟的一种有效方法.     

拉格朗日插值法在渔业资源评估上的应用

应用拉格朗日插值法拟合鱼类体长与体重之间的函数关系 ,使用计算机VB语言进行编写程序 ,与常用的线性回归法、Ricker法进行比较 ,得出体长与体重的关系式为Pn(x) =∑nk=0yklk(x)=∑nk=0(∏nj=0j≠kx-xjxk-xj)yk 。     

Non-Existence of New Meromorphic First Integrals in the Planar Three-Body Problem

We consider the planar three-body problem and prove that, apart from some exceptional cases, there is no additional first integral meromorphic with respect to positions, mutual distances and momenta.     

Crystal-Melt Separation and the Development of Isotopic Heterogeneities in Hybrid Magmas

If a magma is a hybrid of two (or more) isotopically distinctend-members, at least one of which is partially crystalline,separation of melt and crystals after hybridization will leadto the development of isotopic heterogeneities in the magmaas long as some of the pre-existing crystalline material (antecrysts)retains any of its original isotopic composition. This holdstrue whether the hybridization event is magma mixing as traditionallyconstrued, bulk assimilation, or melt assimilation. Once a magma-scaleisotopic heterogeneity is formed by crystal–melt separation,it is essentially permanent, persisting regardless of subsequentcrystallization, mixing, or equilibration events. The magnitudeof the isotopic variability resulting from crystal–meltseparation can be as large as that resulting from differentialcontamination, multiple isotopically distinct sources, or insitu isotopic evolution. In one model, a redistribution of one-thirdof the antecryst cargo yielded a crystal-enriched sample with⁸⁷Sr/⁸⁶Sr of 0·7058, whereas the complementary crystal-poorsample has ⁸⁷Sr/⁸⁶Sr of 0·7068. In other models, crystal-richsamples are enriched in radiogenic Sr. Isotopic heterogeneitiescan be either continuous (controlled by the modal distributionof crystals and melt) or discontinuous (when there is completeseparation of crystals and liquid). The first case may be exemplifiedby some isotopically zoned large-volume rhyolites, formed bythe eruptive inversion of a modally zoned magma chamber. Inthe latter case, the isotopic composition of any (for example)interstitial liquid will be distinct from the isotopic compositionof the bulk crystal fraction. The separation of such an interstitialliquid may explain the presence of isotopically distinct late-stageaplites in plutons. Crystal–melt separation provides anadditional option for the interpretation of isotopically zonedor heterogeneous magmas. This option is particularly attractivefor systems whose chemical variation is otherwise explicableby fractionation-dominated processes. Non-isotopic chemicalheterogeneities can also develop in this fashion. KEY WORDS: isotopic heterogeneity; zoning; hybrid magma; crystal separation; Sr isotopes; aplite; rhyolite     

Numerical Simulation of Residual Circulation due to Bottom Rouglmess Variability Under Tidal Flows in A Semi-Enclosed Bay

Nowadays there are some chronic serious environmental problems, such as eutrophication, blue tide and so on, in a complicated coastal zone or a semi-enclosed bay, because the water exchanges between an inner bay and an outer sea is weak compared with the supply of contaminant. Under this situation, a method to improve the water quality by 3-dimensional small unsymmetrical structures has been proposed by Komatsu et al. In this paper, several numerical simulations of the tidal current and concentration for various arrangements of bottom roughness in a semi-enclosed model bay are carfled out with a depth-averaged 2-D numerical model. The model is solved by the hybrid finite analytic method with nonstaggered grid. And the SIMPLES algorithm with Rhie and Chow＇ s momentum interpolation technique is used for the simulation. The effect of Komatsu＇ s method for water purification is examined by numerical simulation. The result of numerical experiment indicates that it is possible to generate a new tidal residual current and to activate a tidal exchange by bottom roughness arrangement only.     

Estuary/ocean exchange and tidal mixing in a Gulf of Maine Estuary: A Lagrangian modeling study

A Lagrangian particle method embedded within a 2-D finite element code, is used to study the transport and ocean–estuary exchange processes in the well-mixed Great Bay Estuarine System in New Hampshire, USA. The 2-D finite element model, driven by residual, semi-diurnal and diurnal tidal constituents, includes the effects of wetting and drying of estuarine mud flats through the use of a porous medium transport module. The particle method includes tidal advection, plus a random walk model in the horizontal that simulates sub-grid scale turbulent transport processes. Our approach involves instantaneous, massive [O(500,000)] particle releases that enable the quantification of ocean–estuary and inter-bay exchanges in a Markovian framework. The effects of the release time, spring–neap cycle, riverine discharge and diffusion strength on the intra-estuary and estuary–ocean exchange are also investigated.The results show a rather dynamic interaction between the ocean and the estuary with a fraction of the exiting particles being caught up in the Gulf of Maine Coastal Current and swept away. Three somewhat different estimates of estuarine residence time are calculated to provide complementary views of estuary flushing. Maps of residence time versus release location uncover a strong spatial dependency of residence time within the estuary that has very important ramifications for local water quality. Simulations with and without the turbulent random walk show that the combined effect of advective shear and turbulent diffusion is very effective at spreading particles throughout the estuary relatively quickly, even at low (1 m²/s) diffusivity. The results presented here show that a first-order Markov Chain approach has applicability and a high potential for improving our understanding of the mixing processes in estuaries.     

Numerical Analysis of Wave Forces on Arbitrarily Shaped Multi-Pillars Over Varying Topography

The modified hybrid element method (MHEM) is utilized to predict and analyze wave forces on arbitrarily shaped multiple bodies. This method can be applied to waves of all water depths, i. e. shallow, intermediate, and deep waters, on slowly varying seabed. The MHEM employs the ICCG method to save CPU and storage, thus the computation of wave forces for large multi-body systems can be carried out on microcomputers. Numerical results of the present method are compared with experimental data and other solutions. It is shown that the MHEM provides more accurate solutions of the wave forces than other numerical methods do. Therefore, the methodology presented herein can be used in the design of coastal and ocean structures.     

现场海域水平扩散系数的探讨 Ⅱ．黄河三角洲埕北海域的计算

根据Ｔａｙｌｏｒ的湍流扩散理论和Ｈａｙ、Ｐａｓｑｕｉｌｌ提出的拉格朗日流和欧拉流的相关函数的相似性，用黄河三角洲埕北海域的实测资料，按四种不同方案计算出该海域的水平扩散系数和谱密度分布。指出：该海区潮流扩散系数主要受潮流的影响，其分布特征与潮流的分布基本一致，所求的水平扩散系数主要是由周期性的潮运动给予潮流扰动所致。方案４是通过富氏变换滤波消去了周期性潮流扰动，本方案计算的水平扩散系数是由周期为６ｈ以下的小尺度涡提供能量。