一个高分辨对流层物质交换模式

在区域酸沉积模式(RADM)基础上,设计了一个高分辨对流层(TL)物质交换模式(EM3)。EM3在以下五个方面明显不同于RADM:第一,增加行星边界层(PBL)的垂直网格层次;第二、TL内考虑分层的垂直涡旋扩散系数型;第三,使用TL内分层的水平涡旋扩散系数型;第四,采用包含有植物冠层垂直结构影响的新干沉积速度型;第五,使用数值扩散小的二阶矩守恒平流格式.用EM3和RADM对OSCAR Ⅳ实例作对比模拟,结果表明EM3能更好地预测出TL内的三维物质分布规律,在同样条件下,两模式预测的PBL低层中心浓度比值均大于1.8,最大值可大于3。     

A SEARCH TECHNIQUE FOR CHOOSING OPTIMUM PARAMETERS IN NUMERICAL WEATHER PREDICTION

The steepest descent(or ascent)search is employed for finding optimum diffusion coefficients in T42L9Gmodel,with a view to improving the model's computational stability or prediction accuracy.The method of thesteepest descent search is first described,in which the golden section search is chosen as the fundamental one-dimensional search used in the multi-dimentional steepest descent search,and then the optimization of the dif-fusion coefficients is described.     

闾阳－海城－东沟剖面深地震波动力学特征初探

本文利用国内外有关深地震波动力学特征的研究方法，结合现有原始资料的实际条件，对闾阳－海城－东沟深地震测深剖面某些地段的主要波组的动力学特征进行了定性或半定量的分析与研究。在频谱分析，速度振幅的计算与相对理论振幅的对比等方面，获得了一些有益信息。为进一步合理解释深地震测深资料提供依据。     

A Eulerian air pollution model for Europe with nonlinear chemistry

A long-range transport model with nonlinear chemical reactions is described. The model contains 35 pollutants and 70 chemical reactions. This is a Eulerian model defined on a space domain containing the whole of Europe. The spherical space domain (corresponding to the Earth's surface covered by the model) is mapped into a square plane domain and discretized by using a 32×32 grid. The grid increments are equidistant (both along the Ox axis and along the Oy axis). The choice of values of the physical parameters involved in the model and the numerical treatment of the model are shortly discussed. The model is tested with meteorological data for 1985 and 1989. The numerical results are compared with measurements at stations located in different European countries. Extensive comparisons of ozone concentrations for July 1985 with measurements taken at 24 European stations are also carried out. Results concerning three episodes in July 1985 as well as results obtained in the study of the sensitivity of the ozone concentrations to variations of NO _x and/or anthropogenic VOC emissions are presented. The advantages and the limitations of such a model are discussed. The model is continuously improved by adding new modules to it. The plans for improvements in the near future are outlined.     

Accuracy of kinematic wave and diffusion wave approximations for space-independent flows on infiltrating surfaces with lateral inflow neglected in the momentum equation

Error equations for the kinematic wave and diffusion wave approximations with lateral inflow neglected in the momentum equation are derived under simplified conditions for space-independent flows. These equations specify error as a function of time in the flow hydrograph. The kinematic wave, diffusion wave and dynamic wave solutions are parameterized through a dimensionless parameter γ which is dependent on the initial conditions. This parameter reflects the effect of initial flow depth, channel-bed slope, lateral inflow, infiltration and channel roughness when the initial condition is non-vanishing; it reflects the effect of bed slope, channel roughness and acceleration due to gravity when the initial condition is vanishing. The error equations are found to be the Riccati equation. The structure of the error equations in the case when the momentum equation neglects lateral inflow is different from that when the lateral inflow is included.     

Importance of advective zone in longitudinal mixing experiments

One-dimensional Fickian dispersion models such as the advection diffusion equation (ADE) are commonly used to analyse and predict concentration distributions downstream of contamination events in watercourses. Such models are only valid once the tracer had entered the equilibrium zone. This paper compares previous theoretical, experimental and numerical estimates of the distance to reach the equilibrium zone with new experimental values, obtained by examining the change of skewness in a tracer profile, downstream of a cross-sectionally well mixed source. Closer agreement was found with Fischers’ theoretical estimate than prior experimental and numerical studies.     

A filtering approach for estimating lake water quality from remote sensing data

In this paper we consider the estimation of lake water quality constituent distributions from hyperspectral remote sensing data. We present a computational approach that can be used to assimilate information from mathematical evolution models into data processing. The method is based on a reduced order iterated extended Kalman filter, and a convection–diffusion model is used to describe the movement of the water quality constituents. The performance of the technique is evaluated in a simulation study. The results show that the filter approach with an appropriate evolution model yields estimates that have better spatial and temporal resolutions than those obtained with conventional methods. Furthermore, the use of a feasible evolution model may make it possible to obtain information also on the concentrations in the lower parts of the lake.     

基于均匀混合假设的天然河道水质模拟方法

扩散质进入窄浅型河道后,通常采用有限差分法求解一维对流分散方程进行水质模拟,计算相对复杂。为简化计算,提出了根据河道断面特征等水力特性,将河道划分为若干单元,忽略不灵敏因子分散项对单元间扩散质交换的影响,并引入单元均匀混合假定,提出了基于均匀混合假设的天然河道水质模拟方法,建立相应的水质数学模型。算例表明,该模型运算稳定,与有限差分法模拟精度相当,但可大大简化并节省计算工作量。     

扩展卡尔曼滤波数据同化试验--基于Lorenz(1960)系统的研究

本文采用Lorenz(1960)系统，在只考虑初始误差及观测误差而不考虑模式误差的情况下，利用扩展卡尔曼滤波(Extended Kalman Filter)数据同化方法进行了数值模拟试验。数值试验的结果表明：扩展卡尔曼滤波数据同化方法对系统状态的估计有较好的改善作用，能有效的抑制估计误差的增长；加大观测频率，可以进一步改善数据同化的效果，使估计误差进一步减小；由于模式误差的存在，系统的不稳定能量会不断的累积，出现了估计误差的异常增长和计算的不连续现象，在模式预报方程中的均值演变方程加人二阶偏差纠错项，对控制估计误差的异常增长，进一步改善数据同化的效果有较明显作用。     

Sediment transport in analogue flume models compared with real‐world sedimentary systems: a new look at scaling evolution of sedimentary systems in a flume

Evolution of sedimentary systems at large temporal and spatial scales cannot be scaled down to laboratory dimensions by conventional hydraulic Froude scaling. Therefore, many researchers question the validity of experiments aiming to simulate this evolution. Yet, it has been shown that laboratory experiments yield stratigraphic responses to allocyclic forcing that are remarkably similar to those in real‐world prototypes, hinting at scale independency with strong dependence on boundary conditions but weak dependence on the actual sediment transport dynamics. This paper addresses the dilemma by contrasting sediment transport rules that apply in the laboratory with those that apply in real‐world geological systems. It is demonstrated that the generation of two‐dimensional stratigraphy in a flume can be simulated numerically by the non‐linear diffusion equation. Sediment transport theory is used to demonstrate that only suspension‐dominated meandering rivers should be simulated with linear diffusion. With increasing grain‐size (coarse sand to gravel) and shallowness of river systems, the prediction of long‐term transport must be simulated by non‐linear, slope‐dependent diffusion to allow for increasing transport rates and thus change in stratigraphic style. To point out these differences in stratigraphic style, three stages in infill of accommodation have been defined here: (i) a start‐up stage, when the system is prograding to base level (e.g. the shelf edge) with no sediment flux beyond the base‐level point; (ii) a fill‐up stage, when the system is further aggrading while progressively more sediment is bypassing base level with the progression of the infill; and (iii) a keep‐up stage, when more than 90% of the input is bypassing the base level and less than 10% is used for filling the accommodation. By plotting the rate of change in flux for various degrees of non‐linearity (varying the exponent in the diffusion equation) it was found that the error between model and real‐world prototype is largest for the suspension‐dominated prototypes, although never more than 30% and only at the beginning of the fill‐up stage. The error reduces to only 10% for the non‐linear sandy‐gravelly and gravelly systems. These results are very encouraging and open up ways to calibrate numerical models of sedimentary system evolution by such experiments.