Coupled inverse modelling of groundwater flow and mass transport and the worth of concentration data

This paper presents the extension of the self-calibrating method to the coupled inverse modelling of groundwater flow and mass transport. The method generates equally likely solutions to the inverse problem that display the variability as observed in the field and are not affected by a linearisation of the state equations. Conditioning to the state variables is measured by an objective function including, among others, the mismatch between the simulated and measured concentrations. Conditioning is achieved by minimising the objective function by gradient-based methods. The gradient contains the partial derivatives of the objective function with respect to: log conductivities, log storativities, prescribed heads at boundaries, retardation coefficients and mass sources. The derivatives of the objective function with respect to log conductivity are the most cumbersome and need the most CPU-time to be evaluated. For this reason, to compute this derivative only advective transport is considered. The gradient is calculated by the adjoint-state method. The method is demonstrated in a controlled, synthetic study, in which the worth of concentration data is analysed. It is shown that concentration data are essential to improve transport predictions and also help to improve aquifer characterisation and flow predictions, especially in the upstream part of the aquifer, even in the case that a considerable amount of other experimental data like conductivities and heads are available. Besides, conditioning to concentration data reduces the ensemble variances of estimated transmissivity, hydraulic head and concentration.     

Recent ventifact development on the central Oregon coast,western USA

The unusual location of ventifacts, on a boulder‐built jetty at the mouth of the Siuslaw River, Oregon coast, western USA, allows ventifact age and wind abrasion rates to be estimated with some precision. The jetty was built mainly between 1892–1901 and extended throughout the twentieth century. Consideration of historical shoreline position and the history of jetty construction and repair suggests the ventifacts have formed since about 1930. Morphologically the ventifacts are aligned south‐to‐north reflecting winter winds and sediment transport from the adjacent beach. Wind‐parallel grooves and ridges with sharp, sinuous crests are developed on inclined boulder surfaces on top of the jetty and reflect suspended sand transport in wind vortices. Deeply pitted surfaces on steep boulder surfaces nearest the beach reflect impact by saltating sand grains. Based on present wind regimes (1992–2000) from three regional weather stations, southerly winds above the sand transport threshold occur for 21·9–29·6 per cent of the time. Based on estimated depth of loss from boulder surfaces, wind abrasion rates are calculated to be on the order of 0·24–1·63 mm a^?1. This is the first well‐constrained field estimate of ventifact age and ventifaction rate from a modern coastal environment. Copyright © 2003 John Wiley & Sons, Ltd.     

Simulation and measurement of surface shear stress over isolated and closely spaced transverse dunes in a wind tunnel

Topographic interactions generate multidirectional and unsteady air?ow that limits the application of velocity pro?le approaches for estimating sediment transport over dunes. Results are presented from a series of wind tunnel simulations using Irwin‐type surface‐mounted pressure sensors to measure shear stress variability directly at the surface over both isolated and closely spaced sharp‐crested model dunes. Findings complement existing theories on secondary air?ow effects on stoss transport dynamics and provide new information on the in?uence of lee‐side air?ow patterns on dune morphodynamics. For all speeds investigated, turbulent unsteadiness at the dune toe indicates a greater, more variable surface shear, despite a signi?cant drop in time‐averaged measurements of streamwise shear stress at this location. This effect is believed suf?cient to inhibit sediment deposition at the toe and may be responsible for documented intermittency in sand transport in the toe region. On the stoss slope, streamline compression and ?ow acceleration cause an increase in ?ow steadiness and shear stress to a maximum at the crest that is double that at the toe of the isolated dune and 60–70 per cent greater than at ?ow reattachment on the lower stoss of closely spaced dunes. Streamwise ?ow accelerations, rather than turbulence, have greater in?uence on stress generation on the stoss and this effect increases with stoss slope distance and with incident wind speed. Reversed ?ow within the separation cell generates signi?cant surface shear (30–40 per cent of maximum values) for both spacings. This supports ?eld studies that suggest reversed ?ow is competent enough to return sediment to the dune directly or in a de?ected direction. High variability in shear at reattachment indicates impact of a turbulent shear layer that, despite low values of time‐averaged streamwise stress in this region, would inhibit sediment accumulation. Downwind of reattachment, shear stress and ?ow steadiness increase within 6 h (h = dune height) of reattachment and approach upwind values by 25 h. A distance of at least 30 h is suggested for full boundary layer recovery, which is comparable to ?uvial estimates. The Irwin sensor used in this study provides a reliable means to measure skin friction force responsible for sand transport and its robust, simple, and cost‐effective design shows promise for validating these ?ndings in natural dune settings. Copyright © 2003 John Wiley & Sons, Ltd.     

The response of subaqueous dunes to floods in sand and gravel bed reaches of the Dutch Rhine

Abstract The branches of the River Rhine in the Netherlands, characterized by a sand–gravel bed in the upstream part and a sand bed in the downstream part of the river system, show migrating dunes, especially during floods. In the last 20 years, these dunes have been studied extensively. High-resolution echo-sounding measurements of these dunes, made with single and multibeam equipment, were analysed for three different sections of the Rhine river system during several floods. This analysis was done to quantify the growth, decay and migration rates of the dunes during floods. In addition, the migrating dunes were used to calculate bedload transport rates with dune tracking. The results of dune growth and decay and migration rate are shown to be very different for the various sections during the various floods, and these differences are related to differences in grain size of the bed and to differences in the distribution of discharge over the main channel and the floodplain. The relations are used to show that the growth and migration rate of dunes, and the calculated bedload transport rates during the rising stage of a flood wave can be predicted from the mobility of the bed material with simple power relations.     

索马里急流的年际变化及其对半球间水汽输送和东亚夏季降水的影响

利用美国国家环境预报中心和美国国家大气科学研究中 心（NCEP/NCAR）再分析月平均气候资料以及Xie和Arkin分析的月平均降水资料（1968～199 8年），针对索马里低空急流（SMJ）的年际变化及其对东亚夏季降水的影响问题展开了分析 研究. 结果揭示，SMJ作为最主要的越赤道气流，对两个半球间水汽输送起最关键的作用， 它把水汽从冬半球输送到夏半球. 夏季SMJ的年际变化有全球范围内的环流异常与之相联系 ，特别是东亚沿岸的波列状异常分布、南亚高压以及澳大利亚以南的偶极型异常分布；它 也同春季的北印度洋等海区的海温异常有密切关系. 研究还表明，春季SMJ的年际变化对东 亚夏季降水和大气环流有显著影响，由于SMJ影响的超前性，因此它在东亚夏季气候预测上 有重要意义.     

1991年江淮暴雨时期的能量和水汽循环研究

通过对 1991年 5～ 7月江淮暴雨期全球范围的水汽输送和不同降水过程中江淮暴雨区及其临近区域的水汽收支和视热源和视水汽汇的计算分析得到了以下结论 :( 1)从水汽输送的机制来看 ,一方面 ,有大量的水汽以定常涡动的方式从孟加拉湾及南海输送到中国江淮地区 ;另一方面 ,江淮地区的瞬变涡动水汽向北输出 ,这可能与江淮地区频繁活动的α和β中尺度系统有关 ,它们将江淮地区汇集的充沛的水汽除了大部分以降水形式降下外 ,剩余部分继续向水汽较少的高纬地区输出 ,以维持全球水汽的平衡。( 2 )在降水过程中 ,局地蒸发项在水汽的供应中或再循环中十分重要 ,其数值一般为降水量的 13 ～ 12 ,这与1998年的降水情况相似。( 3 )在降水过程中 ,暴雨区的水汽主要是从南边界和西边界流入的 ,东边界和北边界则流出 ,并且水汽的流入、流出主要在中低层进行。 (a)在流入边界上 ,水汽通量的垂直分布存在差异 ,暴雨区西边界和南边界的水汽流入的垂直差异可能与其所在的地理位置有密切关系。 (b)对于暴雨区 ,不同强度的降水过程水汽的主要来源有所不同。( 4)在 5次降水过程中 ,视热源和视水汽汇的较大值对应降水的大值区 ,表明了水汽凝结加热对大气加热所起的主要作用。梅雨期降水 ,以对流性降水为主 ,对流活动随季节变化     

一次沙尘暴中尺度低压和沙尘输送的数值模拟

利用ＰＳＵ／ＮＣＡＲ中尺度数值模式（ＭＭ４），以常规观测资料为初始场，成功地模拟了造成“９４．４．５”沙尘暴天气的主要天气系统张掖热低压。并通过给原有的中尺度模式系统加入沙尘输送方程和驱动地面沙尘抬升参数化方案，用改进后的中尺度数值模式系统，成功地模拟了沙尘抬升和输送，与实况比较一致。     

辽东湾北部浅海区泥沙输送及其沉积特征

根据实测资料，本文定量分析了辽东湾北部泥沙输送及其分布，并对辽河三角洲沉积区划作了初步讨论。研究表明本区泥沙以纵向搬运为特征。双台子河以西来沙和辽河西水道入海泥沙是区域东部拦门沙体和浅滩发育的主导因素。汛季大潮期，泥沙自西向东运移落淤在河口及毗邻浅水域；小潮期，泥沙除向东扩散外，大部泥沙向海方搬运。调查区可划分六个现代沉积作用区，即潮坪沉积区、辽河水下三角洲细粒沉积区、河口沙洲沉积区、波浪潮流冲蚀沉积区、河口冲积沉积区以及潮汐水道沉积区。     

赤峰—朝阳地区金矿成矿物理化学条件研究

系统论述了区内各大类金矿的成矿温度，成矿压力，各种气体的逸度、ｐＨ值，金离子及其配合物的活度，讨论了金的迁移形式，推测了区内各类矿液中金的沉淀机制。     

Transport in a 2-D saturated porous medium: A new method for particle tracking

A new method for solving the transport equation based on the management of a large numbe of particles in a discretized 2-D domain is presented. The method uses numerical variables to represent the number of particles in a given mesh and is more complex than the 1-D problem. The first part of the paper focuses on the specific management of particles in a 2-D problem. The method also would be valid for three dimensions as long as the medium can be modeled similar to a layered system. As the particles are no longer tracked individually, the algorithm is fast and does not depend on the number of particles present. The numerical tests show that the method is nearly numerical dispersion free and permits accurate calculations even for simulations of low-concentration transport. Because each mesh is considered as a closed system between two successive time steps, it is easy to add adsorption phenomenon without any problem of numerical stability. The model is tested under conditions that are extremely demanding for its operating mode and gives a good fit to analytical solutions. The conditions in which it can be used to best advantage are discussed.