青岛上空中层大气密度和温度的激光雷达探测

介绍了中国电波传播研究所瑞利散射激光雷达系统的结构和性能, 阐述了激光雷达探测中层大气密度和温度的工作原理, 给出了青岛地区中层大气密度和温度的初步探测结果. 通过与卫星、探空气球和大气模式数据的结果对比, 验证了激光雷达探测大气温度的可靠性. 基于2008-2009两年的观测, 获得了青岛地区上空中层大气温度的季节变化和平均分布. 激光雷达观测结果表明, 青岛地区平流层温度比CIRA86模式结果高, 且二者偏差呈夏秋季小、冬春季大的特点, 中间层温度则正好相反.      

瑞利激光雷达反演中层大气温度算法

瑞利激光雷达是探测中层大气温度的重要设备,具有高时空分辨率且不存在探测盲区的优点.为充分利用激光雷达探测到的原始数据,改进了传统的Chanin-Haunchecorne方法,采用由均匀搜索生成温度初值的方式反演大气温度.位于北京延庆（40.3°N,116.2°E）的瑞利激光雷达具有589nm和532nm两个通道,将589nm通道用于计算,532nm通道作为参考.在150~250K区间等间隔选取多个温度初值,利用589nm通道反演60~70km高度范围内的大气温度廓线和大气密度廓线,利用参考密度廓线选取准确的温度初值,反演得到准确的大气温度廓线.将589nm通道和532nm通道反演的温度廓线相比较,发现二者具有较高的一致性.改进后的方法有效利用了信噪比较差的光信号,使589nm通道的温度探测上限从60km提高到70km.      

无线电掩星和激光雷达观测结果比较

简要阐述了激光雷达和ＧＰＳ掩星探测大气的基本原理，并利用１９９６年１０月１４日１６２５ ＵＴ时 Ｂ本 ＴＭＵ钠和瑞利散射激光雷Ｂ（３９．２３°Ｅ、 ３５．３７°Ｎ）资＃，反演得到 ３０－５５ ｋｍ高度范围内大气密度和温度的观测结果，与此同时通过激光雷达上空的 ＧＰＳ掩星观测到的掩星切点处的观测结果进行了比较．结果表明：两种方法获得测出的密度廓线符合得很好；温度廓线的变化趋势也基本一致，但在某些高度上有较大的偏离．文中最后对比较结果作了简要分析，讨论了重力波对大气密度和温度廓线的影响．     

基于最优估计法的瑞利激光雷达反演大气温度研究

基于瑞利激光雷达的回波光子信号对中层大气进行探测，结合最优估计法，对大气温度进行反演。本文基于瑞利散射激光雷达方程建立正向模型，选择大气模型的温度廓线作为先验状态信息，构建用于最优化处理的成本函数，利用Levenberg-Marquardt最优化算法对成本函数执行最优化处理，得到大气温度的反演结果，对反演结果的不确定度分析的同时利用平均核矩阵对反演结果中真实信息的贡献进行评估。利用瑞利激光雷达方程产生的模拟回波信号进行了大气温度的反演处理与分析，对中国科学院国家空间中心提供的瑞利激光雷达实测数据进行大气温度的最优估计反演。结果表明，90 km以下的反演不确定度在10 K以内，且相较于CH方法，最优估计法具有反演有效范围高的优势；在回波光子信噪比较高的区域，反演不确定度较小，且真实信息对反演结果的贡献占主导地位。     

我国火箭探测资料与CIRA 1986稿的比较

本文用我国气象火箭探测获得的中层大气风、温度、密度和气压等数据,与国际参考大气模式CIRA 1986稿作了比较.结果看到,我国火箭探测温度数据值在模式给出的平均值剖面上下波动,符合的程度在国际比对偏差范围之内;用探测数据分析得到的风的季节性和声重波强度等与模式给出的描述都很符合.从而说明,用CIRA 1986稿来表征我国上空的中层大气结构是合适的.      

TIMED卫星探测的全球大气温度分布及其与经验模式的比较

利用TIMED卫星遥感探测的全球温度分布与NRLMSISE-00大气经验模式进行了对比研究.研究表明,在中间层下部以下的高度范围内,经验模式与卫星探测的大气温度分布有很好的一致性.但是比较发现,在中层顶区域,经验模式的计算结果与实测结果有较大的差异.卫星探测表明,在春分季节的低纬地区中层顶区存在稳定的逆温层,但是经验模式不能给出低纬地区春分季节中间层逆温层的分布特征.卫星观测表明在全球范围内中层顶有两个非常不同的高度,一个处于100km附近,另一个处于85km附近,但是经验模式不能给出这一中层顶高度的分布特征.同时在低热层,经验模式计算的温度分布与卫星遥感的探测结果有很大的差异.      

武汉中层、低热层大气角谱中频雷达观测

武汉中频雷达是利用中层、低热层中电子密度不均匀体的散射来测量大气的水平风场和电子密度剖面。雷达在计算风场的过程中可得到一些该层大气中电子密度不均匀体的寿命和空间尺度等参量值。结合这些参量和大气风场值可计算得到大气的角谱。其计算方法包括全相关分析技术的谱宽法和空间相关法。前者计算的值被认为是大气角谱的上限值。应用这两种方法,利用2001年2月9日武汉中频雷达的观测数据,对武昌上空中层、低热层大气的角谱进行了计算。结果得到大气角谱随高度增加略微增加,如在68km为6°,90km达最大为11°,其平均值为9.2°。利用2001年2月4—10日86km高度上的数据,得到一个7天平均的日变化曲线,发现大气角谱值在本地夜晚最小。     

廊坊中频雷达流星观测及初步结果

中频雷达用来开展夜间100km高度以上的流星观测,获得流星随时间、高度、方位的分布情况及流星体速度、流星辐射点、流星余迹径向速度等参数,其探测数据可用于流星天文学、中层大气动力学等领域的研究.利用2017年11月16日12:00UT-22:00UT期间廊坊观测站（39.4°N,116.7°E）的中频雷达数据,首次开展了中国中纬度地区夜间流星观测实验,共检测到94个流星回波信号,集中分布在97～115km高度范围内,平均高度为106.5km,计算得到了流星回波的双极扩散系数、方位分布等相关参数,并与国外中频雷达流星探测结果进行了初步比较.      

基于HITRAN数据库的大气激光雷达信号仿真

信号仿真研究是大气激光雷达研究的重要环节，结果可为激光雷达系统的设计和研制提供基础.本文仿真系统可用于激光雷达回波信号模拟.该仿真系统利用HITRAN数据库中的吸收谱线计算大气分子对激光的吸收光谱，并采用大气辐射传输模型中的气溶胶模式模拟气溶胶对激光的衰减.利用激光雷达方程，数值模拟了355nm，532nm以及1064nm的回波信号，并将532nm的数值模拟结果与实测结果进行对比，评估了激光雷达系统的光学效率.      

瑞利激光雷达探测北京上空中间层低逆温层

利用延庆瑞利激光雷达（40.47°N,115.97°E）2012年1-2月及2012年5月至2013年4月的探测数据,分析得到北京上空60～80km高度140个晚间的温度廓线,对这一区域内的低中间层逆温层现象（Lower Mesospheric Inversion Layer,Lower MIL）进行了统计分析,发现平均逆温幅度为23.4K,平均垂直尺度为4.78km,逆温层底部平均高度为68.2km.约有2/3的逆温层存在随时间垂直传播现象,且大多为向下传播.此外,还观测到三个垂直传播速度相差近一倍的特殊双低层MIL演化现象.      

New input to IRI from the worldwide digisonde network

Electron density profiles derived from Digisonde ionograms at Argentia, Millstone Hill, Wallops Island, Bermuda, Dourbes and Karachi are compared with IRI model prediction. Four months of data for 1989/90 were analyzed. For a number of station/months N(h) profiles were available every 15 or 30 minutes providing a good statistical database for the evaluation of the IRI model in terms of diurnal and seasonal variations. The data presented here are part of the VIM study (Validation of Ionospheric Models) initiated by the URSI Working Group G3 on Ionospheric Informatics.     

Modelling bottom and topside electron density and TEC with profile data from topside ionograms

The paper describes the technique that has been implemented to model the electron density distribution above and below the F2 peak making use of only the profiles obtained from the INTERCOSMOS-19 topside ionograms. Each single profile from the satellite height to the ionosphere peak has been fitted by a semi-Epstein layer function of the type used in the DGR model with shape factor variable with altitude. The topside above the satellite height has been extrapolated to match given values of plasmaspheric electron densities to obtain the full topside profile. The bottomside electron density has been calculated by using the maximum electron density and its altitude estimated from the topside ionogram as input for a modified version of the DGR derived profiler that uses model values for the foF1 and foE layers of the ionosphere. Total electron content has also been calculated. Longitudinal cross sections of vertical profiles from latitudes 50° N to 50° S latitude are shown for low and high geomagnetic activity. These cross sections indicate the equatorial anomaly effect and the changes of the shape of low latitude topside ionosphere during geomagnetic active periods. These results and the potentiality of the technique are discussed.     

冕洞中的温度分布与阿尔芬涨落能流

基于两个假设即冕洞中的阿尔芬涨落无阻尼传播和日冕温度范围为9×103K—2.5×106K, 首先导出了快发散流管中阿尔芬涨落无阻尼传播特征;然后利用Munro-Jackson观测结果, 通过数值计算, 确定了可接受的冕洞温度分布;同时确定了冕洞中可能存在的冕底阿尔芬涨落能流为5×105—1×106ergs cm-2sec-1。分析表明, 这样大小的能流对加速冕洞等离子体成为高速风流是有效的和足够的。      

Tidal associated temperature disturbances observed at the middle atmosphere (30–65 km) by a Rayleigh lidar at 23°S

A lidar has been operated in São José dos Campos, Brazil (23.2°S, 45.8°W) since 1972, mainly dedicated to the study of mesospheric sodium at the 589 nm resonant line. The molecular Rayleigh scattering can also be used provided we limit the height to ∼75 km where the sodium scattering begins. Nevertheless, the weak signal obtained only permits the determination of density and temperature profiles by accumulating a large number of shots giving only nocturnal average profiles. Temporal variations in density and temperature on the scale of hours can however, be obtained by performing a superposed epoch analysis for a given time interval and covering a period of several days. In this way we obtained hourly mean profiles grouped by months, seasons and overall, with data acquired from 1993 to 2004. The difference between the hourly temperatures and the nocturnal means shows for some months, with enough data coverage, downward propagating structures that apparently have tidal origin. The seasonal averages show a recurrent feature with high temperatures before and low temperatures after midnight above 50 km. Some similarity is found with the GSWM model, but the observed temperature amplitudes are twice of that for the model.     

Comparing the improved Di Giovanni/Radicella model with sounding-based electron density profiles and with the IRI model

The Di Giovanni/Radicella model (DGR) /1/ determines a bottom side electron densty profile alone from the set of routinely scaled ionogram parameters foE, foF1, foF2 and M(3000)F2 and the total electron content; the smoothed sunspot number R12 appears in the calculation. Present designations are DGR2/2/ and DRR3 /3/ [see Appendix]; they are valid in the northern hemisphere. DGR is compared with electron density profiles derived from ionograms obtained at Juliusruh (54.6°N, 13.4°E), and with the (URSI-based) IRI90 at different conditiones. Experimental total electron content (TEC) data are compared to both models. At the considered station, the profiles obtained by both models are reasonably in agreement amongst themselves and with the experimental data.

The TEC derived from the DGR3 model is in good agreement with experimental TEC, whereas, at high solar activity, IRI90 gives too high TEC values, especially during daytime.     

Accuracy of IRI profiles of ionospheric density and temperatures derived from comparisons to Kharkov incoherent scatter radar measurements

The incoherent scatter radar (ISR) facility in Kharkov, Ukraine (49.6°N, 36.3°E) measures vertical profiles of electron density, electron and ion temperature, and ion composition of the ionospheric plasma up to 1100 km altitude. Acquired measurements constitute an accurate ionospheric reference dataset for validation of the variety of models and alternative measurement techniques. We describe preliminary results of comparing the Kharkov ISR profiles to the international reference ionosphere (IRI), an empirical model recognized for its reliable representation of the monthly-median climatology of the density and temperature profiles during quiet-time conditions, with certain extensions to the storm times. We limited our comparison to only quiet geomagnetic conditions during the autumnal equinoxes of 2007 and 2008. Overall, we observe good qualitative agreement between model and data both in time and with altitude. Magnitude-wise, the measured and modeled electron density and plasma temperatures profiles appear different. We discovered that representation accuracy improves significantly when IRI is driven by observed-averaged values of the solar activity index rather than their predictions. This result motivated us to study IRI performance throughout protracted solar minimum of the 24th cycle. The paper summarizes our observations and recommendations for optimal use of the IRI.     

多种观测数据驱动的三维行星际太阳风MHD模拟

三维磁流体力学（MHD）数值模拟是行星际太阳风研究的重要手段.本文发展了一种由多种观测数据驱动的三维行星际太阳风MHD数值模型.模型的计算区域为0.1AU到1AU附近，使用Lax-Friedrich差分格式在六片网格系统中进行数值求解.边界条件中磁场使用GONG台站观测的光球磁图外推获得，密度通过LASCO观测的白光偏振亮度反演得到，速度根据以上两种观测数据并利用一种基于人工神经网络技术（ANN）的方法得到，温度通过自洽方法根据磁场和密度导出.利用该模型模拟了第2062卡灵顿周（CR2062）时期的行星际太阳风，模拟结果显示出丰富的观测特征，并与OMNI以及Ulysses的实际观测值符合得较好.该模型可用于提供接近真实的行星际太阳风，有助于提高空间天气预报的精度.      

不同上边界条件下的极区电离层数值模拟

利用一维自洽的极区电离层模型,研究了沿磁力线方向不同电离层-磁层耦合条件下极区电离层的响应.此模型在110-610km的电离层空间区域内,综合求解描述极区电离层的连续性方程、动量方程和能量方程,以得到电离层数值解.研究发现,上边界条件在200 km以上的高度能显著地影响电离层参量的形态.较高的O+上行速度对应较低的F层峰值和较高的电子温度.不同边界O+上行速度对应的温度高度剖面完全不同.200km以上电子温度高度剖面不但由来自磁层的热流通量所控制,同时还受到场向O+速度的影响.对利用电离层模型研究电离层内部物理过程提出了建议.      

Reconstruction of topside density profile by using the topside sounder model profiler and digisonde data

To improve the accuracy of the real time topside electron density profiles given by the Digisonde software a new model-assisted technique is used. This technique uses the Topside Sounder Model (TSM), which provides the plasma scale height (H_s), O⁺–H⁺ transition height (H_T), and their ratio Rt = H_s/H_T, derived from topside sounder data of Alouette and ISIS satellites. The Topside Sounder Model Profiler (TSMP) incorporates TSM and uses the model quantities as anchor points in construction of topside density (Ne) profiles. For any particular location, TSMP calculates topside Ne profiles by specifying the values of foF2 and hmF2. In the present version, TSMP takes the F2 peak characteristics – foF2, hmF2, and the scale height at hmF2 – from the Digisonde measurements. The paper shows results for the Digisonde stations Athens and Juliusruh. It is found that the topside scale height used in Digisonde reconstruction is less than that extracted from topside sounder profiles. Rough comparison of their bulk distributions showed that they differ by an average factor of 1.25 for locations of Athens and Juliusruh. When the Digisonde scale heights are adjusted by this factor, the reconstructed topside profiles are close to those provided by TSM. Compared with CHAMP reconstruction profiles in two cases, TSMP/Digisonde profiles show lower density between 400 and 2000 km.