周建印 項杰 黃思訓(xùn)2)?

1) (國防科技大學(xué)氣象海洋學(xué)院, 南京 211101)

2) (國家海洋局第二海洋研究所, 衛(wèi)星海洋環(huán)境動力學(xué)國家重點實驗室, 杭州 310012)

提出了一種確定大氣邊界層頂高度的數(shù)值微分新方法, 該方法使用了正則化技術(shù), 把對彎角廓線求導(dǎo)數(shù)的數(shù)值微分問題轉(zhuǎn)化為求目標(biāo)泛函極小值的問題, 采用雙參數(shù)模型函數(shù)方法來選擇正則化參數(shù), 最后利用最大梯度法確定邊界層頂高度.首先通過兩個數(shù)值實驗驗證了新方法的有效性, 實驗結(jié)果顯示, 隨著掩星資料噪音的增多, 由差分法和結(jié)合L曲線方案的數(shù)值微分方法得到的邊界層頂高度波動增大, 而通過雙參數(shù)模型函數(shù)方法得到的高度很穩(wěn)定, 這說明新方法能夠很好地過濾噪音, 從而保留廓線的主要信息.隨后基于2007—2011年1, 4, 7, 10月的COSMIC彎角數(shù)據(jù), 利用新方法分析了全球海洋大氣邊界層頂高度的季節(jié)特征, 并與用掩星資料自帶的大氣邊界層頂高度數(shù)據(jù)zbalmax得到的季節(jié)分布進(jìn)行對比.結(jié)果表明, 兩者的季節(jié)分布特征十分一致: 海溫相對周圍海域高的區(qū)域, 邊界層頂高度較高, 反之, 邊界層頂高度較低; 在暖流經(jīng)過的海域, 邊界層頂高度較高, 在寒流經(jīng)過的海域, 邊界層頂?shù)母叨认鄬^低.

1 引 言

行星邊界層是對流層底層的大氣, 下墊面主要通過熱量輻射傳輸、摩擦阻力、水汽通量、污染物排放以及地形對邊界層發(fā)生作用[1].邊界層的結(jié)構(gòu)復(fù)雜多變, 邊界層頂高度隨時空變化而變化, 通常定義為混合層(常值通量層)與自由大氣過渡的高度[2], 高度變化范圍在幾百米到幾千米之內(nèi).邊界層頂高度在數(shù)值預(yù)報模式中的邊界層參數(shù)化方案以及氣溶膠的反演中有著重要的作用, 具體表現(xiàn)為: 如果邊界層頂高度太低, 則邊界層與云層解耦,從而抑制了熱量、水分和湍流動能從海洋表面到云層的垂直傳輸, 并可能加速云的消散; 如果邊界層頂高度太高, 則會形成積云而不是層狀云[3], 所以確定邊界層頂高度是一項十分重要的工作.

邊界層頂往往存在濕度、溫度、折射率、氣溶膠粒子以及云滴的急劇變化[4], 無線電掩星探測獲得的彎角和折射率廓線在邊界層頂處往往具有大的垂直梯度, 因此, 最大梯度所在的高度被認(rèn)為是邊界層頂高度[5].在計算邊界層頂高度的觀測數(shù)據(jù)方面, 由于邊界層頂高度的變化在幾百米左右, 而被動紅外和微波探測的分辨率在1—2 km, 故使用分辨率大約在20 m的掩星資料計算邊界層頂高度更加精確, 而且全球定位系統(tǒng)掩星探測(global positioning system radio occultation, GPS-RO)技術(shù)具有全天候、全球覆蓋、長期穩(wěn)定等優(yōu)點, 使得在分析全球邊界層頂高度時更具有優(yōu)勢.

前人在計算邊界層頂高度上, 提出了許多方法.Holzworth[6]通過從最高表面溫度沿干絕熱線隨高度上升到與最近觀測到的溫度廓線的交點, 來得到最大混合層高度; Coulter[7]比較了由激光雷達(dá)、聲雷達(dá)以及溫度廓線得到的混合層頂高度, 發(fā)現(xiàn)激光雷達(dá)與聲雷達(dá)探測的高度基本一致.Van Pul等[8]由無線電探空資料通過氣塊法(parcel method)得到了邊界層頂高度.在 20 世紀(jì), Bianco和Wilczak[9]基于模糊邏輯方法的混合深度識別算法來確定邊界層深度; Lokoshchenko[10]提出了一種方法, 即把參數(shù)積分(如不穩(wěn)定能量)最大值所在的高度作為混合層高度; Balsley等[11]使用最小風(fēng)切變和最大Brunt-V?is?l?頻率對夜晚穩(wěn)定邊界層頂高度進(jìn)行計算; Sokolovskiy等[12]把COSMIC掩星探測折射率梯度的最大下降點的高度確定為邊界層頂高度, 2007年Sokolovskiy等[13]又利用COSMIC掩星彎角資料對邊界層頂高度進(jìn)行了計算; Baars等[14]修正了小波協(xié)方差方法, 利用雷達(dá)資料得出邊界層頂高度; Seidel等[15]用無線電掩星資料比較了計算邊界層頂高度的6個氣象要素(如溫度、濕度、折射率等), 結(jié)果顯示折射率和位溫的梯度最小值得出的邊界層頂高度較高;Dai等[16]使用多個外場觀測項目中飛機(jī)測量的廓線資料, 重新計算了邊界層頂高度; Yan等[17]把帶有正則化技術(shù)的數(shù)值微分方法應(yīng)用于COSMIC掩星彎角資料來確定大氣邊界層頂高度.

正則化技術(shù)是蘇聯(lián)科學(xué)院院士吉洪諾夫(Tikhonov A N)[18]為了獲得線性不適定問題穩(wěn)定的近似解而于1963年提出的.隨著研究的深入,正則化的方法得到了蓬勃的發(fā)展, 目前常見的正則化方法主要有Tikhonov正則化方法[19], Lavrentiev正則化方法[20], 迭代正則化方法(如Tikhonov迭代方法[21]、Landweber迭代方法[22]), 以及基于譜理論的正則化方法(如TSVD (譜截斷)正則化方法[23])等.前人在大氣科學(xué)領(lǐng)域中應(yīng)用正則化方法做了很多有意義的工作[24?28], 而其中正則化方法的一個關(guān)鍵環(huán)節(jié)是正則化參數(shù)的選擇, 常用的正則化參數(shù)選擇方法有L曲線準(zhǔn)則[29?31]、廣義交叉驗證準(zhǔn)則[32,33]、Morozov 偏差準(zhǔn)則[34]等.其中,Morozov偏差準(zhǔn)則是一種應(yīng)用非常廣泛的后驗選取策略, 但是此原則仍然不夠理想, 于是1984年Morozov[35]提出了吸收的Morozov準(zhǔn)則(damped Morozov principle), 它是通過求解吸收的 Morozov偏差方程來得到正則化參數(shù), 但該方法在進(jìn)行數(shù)值求解時, 收斂速度慢且計算量巨大.于是在1998年, Kunisch和Zou[36]提出單參數(shù)模型函數(shù)法來求正則化參數(shù), 從此模型函數(shù)法為減少計算量指出了一條明確的道路.

在Yan 等[17]提出的帶有正則化技術(shù)的數(shù)值微分中, 正則化參數(shù)的選擇采用了L曲線方法(因此, 該方法在下文中簡稱為L曲線方法), 但是此方法計算量大, 且L曲線方法的收斂性問題在理論上并沒有得到解決.因此, 本文提出一種計算大氣邊界層頂高度的新的數(shù)值微分方法, 其中仍然應(yīng)用了正則化技術(shù), 但是正則化參數(shù)的選擇采用了雙參數(shù)模型函數(shù)方法.本文第2節(jié)首先介紹了研究所用到的數(shù)據(jù), 然后介紹雙參數(shù)模型函數(shù)方法的原理;第3節(jié)進(jìn)行數(shù)值實驗, 并與另外兩種方法的結(jié)果作比較, 以驗證方法的合理性; 第4節(jié)應(yīng)用雙參數(shù)模型函數(shù)方法得到全球海洋邊界層頂高度的季節(jié)分布特征, 并與利用掩星資料自帶的邊界層頂高度得到的結(jié)果作對比; 第5節(jié)是總結(jié).

2 數(shù)據(jù)與方法

2.1 數(shù) 據(jù)

FORMOSAT-3/COSMIC是美國與中國臺灣合作的一個GPS掩星探測項目, 本文使用的彎角廓線數(shù)據(jù)來自COSMIC項目, 包括2007—2011年的 1, 4, 7, 10月份的全球的彎角數(shù)據(jù) (https://cdaac-www.cosmic.ucar.edu/cdaac/products.html),總樣本為1074349條廓線, 其中平均每個月的數(shù)量在50000—70000條, 4和7月份比1和10月份數(shù)量多 30000—40000 條 (圖1); 同時, atmPrf文件本身還自帶了利用彎角數(shù)據(jù)得到的大氣邊界層頂高度數(shù)據(jù)zbalmax.考慮到數(shù)據(jù)不可避免會帶有缺省缺測等不足, 所以從以下兩方面進(jìn)行質(zhì)量控制.1)篩選數(shù)據(jù).考慮到邊界層頂基本都在500 m以上, 5 km 以下, 所以本文采用這樣的彎角數(shù)據(jù),即廓線的最低高度在0.5 km以上, 最高高度在5 km 以下.2)控制數(shù)據(jù)長度.選擇有效數(shù)據(jù)的長度至少在5個以上, 以減小由于數(shù)據(jù)缺少而帶來的誤差.

圖1 2007?2011年1, 4, 7, 10月份掩星廓線數(shù)量的時間分布Fig.1.Temporal distribution of the number of occultation profiles in January, April, July and October, 2007?2011.

2.2 雙參數(shù)模型函數(shù)方法

大氣邊界層頂常常會有溫度或水汽壓的急劇變化(溫度急劇上升或水汽含量急劇下降), 因此,可以采用最大梯度法來確定邊界層頂高度.近些年來, 使用“開環(huán)”跟蹤技術(shù)之后, GPS 掩星資料在低層大氣的準(zhǔn)確度有了顯著提高, 使得GPS掩星資料可以提供更多有價值的低對流層(尤其是離地面2 km以下)大氣的信息, 但是由于水汽的影響,低對流層(尤其是離地面2 km以下)大氣廓線數(shù)據(jù)仍然含有負(fù)偏差.Yan 等[17]提出的數(shù)值微分方法中, 正則化參數(shù)的選擇采用了L曲線法, 由于L曲線法具有前面提到的兩個缺陷, 所以本文采用雙參數(shù)模型函數(shù)方法來確定彎角梯度.

2.2.1 模型的構(gòu)建

模型構(gòu)建的基本思路與Yan 等[17]相同, 為了文中敘述的方便, 下面對模型的構(gòu)建作簡要介紹.設(shè)彎角為連續(xù)可微的函數(shù) α (z) , 彎角廓線的梯度為 φ (z)= α′(z) , z 為高度, 范圍為 z0?z?zn.

在高度范圍 I =[z0,zn]中取i=1,2,···,n, 根據(jù)牛頓萊布尼茲公式可以得到

對于(1)式右側(cè), 可以用辛普森法則近似為

其中, 2h=zi+1?zi?1=2(zn?z0)/n, 本文將廓線等分為 n 份, 每一份高度 h 為 10 m, 誤差為error= ?h5/90φ(4)(ξ),ξ∈ [zi?1,zi+1], 如 果 彎 角梯度變化劇烈, 比如存在鋒面過程等劇烈變化的天氣時, 則誤差不容忽視.

將(2)式代入(1)式可以得到

經(jīng)整理后可得到彎角梯度應(yīng)滿足的線性方程:

其中

此時, 問題歸結(jié)為求線性方程(4)的解.

2.2.2 雙參數(shù)模型函數(shù)方法的應(yīng)用

考慮模型(4), 如果觀測資料b帶有高頻成分,則在反演大氣彎角垂直梯度時, 結(jié)果會產(chǎn)生較大的誤差.所以, 為了克服 (4) 式的不適定性, 本文把(4)式轉(zhuǎn)化為求解如下目標(biāo)泛函極小值的問題:

其中X是滿足泛函極小值的解, α ,β 是待求的雙正則化參數(shù), 而矩陣L為一階導(dǎo)數(shù)算子, 表示對x進(jìn)行約束, 使得x的梯度變化不至于過于劇烈,起到平滑的效果.

本文與Yan 等[17]文章中的不同之處是采用雙參數(shù)的模型函數(shù)方法[36?39]來確定目標(biāo)泛函中兩個參數(shù)的最優(yōu)解, 進(jìn)而求解出彎角梯度, 其基本技術(shù)路線如下:

令目標(biāo)函數(shù) F (α,β) 為

此時吸收的Morozov偏差方程為

式中 γ >1,μ>1 為吸收系數(shù), δ 為誤差水平.因為(8)式是關(guān)于 α ,β 的高度非線性方程, 用通常的迭代方法求解結(jié)果并不理想(比如Newton法和擬Newton法), 原因是他們只具有局部收斂的性質(zhì),而且對初值的要求比較高.因此, 我們采用近年來得到廣泛應(yīng)用的模型函數(shù)法來確定正則化參數(shù), 模型函數(shù)法的好處是不但計算量大大減少, 而且還保證了收斂性.模型函數(shù)有各種選擇的方法[38,39], 為簡單起見采用線性模型函數(shù)mk(α,β):=Tk+αCk+βDk表示第k次迭代的結(jié)果以逼近 F (α,β) , 其中

圖2 雙參數(shù)模型函數(shù)法流程圖Fig.2.Flow chart of two-parameter model function method.

將(9)式代入(8)式可得

確定正則化雙參數(shù)的計算流程如圖2所示.因為在正則化方法中物理量往往是無量綱的, 為了使正則化中的量綱一致, 將 α ,β 做了以下處理: 將 α 擴(kuò)大500×std(b)倍, 將 β 擴(kuò)大 s td(b) 倍, 這里 std(·) 表示取標(biāo)準(zhǔn)差.

3 數(shù)值實驗

本節(jié)比較了三種確定邊界層頂高度的方法, 分別是模型函數(shù)法、差分法和L曲線方法.所用數(shù)據(jù)是 2011年1, 4, 7, 10月份的 COSMIC 彎角資料.在處理之前先對數(shù)據(jù)進(jìn)行質(zhì)量控制, 然后借助spline插值法將資料插值到間隔為10 m的網(wǎng)格點上.雙參數(shù)模型函數(shù)法的兩個初始參數(shù) α ,β 都取為 1, 根據(jù)彎角的量級為 10?2, 故取 δ =10?2, 收斂參數(shù) ε =10?4.由差分法、L曲線法以及模型函數(shù)法計算出來的邊界層頂高度分別記為 Hd,HL,HM.

掩星探測資料不可避免地會帶有誤差, 為了比較在存在噪聲情況下幾種方法的表現(xiàn), 本節(jié)將[?δ,+δ]上的隨機(jī)誤差添加到彎角資料, 由于在邊界層中彎角的量級為 10?2, 故令 δ ∈[0.0001,0.01],在帶有噪聲的彎角資料的基礎(chǔ)上, 通過與差分法和L曲線法相比較, 檢驗?zāi)Ｐ秃瘮?shù)法在確定邊界層頂高度上是否具有準(zhǔn)確性.

圖3 由差分法和模型函數(shù)法得到的彎角梯度廓線(BA表示彎角) (a)彎角廓線添加隨機(jī)誤差 δ =0.0025 , 邊界層頂高度HM=1.2km,Hd=1.4km ; (b)彎角廓線添加隨機(jī)誤差 δ =0.005 , 邊界層頂高度 HM=1.1km,Hd=1.4km ; (c)彎角廓線添加隨機(jī)誤 差 δ =0.0075 , 邊界層 頂高 度 HM=1.2km,Hd=1.8km ; (d)彎角 廓線 添加隨 機(jī)誤差 δ =0.01 , 邊界 層頂高度 HM =1.2 km, Hd = 3.7 kmFig.3.Angle gradient profile obtained by the difference method and the model function method using the bending angle gradient profile (BA represents the bending angle): (a) Bending angle profile with uniform random error δ =0.0025 , boundary layer top height HM=1.2km,Hd=1.4km ; (b) bending angle profile with uniform random error δ =0.005 , boundary layer top height HM=1.1km,Hd=1.4km ; (c) bending angle profile with uniform random error δ =0.0075 , boundary layer top height HM=1.2km,Hd=1.8km ; (d) bending angle profile with uniform random error δ =0.01 , boundary layer top height HM=1.2km,Hd=3.7km.

圖3選取了廓線atmPrf_C001.2011.182.05.19.G18_2013.3520_nc, 并且分別添加隨機(jī)誤差δ=0.0025,0.005,0.0075,0.01, 比較了差分法和數(shù)值微分結(jié)合模型函數(shù)方法得到的邊界層頂高度.從結(jié)果可以看出, 圖3(a)是由帶有誤差 δ =0.0025 的彎角得到的彎角梯度, 差分法得到的彎角梯度廓線在一定水平范圍內(nèi)呈密集的鋸齒狀, 其最小梯度在1.4 km左右, 數(shù)值微分方法結(jié)合模型函數(shù)法得到的彎角廓線比較光滑, 其最小梯度在1.2 km左右, 由圖3(b)—(d)可知, 隨著彎角資料誤差的增大, HM始終穩(wěn)定在 1.2 km 左右, 而 Hd由 1.4 km(圖3(a)和圖3(b))變?yōu)?3.7 km(圖3(d)), 變化幅度超過2 km, 說明當(dāng)掩星資料中存在觀測誤差時,模型函數(shù)法能夠很好地過濾噪聲, 保留主要信息,得到的邊界層頂高度的結(jié)果更加穩(wěn)定.

圖4 由差分法、L曲線法和模型函數(shù)法得到的邊界層頂高度 (a) 三種方法基于廓線1得到的邊界層頂高度, Htrue =3.15 km, std(HM) = 0.013, std(HL) = 0.44, std(HM) = 0.61;(b) 三種方法基于廓線2 得到的邊界層頂高度, Htrue = 4.55 km,std(HM) = 0.020, std(HL) = 0.89, std(HM) = 1.19Fig.4.Height of the boundary layer obtained by the three methods: (a) Three methods to get the height of the boundary layer top based on profile 1, Htrue = 3.15 km, std(HM) =0.013, std(HL) = 0.44, std(HM) = 0.61; (b) three methods to get the height of the boundary layer top based on profile 2,Htrue = 4.55 km, std(HM) = 0.020, std(HL) = 0.89, std(HM) =1.19.

圖4給出了三種方法(模型函數(shù)法、差分法、L曲線法)分別在兩條添加誤差后的彎角廓線(atmprf_C001.2001.182.00.22.G23_2013.3520_n c, atmPrf_C001.2011.182.00.39.G20_2013.3520_nc)上的表現(xiàn), 以結(jié)果的標(biāo)準(zhǔn)差衡量方法的穩(wěn)定度.由圖4(a)可以看出, HM比較穩(wěn)定, 保持在 3.15 km左右, 而 Hd和 HL隨著誤差的增大, 波動越來越大,其中, HM,Hd,HL的穩(wěn)定度分別為 0.013, 0.44 和0.61; 同樣由圖4(b) 也可以看出, HM依然是比較穩(wěn)定, 保持在 4.55 km 左右, 而 Hd和 HL隨著誤差的增大, 波動也越來越大, 其中 HM,Hd,HL的穩(wěn)定度分別為0.02, 0.89和1.19.由以上分析可以看出,彎角資料的噪聲越大, Hd振蕩程度就越強(qiáng), HL振蕩的程度次之, 而 HM的穩(wěn)定性則比較好, 說明對于含有噪聲的彎角資料, 模型函數(shù)法也能穩(wěn)定地得到邊界層頂高度, 所以模型函數(shù)法在確定邊界層頂高度時, 準(zhǔn)確性能得到保證.

4 全球邊界層頂高度的季節(jié)特征

本節(jié)應(yīng)用模型函數(shù)法, 對海洋上邊界層頂高度的季節(jié)特征進(jìn)行分析.圖5展示了2007—2011年的 1, 4, 7, 10 四個月份年平均的邊界層頂高度, 可以看出, 1月份南半球的邊界層頂達(dá)到最高, 尤其是中高緯度海洋上的邊界層頂高度達(dá)到2 km以上;4月份南半球的邊界層頂高度開始降低, 南美洲與南極大陸南設(shè)得蘭群島之間的德雷克海峽附近降低的程度最為明顯, 而北極附近海域的邊界層頂高度也降低, 原因在于北極的海冰覆蓋面積在3月份達(dá)到高峰, 4月份海冰的面積比1月份多[40], 所以4月份下墊面溫度更低, 相應(yīng)的邊界層頂高度也在4月份比1月份低; 七月份南半球中高緯度的邊界層頂高度在四個月份中最低, 北半球的低緯度地區(qū)邊界層頂高度明顯增高, 如墨西哥灣地區(qū)、阿拉伯地區(qū)以及日本海地區(qū), 而受西太平洋副熱帶高壓下沉氣流影響, 我國東南沿海、孟加拉灣和菲律賓附近海域邊界層頂高度明顯降低; 10月份南半球中高緯邊界層頂高度略有回升, 北極海域邊界層頂高度開始降低, 而此時我國東南沿海、孟加拉灣和菲律賓附近海域邊界層頂高度明顯升高.

從以上現(xiàn)象分析可以得知, 隨著太陽直射點的南北移動, 海洋上的邊界層頂高度呈現(xiàn)明顯的季節(jié)變化特征, 邊界層頂高度與海面的溫度密切相關(guān),當(dāng)海面的溫度相對周圍平均溫度越低時, 邊界層頂高度越低, 反之, 邊界層頂高度就越高.但是四個月份共同的特征是, 在寒流流經(jīng)的區(qū)域, 邊界層頂高度普遍偏低, 比如巴西寒流和加利福尼亞寒流,暖流流經(jīng)的區(qū)域邊界層頂普遍較高, 比如赤道暖流.圖6展示了利用掩星自帶的zbalmax的資料求得的邊界層頂高度的季節(jié)分布, 通過與圖5相比較, 可以看出模型函數(shù)法得到的邊界層頂高度略高, 但兩者邊界層頂高度的時空分布十分一致, 說明模型函數(shù)法得到的季節(jié)特征具有準(zhǔn)確性.

圖5 用模型函數(shù)法得到的海洋5年平均的邊界層頂高度(所用資料為2007—2011年1, 4, 7, 10四個月份的掩星彎角的資料)(a) 1月份平均邊界層頂高度; (b) 4月份平均邊界層頂高度; (c) 7月份平均邊界層頂高度; (d) 10月份平均邊界層頂高度Fig.5.The 5-year average boundary layer height of the ocean obtained by the model function method, the data used is the bending angle profile of the four months of 2007?2011 in January, April, July, and October: (a) The average height of the boundary layer in January; (b) the average height of the boundary layer in April; (c) the average height of the boundary layer in July; (d) the average height of the boundary layer in October.

圖6 海洋 5年平均的邊界層頂高度 (所用資料為 2007—2011年1, 4, 7, 10 四個月份的掩星自帶 zbalmax 的資料) (a) 1月份平均邊界層頂高度; (b) 4月份平均邊界層頂高度; (c) 7月份平均邊界層頂高度; (d) 10月份平均邊界層頂高度Fig.6.The 5-year average boundary layer height of the ocean obtained by the model function method, the data used is the zbalmax provided by CDAAC of the four months of 2007?2011 in January, April, July, and October: (a) The average height of the boundary layer in January; (b) the average height of the boundary layer in April; (c) the average height of the boundary layer in July; (d) the average height of the boundary layer in October.

5 總 結(jié)

本文提出了一種確定邊界層頂高度的新的數(shù)值微分方法, 通過模型函數(shù)法確定正則化項中的雙參數(shù), 從而得到彎角梯度, 進(jìn)而利用最大梯度法確定邊界層頂高度.文中展示了兩個數(shù)值實驗的結(jié)果, 從中可知, 此方法最主要的優(yōu)點在于, 盡管掩星資料中帶有誤差, 但其基本不受噪聲的影響, 能夠很好地提取主要的信息, 具有很好的穩(wěn)定性和準(zhǔn)確性.

應(yīng)用本文提出的方法, 基于2007—2011年的COSMIC掩星彎角資料, 研究了全球海洋邊界層頂高度的季節(jié)變化特征, 著重分析了 1, 4, 7, 10 四個月份的年平均邊界層頂高度, 結(jié)果顯示, 隨著太陽直射點的南北移動, 海洋上的邊界層頂高度具有明顯的季節(jié)變化特征, 其與海溫密切相關(guān), 具體來說, 相對于周圍海域的海溫越低, 邊界層頂高度越低, 反之, 邊界層頂高度越高, 延伸這一結(jié)論可以得到, 冷洋流對應(yīng)著低邊界層頂, 暖洋流對應(yīng)著高邊界層頂.為了驗證結(jié)論的準(zhǔn)確性, 將其與利用掩星資料自帶的zbalmax得出來的邊界層頂高度的季節(jié)特征作對比, 結(jié)果發(fā)現(xiàn), 兩者邊界層頂高度的時空分布特征十分一致, 說明通過雙參數(shù)模型函數(shù)法得到的季節(jié)特征具有準(zhǔn)確性.

總之, 經(jīng)過數(shù)值實驗中與差分法和L曲線法的比較, 發(fā)現(xiàn)雙參數(shù)模型函數(shù)法在求邊界層頂高度時具有不可替代的優(yōu)勢, 而且經(jīng)過掩星資料自帶的邊界層頂高度數(shù)據(jù)zbalmax的驗證, 表明文中表述的海上邊界層頂高度的季節(jié)特征具有準(zhǔn)確性.我們下一步需要繼續(xù)做的工作是進(jìn)一步改進(jìn)方法, 使其適用于復(fù)雜天氣條件下的大氣波導(dǎo)的研究, 以期得到更好的結(jié)果.

本文利用的數(shù)據(jù)來自CDAAC提供的GPS掩星資料,在此謹(jǐn)表感謝! 感謝國防科技大學(xué)氣象海洋學(xué)院周則明教授對文章行文結(jié)構(gòu)的建議, 感謝國防科技大學(xué)氣象海洋學(xué)院杜華棟副教授和閆申博士在編程方面的幫助, 感謝國防科技大學(xué)氣象海洋學(xué)院陳毓敏、徐海波和北京師范大學(xué)地理學(xué)院孫銘揚(yáng)對文章提供的修改建議.