Hydrodynamic modeling of flushing time in a small estuary of North Bay, Florida, USA

Freshwater fraction method is popular for cost-effective estimations of estuarine flushing time in response to freshwater inputs. However, due to the spatial variations of salinity, it is usually expensive to directly estimate the long-term freshwater fraction in the estuary from field observations. This paper presents the application of the 3D hydrodynamic model to estimate the distributions of salinity and thus the freshwater fractions for flushing time estimation. For a case study in a small estuary of the North Bay in Florida, USA, the hydrodynamic model was calibrated and verified using available field observations. Freshwater fractions in the estuary were determined by integrating freshwater fractions in model grids for the calculation of flushing time. The flushing time in the North Bay is calculated by the volume of freshwater fraction divided by the freshwater inflow, which is about 2.2 days under averaged flow conditions. Based on model simulations for a time series of freshwater inputs over a 2-year period, a power regression equation has been derived from model simulations to correlate estuarine flushing time to freshwater inputs. For freshwater input varying from 12 m³/s to 50 m³/s, flushing time in this small estuary of North Bay changes from 3.7 days to 1.8 days. In supporting estuarine management, the model can be used to examine the effects of upstream freshwater withdraw on estuarine salinity and flushing time.     

Bifurcation modelling in a meandering gravel–sand bed river

The Rhine bifurcation at Pannerden forms the major distribution point for water supply in the Netherlands, distributing not only water and sediment but also flooding risks and navigability. Its morphological stability has been a concern for centuries. We present experiences from more than two decades of numerical morphological modelling of this bifurcation with a gravel–sand bed and a meandering planform. Successive computations have shown the importance of upstream approach conditions, the necessity to include physical mechanisms for grain sorting and alluvial roughness, and the need to assume a thicker active layer of the river bed than is suggested by laboratory flume experiments using a constant discharge. The active layer must be thicker in the model to account for river bed variations due to higher‐frequency discharge variations that are filtered out in morphological modelling. We discuss limitations in calibration and verification, but argue that, notwithstanding these limitations, 2D and 3D morphological models are valuable tools, not only for pragmatic applications to engineering problems, but also for revealing the limitations of established knowledge and understanding of the relevant physical processes. The application of numerical models to the Pannerden bifurcation appeared to reveal shortcomings in established model formulations that do not pose particular problems in other cases. This application is therefore particularly useful for setting the agenda for further research. Copyright © 2012 John Wiley & Sons, Ltd.     

FY-3微波成像仪遥感图像地理定位方法研究

MWRI (MicroWave Radiation Imager) is one of the payloads on our next generation polar meteorological satellite FY-3. MWRI conically scans with a fixed incident angle on the earth surface. It is the first time for Chinese remote sensor to use this scan mode. In this work, we present a geolocation method for FY-3 MWRI’s remote sensing image based on its special scan geometry. The integrated coordinate systems and the specific relationships with these coordinate systems are defined. A spatial relationship model between the remote sensing data and the earth-based coordinate system is established. This method also includes an algorithm of satellite orbit computation, which is used to get the satellite’s instantaneous velocity vector from its position. This method has been applied to MWRI’s remote sensing image geolocation. The results show that the accuracy of this method can achieve 1 pixel. The 33 GCPs (Ground Control Points) which are in the regiones of FY-3 MWRI’s observation have been collected and used to analyze the precision of the geolocation. By statistical analysis, the error along-track is about 1.5km, and the error along-scan is about 3.0km. It is obvious that this method fulfills the requirement of precision for FY-3 MWRI whose space resolution exceeds 5km.     

基于Vulcan的西岭金矿床三维地质模型构建及资源量估算

通过对西岭金矿区域地质和矿区地质的研究,搜集和整理了矿区的地质资料.在充分学习了澳大利亚Vul-can软件的技术后,基于该软件首次建立了矿区地质钻孔数据库,构建了矿区的多个三维可视化地质实体模型和块体模型,并利用距离幂次反比法和普通克吕格法2种品位估值方法分别进行了资源量估算,并与传统的地质块段法的计算结果进行了对比,...     

A bounding surface model for anisotropically overconsolidated clay incorporating thermodynamics admissible rotational hardening rule

A bounding surface model is formulated to simulate the behavior of clays that are subject to an anisotropic consolidation stress history. Conventional rotational hardening is revisited from the perspective of thermodynamics. As the free energy cannot be accumulated infinitely upon critical state failure, the deviatoric back stress must vanish. This requires the rotated yield surface to be turned back to eventually align on the hydrostatic axis in the stress plane. Noting that most of the previous propositions violate this restriction, an innovative rotational hardening rule is formulated that is thermodynamically admissible. The bounding surface framework that employs the modified yield surface is applied to simulate elastoplastic deformations for overconsolidated clays, with which the overprediction of strength on the “dry” side can be greatly improved with reasonable results. Other important features, including contractive or dilative response and hardening or softening behavior, can also be well-captured. It has been shown that the model can simulate three types of reconstituted clays that are sheared with initial conditions over a wide range of anisotropic consolidation stress ratios and overconsolidation ratios under both triaxial undrained and drained conditions. Limitations and potential improvement of the model regarding the fabric anisotropy at critical state have been discussed.     

Numerical simulation of agricultural sediment and pesticide runoff: RZWQM and PRZM comparison

Agricultural sediment and pesticide runoff is a widespread ecological and human health concern. Numerical simulation models, such as Root Zone Water Quality Model (RZWQM) and Pesticide Root Zone Model (PRZM), have been increasingly used to quantify off‐site agricultural pollutant movement. However, RZWQM has been criticized for its inability to simulate sedimentation processes. The recent incorporation of the sedimentation module of Groundwater Loading Effects of Agricultural Management Systems has enabled RZWQM to simulate sediment and sediment‐associated pesticides. This study compares the sediment and pesticide transport simulation performance of the newly released RZWQM and PRZM using runoff data from 2 alfalfa fields in Davis, California. A composite metric (based on coefficient of determination, Nash–Sutcliffe efficiency, index of agreement, and percent bias) was developed and employed to ensure robust, comprehensive assessment of model performance. Results showed that surface water runoff was predicted reasonably well (absolute percent bias <31%) by RZWQM and PRZM after adjusting important hydrologic parameters. Even after calibration, underestimation bias (?89% ≤ PBIAS  ≤ ?36%) for sediment yield was observed in both models. This might be attributed to PRZM's incorrect distribution of input water and uncertainty in RZWQM's runoff erosivity coefficient. Moreover, the underestimation of sediment might be less if the origin of measured sediment was considered. Chlorpyrifos losses were simulated with reasonable accuracy especially for Field A (absolute PBIAS  ≤ 22%), whereas diuron losses were underestimated to a great extent (?98% ≤ PBIAS  ≤ ?65%) in both models. This could be attributed to the underprediction of herbicide concentration in the top soil due to the limitations of the instantaneous equilibrium sorption model as well as the high runoff potential of herbicide formulated as water‐dispersible granules. RZWQM and PRZM partitioned pesticides into the water and sediment phases similarly. According to model predictions, the majority of pesticide loads were carried via the water phase. On the basis of this study, both RZWQM and PRZM performed well in predicting runoff that carried highly adsorptive pesticides on an event basis, although the more physically based RZWQM is recommended when field‐measured soil hydraulic properties are available.     

空间尺度转换与跨尺度信息链接:区域生态水文模拟研究空间尺度转换方法综述

空间尺度转换是近年来区域生态水文研究领域的一个基本研究问题。其需要主要是源于模型的输入数据与所能提供的数据空间尺度不一致以及模型所代表的地表过程空间尺度与所观测的地表过程空间尺度不吻合。综述了目前区域生态水文模拟研究中常用的空间尺度转换研究方法,包括向上尺度转换和向下尺度转换。详细论述了2种向下尺度转换方法: 统计学经验模型和动态模型。前者是通过将GCM大尺度数据与长期的历史观测数据比较从而建立统计学相关模型, 然后利用这个统计学经验模型进行向下的空间尺度转换. 然而动态模型并不直接对GCM数据进行向下尺度的转换,而是对与GCM进行动态耦合的区域气候模型(RCM) 的输出数据进行空间尺度转换. 通常后者所获得的数据精度要比前者高,但是一个主要缺点就是并不是全球所有的研究区域都有对应的RCM。还详细论述了2种向上尺度转换方法: 统计学经验模型和斑块模型。前者是建立一个能代表小尺度信息在大尺度上分布的密度分布概率函数, 然后利用这个函数在所需的大尺度上进行积分而求得大尺度所需的信息。而后者是根据相似性最大化原则将大尺度划分为若干个可操作的小尺度斑块,然后将计算的每个小尺度斑块的信息平均化得到大尺度所需的信息。通常在计算这种斑块化的小尺度信息的时候,对每个小尺度也会采用统计学经验模型来计算代表整个斑块小尺度的信息。建议用斑块模型与统计学经验模型相集合的方法来实现向上的空间尺度转换     

Numerical modeling of 3-D terrain effect on MT field

Using the boundary element method, the numerical modeling problem of three-dimensional terrain effect on magnetotelluric (MT) field is solved. This modeling technique can be run on PC in the case of adopting special net division. The result of modeling test for 2-D terrain by this modeling technique is basically coincident with that by 2-D modeling technique, but there is a great difference between the results of 3-D and 2-D modeling for 3-D terrain. Project supported by the National Natural Science Foundation of China.     

道路网多尺度数据建模的初步研究

本文以路网要素的多尺度表达为例,探讨了路网要素多尺度表达在几何层次、实体层次等方面的多尺度现象,研究了多尺度数据建模的基本特性,提出了构建多尺度空间数据模型的关键策略,详细剖析了路网多尺度建模中不同尺度下的实体、实体间的联系,以及多尺度表达间要素类和要素表达层次上的关系,实现了对路网要素多尺度抽象表达的模型构筑。     

非均匀介质中交错网格高阶有限差分数值模拟

地震波场的数值模拟一直是地球物理学的一个重要的研究领域,而在数值正演模拟方法的研究中,计算精度和计算效率是评价该方法有效性及优越性的二个关键问题。这里从一阶速度—应力弹性波动方程出发,着重介绍如何构造离散化模型的网格,如何求解空间导数,如何选取边界条件等内容,从而更有效地提高数值计算的精度与计算效率。文中构造了不同类型的介质模型,并在交错网格中,利用高阶有限差分模拟非均匀介质的波场传播。模拟结果表明,该方法实现简单,具有很好地稳定性和较高的精度,能够直观、高效地反映出介质中波场的传播规律。