基于COSMIC的平流层重力波参数分析

基于COSMIC卫星观测的2006年12月29日到2008年1月3日30°-40°N纬度内的温度剖面,分别利用垂直滑动窗、双滤波器和单滤波器三种方法计算低平流层重力波的扰动和势能,获得重力波扰动和势能随高度、经度的分布以及多时间尺度变化特性,分析重力波扰动势能与背景温度及风场的变化趋势和特点.比较三种方法得到的结果发现:垂直滑动窗方法只能去除大垂直尺度的背景,无法抑制小尺度的扰动,其得到的结果误差较大;双滤波器法对温度剖面中的大尺度背景和小尺度扰动都能很好地抑制;单滤波器法得到的重力波扰动中基本不包含垂直方向的大尺度背景,但是包含一些小垂直尺度的扰动.因此,对于垂直波长为10km左右的重力波,采用双滤波器法合适;如果需要得到小尺度重力波的变化特性,采用单滤波器法合适.采用双滤波器法无法得到势能随高度的变化,而采用单滤波器法能够给出每月势能随高度的分布.对30°-40°N纬度内的重力波参数进行统计分析得到重力波扰动、势能与背景温度和水平风场的关系.      

Seasonal variation of gravity wave sources from satellite observation

Recent high-resolution satellite observations of gravity waves in the middle atmosphere have shown correlations with the strength of the stratospheric jet stream, surface topography, and tropical convection. Seasonal variations of wave-induced stratospheric radiance variances are often the manifestations of modulations of these sources and refractive influences. In this paper, we focus on the seasonal climatology of gravity waves observed by the UARS MLS, while also showing some new results from GPS and AMSU instruments. Our analysis is aided by MWFM modeling of mountain waves at high latitudes and CMAP precipitation indices in the subtropics to provide a clearer picture of global gravity wave dynamics.     

Dark energy,co-evolution of massive black holes with galaxies,and ASTROD-GW

The detection of low frequency band (100 nHz–100 mHz) and very low frequency band (300 pHz–100 nHz) gravitational waves (GWs) is important for exploration of the equation of state of dark energy and the co-evolution of massive black holes (MBHs) with galaxies. Most galaxies are believed to have a massive black hole in the galactic core. In the formation of these black holes, merging and accretion are the two main processes. Merging of massive black holes generate GWs which could be detected by space GW detectors and Pulsar Timing Arrays (PTAs) to cosmological distances. LISA (Laser-Interferometric Space Antenna) is most sensitive to the frequency band 1 mHz–100 mHz, ASTROD-GW (ASTROD [Astrodynamical Space Test of Relativity using Optical Devices] optimized for Gravitational Wave detection) is most sensitive to the frequency band 100 nHz–1 mHz and PTAs are most sensitive to the frequency band 300 pHz–100 nHz. In this paper, we discuss the sensitivities and outlooks of detection of GWs from binary massive black holes in these frequency bands with an emphasis on ASTROD-GW. The GWs generated by the inspirals, merging and subsequent ringdowns of binary black holes are standard sirens to the cosmological distance. Using GW observations, we discuss the methods for determining the equation of state of dark energy and for testing the co-evolution models of massive black holes. ASTROD-GW is an optimization of ASTROD to focus on the goal of detection of GWs. The mission orbits of the 3 spacecraft forming a nearly equilateral triangular array are chosen to be near the Sun-Earth Lagrange points L3, L4 and L5. The 3 spacecraft range interferometrically with one another with arm length about 260 million kilometers. With 52 times longer in arm length compared to that of LISA, the strain detection sensitivity is 52 times better toward larger wavelength. The scientific aim is focused for gravitational wave detection at low frequency. The science goals include detection of GWs from MBHs, and Extreme-Mass-Ratio Black Hole Inspirals (EMRI), and using these observations to find the evolution of the equation of state of dark energy and to explore the co-evolution of massive black holes with galaxies.     

Gravity wave mixing effects on the OH*-layer

Based on an advanced numerical model for excited hydroxyl (OH*) we simulate the effects of gravity waves (GWs) on the OH*-layer in the upper mesosphere. The OH* model takes into account (1) production by the reaction of atomic hydrogen (H) with ozone (O₃), (2) deactivation by atomic oxygen (O), molecular oxygen (O₂), and molecular nitrogen (N₂), (3) spontaneous emission, and (4) loss due to chemical reaction with O. This OH* model is part of a chemistry-transport model (CTM) which is driven by the high-resolution dynamics from the KMCM (Kühlungsborn Mechanistic general Circulation Model) which simulates mid-frequency GWs and their effects on the mean flow in the MLT explicitly. We find that the maximum number density and the height of the OH*-layer peak are strongly determined by the distribution of atomic oxygen and by the temperature. As a results, there are two ways how GWs influence the OH*-layer: (1) through the instantaneous modulation by O and T on short time scales (a few hours), and (2) through vertical mixing of O (days to weeks). The instantaneous variations of the OH*-layer peak altitude due to GWs amount to 5–10 km. Such variations would introduce significant biases in the GW parameters derived from airglow when assuming a constant pressure level of the emission height. Performing a sensitivity experiment we find that on average, the vertical mixing by GWs moves the OH*-layer down by ~2 to 7 km and increases its number density by more than 50%. This effect is strongest at middle and high latitudes during winter where secondary GWs generated in the stratopause region account for large GW amplitudes.     

Multi-GNSS contributions to differential code biases determination and regional ionospheric modeling in China

Due to the limited number and uneven distribution globally of Beidou Satellite System (BDS) stations, the contributions of BDS to global ionosphere modeling is still not significant. In order to give a more realistic evaluation of the ability for BDS in ionosphere monitoring and multi-GNSS contributions to the performance of Differential Code Biases (DCBs) determination and ionosphere modeling, we select 22 stations from Crustal Movement Observation Network of China (CMONOC) to assess the result of regional ionospheric model and DCBs estimates over China where the visible satellites and monitoring stations for BDS are comparable to those of GPS/GLONASS. Note that all the 22 stations can track the dual- and triple-frequency GPS, GLONASS, and BDS observations. In this study, seven solutions, i.e., GPS-only (G), GLONASS-only (R), BDS-only (C), GPS + BDS (GC), GPS + GLONASS (GR), GLONASS + BDS (RC), GPS + GLONASS + BDS (GRC), are used to test the regional ionosphere modeling over the experimental area. Moreover, the performances of them using single-frequency precise point positioning (SF-PPP) method are presented. The experimental results indicate that BDS has the same ionospheric monitoring capability as GPS and GLONASS. Meanwhile, multi-GNSS observations can significantly improve the accuracy of the regional ionospheric models compared with that of GPS-only or GLONASS-only or BDS-only, especially over the edge of the tested region which the accuracy of the model is improved by reducing the RMS of the maximum differences from 5–15 to 2–3 TECu. For satellite DCBs estimates of different systems, the accuracy of them can be improved significantly after combining different system observations, which is improved by reducing the STD of GPS satellite DCB from 0.243 to 0.213, 0.172, and 0.165 ns after adding R, C, and RC observations respectively, with an increment of about 12.3%, 29.4%, and 32.2%. The STD of GLONASS satellite DCB improved from 0.353 to 0.304, 0.271, and 0.243 ns after adding G, C, and GC observations, respectively. The STD of BDS satellite DCB reduced from 0.265 to 0.237, 0.237 and 0.229 ns with the addition of G, R and GR systems respectively, and increased by 10.6%, 10.4%, and 13.6%. From the experimental positioning result, it can be seen that the regional ionospheric models with multi-GNSS observations are better than that with a single satellite system model.     

Indication of convectively generated gravity waves observed by CLAES

Cryogenic Limb Array Etalon Spectrometer temperature data are analyzed for gravity waves (GWs) using the same methods previously employed for CRISTA data. We obtain nearly 1.5 years of continuous GW data between 34°S and 34°N and good coverage at higher latitudes depending on UARS yaw maneuvers. Correlations of GW variances with sea surface temperature indicate strong convective forcing above the Gulf of Mexico and the Kuro-Shio stream as well as above source regions in the southern subtropics during the summer months of the respective hemisphere. However, it appears difficult to clearly separate between convective forcing and the effects of the background winds, for which similar correlations are found. Further progress in accessing GW sources therefore can be expected mainly from improved measurements, which could provide also information about the GW horizontal wavelength as well as from quantitative comparisons to model predictions.     

Recent advances in observations and modeling of the solar ultraviolet and X-ray spectral irradiance

There have been significant, recent advances in understanding the solar ultraviolet (UV) and X-ray spectral irradiance from several different satellite missions and from new efforts in modeling the variations of the solar spectral irradiance. The recent satellite missions with solar UV and X-ray spectral irradiance observations include the X-ray Sensor (XRS) aboard the series of NOAA GOES spacecraft, the Upper Atmosphere Research Satellite (UARS), the SOHO Solar EUV Monitor (SEM), the Solar XUV Photometers (SXP) on the Student Nitric Oxide Explorer (SNOE), the Solar EUV Experiment (SEE) aboard the Thermosphere, Ionosphere, Mesosphere, Dynamics, and Energetics (TIMED) satellite, and the Solar Radiation and Climate Experiment (SORCE) satellite. The combination of these measurements is providing new results on the variability of the solar ultraviolet irradiance throughout the ultraviolet range shortward of 200 nm and over a wide range of time scales ranging from years to seconds. The solar UV variations of flares are especially important for space weather applications and upper atmosphere research, and the period of intense solar storms in October–November 2003 has provided a wealth of new information about solar flares. The new efforts in modeling these solar UV spectral irradiance variations range from simple empirical models that use solar proxies to more complicated physics-based models that use emission measure techniques. These new models provide better understanding and insight into why the solar UV irradiance varies, and they can be used at times when solar observations are not available for atmospheric studies.     

Impact of four dimensional assimilation of satellite data on long-range simulations over the Indian region during monsoon 2006

A number of experiments were conducted to study the impact of updating model basic fields by satellite data (Quick Scatterometer (QSCAT) surface winds and Atmospheric Infrared Sounder (AIRS) temperature and humidity profiles) on long-range simulation during the Indian summer monsoon 2006. The Pennsylvania State University/National Center for Atmospheric Research (PSU/NCAR) mesoscale model version5 (MM5) and its four dimensional data assimilation (FDDA) technique was used for the numerical simulations. The spatial distribution and temporal variation in model simulated basic meteorological parameters and rainfall were verified against the observed fields from National Center for Environmental Prediction (NCEP) analysis and Tropical Rainfall Measuring Mission (TRMM), respectively. The overall analysis of the results from QSCAT surface wind assimilation as compared to control simulation (CNT; without the satellite data assimilation) suggest that a better representation of a single level wind field during model integration fail to make significant improvement in the model simulation both in the basic meteorological parameters and rainfall. The assimilation of temperature and humidity profiles from the AIRS during model integration significantly improved the rainfall prediction during monsoon period. It is found that the improvement in rainfall prediction is attributed to improved thermodynamics structure due to AIRS profile assimilation and the degree of improvement is more in temperature prediction as compared to humidity prediction. It is also found that the prediction over the regions, such as south west part of India and foothills of Himalaya, where a complex orography exists, is not significantly benefited from satellite data assimilation which highlights the need of improvement in the model in addition to a better representation of atmospheric state.     

卫星观测台风重力波数值模拟与直接对比验证

以ECMWF的T799资料作为WRF-ARW（V3.5）初始场,对卫星高光谱红外大气垂直探测器AIRS观测的2013年超强台风苏力激发平流层重力波过程进行数值模拟,并利用卫星观测对数值模拟结果进行了直接对比验证.数值模拟表明,该台风诱发的重力波在20~40km高度逆着东风背景流向东向上传播,在水平方向呈半圆弧状;大气的垂直扰动随着高度的增加而增强,在40km高度上达到0.5m·s-1.基于三维傅里叶变换的波谱分析表明,平流层重力波水平波长中心值在500km附近,周期为3~5h,垂直波长主要为10~26km.分析表明,在18~40km高度的净纬向动量通量为6.7×10-4~1.89×10-3Pa,背景流强迫计算值为-0.23~1.21m·s-1·d-1,且在18km和40km高度的数值较大.最后,基于辐射传输模式计算的直接对比表明,卫星观测与数值模拟同时揭示了激发的平流层波动可传至40km以上高度及距台风中心2000km以外的区域,且不同资料得到的波动形态、方位以及水平尺度具有较好的一致性.      

Rapid local ionosphere modeling based on Precise Point Positioning over Iran: A case study under 2014 solar maximum

Presently, the ionosphere effect is the main source of the error in the Global Positioning System (GPS) observations. This effect can largely be removed by using the two-frequency measurements, while to obtain the reasonable results in the single-frequency applications, an accurate ionosphere model is required. Since the global ionosphere models do not meet our needs everywhere, the local ionosphere models are developed. In this paper, a rapid local ionosphere model over Iran is presented. For this purpose, the GPS observations obtained from 40 GPS stations of the Iranian Permanent GPS Network (IPGN) and 16 other GPS stations around Iran have been used. The observations have been selected under 2014 solar maximum, from the days 058, 107, 188 and 271 of the year 2014 with different geomagnetic activities. Moreover, ionospheric observables based on the precise point positioning (PPP) have been applied to model the ionosphere. To represent our ionosphere model, the B-spline basis functions have been employed and the variance component estimation (VCE) method has been used to regularize the problem.To show the efficiency our PPP-derived local ionosphere model with respect to the International GNSS Service (IGS) global models, these models are applied on the single point positioning using single-frequency observations and their results are compared with the precise coordinates obtained from the double-differenced solution using dual-frequency observations. The results show that the 95th percentile of horizontal and vertical positioning errors of the single-frequency point positioning are about 3.1 and 13.6?m, respectively, when any ionosphere model are not applied. These values significantly improve when the ionosphere models are applied in the solutions. Applying CODE’s Rapid Global ionosphere map (CORG), improvements of 59% and 81% in horizontal and vertical components are observed. These values for the IGS Global ionosphere map (IGSG) are 70% and 82%, respectively. The best results are obtained from our local ionosphere model, where 84% and 87% improvements in horizontal and vertical components are observed. These results confirm the efficiency of our local ionosphere model over Iran with respect to the global models. As a by-product, the Differential Code Biases (DCBs) of the receivers are also estimated. In this line, we found that the intra-day variations of the receiver DCBs could be significant. Therefore, these variations must be taken into account for the precise ionosphere modeling.     

Recent advances in observing mesoscale ocean dynamics with satellite altimetry

The paper provides a review of recent results on mesoscale ocean dynamics derived from satellite altimetry. Since 1992, we have had an unprecedented 18 year high-resolution monitoring of the ocean mesoscale field. Altimetry is often used in mesoscale studies in synergy with other remote sensing techniques and in situ data. This global, high-resolution data set has allowed oceanographers to quantify the previously unknown seasonal and interannual variations in eddy kinetic energy and eddy heat and salt transports, and investigate their causes. Eddy tracking techniques have allowed us to monitor their propagation pathways, to bring to light the meridional divergence of cyclones and anticyclones, and to question the role of Rossby waves versus non-linear eddies in the mid to high latitude bands. Altimetry has also revealed the presence of zonal fronts and jets everywhere in the ocean, and brought to light how mesoscale eddies can impact back onto the atmospheric circulation. Finally, altimetry, in synergy with other observations and high-resolution numerical models, has helped reveal the complexity of the sub-mesoscale features, associated with stirring and mixing around the mesoscale eddies, of great importance for the vertical exchange of oceanic tracers. Altimetry has revealed the complexity of the mesoscale system… the scientific community is now working to understand the interplay between these mesoscale eddies, the ocean interior and its impact on the overlying atmosphere.     

Precise orbit determination and precision comparison for FY3C and FY3D with receiver antenna GPS and BDS PCV using spaceborne BDS and GPS observation data

The FY3C and FY3D satellites were equipped with global navigation satellite occultation detector (GNOS) receivers that received both GPS and BDS-2 signals. For further improving precise orbit determination (POD) precisions, we estimated receiver GPS and BDS signal phase center variations (PCV) models with 2° and 5° resolutions and set the different weights for GPS and BDS-2 observations in the combined POD. The BDS-based POD precision using BDS-2 satellite antenna phase center offset (PCO) values from the China Satellite Navigation Office (CSNO) are not as accurate as those obtained from the International GNSS Service (IGS) Multi-GNSS experiments project (MGEX). The estimated receiver GPS and BDS PCV models with 2° and 5° resolutions were estimated from the GPS phase residuals of GPS-based POD and BDS phase residuals of combined POD, respectively. In most cases, the POD precisions using the estimated PCVs with 2° resolution are superior to those with 5° resolution. The precisions of the BDS-based POD and combined POD were both improved by introducing the receiver BDS PCV models. The weighting for GPS and BDS-2 observations can further improve the precision of the combined POD. The tested results of selected weights are better than those with equal weight in the combined POD. The experiment results show that orbital precisions of FY3C are worse than those of FY3D.     

全球海表流场多尺度结构观测卫星计划

全球海表流场多尺度结构观测卫星计划（Ocean Surface Current multiscale Observation Mission, OSCOM）首次提出海表流场、海面风场和海浪谱（简称 “流–风–浪”）一体化探测的多普勒散射计（Doppler Scatterometer, DOPS）测量原理和系统体制。OSCOM采用Ka-Ku双频多波束圆锥扫描体制的真实孔径雷达,将实现超过1000 km观测刈幅、公里级分辨率的“流–风–浪”一体化卫星直接观测。OSCOM将突破海洋亚中尺度非平衡态动力学、海洋多尺度相互作用、海气耦合的研究瓶颈,支撑实现海洋系统科学、气候变化等理论研究的重大突破。未来,应用OSCOM海表流速观测的模式改进,将奠定海洋非平衡态过程数值模拟、同化和预报的动力学基础,实现海洋和海气耦合模式的重大改进。通过与多源数据融合,OSCOM海流观测的应用将为海洋生物地球化学循环、碳收支研究和国家重大任务提供支撑。OSCOM科学卫星的实施对于我国地球系统科学和卫星对地观测重大应用的突破有至关重要的意义,有望带动我国应用卫星的发展从追赶、并行走向领跑。      

Estimation of differential code biases with Jason-2/3 onboard GPS observations

With the continuous deployment of Low Earth Orbit (LEO) satellites, the estimation of differential code biases (DCBs) based on GNSS observations from LEO has gained increasing attention. Previous studies on LEO-based DCB estimation are usually using the spherical symmetry ionosphere assumption (SSIA), in which a uniform electron density is assumed in a thick shell. In this study, we propose an approach (named the SHLEO method) to simultaneously estimate the satellite and LEO onboard receiver DCBs by modeling the distribution of the global plasmaspheric total electron content (PTEC) above the satellite orbit with a spherical harmonic (SH) function. Compared to the commonly used SSIA method, the SHLEO model improves the GPS satellite DCB estimation accuracy by 13.46% and the stability by 22.34%, respectively. Compared to the GPS satellite DCBs estimated based on the Jason-3-only observations, the accuracy and monthly stability of the satellite DCBs can be improved by 14.42% and 26.8% when both Jason-2 and Jason-3 onboard observations are jointly processed. Compared with the Jason-2 solutions, the GPS satellite DCB estimates based on the fusion of Jason-2 and Jason-3 observations have an improved consistency of better than 18.26% and 9.71% with the products provided by the Center for Orbit Determination in Europe (CODE) and Chinese Academy of Sciences (CAS). Taking the DCB products provided by the German Aerospace Center (DLR) as references, there is no improvement in accuracy of the GPS satellite DCB estimates based on the fusion of Jason-2 and Jason-3 observations than the Jason-2 solutions alone. A periodic variation is found in the time series of both the Jason-3 and Jason-2 onboard receiver DCB estimates. Preliminary analysis of the PTEC distribution based on the estimated SH coefficients are also presented.     

An Integer Precise Point Positioning technique for sea surface observations using a GPS buoy

GPS data dedicated to sea surface observation are usually processed using differential techniques. Unfortunately, the precision of resulting kinematic positions is baseline-length dependent. So, high precision sea surface observations using differential GPS techniques are limited to coasts, lakes, and rivers. Recent improvements in GPS satellite products (orbits, clocks, and phase biases) make phase ambiguity fixing at the zero difference level achievable and opens up the observation of the sea surface without geographical constraints. This paper recalls the concept of the Integer Precise Point Positioning technique and discusses the precision of GPS buoy positioning. A sequential version of the GINS software has been implemented to achieve single epoch GPS positioning. We used 1 Hz data from a two week GPS campaign conducted in the Kerguelen Islands. A GPS buoy has been moored close to a radar gauge and 90 m away from a permanent GPS station. This infrastructure offers the opportunity to compare both kinematic Integer Precise Point Positioning and classical differential GPS positioning techniques to in situ radar gauge data. We found that Precise Point Positioning results are not significantly biased with respect to radar gauge data and that horizontal time series are consistent with differential processing at the sub-centimetre precision level. Nevertheless, standard deviations of height time series with respect to radar gauge data are typically [4–5] cm. The dominant driver for noise at this level is attributed to errors in tropospheric estimates which propagate into position solutions.     

DORIS and GPS monitoring of the Gavdos calibration site in Crete

Due to its specific geographical location as well as its geodetic equipment (DORIS, GNSS, microwave transponder and tide gauges), the Gavdos station in Crete, Greece is one of the very few sites around the world used for satellite altimetry calibration. To investigate the quality of the Gavdos geodetic coordinates and velocities, we analyzed and compared here DORIS and GPS-derived results obtained during several years of observations. The DORIS solution is the latest ignwd11 solution at IGN, expressed in ITRF2008, while the GPS solution was obtained using the GAMIT software package. Current results show that 1–2 mm/yr agreement can be obtained for 3-D velocity, showing a good agreement with current geophysical models. In particular, the agreement obtained for the vertical velocity is around 0.3–0.4 mm/yr, depending on the terrestrial reference frame. As a by-product of these geodetic GPS and DORIS results, Zenith Tropospheric Delays (ZTDs) estimations were also compared in 2010 between these two techniques, and compared to ECMWF values, showing a 6.6 mm agreement in dispersion without any significant difference between GPS and DORIS (with a 97.6% correlation), but with a 13–14 mm agreement in dispersion when comparing to ECMWF model (with only about 90% correlation for both techniques). These tropospheric delay estimations could also provide an external calibration of the tropospheric correction used for the geophysical data of satellite altimetry missions.     

Development of a GNSS water vapour tomography system using algebraic reconstruction techniques

A GNSS water vapour tomography system developed to reconstruct spatially resolved humidity fields in the troposphere is described. The tomography system was designed to process the slant path delays of about 270 German GNSS stations in near real-time with a temporal resolution of 30 min, a horizontal resolution of 40 km and a vertical resolution of 500 m or better. After a short introduction to the GPS slant delay processing the framework of the GNSS tomography is described in detail. Different implementations of the iterative algebraic reconstruction techniques (ART) used to invert the linear inverse problem are discussed. It was found that the multiplicative techniques (MART) provide the best results with least processing time, i.e., a tomographic reconstruction of about 26,000 slant delays on a 8280 cell grid can be obtained in less than 10 min. Different iterative reconstruction techniques are compared with respect to their convergence behaviour and some numerical parameters. The inversion can be considerably stabilized by using additional non-GNSS observations and implementing various constraints. Different strategies for initialising the tomography and utilizing extra information are discussed. At last an example of a reconstructed field of the wet refractivity is presented and compared to the corresponding distribution of the integrated water vapour, an analysis of a numerical weather model (COSMO-DE) and some radiosonde profiles.     

基于多种观测手段的东亚低纬电离层不规则体事件分析

利用海南台站（19.5°N,109.1°E,dip:13.6°N）和磁赤道区的多种地基和天基观测数据,对2011年11月20日观测到的电离层不规则体事件进行了分析.海南台站VHF雷达、电离层闪烁和数字测高仪的综合观测结果表明,当天日落附近发生了强的电离层不规则体事件,主要表现为雷达羽和强闪烁的形态.结合磁赤道区GPS和C/NOFS卫星观测结果进行分析可知,海南台站日落附近出现的雷达羽和强闪烁与南海磁赤道区产生的主等离子体泡存在明显联系.      

重力波非线性传播过程中的饱和与破碎

采用水平方向的显式算法与垂直方向的隐式算法相结合的时间分裂法,建立了二维可压缩大气中重力波非线性传播的数值模式.用本模式对小振幅重力波传播过程的模拟结果与线性重力波理论预测的结果吻合很好,从而验证了本模式的正确性.我们用此模式模拟了有限振幅重力波在非线性传播过程中的饱和与破碎,结果表明,(1)翻转出现在饱和之前,但向破碎演化仍需要一段时间,由于非线性波-波和波-流相互作用使得非线性数值模拟的饱和高度(出现时间)高(早)于线性饱和理论预测的结果;(2)重力波在不稳定之前已经有能量向背景场中转移,破碎直接导致非线性波-波相互作用,造成能量向小尺度短波上转移;(3)背景风场的加速方向,形成射流的方向与重力波的水平传播方向一致,表明重力波与背景流的非线性相互作用加剧了背景风剪切和不稳定性的发展.      

一次暴雨激发平流层重力波的卫星观测与数值模拟

针对卫星Aqua/AIRS观测到的与2011年7月25日山东省乳山市特大暴雨相伴的一次平流层重力波过程,利用中尺度数值模式WRF进行暴雨诱发平流层重力波的数值模拟.对模式输出的垂直速度场和温度扰动场的分析表明,暴雨在平流层内的弧状波结构主要集中在降水云系东侧,水平影响范围大于1000km,且随着高度的增加,圆弧状结构趋于闭合,波动能量显著增强.此外,对垂直速度剖面结构分析表明,受高空东风和风切变的影响,重力波在上传过程中逆着背景风场向东传输,不同高度波动形态各异.基于快速傅里叶变换（FFT）的功率谱分析结果表明,此次暴雨激发的平流层重力波在35km高度的周期为7~20h,水平波长约为1000km,垂直波长为5~10km.通过分析动量通量的垂直输送,定量反映出重力波上传过程中的动力学变化特征.