2002年7～8月海河流域面雨量预报的误差分析

对天津市气象台和河北省气象台在2002年7—8月的海河流域面雨量24h、48h预报进行了误差对比分析,列举了几种预报误差,并利用模糊数学中的综合评判法,计算了它们的模糊评分。结果表明：海河流域的面雨量预报有一定的可信度,把其作为水文模式的初值,有较高的利用价值;对于各气象台各自辖区内的预报,范围越小,预报越准确;在对面雨量预报效果的评估上,80％(或80)这一数值是客观、可信的。     

非线性大气动力学的进展

总结了自新中国成立以来中国科学院大气物理研究所在非线性大气动力学领域的进展,其中主要介绍了在非线性适应过程、非线性稳定与不稳定、全局分析理论、可预报性、低阶谱方法与多平衡态动力学、非线性波动、波流相互作用、阻塞高压非线性动力学、中小尺度非线性动力学等几个方面的研究成果及最新进展.     

流域坡面汇流研究现状述评

对当前坡面汇流计算方法的研究进展进行了较为系统的总结与分析,并对坡面汇流的非线性效应以及城市低影响开发中的雨水入渗与蓄集对坡面汇流的控制作用进行了简要分析。从模型简单实用的角度出发,认为以流域时间-面积关系与线性水库相串联的ModClark法等为代表的概念性分布式坡面汇流模型具有良好的发展前景;考虑到基于等流时单元的变动等流时线法在反映雨强非线性影响中存在的问题,认为根据水文响应单元在不同雨强条件下汇流时间的变化,调整其汇流参数以反映坡面汇流的非线性效应,对于流域坡面汇流的分布式模拟更具有实际意义;针对目前低影响开发设施长时间序列大空间尺度的室外降雨径流监测资料普遍较为缺乏的现状,给出了后期应积极选择合适的技术以加强低影响开发性能监测工作的建议。     

n维非线性波动方程对应的线性方程Cauchy问题解的一个估计式

对于地震波在地球介质中传播的传播非线性而言,在求解n维非线性波动方程Canchy问题时,需要求解其对应的线性波动方程Cauchy问题和建立一些解的估计式。我们已经求得了其对应的线性波动方程Canchy问题解的表达式。在此基础上,本文应用函数空间L^1(R^n)、L^∞（R）的范数,建立了n维非线性波动方程对应的线性波动方程Cancny问题解的一个估计式,为求解n维非线性波动方程莫定了基础。     

基于非线性动力学模型的旅游生态容量分析

近年来,大量增加的游客给景区的旅游生态容量带来了较大的压力,矛盾日渐突出。为了实现景区的环境保护与可持续发展,研究游客活动量与景区旅游生态容量之间的关系是十分必要的。本文从非线性动力学观点出发,构建二者非线性动力学模型,并进行平衡态分析。研究结果表明:(1)系统要达到可持续发展,必须建立在景区环境的自我更新能力及人工恢复程度之和大于外界污染水平,且游客的增长率必须保持在较低水平的基础上;(2)景区旅游生态容量的承载力是有限的,所能容纳的游客也是有限的,超过一定的限度就会引起生态系统失衡。采取严格控制游客量,加强对景区内外环境的保护和改善、提升游客素养等措施,对于系统的协调发展至关重要。     

2007/2008年冬季平流层环流异常及平流层—对流层耦合特征

利用NCEP/NCAR再分析资料,分析了2007/2008年冬季平流层环流振荡的过程和异常特征,以及平流层—对流层的环流耦合特征。结果发现：2007/2008年冬季平流层极涡持续偏强,且极涡从12月下旬开始向东亚—大西洋地区偏心发展,同时在欧洲—地中海地区减弱向极区收缩,从而在欧亚地区的平流层高层最早形成“一脊一槽”1波型的环流形势;该平流层的环流异常信号逐步向低层传播;从平流层10 hPa向下传播到对流层中、低层需要约15~20天的时间;此时对流层500 hPa上,与平流层向下传播的信号相对应,呈现为高纬度“西高东低”和副热带“东高西低”的环流异常型,即分别对应着乌拉尔山阻塞高压、鄂霍茨克海低槽以及副热带高压异常偏强和副热带欧亚大陆中西部低槽活动频繁。这种环流异常配置即形成了有利于我国南方降雪持续的大尺度背景条件。进一步研究还发现,平流层异常信号存在由热带向极传播的趋势。2008年1月的平流层极涡振荡正异常所对应的环流异常信号可以追溯到2007年秋季10月份最早出现在热带地区的环流异常信号,该异常信号在向极传播过程中同时向东移动,2008年1月传播到极区,对应极涡振荡的正异常。说明前期平流层的异常信号及其演变,对2008年1月的气候异常具有重要指示意义。     

Upper bound limit analysis using linear finite elements and non‐linear programming

A new method for computing rigorous upper bounds on the limit loads for one‐, two‐ and three‐dimensional continua is described. The formulation is based on linear finite elements, permits kinematically admissible velocity discontinuities at all interelement boundaries, and furnishes a kinematically admissible velocity field by solving a non‐linear programming problem. In the latter, the objective function corresponds to the dissipated power (which is minimized) and the unknowns are subject to linear equality constraints as well as linear and non‐linear inequality constraints. Provided the yield surface is convex, the optimization problem generated by the upper bound method is also convex and can be solved efficiently by applying a two‐stage, quasi‐Newton scheme to the corresponding Kuhn–Tucker optimality conditions. A key advantage of this strategy is that its iteration count is largely independent of the mesh size. Since the formulation permits non‐linear constraints on the unknowns, no linearization of the yield surface is necessary and the modelling of three‐dimensional geometries presents no special difficulties. The utility of the proposed upper bound method is illustrated by applying it to a number of two‐ and three‐dimensional boundary value problems. For a variety of two‐dimensional cases, the new scheme is up to two orders of magnitude faster than an equivalent linear programming scheme which uses yield surface linearization. Copyright © 2001 John Wiley & Sons, Ltd.     

Evaluation and forecasting of daily groundwater outflow in a small chalky watershed

The purpose of this research is to include expert knowledge as one part of the modelling system and therefore offer the chance to create a productive interaction system between expert, mathematical model (MMO8) and artificial neural networks (ANNs). In the present project, the first objective is to determine some parameters by the MMO8 model, introduced as ANN input parameters to forecast spring outflow. The second objective is first to investigate the effect of temporal information by taking current and past data sets and then to forecast spring outflow. The good results obtained reveal the merit of the ANNs–MMO8 combination, and specifically multilayer perceptron (MLP) models. This methodology, for a network with lower, lag and number hidden layer, consistently produced better performance. Copyright © 2003 John Wiley & Sons, Ltd.     

ANALYSES OF THE DYNAMIC EFFECTS ON WINTER CIRCULATION OF THE TWO MAIN MOUNTAINS IN THE NORTHERN HEMISPHERE——Ⅱ.VERTICAL PROPAGATION OF PLANETARY WAVES

A linear,hemispheric and stationary spectral model with multilayers in the vertical was employed to simulate thevertical propagation of waves triggered by mountains.Results show that,in cooperation with the East Asia zonal meanflow,Tibetan Plateau can excite a strong wavenumber 1 perturbation in the stratosphere with its ridge and trough lo-cated over the Pacific and Atlantic Oceans respectively.On the other hand,the stratospheric wavenumber 1 perturbationcaused by the mechanical forcing of the Rocky Mountains in cooperation with the North America zonal mean flow isvery weak.Calculations from observational data of the vertical profile of critical wavenumber for vertically propagatingwaves imply that the tropospheric wavenumber 1 perturbation can hardly penetrate the North America tropopause up-wards,whereas it can freely propagate through the East Asia tropopause into the stratosphere.Two-dimensional E-Pcross-sections obtained from both observational data and simulated results also demonstrate that waves excited by theRocky Mountains are refracted towards low latitudes in the troposphere during their upward propagation:whereas,inaddition to the above mentioned equatorward leaning branch,the wavenumber 1 and 2 planetary waves excited by theTibetan Plateau possess another branch which is refracted to high latitudes during upward propagation and penetratesthe tropopause into the stratosphere.It is therefore concluded that the difference in the horizontal and vertical wavepropagations in the two hemispheres is a result of the different dynamical forcing induced by the two main mountains inthe Northern Hemisphere.     

Spatial pattern recognition by decomposition

The decomposition problem—the assignment of sample observations to component populations—is studied in a spatial context. The observations are spatially located and the assignment to component populations takes into account the value of each observation as well as the values of neighboring observations. Both parameter estimation and assignment rules use a new method that integrates a standard multivariate decomposition algorithm with nonlinear regression. The method is illustrated and tested with artificial data. The distribution of the trace component Cr₂O₃ in recent Lake Michigan sediments is, then, analyzed by the method. It yields a pattern of component populations that is correlated with the Lake's bottom structure and depositional environments. Study carried out as visiting scholar at the Department of Geological Sciences, Northwestern University, Evanston, Illinois U.S.A.