• <tr id="yyy80"></tr>
  • <sup id="yyy80"></sup>
  • <tfoot id="yyy80"><noscript id="yyy80"></noscript></tfoot>
  • 99热精品在线国产_美女午夜性视频免费_国产精品国产高清国产av_av欧美777_自拍偷自拍亚洲精品老妇_亚洲熟女精品中文字幕_www日本黄色视频网_国产精品野战在线观看 ?

    人工凍土觀測的凍點對應(yīng)的土壤溫度研究

    2021-07-14 12:41:56孫艷云苗傳海叢郁王若男侯暢韓瑩李琳肖楠舒
    氣象與環(huán)境學(xué)報 2021年3期
    關(guān)鍵詞:凍土土壤溫度淺層

    孫艷云 苗傳海 叢郁 王若男 侯暢 韓瑩 李琳 肖楠舒

    (遼寧省氣象裝備保障中心,遼寧 沈陽 110166)

    引言

    《地面氣象觀測規(guī)范》[1](以下,簡稱《規(guī)范》)指出,凍土是指含有水分的土壤因溫度下降到0 ℃或0 ℃以下而呈凍結(jié)的狀態(tài)。凍土是反映土壤熱狀況的一項指標(biāo),在研究地表和大氣的熱量交換上有重要意義。土壤凍結(jié)深度與農(nóng)事活動、建筑、道路橋梁、鐵路設(shè)計等關(guān)系密切,凍土對氣候變化具有敏感性,是氣候變化的重要指標(biāo)[2-5]。為了研究方便,將人工凍土器內(nèi)管中水相態(tài)發(fā)生變化的點所在的深度稱為凍點。凍點有兩種情況,一種情況是由液態(tài)轉(zhuǎn)化為固態(tài)的點,另一種情況是由固態(tài)轉(zhuǎn)化為液態(tài)的點。凍點表示的是土壤凍結(jié)的上限或下限。

    地面氣象觀測中的凍土觀測始于1954年,至2019年凍土觀測全部為人工觀測,所用的設(shè)備為人工凍土觀測器(人工凍土器),主要設(shè)備型號為TB1-1。自2020年起,凍土自動觀測儀開始在氣象觀測業(yè)務(wù)中應(yīng)用。中國氣象局考核通過的凍土自動觀測儀有測溫式、凍阻式兩種類型[6]。在溫度需要通過玻璃制成的溫度計進(jìn)行人工觀測的年代,人工凍土器非常好的實現(xiàn)了人們動態(tài)掌握土壤熱狀況的需求。研究表明,凍點的溫度與0 ℃相關(guān)但又存在差別[7-9]。王吉興[7]通過溫度插值法對達(dá)尼林凍土器獲得的凍土數(shù)據(jù)進(jìn)行研究指出,凍點對應(yīng)的溫度基本為0 ℃。孟慶勇等[8]結(jié)合地表溫度對凍土數(shù)據(jù)分析指出,08時土壤溫度達(dá)到了-2.4 ~-8.3 ℃而不出現(xiàn)凍結(jié)的情況。陳志超等[9]指出,地面溫度在-2 ℃左右就可能出現(xiàn)凍土現(xiàn)象,而當(dāng)?shù)孛鏈囟鹊陀?5 ℃,5 cm、10 cm 土壤溫度已經(jīng)低于0 ℃,表面土層開始凍結(jié),儀器觀測卻無凍土層。人工凍土觀測的實質(zhì)是什么,凍點對應(yīng)的土壤溫度是多少?是應(yīng)用凍土數(shù)據(jù)時需要解決的問題。由于此前對土壤溫度的測量數(shù)據(jù)少,無法開展精細(xì)的研究。本文利用測溫式凍土自動觀測儀(DTD4,簡稱:測溫式凍土儀)獲得的高時空密度溫度數(shù)據(jù),輔以實驗室測量及氣溫資料等,闡述了地面氣象觀測凍土觀測的實質(zhì),分析了人工凍土觀測的凍點對應(yīng)的溫度,并研究了土壤溫度變化特征。

    1 資料與方法

    1.1 儀器設(shè)備及相關(guān)實驗

    1.1.1 人工凍土器

    人工凍土器由外管和內(nèi)管組成(圖1)。外管為一標(biāo)有0 cm刻度線的硬管,最早為橡膠管,后大多使用PP-P管材,長度為50—450 cm(根據(jù)觀測需要確定長度)。內(nèi)管為一根帶有厘米刻度、壁厚2 mm的橡皮管,長度與外管相匹配,管內(nèi)有固定冰用的鏈子或銅絲、線繩,其結(jié)構(gòu)如圖1。

    圖1 人工凍土器結(jié)構(gòu)示意(a)和橫截面(b)Fig.1 Structure (a) and cross-section (b) of the manual frozen soil observation apparatus

    《氣象觀測規(guī)范》[6]要求,人工凍土器內(nèi)管灌注干凈的河水、井水或自來水至0刻度線處。遼寧省承擔(dān)凍土觀測業(yè)務(wù)的臺站人工凍土器內(nèi)管灌注用水均為臺站自來水。觀測時,一手把凍土器的鐵蓋連同注水軟膠管提到地面之上,用另一只手摸測內(nèi)管冰(包括凍結(jié)得不夠堅實的冰柱)所在位置,從管壁刻度線上讀出冰柱上下兩端的相應(yīng)刻度數(shù)??梢?,人工凍土觀測根據(jù)內(nèi)管水的相態(tài)來判斷土壤的狀態(tài),沒有考慮土壤水分、酸堿度、土質(zhì)等因素對土壤凍結(jié)的影響。而水凍結(jié)的溫度為0 ℃,所以說人工凍土觀測的本質(zhì)是觀測土壤溫度0 ℃點所在的位置。人工凍土器的感應(yīng)部位是裝在內(nèi)管中的水,是一種接觸式溫度感應(yīng)設(shè)備。

    1.1.2 測溫式凍土儀

    測溫式凍土儀通過對垂直方向每1 cm間隔設(shè)置數(shù)字溫度傳感器,每分鐘獲取各溫度數(shù)據(jù),再由內(nèi)置程序提取“凍點”所在的垂直位置,從而得到凍土上、下限及深度數(shù)據(jù)。所用溫度傳感器將溫度與晶振頻率建立線性關(guān)系,從而將溫度轉(zhuǎn)變成數(shù)字信號[10]。測試所用測溫式凍土儀外管管壁為PP-P管,比人工凍土器的外管細(xì)、薄。管內(nèi)放置傳感器電路板,傳感器固化在電路板上,其結(jié)構(gòu)和參數(shù)如圖2和表1所示。

    圖2 測溫式凍土儀外觀(a)和橫截面(b)Fig.2 Appearance (a) and cross-section (b) of a temperature-measuring frozen soil observation apparatus

    1.2 資料來源

    本文選取喀左站(119.77°E、41.11°N)、沈陽站(123.51°E、41.73°N)、遼陽站(123.17°E、41.24°N)、滿洲里站(117.43°E、49.62°N)4個國家級地面氣象站開展對比觀測,觀測時段為2018年12月至2020年6月。滿洲里站安裝了兩套設(shè)備(001設(shè)備和002設(shè)備),為獲取更多數(shù)據(jù),除常規(guī)08:00人工觀測外,部分觀測時段增加了14:00和20:00兩個時次的人工觀測。研究所用設(shè)備主要有人工凍土器和測溫式凍土儀。為保證觀測數(shù)據(jù)的可比性,兩種設(shè)備水平安裝位置直線距離為50 cm。同時結(jié)合國家級地面氣象站觀測數(shù)據(jù)進(jìn)行分析。

    表1 設(shè)備結(jié)構(gòu)參數(shù)比較Table 1 Comparison of equipment structure parameters

    1.3 研究方法

    將凍土凍點區(qū)分為淺層凍點和深層凍點以及凍結(jié)階段凍點和融化階段凍點,將各凍點對應(yīng)的溫度數(shù)據(jù)按不同分類方法分別統(tǒng)計分析。數(shù)據(jù)分析主要運用Excel進(jìn)行統(tǒng)計和對比分析,并繪制圖表。凍點對應(yīng)深度土壤溫度采用自編程序從測溫式凍土自動觀測儀獲取的原始溫度中提取。通過Matlab繪制土壤溫度日等溫線圖及冬季等溫線圖,分析土壤溫度變化特征及其與人工凍土凍點關(guān)系。針對需要進(jìn)一步驗證的內(nèi)容,設(shè)計了實驗室測試。

    2 結(jié)果分析

    2.1 凍點對應(yīng)土壤溫度的分析

    對2018—2019年滿洲里站的融化階段每天08:00、14:00、20:00三個時次和喀左站2018—2019年和2019—2020年冬季每天08:00人工觀測得到的凍點對應(yīng)的土壤溫度數(shù)據(jù)統(tǒng)計分析,并對滿洲里站進(jìn)行了分時次統(tǒng)計,對喀左站2019—2020年冬季凍點區(qū)分凍結(jié)過程和融化過程統(tǒng)計(表2和表3)。其中,凍結(jié)過程是指凍土器內(nèi)管的水由液態(tài)變?yōu)楣虘B(tài)的過程,主要包括第一層、第二層下限加深階段;融化過程是指凍土器內(nèi)管的水由固態(tài)變?yōu)橐簯B(tài)過程,主要包括第一層上限加深,第一層、第二層下限減小階段。通過散點圖(圖3)顯示凍土凍點分布情況,通過3個時次凍點對應(yīng)的溫度曲線圖(圖4)揭示深層凍土凍點對應(yīng)的溫度與觀測時次的關(guān)系。

    表2 滿洲里站凍點對應(yīng)的溫度統(tǒng)計Table 2 Statistics of temperature corresponding to freezing point at Manzhouli station

    表3 喀左站凍點對應(yīng)的溫度統(tǒng)計Table 3 Statistics of temperature corresponding to freezing point at Kazuo station

    圖3 2019年滿洲里春季深層(>35cm)融化階段凍點對應(yīng)的溫度Fig.3 The temperature corresponding to the freezing point when the soil in the deep layer (>35 cm) melted in the spring of 2019 in Manzhouli

    圖4 2019年滿洲里春季融化階段深層分時次凍點對應(yīng)的溫度Fig.4 The temperatures at different times corresponding to freezing point when soil in deep layer melted in spring of 2019 in Manzhouli

    表2和表3淺層數(shù)據(jù)顯示滿洲里三個不同的觀測時次得到的淺層凍土凍點對應(yīng)的溫度差異較大,同一觀測時次獲得的凍土凍點對應(yīng)的溫度極差高達(dá)6.05 ℃;14:00凍點對應(yīng)的溫度最高,平均溫度為1.63 ℃,最大值達(dá)到5.25 ℃,極差也最大,跳躍性變化明顯。淺層凍融交替階段,人工凍土觀測到的凍點對應(yīng)的溫度與0 ℃存在較大差異。淺層土壤溫度變化快,人工凍土器無法快速跟隨土壤溫度變化做出相態(tài)變化,致使淺層出現(xiàn)觀測結(jié)果與實際狀況不相符的情況。

    表2和表3深層數(shù)據(jù)顯示喀左、滿洲里深層凍點對應(yīng)的溫度平均值較淺層更接近于0 ℃,極差較小。從圖3和圖4可知,一定時間段凍土凍點對應(yīng)的溫度與0 ℃同向偏離,凍點溫度跳躍性變化減少,凍點溫度最大值和最小值均低于0—35 cm的凍點溫度最大值和最小值;分三個時次的統(tǒng)計結(jié)果表明,三個時次凍土凍點平均極差均小于淺層凍土凍點的極差,且圖4三個時次凍點對應(yīng)的溫度曲線也無明顯差別,說明三個時次數(shù)據(jù)一致性較好。深層土壤融化的凍點溫度稍高于但更接近于0 ℃;深層的凍結(jié)與融化無日周期變化特征,其凍點可以表征土壤熱狀況??ψ笳緟^(qū)分凍點的凍結(jié)和融化過程統(tǒng)計得到的數(shù)據(jù)顯示,融化過程凍點對應(yīng)的溫度高于凍結(jié)過程凍點對應(yīng)的溫度。

    2.2 測溫式凍土儀0 ℃線與人工凍土觀測的凍結(jié)線的比較

    為了比較測溫式凍土儀的0 ℃點與人工凍土觀測的凍點的關(guān)系,以喀左2018—2019年和2019—2020年兩個冬季數(shù)據(jù)為基礎(chǔ),在測溫式凍土儀獲得的冬季08:00土壤溫度圖上疊加人工凍土凍點(圖5)。以遼陽和沈陽數(shù)據(jù)為基礎(chǔ),將測溫式凍土儀0 ℃所在深度,與人工凍土凍點深度在同一坐標(biāo)中顯示,如圖6。

    從圖5可知,凍結(jié)階段人工凍土凍點落于0 ℃線附近,融化階段,凍點落在0.3 ℃附近。圖5a中A—B階段,0 ℃線與凍點出現(xiàn)較大偏差,對應(yīng)的觀測時間是2018年12月15—22日,測溫式凍土儀探測到0 ℃線上移,而人工凍土凍點持續(xù)。結(jié)合該站其他氣象數(shù)據(jù)(表4),可以確定有氣溫持續(xù)回升過程。表4是喀左站2018年相應(yīng)時段國家級地面氣象站觀測數(shù)據(jù),平均氣溫均高于-6 ℃,最高氣溫達(dá)到了7.2 ℃,持續(xù)天數(shù)為8 d。類似的情況出現(xiàn)在圖5b中,對應(yīng)時間為2019年12月10—20日;同樣出現(xiàn)在圖6a中,對應(yīng)的時間為2018年12月19—22日。在圖5a中C點位置,測溫式凍土儀0 ℃線于1月21日出現(xiàn)最大值,然后減小,而人工凍土觀測的凍點在最大值處持續(xù)9 d,于1月31日開始減小。在圖5a中E點,測溫式凍土儀0 ℃線于2月19日出現(xiàn)最大值,然后減小,而人工凍土觀測的凍點在最大值處持續(xù)6 d,于2月25日開始減小。同樣的情況,在圖5b、圖6a和圖6b中均有體現(xiàn)。

    圖中紅色線為溫度的0 ℃線,藍(lán)色線為人工凍土觀測的凍結(jié)線;藍(lán)色代表低溫區(qū);黃色代表高溫區(qū)圖5 2018—2019年(a)和2019—2020年(b)喀左站冬季土壤溫度與凍點疊加圖Fig.5 Superposition maps of ground temperature and freezing point at Kazuo station in winter during 2018-2019(a) and 2019-2020(b)

    圖6 2018—2019年遼陽(a)、沈陽(b)冬季人工凍土與0 ℃線變化Fig.6 Variation of manual frozen soil and zero degree line at Liaoyang(a) and Kazuo(b) stations in winter from 2018 to 2019

    表4 喀左國家級地面氣象站數(shù)據(jù)Table 4 Data of Kazuo national automatic weather station

    分析得到,凍結(jié)階段人工凍土凍點對應(yīng)的溫度為0 ℃點附近,融化階段人工凍點對應(yīng)的溫度稍高于0 ℃;人工觀測的最大凍點與測溫式凍土儀0 ℃線最深點相近但高于后者;人工凍土觀測的“拐點”出現(xiàn)的時間晚于測溫式凍土儀的“拐點”;從完全融化時間上來看,人工凍土器完全融化的時間晚于測溫式凍土儀完全融化的時間。以上特征出現(xiàn)的根源是人工凍土器的滯后性。

    2.3 人工凍土器內(nèi)管凍結(jié)特征分析

    為探尋人工凍土器內(nèi)管凍結(jié)與所灌注的水是否有關(guān)系,設(shè)計了實驗室測試。所用的設(shè)備主要是恒溫槽,它以酒精為介質(zhì),具備保持槽內(nèi)溫度恒定的能力;溫度標(biāo)準(zhǔn)器為二等標(biāo)準(zhǔn)鉑電阻溫度計。被測物體為灌注來自不同臺站自來水的人工凍土器內(nèi)管。被測物起始狀態(tài)有室溫液態(tài)和一定溫度下的凍結(jié)狀態(tài)。測試時,控制恒溫槽的溫度到預(yù)定的測量點,加入被測物和溫度標(biāo)準(zhǔn)器,持續(xù)一定時間后,識別被測物狀態(tài)并讀取標(biāo)準(zhǔn)器的溫度,各測試點的溫度以標(biāo)準(zhǔn)器讀取的溫度為準(zhǔn)。

    人工凍土器內(nèi)管的實驗室測試獲得19組有效數(shù)據(jù),如表5。由表5可知,灌注不同臺站水的凍土器內(nèi)管在相同的溫度環(huán)境下,凍結(jié)與融化狀態(tài)無明顯區(qū)別;在標(biāo)準(zhǔn)器讀數(shù)為-0.73 ℃時,經(jīng)過12 h內(nèi)管依然沒有凍結(jié),凍結(jié)需要時間長度超乎想象。可以從人工凍土器和測溫式凍土儀結(jié)構(gòu)及原理上分析,如1.1所述,二者本質(zhì)上都是接觸式溫度感應(yīng)設(shè)備,都會存在測溫設(shè)備普遍存在的熱滯效應(yīng)[11],因各自結(jié)構(gòu)的不同,滯后的程度不同。人工凍土觀測器外管管壁更厚,又比測溫式凍土儀多內(nèi)管,這種結(jié)構(gòu)導(dǎo)致了其滯后性大于測溫式凍土儀。同時人工凍土器內(nèi)管需要灌注水,需要水或冰的狀態(tài)發(fā)生變化后,才能被觀測到。水的比熱容高于土壤,而且水是有固定相變溫度的物質(zhì),在發(fā)生相變的時候會釋放(由液態(tài)變?yōu)楣虘B(tài))或者吸收熱量(由固態(tài)變?yōu)橐簯B(tài)),也就是存在相變焓。所以人工凍土器在溫度達(dá)到0 ℃時,還需要持續(xù)的同向溫差將內(nèi)管中水變成冰(或冰變成水)過程所釋放出(或吸收)的熱量抵消,才能完成水(或冰)相態(tài)變化,這是人工凍土器滯后的另一個因素。有研究表明[12],在時間足夠長的情況下,飽水土的凍結(jié)溫度為-0.06 ℃,飽水砂的凍結(jié)溫度為-0.02 ℃。

    表5 人工凍土器內(nèi)管實驗室測量數(shù)據(jù)Table 5 Measurement data of inner pipe of manual frozen soil observation apparatus in the laboratory

    2.4 從土壤溫度日變化特征看淺層、深層凍點對應(yīng)溫度的差異

    2.1中分析得到淺層凍點對應(yīng)的溫度不穩(wěn)定,呈跳躍性變化特征,而深層凍點對應(yīng)的溫度較穩(wěn)定,為了探究其原因,對土壤溫度變化特征進(jìn)行分析。利用測溫式凍土儀獲取的日溫度數(shù)據(jù),分析土壤溫度日變化特征。為使數(shù)據(jù)更全面,選取四個站同一日(2019年4月5日)和沈陽站不同日的兩組土壤溫度數(shù)據(jù)進(jìn)行分析。土壤日溫度數(shù)據(jù)由每10 min間隔的土壤溫度組成。因土壤深度大于60 cm的區(qū)域土壤溫度日變化很小,形成略有傾斜的豎線,為突出淺層變化的特征,截取0—60 cm深度范圍形成圖7和圖8。

    圖7和圖8顯示,淺層土壤受太陽輻射影響,白天溫度升高,夜間溫度降低。在0—25 cm,日土壤溫度周期性變化特征明顯;25—35 cm,日土壤溫度有周期性但幅度減小,而35 cm深度以下,無明顯周期變化特征。龔強等[13]指出,40 cm以下深度土壤溫度幾乎不受地表日變化的影響。本文將淺層與深層的分界點確定為35 cm。在2.3中分析了人工凍土器的滯后性,淺層土壤的周期性變化,使得人工凍土器不能隨著溫度的變化而及時做出相態(tài)的改變,形成了淺層人工凍土凍點對應(yīng)的溫度變化幅度大的情況。而在深層,土壤溫度無周期性且日變化很小,深層凍土凍點對應(yīng)的溫度與0 ℃更接近。

    圖7 沈陽站1月5日(a)、2月5日(b)、3月5日(c)和4月5日(d)土壤溫度變化Fig.7 Variations of soil temperature at Shenyang station on January 5,February 5,March 5 and April 5,2019

    圖8 2019年4月5日喀左(a)、遼陽(b)、沈陽(c)、滿洲里(d)站日土壤溫度變化Fig.8 Variations of soil temperature at Kazuo,Shenyang,Liaoyang,Manzhouli stations on April 5,2019

    由圖7和圖8可知,在土壤溫度下降階段,夜間溫度降低的影響大于白天溫度升高的影響,如圖7a和圖7b。而在氣溫升高階段,白天溫度升高的影響大于夜間溫度降低的影響,如圖7c、圖7d、圖8a、圖8b、圖8c和圖8d。圖7中,不同時間的日數(shù)據(jù)等溫線密度不同,也就是土壤溫度變化快慢不同,圖7a和圖7c等溫線密度大,溫度變化快,而圖7b和圖7d等溫線密度小,溫度較穩(wěn)定。由日平均土壤溫度形成的圖可以展示年平均土壤溫度情況,以此為基礎(chǔ)可分析氣候變化等。

    3 結(jié)論與討論

    (1)淺層凍融交替階段,人工凍土觀測到的凍點對應(yīng)的溫度與0 ℃存在較大差異且跳躍性變化明顯。深層凍點對應(yīng)的溫度平均值較淺層更接近于0 ℃,差值較小;一定時間段凍土凍點對應(yīng)的溫度與0 ℃同向偏離。人工觀測的最深凍點與測溫式凍土儀0 ℃線最深點相近但大于后者;人工凍土觀測的“拐點”出現(xiàn)的時間晚于測溫式凍土儀的“拐點”;從完全融化時間上來看,人工凍土器完全融化的時間晚于測溫式凍土儀完全融化的時間。

    (2)人工凍土觀測的實質(zhì)是獲得土壤溫度0 ℃點所在位置。人工凍土器與測溫式凍土儀相比,外管管壁更厚,有2 mm厚的內(nèi)管,內(nèi)管需要灌注水,水的比熱大于土壤,且水在發(fā)生相變時需持續(xù)的同向溫差消除相變焓,導(dǎo)致了其滯后性大于測溫式凍土儀。土壤溫度0—25cm深度日周期性變化特征明顯,25—35 cm有日周期性但變化幅度減小,而35 cm深度以下,無周期性變化特征。淺層凍點對應(yīng)的溫度范圍大,主要是人工凍土器的滯后性和淺層土壤(0—35 cm深度范圍)溫度的日周期變化性共同作用使得人工凍土器內(nèi)管中的水無法快速跟隨土壤溫度的周期性變化及時改變相態(tài)。

    (3)測溫式凍土儀除了獲取凍土觀測數(shù)據(jù)外,其獲得的高時空密度土壤溫度數(shù)據(jù)具有重要的科學(xué)研究價值。地球系統(tǒng)各圈層間的作用和制約過程都伴隨著能量傳遞,因而與溫度有著密切聯(lián)系,長期地溫監(jiān)測在地球科學(xué)研究中有重要意義[17]。溫度是土壤呼吸速率時間動態(tài)的主要影響因素[18]。測溫式凍土自動觀測儀首次實現(xiàn)了對每厘米深度土壤溫度的精密探測,精細(xì)展現(xiàn)了土壤溫度變化過程,即使在沒有凍土的低緯度地區(qū),高時空密度的土壤溫度數(shù)據(jù)也是氣候研究、農(nóng)業(yè)服務(wù)急需的資料,可為研究土壤溫度變化特征及土壤與大氣熱交換提供充足的數(shù)據(jù)支撐。因此,有必要將測溫式凍土自動觀測儀獲得的高時空密度土壤溫度業(yè)務(wù)化應(yīng)用。

    (4)測溫式凍土儀在正式生產(chǎn)時,增加了外徑40 mm,內(nèi)徑27 mm,壁厚6.5 mm的PPR管作為外管,其對觀測結(jié)果的影響有待于進(jìn)一步研究。但作為接觸式測溫設(shè)備,減小外因產(chǎn)生的時滯是提高其靈敏度的重要環(huán)節(jié)。建議測溫式凍土儀的外管壁使用溫度滯后效應(yīng)更小的金屬外管。

    猜你喜歡
    凍土土壤溫度淺層
    淺層換填技術(shù)在深厚軟土路基中的應(yīng)用
    基于淺層曝氣原理的好氧顆粒污泥的快速培養(yǎng)
    北極凍土在求救
    凍土下的猛犸墳場
    遼東山區(qū)3種人工林土壤呼吸對土壤溫度和土壤水分的響應(yīng)
    淺層地下水超采區(qū)劃分探究
    管群間歇散熱的土壤溫度響應(yīng)與恢復(fù)特性
    華北地質(zhì)(2015年3期)2015-12-04 06:13:29
    土壤濕度和土壤溫度模擬中的參數(shù)敏感性分析和優(yōu)化
    26
    欧美日韩国产亚洲二区| 国产免费av片在线观看野外av| 黄色视频,在线免费观看| 午夜福利视频1000在线观看| 午夜福利成人在线免费观看| 日韩精品中文字幕看吧| 美女大奶头视频| 国产成年人精品一区二区| 亚洲中文av在线| 亚洲精品久久国产高清桃花| 法律面前人人平等表现在哪些方面| 又黄又粗又硬又大视频| 亚洲国产精品999在线| 久久久久久久久久黄片| 欧美在线一区亚洲| 国内久久婷婷六月综合欲色啪| 国产午夜精品论理片| 又粗又爽又猛毛片免费看| www.999成人在线观看| 亚洲七黄色美女视频| 午夜精品久久久久久毛片777| 国产精品一区二区免费欧美| 国产一级毛片七仙女欲春2| 欧美午夜高清在线| 久久性视频一级片| 亚洲av成人精品一区久久| 国产精品1区2区在线观看.| 亚洲欧美日韩无卡精品| 精品免费久久久久久久清纯| 欧美另类亚洲清纯唯美| 成年版毛片免费区| 久久天堂一区二区三区四区| 国产亚洲av嫩草精品影院| 久久久久国产一级毛片高清牌| 精品久久久久久久人妻蜜臀av| 99国产精品一区二区三区| 嫩草影视91久久| 日韩三级视频一区二区三区| 久久久精品国产亚洲av高清涩受| 丁香六月欧美| 日韩欧美 国产精品| 天堂√8在线中文| 青草久久国产| 窝窝影院91人妻| 人妻夜夜爽99麻豆av| 国产成年人精品一区二区| 男人的好看免费观看在线视频 | 老熟妇仑乱视频hdxx| 亚洲男人天堂网一区| 国产成人精品久久二区二区91| 女人被狂操c到高潮| 国产激情欧美一区二区| 男人舔奶头视频| 久久久久久久午夜电影| 国产区一区二久久| 国产精品亚洲美女久久久| 亚洲乱码一区二区免费版| 欧美日韩亚洲综合一区二区三区_| 精品一区二区三区四区五区乱码| 欧美乱色亚洲激情| 久久精品夜夜夜夜夜久久蜜豆 | 一区二区三区国产精品乱码| 日本免费一区二区三区高清不卡| 久久香蕉国产精品| 久久这里只有精品19| 一级毛片女人18水好多| 亚洲av成人不卡在线观看播放网| 黄频高清免费视频| 亚洲国产中文字幕在线视频| 亚洲午夜精品一区,二区,三区| 国产又色又爽无遮挡免费看| 一边摸一边抽搐一进一小说| 日韩大尺度精品在线看网址| 悠悠久久av| 变态另类丝袜制服| 日本一区二区免费在线视频| 中文字幕高清在线视频| 少妇人妻一区二区三区视频| 午夜亚洲福利在线播放| 级片在线观看| 国产精品综合久久久久久久免费| 日本免费a在线| 久久这里只有精品19| 国产成人精品久久二区二区免费| 亚洲专区中文字幕在线| 九九热线精品视视频播放| 亚洲欧美激情综合另类| 欧美激情久久久久久爽电影| 国产精品,欧美在线| 男插女下体视频免费在线播放| 99热这里只有是精品50| 久久婷婷成人综合色麻豆| 女生性感内裤真人,穿戴方法视频| 精品熟女少妇八av免费久了| 熟妇人妻久久中文字幕3abv| 他把我摸到了高潮在线观看| 国产精品美女特级片免费视频播放器 | 黄片大片在线免费观看| 久久久久久久久免费视频了| 身体一侧抽搐| 国产精品精品国产色婷婷| 免费看日本二区| 亚洲国产精品合色在线| 亚洲激情在线av| 禁无遮挡网站| 脱女人内裤的视频| 亚洲国产精品999在线| 69av精品久久久久久| 国产久久久一区二区三区| 午夜两性在线视频| 亚洲熟妇中文字幕五十中出| 久久精品国产99精品国产亚洲性色| 久久久久国产精品人妻aⅴ院| 草草在线视频免费看| 在线观看美女被高潮喷水网站 | 亚洲真实伦在线观看| 搡老妇女老女人老熟妇| 久久国产精品人妻蜜桃| 人妻夜夜爽99麻豆av| 90打野战视频偷拍视频| 国产69精品久久久久777片 | 18禁黄网站禁片免费观看直播| 国产成人精品久久二区二区91| 日本在线视频免费播放| 日本黄色视频三级网站网址| 最近视频中文字幕2019在线8| 麻豆国产av国片精品| 麻豆一二三区av精品| av天堂在线播放| 亚洲七黄色美女视频| 黄色丝袜av网址大全| 老汉色∧v一级毛片| 欧美性猛交黑人性爽| 久久久国产成人精品二区| 亚洲一区中文字幕在线| 国产欧美日韩精品亚洲av| 亚洲欧美激情综合另类| 成人高潮视频无遮挡免费网站| 国产精品久久久久久精品电影| 两性夫妻黄色片| 亚洲av中文字字幕乱码综合| 国产在线观看jvid| 久久久国产成人免费| 中文资源天堂在线| 国产精品一及| 国产视频内射| 日本黄色视频三级网站网址| 性色av乱码一区二区三区2| 日本在线视频免费播放| 神马国产精品三级电影在线观看 | 国产精品久久久久久久电影 | 美女高潮喷水抽搐中文字幕| 久久久国产成人精品二区| 国产精品99久久99久久久不卡| 日本精品一区二区三区蜜桃| 人妻久久中文字幕网| 亚洲 欧美 日韩 在线 免费| 日韩欧美免费精品| 国产成人av激情在线播放| 999久久久国产精品视频| 免费看美女性在线毛片视频| 亚洲天堂国产精品一区在线| 亚洲美女黄片视频| 一本久久中文字幕| 亚洲成av人片在线播放无| avwww免费| 一本一本综合久久| 欧美日韩瑟瑟在线播放| 国产高清视频在线观看网站| 欧美色欧美亚洲另类二区| 欧美激情久久久久久爽电影| 九九热线精品视视频播放| 欧美极品一区二区三区四区| 午夜福利高清视频| 一个人免费在线观看电影 | 亚洲一区中文字幕在线| 午夜福利免费观看在线| 亚洲第一电影网av| 欧美日本亚洲视频在线播放| 亚洲精品久久国产高清桃花| xxx96com| 国产精品香港三级国产av潘金莲| 国内精品久久久久久久电影| 国产精品98久久久久久宅男小说| videosex国产| 国产aⅴ精品一区二区三区波| 亚洲av成人不卡在线观看播放网| 禁无遮挡网站| 国产精品一区二区三区四区免费观看 | 在线观看www视频免费| 中文字幕熟女人妻在线| 淫妇啪啪啪对白视频| 精品人妻1区二区| a级毛片在线看网站| 看免费av毛片| 他把我摸到了高潮在线观看| 成年女人毛片免费观看观看9| 欧美三级亚洲精品| 啪啪无遮挡十八禁网站| 亚洲成av人片在线播放无| 黄色女人牲交| 亚洲一区高清亚洲精品| 午夜视频精品福利| 一级片免费观看大全| 啦啦啦韩国在线观看视频| 一级a爱片免费观看的视频| 午夜精品在线福利| 成年版毛片免费区| 亚洲中文av在线| 亚洲国产精品久久男人天堂| 此物有八面人人有两片| 最近最新免费中文字幕在线| 在线观看舔阴道视频| 久久99热这里只有精品18| 欧美日韩福利视频一区二区| 亚洲人成网站在线播放欧美日韩| 搡老妇女老女人老熟妇| 亚洲中文日韩欧美视频| 岛国视频午夜一区免费看| 国产69精品久久久久777片 | 舔av片在线| 中文在线观看免费www的网站 | svipshipincom国产片| www.自偷自拍.com| 亚洲欧美日韩高清专用| 国产亚洲av高清不卡| 中文字幕熟女人妻在线| 99精品在免费线老司机午夜| 欧美午夜高清在线| 亚洲欧美激情综合另类| 久久精品91蜜桃| 波多野结衣高清作品| 给我免费播放毛片高清在线观看| 国产伦一二天堂av在线观看| 啦啦啦免费观看视频1| 亚洲国产欧美一区二区综合| 天天一区二区日本电影三级| 最近视频中文字幕2019在线8| 最近最新中文字幕大全电影3| 嫩草影视91久久| 丝袜人妻中文字幕| 曰老女人黄片| 一个人免费在线观看电影 | 18禁观看日本| 91国产中文字幕| 亚洲va日本ⅴa欧美va伊人久久| 亚洲精品美女久久av网站| 色综合欧美亚洲国产小说| 成人av一区二区三区在线看| 99riav亚洲国产免费| 在线国产一区二区在线| 国产av不卡久久| 在线看三级毛片| 国产99白浆流出| 国产成人精品久久二区二区91| 亚洲av成人精品一区久久| 午夜福利欧美成人| 国产精品亚洲av一区麻豆| 我的老师免费观看完整版| 母亲3免费完整高清在线观看| 亚洲av中文字字幕乱码综合| 美女高潮喷水抽搐中文字幕| 国产激情久久老熟女| 午夜免费成人在线视频| 国产成人aa在线观看| 久久国产精品影院| 99久久无色码亚洲精品果冻| 亚洲成av人片在线播放无| 亚洲18禁久久av| 欧美极品一区二区三区四区| 日韩欧美国产一区二区入口| 国产精品美女特级片免费视频播放器 | 午夜视频精品福利| 给我免费播放毛片高清在线观看| 一边摸一边做爽爽视频免费| 免费人成视频x8x8入口观看| 国产免费男女视频| 国产精品自产拍在线观看55亚洲| 中文字幕av在线有码专区| 亚洲激情在线av| 特大巨黑吊av在线直播| 成人三级黄色视频| e午夜精品久久久久久久| 亚洲精品中文字幕在线视频| 亚洲精品久久国产高清桃花| 变态另类成人亚洲欧美熟女| 99久久无色码亚洲精品果冻| 激情在线观看视频在线高清| 国产成人精品无人区| 中文字幕av在线有码专区| 亚洲成人国产一区在线观看| 校园春色视频在线观看| 国产在线观看jvid| 成年免费大片在线观看| 久久欧美精品欧美久久欧美| 国产成人一区二区三区免费视频网站| 国产爱豆传媒在线观看 | 精品人妻1区二区| 午夜精品久久久久久毛片777| 欧美成人午夜精品| 亚洲五月天丁香| 日本熟妇午夜| 国产亚洲精品久久久久5区| 亚洲精品av麻豆狂野| 99热只有精品国产| www.熟女人妻精品国产| 国产又色又爽无遮挡免费看| 亚洲av电影不卡..在线观看| 成人三级黄色视频| 狂野欧美激情性xxxx| 嫩草影视91久久| 久9热在线精品视频| 一级黄色大片毛片| 欧美黑人巨大hd| 少妇熟女aⅴ在线视频| 欧美性猛交黑人性爽| 啦啦啦观看免费观看视频高清| 18禁国产床啪视频网站| 午夜久久久久精精品| 国产不卡一卡二| 欧美中文日本在线观看视频| 少妇人妻一区二区三区视频| www国产在线视频色| 亚洲精品一区av在线观看| 国产成人aa在线观看| 久久婷婷成人综合色麻豆| 午夜免费激情av| 黄色成人免费大全| 一边摸一边抽搐一进一小说| 国产精品亚洲一级av第二区| 嫁个100分男人电影在线观看| 丰满的人妻完整版| 黄色丝袜av网址大全| 国内揄拍国产精品人妻在线| 蜜桃久久精品国产亚洲av| 在线视频色国产色| 少妇的丰满在线观看| 免费看a级黄色片| 久久久精品国产亚洲av高清涩受| 神马国产精品三级电影在线观看 | 亚洲成人精品中文字幕电影| aaaaa片日本免费| 日韩免费av在线播放| 久久香蕉激情| 老熟妇仑乱视频hdxx| 亚洲欧美日韩高清专用| 国产高清视频在线播放一区| 亚洲人与动物交配视频| 欧美+亚洲+日韩+国产| 国产精品美女特级片免费视频播放器 | 国产亚洲精品第一综合不卡| 女人被狂操c到高潮| 亚洲美女黄片视频| 成人国语在线视频| 免费看a级黄色片| 久久久久性生活片| 无遮挡黄片免费观看| 亚洲欧美激情综合另类| 日韩成人在线观看一区二区三区| 亚洲成a人片在线一区二区| 婷婷精品国产亚洲av在线| 亚洲人与动物交配视频| 好男人电影高清在线观看| 一进一出抽搐动态| 欧美又色又爽又黄视频| 亚洲免费av在线视频| 精品人妻1区二区| 亚洲欧美激情综合另类| 欧美av亚洲av综合av国产av| 亚洲av电影不卡..在线观看| 久久精品国产亚洲av香蕉五月| 亚洲av电影不卡..在线观看| 午夜福利欧美成人| www.熟女人妻精品国产| 88av欧美| 久久 成人 亚洲| 久久精品综合一区二区三区| 国产激情偷乱视频一区二区| 最近视频中文字幕2019在线8| 国产亚洲精品av在线| 极品教师在线免费播放| 欧美精品啪啪一区二区三区| 人成视频在线观看免费观看| 黑人巨大精品欧美一区二区mp4| 久久久久国内视频| 身体一侧抽搐| 国产熟女午夜一区二区三区| 日韩欧美国产一区二区入口| 久久精品国产99精品国产亚洲性色| 少妇人妻一区二区三区视频| 国产精品久久电影中文字幕| 国产精品久久久久久精品电影| 级片在线观看| 长腿黑丝高跟| 一区二区三区国产精品乱码| 免费观看人在逋| www日本在线高清视频| 岛国在线免费视频观看| 亚洲欧美精品综合一区二区三区| 成人手机av| 国产99白浆流出| 欧美日韩一级在线毛片| 99国产极品粉嫩在线观看| 欧美色欧美亚洲另类二区| 亚洲自拍偷在线| 精品欧美一区二区三区在线| 国产精品久久久人人做人人爽| 国产高清videossex| 免费电影在线观看免费观看| 国产精品电影一区二区三区| 午夜福利免费观看在线| 不卡av一区二区三区| 人成视频在线观看免费观看| 亚洲成人免费电影在线观看| 久9热在线精品视频| 91字幕亚洲| av片东京热男人的天堂| 欧美在线黄色| 天堂影院成人在线观看| 久久久水蜜桃国产精品网| 日日夜夜操网爽| 亚洲av电影在线进入| 欧美乱色亚洲激情| 中亚洲国语对白在线视频| 欧美日韩一级在线毛片| 成在线人永久免费视频| 亚洲免费av在线视频| 此物有八面人人有两片| 麻豆成人午夜福利视频| 国产99白浆流出| 日本 欧美在线| 欧美 亚洲 国产 日韩一| 视频区欧美日本亚洲| 最近最新中文字幕大全电影3| www日本黄色视频网| 99热6这里只有精品| 三级国产精品欧美在线观看 | 亚洲九九香蕉| 精品久久久久久久人妻蜜臀av| 欧美黑人欧美精品刺激| 蜜桃久久精品国产亚洲av| 亚洲美女黄片视频| 18禁黄网站禁片午夜丰满| 精品一区二区三区视频在线观看免费| 国模一区二区三区四区视频 | 亚洲成人久久性| 又大又爽又粗| 又黄又爽又免费观看的视频| 亚洲九九香蕉| 亚洲中文日韩欧美视频| 久久久水蜜桃国产精品网| 最近最新中文字幕大全电影3| 久久亚洲精品不卡| www国产在线视频色| 日本五十路高清| 99久久综合精品五月天人人| 黄片小视频在线播放| 日日爽夜夜爽网站| 国产伦人伦偷精品视频| 精品不卡国产一区二区三区| 18禁黄网站禁片午夜丰满| 成人三级黄色视频| 最新美女视频免费是黄的| 国产三级黄色录像| 欧美性长视频在线观看| 日韩大码丰满熟妇| 好男人电影高清在线观看| 免费看日本二区| 香蕉av资源在线| 国产av麻豆久久久久久久| 国产亚洲精品久久久久5区| 亚洲专区国产一区二区| 欧美日韩福利视频一区二区| 巨乳人妻的诱惑在线观看| 91大片在线观看| 色综合站精品国产| 成人高潮视频无遮挡免费网站| 18禁裸乳无遮挡免费网站照片| 亚洲人成网站在线播放欧美日韩| 午夜激情福利司机影院| 久久精品夜夜夜夜夜久久蜜豆 | 99久久精品热视频| 免费观看精品视频网站| 免费观看人在逋| 国产视频一区二区在线看| 亚洲人与动物交配视频| av在线播放免费不卡| 精品一区二区三区av网在线观看| 丝袜美腿诱惑在线| 国产亚洲av高清不卡| bbb黄色大片| 啦啦啦免费观看视频1| 男女那种视频在线观看| 久久中文字幕人妻熟女| 在线观看免费视频日本深夜| 一进一出抽搐动态| avwww免费| 精品国内亚洲2022精品成人| 搡老妇女老女人老熟妇| 精品一区二区三区视频在线观看免费| 桃色一区二区三区在线观看| 少妇被粗大的猛进出69影院| 精品高清国产在线一区| 国产精品 国内视频| 国产精品综合久久久久久久免费| 成人三级做爰电影| 五月玫瑰六月丁香| 欧美乱色亚洲激情| 青草久久国产| 少妇人妻一区二区三区视频| 国产在线精品亚洲第一网站| 香蕉丝袜av| 亚洲成人免费电影在线观看| xxxwww97欧美| 免费av毛片视频| 波多野结衣巨乳人妻| av视频在线观看入口| 老汉色∧v一级毛片| www国产在线视频色| 99热只有精品国产| 一本大道久久a久久精品| 亚洲五月婷婷丁香| 亚洲成人精品中文字幕电影| 亚洲av成人精品一区久久| 久久人妻福利社区极品人妻图片| 夜夜看夜夜爽夜夜摸| 桃色一区二区三区在线观看| 国产精品,欧美在线| 久久久久免费精品人妻一区二区| 1024手机看黄色片| 巨乳人妻的诱惑在线观看| 国内精品久久久久精免费| 18禁美女被吸乳视频| 国产午夜福利久久久久久| 久久久国产成人精品二区| 88av欧美| 亚洲欧美日韩东京热| 国产熟女午夜一区二区三区| 一a级毛片在线观看| 91大片在线观看| 欧美日韩亚洲综合一区二区三区_| 亚洲熟女毛片儿| 国产69精品久久久久777片 | 久久香蕉激情| 午夜a级毛片| 一区二区三区国产精品乱码| 国模一区二区三区四区视频 | 黄频高清免费视频| 三级国产精品欧美在线观看 | 黄色成人免费大全| 老司机午夜十八禁免费视频| 国产精品免费一区二区三区在线| 又黄又粗又硬又大视频| 国产欧美日韩一区二区精品| 一a级毛片在线观看| 99国产精品99久久久久| 在线观看午夜福利视频| 国产亚洲精品第一综合不卡| 淫妇啪啪啪对白视频| 国产精品精品国产色婷婷| 亚洲国产高清在线一区二区三| www国产在线视频色| 亚洲精品美女久久av网站| 白带黄色成豆腐渣| 久久久精品国产亚洲av高清涩受| 亚洲激情在线av| 一本大道久久a久久精品| 最近在线观看免费完整版| 国产一区在线观看成人免费| 国产精品野战在线观看| 午夜福利在线观看吧| 婷婷精品国产亚洲av在线| 91大片在线观看| 露出奶头的视频| 成人高潮视频无遮挡免费网站| 久久天躁狠狠躁夜夜2o2o| 中文亚洲av片在线观看爽| 亚洲成人久久爱视频| 一级a爱片免费观看的视频| 蜜桃久久精品国产亚洲av| 国产蜜桃级精品一区二区三区| 国产成人av激情在线播放| 听说在线观看完整版免费高清| 国产精品国产高清国产av| 香蕉国产在线看| 女人高潮潮喷娇喘18禁视频| 色尼玛亚洲综合影院| 国产av一区二区精品久久| 亚洲精品av麻豆狂野| 精品免费久久久久久久清纯| 午夜两性在线视频| 麻豆国产av国片精品| 丁香欧美五月| 国产成人av激情在线播放| 国内毛片毛片毛片毛片毛片| 在线观看免费视频日本深夜| 国产真实乱freesex| 亚洲自拍偷在线| 观看免费一级毛片| 亚洲av电影在线进入| 国产精品一区二区精品视频观看| 后天国语完整版免费观看| 村上凉子中文字幕在线| 国产又黄又爽又无遮挡在线| 露出奶头的视频| 亚洲午夜精品一区,二区,三区| 午夜精品在线福利| 国产亚洲精品第一综合不卡| 正在播放国产对白刺激| 久久久久久大精品| 色老头精品视频在线观看| 久久性视频一级片| 神马国产精品三级电影在线观看 | 国产男靠女视频免费网站| 成熟少妇高潮喷水视频|