Integrated surface–groundwater flow modeling: A free-surface overland flow boundary condition in a parallel groundwater flow model

Interactions between surface and groundwater are a key component of the hydrologic budget on the watershed scale. Models that honor these interactions are commonly based on the conductance concept that presumes a distinct interface at the land surface, separating the surface from the subsurface domain. These types of models link the subsurface and surface domains via an exchange flux that depends upon the magnitude and direction of the hydraulic gradient across the interface and a proportionality constant (a measure of the hydraulic connectivity). Because experimental evidence of such a distinct interface is often lacking in field systems, there is a need for a more general coupled modeling approach.     

Finite difference modeling of contaminant transport using analytic element flow solutions

The two-dimensional implementation of the analytic element method (AEM) is commonly used to simulate steady-state saturated groundwater flow phenomena at regional and local scales. However, unlike alternative groundwater flow simulation methods, AEM results are not ordinarily used as the basis for simulation of reactive solute transport. The use of AEM-simulated flow fields is impeded by the discrepancy between a continuous representation of flow and a typically discrete representation of transport, and requires translation of the flow solution to a discrete analog. This paper presents a variety of methods for analytically calculating conservative discrete water fluxes and integrated components of the dispersion tensor across cell interfaces. An Eulerian finite difference method based on these AEM-derived parameters is implemented for use in simulation of 2D (vertically averaged) solute transport. This implementation is first benchmarked against existing methods that use standard finite difference flow solutions, then used to investigate the effects of an inaccurate discrete water balance. It is shown that improper translation of AEM fluxes leads to significant water balance errors and inaccurate simulation of contaminant transport.     

中东亚干旱半干旱区C3植物δ13 C值的分布 及其对气候的响应*

植物化石和土壤中的有机质碳同位素指标常用来反映古气候的变化,然而碳同位素这个指标在特定地区反映气候的定量关系缺乏检验。研究剖面选择自中国的秦岭(34°14'24″N,106°55'30″E)到蒙古人民共和国北部,接近贝加尔湖地区(51°35'08″N, 100°45'49″E)的研究剖面线,选择了3种C₃植物(Artemisia scoparia, Ajania achilleides 和 Artemisia frigida),在剖面线上沿南北方向上每隔4'到5'采取一个样点,共选取161个C₃植物茎叶样品进行了δ¹³ C值测定。同时收集了剖面线附近气象站的降水、气温等资料,用插值方法得到每个采样点的气温、降水数据。分析表明:C₃植物的δ¹³ C值分布范围为-30 ‰ ~-22 ‰ ,其平均值为-26.81 ‰ ,该平均值较全球C₃植物δ¹³ C平均值偏正。通过对比C₃植物δ¹³ C与年均温、年均降水量、生长季节的干燥度等随纬度的变化规律,发现C₃植物δ¹³ C、年均降水量、生长季节的干燥度有非常一致的变化趋势,而C₃植物δ¹³ C和年均温不具有一致性。通过一元回归分析也同样发现C₃植物δ¹³ C与年均降水量呈线性负相关关系(y=-0.0077x-24.838,n=161,R²=0.4418,p=0.01),与生长季节的干燥度呈线性正相关关系(y=0.7328x-28.806,n=161,R²=0.3685,p=0.01),而与年均温度没有明显的相关关系(y=-0.0461x-26.756,n=161,R²=0.0232,p=0.01)。在本研究区C₃植物δ¹³ C对年均降水量和生长季节的干燥度响应十分显著,而对温度的响应不明显。研究区具有明显的降水和温度的梯度分布特征,是验证植物碳同位素与气候关系的理想场所,而土壤中的有机质碳同位素与其地面上的植物碳同位素息息相关。研究也说明,在本研究区或其他气候植物组合相似的地区可以利用古土壤中的有机质碳同位素来定量或半定量地反映古气候的变化。     

三维GIS技术在露采矿山环境管理中的应用

在简要介绍GIS技术的基础上,对传统平面图与三维模型进行了比较。结合实例介绍MAPGIS三维模型的建立及其在露采矿山环境管理信息系统中的应用,并对露采矿山生产、管理及决策上提供更加直观、有效的支持。     

Global force balance of region 1 current system

At times of strong solar wind forcing such as those that produce major magnetic storms, the region 1 current system dominates over the Chapman–Ferraro current system in mediating the transfer of force between the solar wind and the terrestrial system. The global force balance can be broken into two components, one involving the high-altitude part of the region 1 current system that is in contact with the solar wind (labeled here the HRS) and the other involving the low-altitude part of the region 1 current system that lies in the ionosphere (the LRS). Both communicate their J×B force to the geomagnetic dipole via a gradient in the magnetic field that they generate. In the HRS case the force acts to push the dipole away from the sun. This is the region 1 analog of the Chapman–Ferraro mechanism for transmitting the solar wind's force to the Earth. However, in the LRS case, the force (which is much stronger than in the HRS case) acts to push the dipole toward the sun, seemingly paradoxically. The LRS balances the ‘paradoxical’ sunward force on the dipole with an opposite force on the atmosphere. This paper uses MHD simulations to demonstrate the presence of both the normal force-transmitting gradient generated by the Chapman–Ferraro and the counter-Chapman–Ferraro gradient in the magnetic field generated by the region 1 current system.     

Static and dynamic response of massive caisson foundations with soil and interface nonlinearities—validation and results

The static, cyclic, and dynamic response of a massive caisson foundation embedded in nonlinear layered or inhomogeneous soil and loaded at its top is investigated. The caisson is supported against horizontal displacement and rotation by four types of inelastic springs and dashpots, described with the BWGG model that was developed in the preceding companion paper [Gerolymos N, Gazetas G. Development of winkler model for static and dynamic response of caisson foundations with soil and interface nonlinearities. Soil Dyn Earthq Eng, submitted companion paper]. The prediction of the model is satisfactorily compared with results from 3D-finite element analysis. Some experimental corroboration of the method is provided with the help of a 1/3-scale lateral load test that had been conducted in the field by EPRI. An illustrative example of a caisson embedded in linearly-inhomogeneous clay and subjected to static and dynamic loading is analysed. Characteristic results are presented highlighting the role of soil inelasticity and its interplay with the two dominant interface nonlinearities: separation (gapping) of the caisson shaft from the surrounding soil, and uplifting of the base from the underlying soil.     

1-（4-硝基苯基）-3-（3，5-二溴吡啶）三氮烯的合成及其与镍的显色反应

合成了1-(4-硝基苯基)-3-(3,5-二溴吡啶)三氮烯(NPDBPDT),并研究了NPD-BPDT与镍的显色反应。在pH 10.0的Na2B4O7-NaOH缓冲溶液介质中,表面活性剂Triton X-100存在下,镍与试剂生成摩尔比为1∶4的络合物在470 nm处有最大正吸收峰,表观摩尔吸光系数为1.02×105L.mol-1.cm-1;在540 nm处有最大负吸收峰,表观摩尔吸光系数为1.15×105L.mol-1.cm-1;以540 nm为参比波长、470 nm为测量波长进行双波长测定,表观摩尔吸光系数为2.17×105L.mol-1.cm-1。镍含量在0~480μg/L内符合比尔定律。方法用于合金样品中镍的测定,结果与原子吸收法相符,加标回收率为98.0%~102.0%;精密度小于5%(RSD,n=6)。     

地质公园三维电子地图的设计与实现

结合三维建模和电子地图的基本理论,根据地质公园的特点,利用ERDAS实现大面积区域三维可视化,并结合3DS MAX丰富的建模技术及强大的动画制作功能,在地质景观演化方面做了初步探讨。采用“超地图”的形式来组织各种地图数据以及超媒体相关信息,从内容结构设计及特点进行论述,实现了伊通火山群的三维电子地图,展示出“伊通型火山”的特点。为科研、科普和旅游提供了方便、直观的信息表现形式,拓宽了电子地图的应用领域,同时也丰富了电子地图的研究内容。     

利用高分辨率影像建立三维虚拟校园——以"数字成都理工大学"为例

以数字化校园工程———“数字成都理工大学”为例,探讨利用航空影像和数字摄影测量技术,获取基础空间地理信息的方法及途径。并采用三维建模软件Cyber C ity对建筑物进行三维建模。而对于复杂地物,采用三维地理信息系统IMAGIS软件建立单个模型,从而实现了成都理工大学校区的三维可视化虚拟重建和漫游。实践证明,利用数字摄影测量技术和遥感影像,能够有效地获取大范围区域地表地物的三维信息,并快速地建立地面模型和三维景观。     

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

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