遥测地表温度与比辐射率的迭代反演方法──理论推导与数值模拟

首先在地表比辐射率为已知的条件下,提出一个非线性迭代温度反演模型,我们对不同的地表和大气条件进行了模拟计算,结果表明当大气温度廓线误差－２Ｋ－２Ｋ,水汽廓线误差±２０％时的温度均方根误差为０．４８Ｋ。当大气模式误差一个模式时反演的温度均方根误差为０．８５Ｋ。在此基础上,引人相邻像元的概念,相邻像元的大气状况可以认为是相同的,应用两个时相的遥感影像数据,假定在两个相近时相之间地表比辐射率值不变,建立地表比辐射率与温度的反演模型。我们对不同的地表和大气条件进行了模拟计算,结果表明当大气温度廓线误差－２Ｋ－２Ｋ,水汽廓线误差±２０％时地表温度均方根误差小于１．５Ｋ,地表比辐射率均方根误差小于０．０２,地表辐射均方根误差为１％;当大气温度廓线误差－２Ｋ－２Ｋ,水汽廓线误差±１０％时,地表温度均方根误差小于１．０Ｋ,地表辐射均方根误差小于０．６％。     

Inversion, error analysis, and validation of GPS/MET occultation data

The global positioning system meteorology (GPS/MET) experiment was the first practical demonstration of global navigation satellite system (GNSS)-based active limb sounding employing the radio occultation technique. This method measures, as principal observable and with millimetric accuracy, the excess phase path (relative to propagation in vacuum) of GNSS-transmitted radio waves caused by refraction during passage through the Earth’s neutral atmosphere and ionosphere in limb geometry. It shows great potential utility for weather and climate system studies in providing an unique combination of global coverage, high vertical resolution and accuracy, long-term stability, and all-weather capability. We first describe our GPS/MET data processing scheme from excess phases via bending angles to the neutral atmospheric parameters refractivity, density, pressure and temperature. Special emphasis is given to ionospheric correction methodology and the inversion of bending angles to refractivities, where we introduce a matrix inversion technique (instead of the usual integral inversion). The matrix technique is shown to lead to identical results as integral inversion but is more directly extendable to inversion by optimal estimation. The quality of GPS/MET-derived profiles is analyzed with an error estimation analysis employing a Monte Carlo technique. We consider statistical errors together with systematic errors due to upper-boundary initialization of the retrieval by a priori bending angles. Perfect initialization and properly smoothed statistical errors allow for better than 1 K temperature retrieval accuracy up to the stratopause. No initialization and statistical errors yield better than 1 K accuracy up to 30 km but less than 3 K accuracy above 40 km. Given imperfect initialization, biases ≫ 2 K propagate down to below 30 km height in unfavorable realistic cases. Furthermore, results of a statistical validation of GPS/MET profiles through comparison with atmospheric analyses of the European Centre for Medium-range Weather Forecasts (ECMWF) are presented. The comparisons indicate the high utility of the occultation data in that very good agreement of upper troposphere/lower stratosphere temperature (better than 1.5 K rms, ≪ 0.5 K bias) is found for a region (Europe+USA) where the ECMWF analyses are known to be good, but poorer agreement for a region (Southern Pacific) where the analyses are known to be degraded.     

On the conception ‘Fair weather condition’ in atmospheric electricity

Measurements of atmospheric electrical and meteorological parameters during different meteorological conditions indicate that the use of the conception Fair weather condition in atmospheric electricity is discussable. Fair weather contains a very broad stability range, from very unstable to strong stable stability of the atmosphere. For turbulent fluctuations of the electric parameters (the most local variations) the variations are determined by the micrometeorological processes for all stability conditions.These fluctuations represent frequencies greater than one period per four minutes. For lower frequencies (less local variations), however, the stability dependency increases. During stable conditions the electric field and vertical current density were nearly wholly influenced by the charges and their transfer in the nearest layer. During near-neutral and unstable conditions the electrical parameters were influenced by more separated sources. Measurements of how well Ohm's law was fulfilled also indicate the difficulties by using the conception fair weather. The measurements also indicate the importance of taking the convection current density into consideration in studies of the electric charge transfer in the atmosphere.In the more large scale of variations measurements of the electric field by radiosoundings show that 88% of the ionospheric potential is derived from the troposphere, where the meteorological processes are of fundamental character for the atmospheric electrical phenomena. The relatively great stability of the diurnal variation of the large scale or global electric field is also valid for the meteorological processes in this scale of variations.     

Intercomparison of New Zealand filter ozonometer and Dobson spectrophotometer total ozone measurements

A commercial version of the prototype New Zealand narrowband interference filter instrument has been run at Garmisch-Partenkirchen since autumn 1977. The vicinity of the Dobson stations Hohenpeissenberg and Arosa permitted a check of the filter instrument's utility, reliability and long-term stability in total ozone measurement by (i) station to station-or short-range intercomparisons of daily means and (ii) a several days' direct intercomparison at Arosa in spring 1978.The regression analysis with the Hohenpeissenberg data-covering the five month autumn-winter 1977/78 period with highly variable ozone-shows that the filter instrument's direct-sun modeX _AD total ozone values were systematically too low by a conversion factor of 0.93. The excellent proportionality between the instruments is indicated by the standard deviation of only 1,3%. The (–7.0±1.3%) bias perfectly agrees with the result reported byBasher [1] for a similar, yet direct, intercomparison.During the direct Arosa intercomparison the filter instrument also tracked down short-time variations in a remarkable, though biased, parallelism with the Dobson instruments, and theX _AD daily means fit well into the prognostic relationship derived from the preceding autumn-winter short-range intercomparison.A creeping deterioration of the characteristics of the 305.5 nm double filter then prevented any furtherX _AD total ozone intercomparison. No final decision is possible whether this behaviour after a half year's run is typical or just an unusual exception. The otherwise excellent long-term stability is demonstrated by the comparison of theX _CD daily means with the HohenpeissenbergX _AD data for a 1 1/2 year's period (autumn 1977–spring 1979), resulting in the same bias relationship as for theX _AD total ozone measurements.     

Measurement of lower-stratospheric aerosol by a balloon-borne photometer

Height distribution of the stratospheric aerosol extinction coefficient was measured in the altitude range 10 to 20 km by a balloon-borne multi-color sunphotometer in May 1978. It is demonstrated that detailed structures of the distribution of stratospheric aerosol can be remotely measured by the solar occultation method as well as by lidar andin situ particle counter observations. In the aerosol layer appearing at 18 km altitude the extinction coefficient at 800–1000 nm wavelength reached to 3×10^–7 m^–1, which was reasonable compared with lidar observations. Wavelength dependence of the aerosol optical depth was crudely estimated to be proportional to ^–1.5.     

The characteristics of rain electricity in Nigeria. I — Magnitudes and variations

The characteristics of rain and point charges based on routine measurements extending over four rainy seasons are presented. An average rain current density of (1.0±0.1)×10^–10 A m^–2 and charge per unit volume of rain water of (0.43±0.02)×10^–4 C m^–3 for the locality are obtained, which are compared with data obtained elsewhere by other workers. The point-discharge current measurements lead to a revised estimate of (0.86±0.08)×10^–9 A m^–2 for the average point discharge current below storm clouds.     

Equations of kinetic and available potential energy evolution in wave-number frequency space

The equations of kinetic and available potential energy for the analysis of the growth and decay of atmospheric waves in wave-number-frequency space are derived. The effects of linear and nonlinear interactions of waves in the velocity and temperature fields, the conversions between the available potential and kinetic energies, and the Reynolds and molecular stresses on wave growth and decay are discussed.     

Global distribution of gaseous mercury in the troposphere

Atmospheric mercury concentrations were measured during a nautical expedition on the Atlantic Ocean between Hamburg (54°N, 10°E) and Santo Domingo (20°N, 67°W). In addition, samples were taken during flights on a commerical aircraft in the upper and middle troposphere between 60°N and 55°S, mostly over the Pacific Ocean. The data obtained in the lower troposphere over the Northern Atlantic show considerable variation in the Hg concentrations, with values ranging between 1 and 11 ng/m³; the average concentration was found to be 2.8 ng/m³. The upper tropospheric data show an interhemispheric difference with average values of 1.45 ng/m³ and 1.08 ng/m³ in the Northern and Southern Hemisphere, respectively. This suggests that mercury production occurs predominantly over the continents both by natural and anthropogenic processes. The mercury content in aerosols was found to be 0.3 ng/m³, or one-tenth of the atmospheric concentration. The data indicate a mean residence time of mercury in the atmosphere of a few months to one year.     

大气重力信号的理论计算及其检测

基于标准大气定律和大气圆柱体分布模型,本文引进了大气重力格林函数,用离散格积方法求得了大气对重力场观测的影响,对台站高程、周围地形和地表温度变化等因素的影响问题进行了讨论．结果说明台站近区气压变化是大气重力信号的主要贡献者,考虑大气质量负荷引起的弹性地球形变效应后,对距台站0．5°的区域积分获得的大气重力导纳值为-0．3603μGal／hPa,占全球大气变化引起的总信号的90％以上,这一理论模型结果与超导重力仪实测结果相吻合,并能较有效地用于消除重力观测中的气压干扰成分．