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.     

Two different methods of geoidal height determinations using a spherical harmonic representation of the geopotential, topographic corrections and the height anomaly–geoidal height difference

 It is suggested that a spherical harmonic representation of the geoidal heights using global Earth gravity models (EGM) might be accurate enough for many applications, although we know that some short-wavelength signals are missing in a potential coefficient model. A `direct' method of geoidal height determination from a global Earth gravity model coefficient alone and an `indirect' approach of geoidal height determination through height anomaly computed from a global gravity model are investigated. In both methods, suitable correction terms are applied. The results of computations in two test areas show that the direct and indirect approaches of geoid height determination yield good agreement with the classical gravimetric geoidal heights which are determined from Stokes' formula. Surprisingly, the results of the indirect method of geoidal height determination yield better agreement with the global positioning system (GPS)-levelling derived geoid heights, which are used to demonstrate such improvements, than the results of gravimetric geoid heights at to the same GPS stations. It has been demonstrated that the application of correction terms in both methods improves the agreement of geoidal heights at GPS-levelling stations. It is also found that the correction terms in the direct method of geoidal height determination are mostly similar to the correction terms used for the indirect determination of geoidal heights from height anomalies. Received: 26 July 2001 / Accepted: 21 February 2002     

Measured And Simulated Latent And Sensible Heat Fluxes At Two Marine Sites In The Baltic Sea

In this study, turbulent heat flux data from two sites within the Baltic Sea are compared with estimates from two models. The main focus is on the latent heat flux. The measuring sites are located on small islands close to the islands of Bornholm and Gotland. Both sites have a wide wind direction sector with undisturbed over-water fetch. Mean parameters and direct fluxes were measured on masts during May to December 1998.The two models used in this study are the regional-scale atmospheric model HIRLAM and the ocean model PROBE-Baltic. It is shown that both models overestimate the sensible and latent heat fluxes. The overestimation can, to a large extent, be explained by errors in the air-water temperature and humidity differences. From comparing observed and modelled data, the estimated 8-month mean errors in temperature and humidity are up to 1 °C and 1 g kg^-1, respectively. The mean errors in the sensible and latent heat fluxes for the same period are approximately 15 and 30 W m^-2, respectively.Bulk transfer coefficients used for calculating heat and humidity fluxes at the surface were shown to agree rather well with the measurements, at least for the unstable data. For stable stratification, the scatter in data is generally large, and it appears that the bulk formulation chosen overestimates turbulent heat fluxes.     

扁椭球形雨滴和冰粒后向散射实验与理论研究

本文介绍了对扁椭球形雨滴和冰（雹）粒散射能力的实验模拟和测量技术，报告了一个系列扁水(冰)椭球后向散射截面的实验测量结果；并用扩张边界条件法(EBCM)计算出相应的理论值，用Mie公式计算出等容积球形后向散射截面。将实验结果与理论计算作了综合对比分析，得到了若干有用的结论。     

体应变观测中的气压干扰机制和排除方法的研究

建立体应变理论气压干扰模型,通过实际观测资料研究了气压对体应变月变和日变的干扰过程,探讨消除干扰的方法,且和水位资料作了对比分析。     

大型抽水试验最优化设计研究

抽水试验最优化设计,是指针对一个具体的地下水系统优选出一个抽水试验设计,它使抽水试验的费用最小化,而试验所获资料能满足识别参数所要求的可靠性,它可以概化为一个非线性的混合整数规划问题,其决策变量是确定抽水井和观测井的位置及井数、抽水量和观测频率等,采用搜索算法求解了这一优化问题,得到了理论最优解,并进行了实例演算。     

Sensitivity of GPS and GLONASS orbits with respect to resonant geopotential parameters

The Center for Orbit Determination in Europe (CODE) has been involved in the processing of combined GPS/GLONASS data during the International GLONASS Experiment (IGEX). The resulting precise orbits were analyzed using the program SORBDT. Introducing one satellites positions as pseudo-observations, the program is capable of fitting orbital arcs through these positions using an orbit improvement procedure based on the numerical integration of the satellites orbit and its partial derivative with respect to the orbit parameters. For this study, the program was enhanced to estimate selected parameters of the Earths gravity field. The orbital periods of the GPS satellites are —in contrast to those of the GLONASS satellites – 2:1 commensurable (P _Sid:P _GPS) with the rotation period of the Earth. Therefore, resonance effects of the satellite motion with terms of the geopotential occur and they influence the estimation of these parameters. A sensitivity study of the GPS and GLONASS orbits with respect to the geopotential coefficients reveals that the correlations between different geopotential coefficients and the correlations of geopotential coefficients with other orbit parameters, in particular with solar radiation pressure parameters, are the crucial issues in this context. The estimation of the resonant geopotential terms is, in the case of GPS, hindered by correlations with the simultaneously estimated radiation pressure parameters. In the GLONASS case, arc lengths of several days allow the decorrelation of the two parameter types. The formal errors of the estimates based on the GLONASS orbits are a factor of 5 to 10 smaller for all resonant terms. AcknowledgmentsThe authors would like to thank all the organizations involved in the IGS and the IGEX campaign, in particular those operating an IGS or IGEX observation site and providing the indispensable data for precise orbit determination.     

Coda Q in the Kachchh Basin, Western India Using Aftershocks of the Bhuj Earthquake of January 26, 2001

Q_C-estimates of Kachchh Basin in western India have been obtained in a high frequency range from 1.5 to 24.0 Hz using the aftershock data of Bhuj earthquake of January 26, 2001 recorded within an epicentral distance of 80 km. The decay of coda waves of 30 sec window from 186 seismograms has been analysed in four lapse time windows, adopting the single backscattering model. The study shows that Q_c is a function of frequency and increases as frequency increases. The frequency dependent Q_c relations obtained for four lapse-time windows are: Q_c=82 f^1.17 (20–50 sec), Q_c=106 f^1.11 (30–60 sec), Q_c=126f^1.03 (40–70 sec) and Q_c=122f^1.02 (50–80 sec). These empirical relations represent the average attenuation properties of a zone covering the surface area of about 11,000, 20,000, 28,000 and 38,000 square km and a depth extent of about 60, 80, 95, 110 km, respectively. With increasing window length, the degree of frequency dependence, n, decreases marginally from 1.17 to 1.02, whereas Q₀ increases significantly from 82 to 122. At lower frequencies up to 6 Hz, Q_c⁻¹ of Kachchh Basin is in agreement with other regions of the world, whereas at higher frequencies from 12 to 24 Hz it is found to be low.     

基于抽水试验的土体渗透性耦合响应模型研究

土体中大量降排水的过程实际上是土体中有效应力增长的过程。通过对抽水试验中土体渗透系数的回归分析,建立了渗透系数与有效应力增量之间的耦合响应模型关系,说明降排水过程中土体渗透系数是呈非线性变化的。