吉林市第二松花江江段污染物-NH4+分布特征研究

针对吉林市铁合金厂排污口下游江水的污染情况，设计研究江段，分段、分层取样。通过野外与室内实验，分析研宄了江水污染物的主要成分，并对NH4^ 的纵向、横向、垂向上的分布特征及其影响因素进行了详细研究。结果表明，垂直方向污染物扩散可以瞬时完成，在极短的时间内基本能达到完全混合；在横向和纵向上，岸边排放污染物一般只在近岸扩散，污染范围一般也只在靠近排污口一带；随着离排污口距离增大，污染物浓度降低，横向变化率较纵向大。本研究为排污口的污染物控制与治理、水资源保护提供了科学依据。     

水质大涡模拟数学模型研究

建立了基于弱可压缩流基础上的水质大涡模拟数学模型及污染物质扩散系数与平均运动的关系。数值计算采用有限体积法和预测校正法,固壁边界条件采用"部分滑移条件"。并应用该模型对重庆嘉陵江、长江两江交汇段的流场,污染物质浓度分布进行了计算,计算结果与实测值比较吻合。研究表明,横向扩散与流动状态有关,当支流流量相对较小时,横向扩散较弱,形成岸边污染带,当支流流量相对较大时,横向扩散较强,污染带向江心移动。     

Postglacial topographic evolution of glaciated valleys: a stochastic landscape evolution model

The retreat of valley glaciers has a dramatic effect on the stability of glaciated valleys and exerts a prolonged influence on the subsequent fluvial sediment transport regime. We have studied the evolution of an idealized glaciated valley during the period following retreat of ice using a numerical model. The model incorporates a stochastic process to represent deep‐seated landsliding, non‐linear diffusion to represent shallow landsliding and an approximation of the Bagnold relation to represent fluvial sediment transport. It was calibrated using field data from several recent surveys within British Columbia, Canada. We present ensemble model results and compare them with results from a deterministic linear‐diffusion model to show that explicit representation of large landslides is necessary to reproduce the morphology and channel network structure of a typical postglacial valley. Our model predicts a rapid rate of fluvial sediment transport following deglaciation with a subsequent gradual decline, similar to that inferred for Holocene time. We also describe how changes in the model parameters affect the estimated magnitude and duration of the paraglacial sediment pulse. Copyright © 2005 John Wiley & Sons, Ltd.     

Intergranular diffusion rates from the analysis of garnet surfaces: Implications for metamorphic equilibration

Novel approaches to garnet analysis have been used to assess rates of intergranular diffusion between different matrix phases and garnet porphyroblasts in a regionally metamorphosed staurolite‐mica‐schist from the Barrovian‐type area in Scotland. X‐ray maps and chemical traverses of planar porphyroblast surfaces reveal chemical heterogeneity of the garnet grain boundary linked to the nature of the adjacent matrix phase. The garnet preserves evidence of low temperature retrograde exchange with matrix minerals and diffusion profiles documenting cation movement along the garnet boundaries. Garnet–quartz and garnet–plagioclase boundaries preserve evidence of sluggish Mg, Mn and Fe diffusion at comparable rates to volume diffusion in garnet, whereas diffusion along garnet–biotite interfaces is much more effective. Evidence of particularly slow Al transport, probably coupled to Fe³⁺ exchange, is locally preserved on garnet surfaces adjacent to Fe‐oxide phases. The Ca distribution on the garnet surface shows the most complex behaviour, with long‐wavelength heterogeneities apparently unrelated to the matrix grain boundaries. This implies that the Ca content of garnet is controlled by local availability and is thought likely to reflect disequilibrium established during garnet growth. Geochemical anomalies on the garnet surfaces are also linked to the location of triple junctions between the porphyroblasts and the matrix phases, and imply enhanced transport along these channels. The slow rates of intergranular diffusion and the characteristics of different boundary types may explain many features associated with the prograde growth of garnet porphyroblasts. Thus, minerals such as quartz, Fe‐oxides and plagioclase whose boundaries with garnet are characterized by slow intergranular diffusion rates appear to be preferentially trapped as inclusions within porphyroblasts. As such grain boundary diffusion rates may be a significant kinetic impediment to metamorphic equilibrium and garnet may struggle to maintain chemical and textural equilibrium during growth in pelites.     

Interpreting zirconium‐in‐rutile thermometric results

At first sight, experimental results and observations on rocks suggest that the Zr content in rutile, where equilibrated with quartz and zircon, should be a useful thermometer for metamorphic rocks. However, diffusion data for Zr in rutile imply that thermometry should not, for plausible rates of cooling, give the high temperatures commonly observed in high‐grade metamorphic rocks. It is suggested here that such observations can be accounted for by high‐T diffusive closure of Si in rutile, causing the interior of rutile grains to become insensitive to the thermometer equilibrium well above the temperature of Zr diffusive closure. Paired with comparatively slow grain boundary diffusion and problematic zircon nucleation, this allows for cases of Zr retention in rutile through temperatures where Zr is still diffusively mobile within rutile grains. Other observations that may be accounted for in this context are large inter‐grain ranges of rutile Zr contents uncorrelated with rutile grain size, and flat Zr profiles across individual rutile grains, counter to what would be expected from diffusive closure. A consequence is that it is unlikely that Zr‐in‐rutile thermometry will be useful for estimating rock cooling rates.     

Controls on metamorphic equilibration: the importance of intergranular solubilities mediated by fluid composition

A dramatic demonstration of the role of intergranular solubility in promoting chemical equilibration during metamorphism is found in the unusual zoning of garnet in pelitic schist exposed at Harpswell Neck, Maine, USA. Many garnet crystals have irregular, patchy distributions of Mn, Cr, Fe and Mg in their inclusion‐rich interiors, transitioning to smooth, concentric zoning in their inclusion‐poor outer rims; in contrast, zoning of Ca and Y is comparatively smooth and concentric throughout. We re‐assess the disputed origin of these zoning features by examining garnet growth in the context of the thermal and structural history of the rocks, and by evaluating the record of fluid–rock interaction revealed in outcrop‐scale veining and fluid‐inclusion assemblages. The transition in the character of garnet zoning correlates with the onset of a synkinematic, simple‐shear‐dominated phase of garnet growth and with a shift in the composition of the intergranular fluid from CO₂‐rich to H₂O‐rich. Compositional variations in garnet are therefore best explained by a two‐stage growth history in which intergranular diffusive fluxes reflect differences in the concentration of dissolved species in these two contrasting fluids. Interiors of garnet crystals grew in the presence of a CO₂‐rich fluid, in which limited solubility for Mn and Cr (and perhaps Fe and Mg) produced patchy disequilibrium overprint zoning, while appreciable solubility for Ca and Y permitted their rock‐wide equilibration. Rims grew in the presence of an H₂O‐rich fluid, in which high intergranular concentrations for all elements except Cr enabled diffusion over length scales sufficient for rock‐wide equilibration. This striking example of partial chemical equilibrium during reaction and porphyroblast growth implies that thermal effects may commonly be subsidiary in importance to solubilities in the intergranular medium as determinants of length scales for metamorphic equilibration.     

Oxygen isotope speedometry in granulite facies garnet recording fluid/melt–rock interaction (Sør Rondane Mountains,East Antarctica)

In situ analysis of a garnet porphyroblast from a granulite facies gneiss from Sør Rondane Mountains, East Antarctica, reveals discontinuous step‐wise zoning in phosphorus and large δ¹⁸O variations from the phosphorus‐rich core to the phosphorus‐poor rim. The gradually decreasing profile of oxygen isotope from the core (δ¹⁸O = ~15‰) to the rim (δ¹⁸O = ~11‰) suggests that the ¹⁸O/¹⁶O zoning was originally step‐wise, and modified by diffusion after the garnet rim formation at ~800°C and 0.8 GPa. Fitting of the ¹⁸O/¹⁶O data to the diffusion equation constrains a duration of the high‐T event (~800°C) to c. 0.5–40 Ma after the garnet rim formation. The low δ¹⁸O value of the garnet rim, together with the previously reported low δ¹⁸O values in metacarbonates, indicates regional infiltration, probably along a detachment fault, of low δ¹⁸O fluid/melt possibly derived from meta‐mafic to ultramafic rocks.     

Diffusive mass exchange of non‐reactive substances in dual‐porosity porous systems – column experiments under saturated conditions

Diffusive mass exchange into immobile water regions within heterogeneous porous aquifers influences the fate of solutes. The percentage of immobile water is often unidentified in natural aquifers though. Hence, the mathematical prediction of solute transport in such heterogeneous aquifers remains challenging. The objective of this study was to find a simple analytical model approach that allows quantifying properties of mobile and immobile water regions and the portion of immobile water in a porous system. Therefore, the Single Fissure Dispersion Model (SFDM), which takes into account diffusive mass exchange between mobile and immobile water zones, was applied to model transport in well‐defined saturated dual‐porosity column experiments. Direct and indirect model validation was performed by running experiments at different flow velocities and using conservative tracer with different molecular diffusion coefficients. In another column setup, immobile water regions were randomly distributed to test the model applicability and to determine the portion of immobile water. In all setups, the tracer concentration curves showed differences in normalized maximum peak concentration, tailing and mass recovery according to their diffusion coefficients. These findings were more pronounced at lower flow rates (larger flow times) indicating the dependency of diffusive mass exchange into immobile water regions on tracers' molecular diffusion coefficients. The SFDM simulated all data with high model efficiency. Successful model validation supported the physical meaning of fitted model parameters. This study showed that the SFDM, developed for fissured aquifers, is applicable in porous media and can be used to determine porosity and volume of regions with immobile water. Copyright © 2015 John Wiley & Sons, Ltd.     

Efficient DEM‐based overland flow routing using integrated recursive algorithms

Two primary concerns in performing watershed overland flow routing are the numerical instability and computational efficiency. The stability of executing an explicit scheme has to be maintained by observing the Courant–Friedrich–Lewy criterion, which is adopted to confirm that the numerical marching speed is larger than the wave celerity. Moreover, there is another criterion of time step devised in previous studies to avoid back‐and‐forth refluxing between adjacent grids. The situation of refluxing usually occurs on flat regions. In light of this, the selection of a small time increment to honor both restrictions simultaneously is believed to decrease the computational efficiency in performing overland flow routing. This study aims at creating a robust algorithm to relax both restrictions. The proposed algorithm was first implemented on a one‐dimensional overland plane to evaluate the accuracy of the numerical result by comparing it with an analytical solution. Then, the algorithm was further applied to a watershed for 2D runoff simulations. The results show that the proposed integrated algorithm can provide an accurate runoff simulation and achieve satisfactory performance in terms of computational speed.     

A discontinuous finite element suspended sediment transport model for water quality assessments in river networks

Suspended sediment plays an important role in the distribution and transport of many pollutants (such as radionuclides) in rivers. Pollutants may adsorb on fine suspended particles (e.g. clay) and spread according to the suspended sediment movement. Hence, the simulation of the suspended sediment mechanism is indispensable for realistic transport modelling. This paper presents and tests a simple mathematical model for predicting the suspended sediment transport in river networks. The model is based on the van Rijn suspended load formula and the advection–diffusion equation with a source or sink term that represents the erosion or deposition fluxes. The transport equation is solved numerically with the discontinuous finite element method. The model evaluation was performed in two steps, first by comparing model simulations with the measured suspended sediment concentrations in the Grote Nete–Molse Nete River in Belgium, and second by a model intercomparison with the sediment transport model NST MIKE 11. The simulations reflect the measurements with a Nash‐Sutcliffe model efficiency of 0.6, while the efficiency between the proposed model and the NST MIKE 11 simulations is 0.96. Both evaluations indicate that the proposed sediment transport model, that is sufficiently simple to be practical, is providing realistic results.