山基GPS掩星观测实验及其反演原理

分别在湖北省药姑山和九宫山的山头上开展GPS掩星观测实验,成功获取山基掩星观测数据,对掩星事件进行了分析和统计.给出利用山基掩星观测数据反演大气折射指数剖面和电波弯曲角的原理和算法.利用山基GPS掩星观测模拟数据,对反演方法进行试算和检验,结果表明反演方法准确可行.将该反演方法应用于观测数据的反演,获得了观测点高度以下的大气折射率剖面,以及电波弯曲角.实验结果和原理研究表明,山基掩星观测技术是一种潜在的低层大气环境监测新技术.     

COSMIC大气掩星开环数据反演方法

COSMIC星座GPS无线电掩星探测利用GPS开环接收技术提高低层大气观测数据的质量和对上升掩星事件的跟踪能力.开环掩星观测数据受到GPS导航数据调制的影响,在其数据后处理中必须消除该影响以获得高质量的科学反演结果.利用GPS导航数据调制码数据和利用开环数据本身内在的关系等两种方法可以消除该影响.将上述方法应用于COSMIC的掩星事件个例反演,获得了修复的大气附加相位数据;并利用几何光学近似反演方法和全谱反演方法,获得了射线弯曲角. 全谱反演方法获得的弯曲角及其温度反演结果与COSMIC数据中心的结果一致,说明我们的算法是有效的.     

地基GNSS VTEC约束的大气掩星电离层误差修正方法

介绍了一种考虑电离层沿GNSS掩星射线路径分布不对称且兼顾一阶和二阶项的大气掩星电离层误差修正新方法,它综合地基GNSS VTEC的水平变化信息和电离层模式的垂直变化信息,在沿"入射线"与"出射线"双边局部球对称假设下,估算GNSS掩星射线路径上的电子密度;进而计算一阶和二阶项弯曲角电离层误差廓线.采用太阳活动低年的2008年7月15日和太阳活动较高年的2013年7月15日两天的MetOp-A和GRACE掩星观测资料和IGS的GNSS VTEC数据,计算了GNSS大气掩星弯曲角一阶和二阶项电离层误差廓线.对比分析表明：新方法经验模型和理论模型的二阶项弯曲角电离层残差,以及经验模型与实测数据的一阶和二阶项弯曲角电离层误差均具有良好的一致性,因此该方法可用于L2信号数据质量较差或L2信号中断的掩星事件,同时修正大气掩星一阶和二阶电离层误差,从而提高弯曲角精度和掩星观测资料的利用率.     

山基GPS掩星折射率与探空折射率比较

本文简要阐述了山基GPS掩星数据处理过程,介绍了雾灵山(40.60.N,117.48.E,～2118 m)山基GPS掩星观测实验的基本情况,将雾灵山山基GPS掩星折射率与探空折射率进行了比较.按照山基掩星和无线电探空在空间(经纬度)上相差在2.以内、时间上相差在1.5 h以内的条件,共选出了206个山基掩星折射率廓线样本用于比较分析.山基掩星折射率与探空折射率的个例比较,表明山基掩星折射率廓线合理,且和探空折射率基本一致.统计结果显示,在2100 m 以下,山基掩星折射率与探空折射率相对偏差的平均值为-2.0%,标准偏差为7.3%.     

山基GPS掩星与自动气象站观测结果比较

山基GPS掩星技术可以精确获得低层大气折射率信息,其观测原理是:在高山的山顶安装GPS接收机,跟踪GPS卫星的低仰角和负仰角信号,通过载波相位的变化求出大气折射造成的弯曲角剖面,从而利用科学反演方法获得观测点高度以下的大气折射率剖面.为了对山基掩星测量数据进行验证,2005年8月1日～29日,在河北雾灵山(40.60°N,117.48°E,2118m)开展了为期一个月的山基GPS掩星观测与其他探测的时比实验.实验所用的JAVAD双频GPS接收机共工作576小时,接收机朝南观测,记录到山基掩星事件1136次,其中上升掩星621次,下降掩星515次.实验期间,在雾灵山顶利用自动气象观测仪每天24小时观测当地的温度、湿度和气压,由此可计算大气折射率.本文利用雾灵山GPS掩星实验所获得的观测数据,成功反演得到大气折射率剖面,并将所得接收机高度的折射率与雾灵山顶自动气象站观测结果进行比对,分析结果表明:山基GPS掩星和自动气象站观测结果是一致的.山基GPS掩星观测可为低层大气环境监测提供大量数据,具有潜在的应用前景.     

地基GPS低高度角观测反演大气折射率廓线的模拟仿真

用价值函数的最优估计方法,发展了地基GPS弯曲角和相位观测联合反演大气折射率廓线的算法. 在价值函数中提出了一个比较合理的权系统;通过四参数形式的折射率模型,采用变分同化技术,反演出大气折射率廓线. 模拟结果显示,地基GPS观测反演得到的低层大气折射率廓线与真值吻合较好,并可以应用到大气波导的探测.     

COSMIC数据验证FY-3 MWTS通道4的观测结果

利用无线电掩星探测在平流层具有较高精度和稳定性的优势,本文采用COSMIC资料验证我国风云3号(FY-3)系列卫星上微波温度探测仪(MWTS)通道4表征的平流层低层测量结果.两年的FY-3A/3B上MWTS通道4观测亮温与掩星数据模拟的亮温相比较分析表明:两颗卫星平台上MWTS通道4观测亮温都偏高,尤其在热带地区和极区夏季;FY3A/3B上MWTS通道4亮温偏差的月均值变化趋势较为一致,两个年份的变化趋势也较为一致.MWTS通道4表征的平流层低层观测亮温在不同纬度带呈现不同偏差分布特征:热带地区明显系统偏高2~4 K,中纬度地区偏高1 K;而高纬度地区则出现随季节显著变化的偏差,尤其在南极地区季节差异达到5 K,这种依赖于环境温度出现的较大亮温偏差不完全是频率偏移引起.热带地区观测与模拟亮温差异显著,并且热带地区样本对于全球总体亮温偏差影响程度达到20%,这表明COSMIC数据在热带地区的验证结果需谨慎使用.     

电离层掩星反演的变分同化方法

电离层掩星数据反演的传统方法是Abel反演法及其改进方法,而实际电离层的非球对称给电离层电子密度反演带来误差.本文研究了在气象领域广泛应用的变分同化方法在掩星数据反演电离层电子密度廓线的可行性,利用IRI和Nequiek模型模拟掩星真值场和背景场进行数值模拟反演,并与Abel反演法反演结果进行对比.结果表明,变分同化反演方法能够有效综合模式和观测数据,使得反演结果精度较高;与Abel反演法相比,反演的电子密度廓线F2层峰值浓度误差在10％以下,而Abel反演法在20％～30％之间;而且变分同化反演法对误差扰动有较好的过滤性,因此实用性较强.     

电离层GPS掩星观测改正TEC反演方法

电离层掩星观测中,当低轨卫星（LEO）轨道高度较低时,轨道以上的电离层电子总含量（TEC）对掩星反演的影响不能忽略.目前,一般采用指数函数等外推方法来处理该问题,对反演结果可能引起较大误差.为提高电离层掩星反演精度,本文研究利用LEO处于非掩星一侧GPS观测数据的改正TEC新反演方法.用三维射线追踪程序计算出电离层掩星观测模拟数据,分别应用改正TEC方法和外推方法进行反演,将反演结果与所用模式值进行比较.结果表明：对于轨道高度约800km的GPS/MET掩星模拟数据,外推方法和改正TEC方法反演结果都与模式值基本一致;对于轨道高度约400km的CHAMP掩星模拟数据,外推方法误差较大,改正TEC方法反演结果与模式值相符得较好.将改正TEC方法应用于GPS/MET实测数据的反演,取得了合理的结果.这些说明,改正TEC算法是一种有效的电离层掩星反演方法,尤其是对于轨道较低的LEO的电离层掩星观测反演特别有用.     

电离层残差对掩星反演温度精度的影响

本文以MSIS90大气模式和3D NeUoG电离层模式为大气背景,用三维射线追踪法模拟研究了太阳活动强度、地方时、掩星平面方位角对弯曲角电离层残差和温度电离层残差的影响,以及电离层残差对全球日平均温度的影响.结果表明：电离层残差是平流层顶部（35~50 km）和中间层底部（50~70 km）掩星大气温度反演的主要误差.在太阳活动活跃期,电离层残差对单一掩星事件的平流层顶部平均温度的影响可达1.8 K,中间层底部平均温度的影响可达7 K;对全球日平均温度的影响在平流层顶可达-0.6 K,在70 km高度处可达1.2 K.发展新的电离层改正方法或电离层残差修正算法对提高掩星大气反演精度和全球气候监测意义重大.     

Typhoon-induced ocean responses off the southwest coast of Taiwan

Typhoons can cause substantial sea surface cooling (typically 2–4 °C), which is usually biased to the right side of the storm track. Under influence of the complex bathymetry of the southern Taiwan Strait (TS), two types of sea surface temperature (SST) response, cooling and warming, each associated with a different type of typhoon track were identified using satellite and mooring observations. When a typhoon moved westward (or northwestward) and passed through the TS (track A), the SST cooling in the TS was biased toward the left of the storm track. Numerical model results indicated that in track A, strong wind stress accelerates the flow east of the Taiwan Banks and drove the bottom flow to uplift due to the topography. Moreover, both wind stress and wind stress curl enhanced the Luzon loop. After the typhoon passed, the mean circulation was modified around strong cooling in the southern TS, causing more South China Sea surface water to be distributed to the Kuroshio region. However, when a typhoon moved westward (or northwestward) and passed south of the TS, SST warming was induced in the southern TS (track B). The model results indicated that when the typhoon passed to the south of the TS, the typhoon-induced horizontal divergent flow travelled to the north, where it encountered the shallower shelf of the TS that was confined to the water, causing warm water transported into this area to accumulate and downwelling to occur. This can be regarded as redistributing the heat content in the shelf area. After the typhoon, the thickened mixed layer resulting from downwelling prevented the formation of near-inertial waves and reduced the vertical mixing.     

The four dimensional variational data assimilation with multiple regularization parameters as a weak constraint (Tikh-4D-Var) and its preliminary application on typhoon initialization

Traditional variational data assimilation (VDA) with only one regularization parameter constraint cannot produce optimal error tuning for all observations. In this paper, a new data assimilation method of “four dimensional variational data assimilation (4D-Var) with multiple regularization parameters as a weak constraint (Tikh-4D-Var)” is proposed by imposing different regularization parameters for different observations. Meanwhile, a new multiple regularization parameters selection method, which is suitable for actual high-dimensional data assimilation system, is proposed based on the posterior information of 4D-Var system. Compared with the traditional single regularization parameter selection method, computation of the proposed multiple regularization parameters selection method is smaller. Based on WRF3.3.1 4D-Var data assimilation system, initialization and simulation of typhoon Chaba (2010) with the new Tikh-4D-Var method are compared with its counterpart 4D-Var to demonstrate the effectiveness of the new method. Results show that the new Tikh-4D-Var method can accelerate the convergence with less iterations. Moreover, compared with 4D-Var method, the typhoon track, intensity (including center surface pressure and maximum wind speed) and structure prediction are obviously improved with Tikh-4D-Var method for 72-h prediction. In addition, the accuracy of the observation error variances can be reflected by the multiple regularization parameters.     

Numerical study of upper ocean response to a typhoon moving zonally across the Luzon Strait

The Luzon Strait (LS) is a wide channel between Taiwan and the Luzon islands. Eastward of the LS, the Kuroshio Current (KC) flows northward along the eastern coasts of Luzon and Taiwan. A typhoon is a strong and localized low-pressure weather system that occurs frequently in the vicinity of the Taiwan area. One typical typhoon track, passing through the seas surrounding Taiwan, is a zonal path across the LS. The satellite measured SST, corresponding to typhoons Pabuk (5 Aug. 2007–12 Aug. 2007) and Dujuan (29 Aug. 2003–5 Sep. 2003), which both moved along this path, demonstrated that a classic right bias cooling occurs to the north of the storm track, during the typhoon forced period. However, some cold anomaly water also present toward left (south) of the storm track east of the LS in the relaxation period. This paper adopted a three-dimensional hydrostatic primitive equation model to study the possible causes of this southward transport of cold waters. Both model results and the observed SST anomaly revealed that the strength of the upper ocean cooling depends on whether a resonant regime between the typhoon winds and the near-inertial currents can be excited. To the east of the LS, the convergence between the warm Kuroshio water and the cold wakes in the poststorm period will enhance the southward spreading of cold anomaly water. The enhanced vertical mixing, induced by the southward propagation of nearly inertial waves associated with the cold wakes, can also produce some cold anomaly to the south of a storm track in the poststorm period. Both mechanisms can contribute to the occurrence of some cold anomaly water to the south of the storm track east of the KC. To the west of the LS, the convergence between the warm Kuroshio water and the cold upwelling water from the northern South China Sea can further strengthen the Kuroshio front in the LS.     

两种边界层参数化方案模拟热带气旋Megi(2010)路径差异的机理分析

本文利用中尺度数值模式WRF,分别采用YSU和MYJ两种边界层参数化方案对2010年超强台风Megi的移动路径进行了模拟,研究了热带气旋(TC)路径模拟对边界层方案的敏感性,并从模拟TC尺度差异所造成的影响角度揭示了模式边界层方案影响TC路径的机理.结果表明:由于两种方案对边界层垂直混合作用过程的描述不同,两个试验模拟的低层水汽垂直输送存在差异.相对于能很好模拟出Megi路径的MYJ方案,YSU方案模拟的TC外围螺旋雨带更活跃,造成TC尺度增大,引起TC中心北侧外围气压梯度和径向风速增加,使得由副高向TC中心输送更多的质量,造成副高异常减弱,从而导致由副高主导的引导气流发生改变,最终使得采用YSU方案模拟的Megi路径出现提前转向.     

Multi‐scale assimilation of root zone soil water predictions

When hydrology model parameters are determined, a traditional data assimilation method (such as Kalman filter) and a hydrology model can estimate the root zone soil water with uncertain state variables (such as initial soil water content). The simulated result can be quite good. However, when a key soil hydraulic property, such as the saturated hydraulic conductivity, is overestimated or underestimated, the traditional soil water assimilation process will produce a persistent bias in its predictions. In this paper, we present and demonstrate a new multi‐scale assimilation method by combining the direct insertion assimilation method, particle swarm optimisation (PSO) algorithm and Richards equation. We study the possibility of estimating root zone soil water with a multi‐scale assimilation method by using observed in situ data from the Wudaogou experiment station, Huaihe River Basin, China. The results indicate there is a persistent bias between simulated and observed values when the direct insertion assimilation surface soil water content is used to estimate root zone soil water contents. Using a multi‐scale assimilation method (PSO algorithm and direct insertion assimilation) and an assumed bottom boundary condition, the results show some obvious improvement, but the root mean square error is still relatively large. When the bottom boundary condition is similar to the actual situation, the multi‐scale assimilation method can well represent the root zone soil water content. The results indicate that the method is useful in estimating root zone soil water when available soil water data are limited to the surface layer and the initial soil water content even when the soil hydraulic conductivities are uncertain. Copyright © 2011 John Wiley & Sons, Ltd.     

斜腿夹角对V形墩连续刚构桥地震响应的影响研究

为研究斜腿夹角对V形墩连续刚构桥地震响应的影响及合理斜腿夹角角度,以一座典型V形墩预应力混凝土连续刚构桥为研究对象,采用有限元分析方法研究了斜腿夹角θ对桥梁内力及位移的影响,得出了θ对结构地震响应的影响规律和变化曲线。研究结果表明:随着斜腿夹角的增加,在纵向地震力作用下,墩底纵向弯矩逐渐减小,墩顶和主梁墩顶支撑处纵向弯矩逐渐增大;在横向地震力作用下,跨中横向弯矩逐渐减小,墩底横向弯矩逐渐增大,墩顶横向弯矩基本不变;在竖向地震力作用下,墩底和墩顶竖向弯矩逐渐增大,主梁支撑处竖向弯矩逐渐减小;斜腿夹角对纵向或横向地震力作用下结构位移影响不大,对竖向地震力作用下的位移影响较大。在满足静力设计的前提下,当两斜腿夹角为90°时,结构地震响应相对较小,受力合理性最优。研究成果可为该类桥梁的抗震设计与斜腿夹角角度选取提供参考和依据。     

The Assimilation of GPS Radio Occultation Data and its Impact on Rainfall Prediction along the West Coast of India during Monsoon 2002

In this study, the Weather Research and Forecasting (WRF-2.0.3.1) model with three-dimensional variational data assimilation (3DVAR) was utilized to study a heavy rainfall event along the west coast of India with and without the assimilation of GPS occultation refractivity soundings in the monsoon period of 2002. The WRF model is a next-generation mesoscale numerical weather prediction system designed to serve both operational forecasting and atmospheric research communities. The Global Positioning System (GPS) radio occultation (RO) refractivity data, processed by UCAR, were obtained from the CHAMP and SAC-C missions. This study investigates the impact of thirteen GPS occultation refractivity soundings only, as assimilated into the WRF model with 3DVAR, on the rainfall prediction over the western coastal mountain of India. The model simulation, with the finest resolution of 10 km, was in good agreement with rainfall observations, up to 72-h forecast. There are some subtle but important differences in predicted rainfalls between the control run CN (without the assimilation of refractivity soundings) and G13 (with the assimilation of thirteen GPS RO soundings). In general, the assimilation run G13 gives a better prediction in terms of both rainfall locations and amounts at later times. The moisture increments were analyzed at the initial and forecast times to assess the impact of GPS RO data assimilation. The results indicate that remote soundings in the forcing region could have significant impacts on distant downstream regions. It is anticipated, based on this study, that considerably occultation soundings available from the six-satellite constellation of FORMOSAT-3/COSMIC would have even more significant impacts on weather prediction in this region.     

强台风对汶川大地震和昆仑山大地震“震前扰动”影响的分析

2008年8.0级汶川大地震发生前约40 h起,中国地震局测震台网中的宽带地震仪、重力仪和倾斜仪都记录到了明显的低频扰动现象.然而,在此期间西太平洋上空产生了一次强台风Rammasun（威马逊）;与此类似,2001年8.1级昆仑山大地震发生前约60 h起,国家地震局测震台网中的宽带地震仪也记录到了明显的低频扰动现象,在此期间中国南海上空产生了一次强台风LingLing（玲玲）.本文针对这两大地震的“震前扰动”现象进行了研究,旨在探讨其“震前扰动”现象是否与台风引起的巨浪产生的地脉动效应有关.我们的研究结果表明：在中国内陆观测到的异常扰动与台风路径以及台风的强度密切相关.近海台风可在中国大陆内地引起明显的异常扰动;而远离大陆海岸位于深水海域上空的台风,很难在中国大陆内地引起明显的异常地脉动变化.台风LingLing的路径距离中国海岸最近处仅有约700 km,而台风Rammasun距离中国大陆海岸最近处却达1300 km.本文的对比分析研究表明：路径和强度类似于LingLing的台风在中国内陆可引起明显的异常扰动.因此,中国内陆和沿海许多地震台记录到的昆仑山大地震的“震前扰动”源于LingLing地脉动效应的可能性较大;但位于新疆的地震台记录到的扰动信号频率却明显偏低,这一现象值得注意.路径和强度类似于Rammasun的台风很难在中国内陆引起明显的异常扰动,因此汶川大地震的“震前扰动”完全源于台风Rammasun地脉动效应的可能性较小.汶川大地震的“震前扰动”产生的原因和机理值得我们去进行深入的探讨和研究.     

A soil moisture assimilation scheme based on the ensemble Kalman filter using microwave brightness temperature

This study presents a soil moisture assimilation scheme, which could assimilate microwave brightness temperature directly, based on the ensemble Kalman filter and the shuffled complex evolution method (SCE-UA). It uses the soil water model of the land surface model CLM3.0 as the forecast operator, and a radiative transfer model (RTM) as the observation operator in the assimilation system. The assimilation scheme is implemented in two phases: the parameter calibration phase and the pure soil moisture assimilation phase. The vegetation optical thickness and surface roughness parameters in the RTM are calibrated by SCE-UA method and the optimal parameters are used as the final model parameters of the observation operator in the assimilation phase. The ideal experiments with synthetic data indicate that this scheme could significantly improve the simulation of soil moisture at the surface layer. Furthermore, the estimation of soil moisture in the deeper layers could also be improved to a certain extent. The real assimilation experiments with AMSR-E brightness temperature at 10.65 GHz (vertical polarization) show that the root mean square error (RMSE) of soil moisture in the top layer (0–10 cm) by assimilation is 0.03355 m³ · m⁻³, which is reduced by 33.6% compared with that by simulation (0.05052 m³ · m⁻³). The mean RMSE by assimilation for the deeper layers (10–50 cm) is also reduced by 20.9%. All these experiments demonstrate the reasonability of the assimilation scheme developed in this study.