VICARIOUS RADIOMETRIC CALIBRATION OF SATELLITE FY-ID SENSORS AT VISIBLE AND NEAR INFRARED CHANNELS

For the problem of in-flight test site radiometric calibration for the FY-1D Meteorological Satellite onboard sensor's visible and near infrared channels,this paper described the calibration method,satellite-ground synchronous observation data acquired at China Dunhuang Calibration Test Site,parameter derivation and radiative transfer computation during the calibration,and the calibration result analysis.At FY-1D seven channels (with central wavelength at Channels 1:630 nm,2:865 nm;6:1610 nm;7:455 nm,8:505 nm;9:555 nm;10:932 nm),calibration coefficients obtained during the ground site calibration were compared with that respectively of prelaunch calibration.It is demonstrated that results of FY-1D onboard two sensors' (A and B) at channels 1,2,6 and 10 were close to pre-calibration,and the absolute difference of Gobi desert reflectance computed using test site and pre-launch calibration coefficients was no more than 2%. At other channels,large pre-launch calibration errors resulted in a poor consistency between the test site and pre-launch calibration.The errors can be corrected by the test site calibration results. Based on a rough estimation,the overall error of the calibration was about 6%. The paper also presented the in-flight vicarious calibration at the visible and near infrared channels of FY-1C sensor A which was launched in 1999 and has been put into operational mode since 2002.The results exhibit that FY-1C sensor's response has 23% attenuation at Channels 7 and 8,while only minor degradation at the other channels was found. During the mission,calibrations were also conducted at NOAA-17's Channels 1 and 2 (1:430-830 nm,2:500-1072 nm).A very good consistency has been achieved between the test site and pre-launched calibration results.     

多窗谱分析方法及其在全球变暖研究中的应用

多窗谱分析方法是一种低方差、高分辨的谱分析方法 ,尤其适合于非线性气候系统中高噪声背景下弱信号、时频演变信号的诊断分析。文中简要地介绍了多窗谱分析方法的基本原理 ,并将其用于近 1 5 0 a来全球及南、北半球温度变化的研究。结果表明 :(1 )在温度呈线性变暖的背景趋势下 ,北半球、全球年平均温度还具有显著的 40～ 70 a的准周期低频振荡 ,它们与变暖过程中的波动性密切相关 ;(2 )年代际振荡信号在北半球温度序列中也相当明显 ,南半球则存在多种低于 5 a尺度的显著振荡周期 ,半球或全球 ENSO振荡分量、QBO周期信号的振幅都具有缓变包络的特性 ;(3 )与传统功率谱估计、最大熵谱估计结果的对比表明 ,多窗谱分析方法方法得到的谱估计分辨力高 ,稳定性强。     

大气谱模式中质量守恒格式的研究

当前在世界上广泛采用的谱模式，是以地面气压ps取对数作为展谱变量，模式大气总质量不守恒。本文研究这样的问题，提出了以热力学变量比容α代替温度变量T的方法，并以地面气压ps本身作为展谱变量，这样可以保持大气总质量守恒。同时，由于以α作展谱变量，可使重要项气压梯度项为二阶变量组成，截谱后，不会出现混淆误差。计算结果表明，这样改变后，大气总质量保持守恒。同时消除了混淆误差的出现。     

实际预报可预报性的时空依赖性分析

本文利用国家气象中心1990年1月1日至1992年2月29日的1～5 d、500 hPa高度场业务预报结果研究了不同球谐函数谱分量和经验正交函数(EOF)谱分量的可预报性。对球谐函数谱分量的研究表明可预报性并不是随着空间尺度的增加呈现单调的递减关系，主要表现为在总波数n等于5附近具有最大可预报性。可预报性主要依赖于总波数n，经卷大圆上的波数和纬向圆周上的波数对可预报性具有完全相同的重要意义。对EOF谱分量的分析表明，随着EOF指数(即序号数)的增大，可预报性依次减小。从EOF和球谐函数的关系发现前几个EOF分量正是对应着最可预报的球谐函数分量。进一步分析表明，最可预报的分量对应大气运动的慢变流型。这些结果对如何使用数值预报产品以及如何进行延伸预报具有重要意义。     

利用不同方法反演冬小麦叶面积指数研究

基于高光谱数据的叶面积指数监测是快速获取冬小麦叶面积指数的重要方法。为了探究回归方法和高光谱数据变换对冬小麦叶面积指数反演精度的影响,采用逐步回归和偏最小二乘回归方法,分别建立基于冬小麦拔节期冠层高光谱数据、一阶导数光谱数据、二阶导数光谱数据和对数光谱数据的叶面积指数多元线性回归模型。结果显示,导数和对数变换能够提高冬小麦LAI反演精度,以蓝紫光、绿光、红光和近红外波段建立的一阶导数光谱数据逐步回归模型最优,建立回归模型的决定系数R2为0.974,交叉验证的RMSE为0.131,可为冬小麦LAI估算的方法选择和数据处理提供依据和参考。     

高光谱大气红外探测仪(AIRS)反演大气不稳定度指数在强对流天气个例中的应用试验

使用探空、NCEP-FNL,数据和高光谱分辨率大气垂直探测仪(AIRS)标准反演数据计算大气不稳定度指数,对2011年6月23日北京强对流天气发生前的本地及上游大气中不稳定能量进行分析研究。分析发现:利用08、14和20时探空数据计算的北京站不稳定度指数显示了在"6·23暴雨"过程发生前后北京上空不稳定能量变化,上游关键区无探空数据;利用NCEP和AIRS数据计算的不稳定度指数显示,强对流天气发生前,在北京的上游关键区大气处于极端不稳定状态(K指数大于40,SI指数小于一5),有利于强对流天气发生。文章的研究结果表明,探空数据时空分辨率较低,不利于监测强对流天气的发生;质量控制后AIRS数据计算的不稳定度指数可以监测对流天气的发生;空间分辨率较低的NCEP数据监测小范围大气不稳定层结能力较低。综上所述,AIRS反演产品具有弥补探空资料时空分辨率不足的优势,利用AIRS L2反演产品计算晴空大气不稳定度指数产品可以监测到"6·23暴雨"天气发生前上游关键区大气层结稳定度状态,为预报员决策提供有效的辅助信息。     

Distribution of L_g coda Q in the Chinese continent and its adjacent regions

In continent Lg is usually one of the predominant phases recorded by short-period or broad-band seismometers. A ray-theoretical approach shows that Lg wave is the superposition of higher-mode surface waves propagating in the continental crust[1—4]. The g…     

Analysis of a large extent aeromagnetic survey near the geomagnetic equator (Minas Gerais,Brazil)

A large extent aeromagnetic survey, covering a 9°×9° degree of a Precambrian shield area in SE Brazil, was used to analyse long wavelength magnetic anomalies with the purpose of determining the magnetic structure of the lithosphere.Data analysis based on a novel approach to the estimation of the parameters in the statistical model ofSpector andGrant (1970), reveals a highly magnetized lower crust, bounded by an essentially non-magnetic upper mantle. The highly magnetized layer can be associated with a higher density lower crust model, derived from gravity and seismic data available in the area, suggesting a compositional transition from upper to lower crust.It was also verified that the operation of reduction to the pole, which is singular at the geomagnetic equator, performs quite satisfactorily at the low geomagnetic inclination of the area, bringing into accordance the long wavelength anomalies with the structural trends.     

Upper mantle seismic anisotropy beneath a convergent boundary: SKS waveform modeling in central Tibet

This study used SKS waveforms from the International Deep Profiling of Tibet and the Himalayas (INDEPTH) III dataset and a new 2D method for modeling seismic waves in anisotropic media to construct an image of anisotropic structures beneath central Tibet. A preferred model revealed three-segment anisotropic structures in the upper mantle beneath the study region. Waveform modeling demonstrated that the anisotropy was mainly generated by the lithosphere but not the asthenosphere, and that an anisotropic model with a flatter axis of symmetry provides a more consistent interpretation of the observations than models having steeply dipping symmetry axes. A relatively low velocity zone may underlie or intermingle with the anisotropic structures in the northern portion of the region. Synthetic tests also indicate that variations in the elastic constants and depth extent of the anisotropy assumed by the calculations do not affect the general conclusions, although trade-offs exist among certain model parameters. The modeling results suggest that the complex seismic structures in central Tibet were associated with underthrusting of the Indian lithosphere beneath the Asian lithosphere; the inferred flat symmetry axis of the anisotropy was likely generated during this collision process. If this were not the case, the inherited anisotropy would exhibit a steeply dipping axis of symmetry, parallel to the direction of underthrusting.