• <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亚洲精品天堂| 免费在线观看黄色视频的| 天天添夜夜摸| 他把我摸到了高潮在线观看 | 考比视频在线观看| 午夜免费成人在线视频| 五月天丁香电影| 啦啦啦在线免费观看视频4| 亚洲第一欧美日韩一区二区三区 | 人人澡人人妻人| 大片电影免费在线观看免费| 国产精品秋霞免费鲁丝片| 国产成人影院久久av| 叶爱在线成人免费视频播放| 老汉色av国产亚洲站长工具| 亚洲精品久久成人aⅴ小说| 美女视频免费永久观看网站| 国产一卡二卡三卡精品| 男男h啪啪无遮挡| 黑人猛操日本美女一级片| 亚洲人成电影免费在线| 中国美女看黄片| 国产亚洲精品久久久久5区| 国产成人精品在线电影| 亚洲成av片中文字幕在线观看| 在线观看一区二区三区激情| 精品一区二区三区视频在线观看免费 | 午夜91福利影院| 怎么达到女性高潮| 久久人妻av系列| 熟女少妇亚洲综合色aaa.| 女人精品久久久久毛片| 久久这里只有精品19| 波多野结衣一区麻豆| 一边摸一边抽搐一进一小说 | 高清黄色对白视频在线免费看| 性高湖久久久久久久久免费观看| 成在线人永久免费视频| 免费观看av网站的网址| 汤姆久久久久久久影院中文字幕| 日本vs欧美在线观看视频| 亚洲欧洲精品一区二区精品久久久| 一级片'在线观看视频| 在线观看免费午夜福利视频| 日韩人妻精品一区2区三区| 一级,二级,三级黄色视频| 男女边摸边吃奶| 中亚洲国语对白在线视频| 中亚洲国语对白在线视频| 国产精品 国内视频| 大型av网站在线播放| 性色av乱码一区二区三区2| 久久人妻福利社区极品人妻图片| 欧美日韩亚洲综合一区二区三区_| 精品久久久久久久毛片微露脸| 日韩熟女老妇一区二区性免费视频| 十八禁网站免费在线| 一进一出好大好爽视频| h视频一区二区三区| 老汉色av国产亚洲站长工具| 久久亚洲真实| 视频区图区小说| 日韩有码中文字幕| 一区在线观看完整版| 国产成+人综合+亚洲专区| 黄色视频在线播放观看不卡| 色婷婷久久久亚洲欧美| 五月天丁香电影| 国产精品国产高清国产av | 国产精品98久久久久久宅男小说| 国产无遮挡羞羞视频在线观看| 久久狼人影院| 狠狠婷婷综合久久久久久88av| 操出白浆在线播放| 黄片播放在线免费| 欧美久久黑人一区二区| 欧美激情久久久久久爽电影 | 成人免费观看视频高清| 丝袜人妻中文字幕| 狂野欧美激情性xxxx| 久久99热这里只频精品6学生| 国产1区2区3区精品| 丰满饥渴人妻一区二区三| 亚洲精品久久成人aⅴ小说| 肉色欧美久久久久久久蜜桃| 亚洲一卡2卡3卡4卡5卡精品中文| 成人三级做爰电影| 久久久久久久久免费视频了| 国产亚洲精品一区二区www | 女性被躁到高潮视频| 老汉色∧v一级毛片| 国产xxxxx性猛交| 欧美黄色片欧美黄色片| 麻豆av在线久日| 一本—道久久a久久精品蜜桃钙片| 91老司机精品| 国产成+人综合+亚洲专区| 自线自在国产av| 99在线人妻在线中文字幕 | 一个人免费在线观看的高清视频| 男女边摸边吃奶| 精品久久久久久电影网| 久久久国产一区二区| 久久国产亚洲av麻豆专区| 9热在线视频观看99| 亚洲精品在线观看二区| 天堂动漫精品| 纵有疾风起免费观看全集完整版| 午夜老司机福利片| 成人影院久久| 精品国产国语对白av| 一边摸一边做爽爽视频免费| 免费看a级黄色片| 日韩三级视频一区二区三区| 国产麻豆69| 丰满少妇做爰视频| 一区二区日韩欧美中文字幕| 色精品久久人妻99蜜桃| 亚洲情色 制服丝袜| 蜜桃国产av成人99| 91成年电影在线观看| 99久久精品国产亚洲精品| 久久精品亚洲熟妇少妇任你| 亚洲第一青青草原| 国产成人影院久久av| 国产99久久九九免费精品| 国产无遮挡羞羞视频在线观看| 五月开心婷婷网| 桃花免费在线播放| 少妇 在线观看| a级片在线免费高清观看视频| 如日韩欧美国产精品一区二区三区| 精品一区二区三区av网在线观看 | 俄罗斯特黄特色一大片| 老司机午夜福利在线观看视频 | 黄色视频,在线免费观看| 久久精品国产99精品国产亚洲性色 | 999久久久国产精品视频| 丝袜美腿诱惑在线| 成人18禁在线播放| 老熟女久久久| 日韩大码丰满熟妇| 亚洲av片天天在线观看| 亚洲视频免费观看视频| 久久这里只有精品19| 免费在线观看视频国产中文字幕亚洲| 国产一卡二卡三卡精品| 三级毛片av免费| 一本大道久久a久久精品| 99re在线观看精品视频| 欧美日韩视频精品一区| 精品卡一卡二卡四卡免费| 欧美在线黄色| 99热国产这里只有精品6| 亚洲欧美一区二区三区久久| 亚洲男人天堂网一区| 国产国语露脸激情在线看| 一区二区三区乱码不卡18| 日韩 欧美 亚洲 中文字幕| 又大又爽又粗| 女人爽到高潮嗷嗷叫在线视频| 国产高清激情床上av| 大香蕉久久成人网| 高潮久久久久久久久久久不卡| 丝瓜视频免费看黄片| 黄频高清免费视频| av一本久久久久| 超碰97精品在线观看| videosex国产| 制服人妻中文乱码| 91成人精品电影| 一级片免费观看大全| 国产成人av激情在线播放| 丝瓜视频免费看黄片| 亚洲av欧美aⅴ国产| 人人妻人人添人人爽欧美一区卜| 老熟女久久久| 青草久久国产| 在线十欧美十亚洲十日本专区| 亚洲第一青青草原| 欧美成人免费av一区二区三区 | 91字幕亚洲| 国产福利在线免费观看视频| 中文字幕av电影在线播放| 色在线成人网| 国产麻豆69| 欧美日韩精品网址| 午夜福利在线免费观看网站| 91老司机精品| 国产主播在线观看一区二区| 两性夫妻黄色片| 精品国产国语对白av| 天天影视国产精品| 国产成人av教育| 国产老妇伦熟女老妇高清| 男女床上黄色一级片免费看| 一本色道久久久久久精品综合| 国产精品免费视频内射| 亚洲av电影在线进入| 亚洲成a人片在线一区二区| 久热爱精品视频在线9| 国产在线免费精品| 自拍欧美九色日韩亚洲蝌蚪91| 成人特级黄色片久久久久久久 | 两个人看的免费小视频| 久久精品国产a三级三级三级| 午夜福利在线免费观看网站| 黄色怎么调成土黄色| av线在线观看网站| 一级片免费观看大全| 成人精品一区二区免费| 少妇被粗大的猛进出69影院| 精品久久蜜臀av无| 亚洲 国产 在线| 天天添夜夜摸| 亚洲国产av新网站| 一本大道久久a久久精品| 免费少妇av软件| 久久国产精品大桥未久av| 久久久国产一区二区| 老熟妇仑乱视频hdxx| 精品免费久久久久久久清纯 | 黄色成人免费大全| 亚洲av电影在线进入| 国产片内射在线| 999久久久国产精品视频| 性色av乱码一区二区三区2| 日本av手机在线免费观看| 搡老岳熟女国产| 别揉我奶头~嗯~啊~动态视频| 亚洲一区中文字幕在线| 午夜福利在线免费观看网站| 黄色a级毛片大全视频| 国产精品麻豆人妻色哟哟久久| 这个男人来自地球电影免费观看| 91九色精品人成在线观看| 亚洲成人国产一区在线观看| 国产成人啪精品午夜网站| 亚洲专区中文字幕在线| 欧美精品一区二区免费开放| 亚洲色图av天堂| 久久久水蜜桃国产精品网| 一本色道久久久久久精品综合| 亚洲欧美日韩另类电影网站| 国产一区二区三区在线臀色熟女 | 成人国产av品久久久| 日韩大码丰满熟妇| 免费在线观看视频国产中文字幕亚洲| 亚洲人成77777在线视频| 欧美久久黑人一区二区| 在线永久观看黄色视频| 18禁裸乳无遮挡动漫免费视频| 悠悠久久av| 亚洲黑人精品在线| 一本久久精品| 国产免费现黄频在线看| 高潮久久久久久久久久久不卡| 精品视频人人做人人爽| aaaaa片日本免费| 一边摸一边做爽爽视频免费| 久久99一区二区三区| 免费一级毛片在线播放高清视频 | 久久 成人 亚洲| 亚洲视频免费观看视频| 精品久久蜜臀av无| 丁香六月天网| 成人黄色视频免费在线看| 午夜激情av网站| 日韩视频一区二区在线观看| 日日夜夜操网爽| 欧美黑人精品巨大| 大码成人一级视频| 国产日韩欧美亚洲二区| 18禁黄网站禁片午夜丰满| 黑人巨大精品欧美一区二区mp4| 久久中文字幕一级| 9191精品国产免费久久| 90打野战视频偷拍视频| 窝窝影院91人妻| 宅男免费午夜| 免费av中文字幕在线| 天天躁狠狠躁夜夜躁狠狠躁| 亚洲人成伊人成综合网2020| www.自偷自拍.com| 欧美人与性动交α欧美精品济南到| av免费在线观看网站| 成在线人永久免费视频| 色播在线永久视频| 国产免费现黄频在线看| 午夜福利视频精品| 精品少妇黑人巨大在线播放| 9色porny在线观看| 蜜桃在线观看..| 精品欧美一区二区三区在线| 亚洲欧美色中文字幕在线| 国产精品秋霞免费鲁丝片| 亚洲精品国产区一区二| 高清毛片免费观看视频网站 | 精品久久蜜臀av无| 欧美精品高潮呻吟av久久| 午夜老司机福利片| 丝瓜视频免费看黄片| 一区在线观看完整版| 国产日韩欧美视频二区| 十八禁人妻一区二区| 色婷婷久久久亚洲欧美| 久久精品国产亚洲av高清一级| 精品国产国语对白av| 国产午夜精品久久久久久| 搡老乐熟女国产| 精品一区二区三卡| 国产亚洲欧美精品永久| 日韩精品免费视频一区二区三区| 亚洲精品国产一区二区精华液| 免费黄频网站在线观看国产| 国产精品亚洲av一区麻豆| 99久久人妻综合| 操出白浆在线播放| 亚洲全国av大片| 亚洲天堂av无毛| 女同久久另类99精品国产91| 日韩大码丰满熟妇| 天堂俺去俺来也www色官网| 午夜福利在线观看吧| 老司机亚洲免费影院| 天天添夜夜摸| 午夜福利影视在线免费观看| 99国产综合亚洲精品| 黄色毛片三级朝国网站| 成人免费观看视频高清| 午夜日韩欧美国产| 人人妻人人澡人人爽人人夜夜| 精品国产乱子伦一区二区三区| 久久久欧美国产精品| 免费av中文字幕在线| 久久精品熟女亚洲av麻豆精品| 成人国产av品久久久| 国产一区二区在线观看av| 999久久久精品免费观看国产| 欧美亚洲 丝袜 人妻 在线| 下体分泌物呈黄色| 午夜福利在线观看吧| 欧美日韩精品网址| 免费看十八禁软件| 日本黄色视频三级网站网址 | 黄片播放在线免费| 高清视频免费观看一区二区| 人人澡人人妻人| a级片在线免费高清观看视频| av一本久久久久| 久久久精品免费免费高清| 成年女人毛片免费观看观看9 | 他把我摸到了高潮在线观看 | 国产亚洲欧美在线一区二区| 黄色毛片三级朝国网站| a级片在线免费高清观看视频| √禁漫天堂资源中文www| 日本欧美视频一区| 91麻豆av在线| 国产真人三级小视频在线观看| 日韩大片免费观看网站| 丁香六月天网| 亚洲综合色网址| 日韩成人在线观看一区二区三区| 少妇的丰满在线观看| 国产在线观看jvid| 777米奇影视久久| 一个人免费看片子| 国产亚洲欧美精品永久| 欧美av亚洲av综合av国产av| 日韩一卡2卡3卡4卡2021年| 99国产极品粉嫩在线观看| 一进一出抽搐动态| 极品少妇高潮喷水抽搐| 亚洲第一欧美日韩一区二区三区 | 国产成人一区二区三区免费视频网站| 精品人妻1区二区| 亚洲成a人片在线一区二区| 精品亚洲成a人片在线观看| 菩萨蛮人人尽说江南好唐韦庄| www.自偷自拍.com| 成人国语在线视频| 成年动漫av网址| 欧美中文综合在线视频| a级片在线免费高清观看视频| 丰满迷人的少妇在线观看| 深夜精品福利| 国产1区2区3区精品| 色94色欧美一区二区| 一区二区三区国产精品乱码| 中文字幕人妻丝袜一区二区| 成人三级做爰电影| 飞空精品影院首页| 欧美在线一区亚洲| 美国免费a级毛片| 日韩欧美三级三区| 精品国内亚洲2022精品成人 | 久久中文字幕一级| 美女国产高潮福利片在线看| a级毛片在线看网站| 亚洲欧美日韩另类电影网站| tocl精华| 91字幕亚洲| 日本撒尿小便嘘嘘汇集6| 亚洲精品在线观看二区| 777米奇影视久久| 国产精品香港三级国产av潘金莲| 成人18禁在线播放| 99re6热这里在线精品视频| 精品国产乱子伦一区二区三区| 黄色视频不卡| av电影中文网址| 国产免费福利视频在线观看| 在线观看66精品国产| 久久久久国内视频| 少妇粗大呻吟视频| 亚洲国产av新网站| 巨乳人妻的诱惑在线观看| 成人三级做爰电影| 1024香蕉在线观看| 啦啦啦视频在线资源免费观看| 亚洲中文av在线| 精品少妇内射三级| 久久精品国产亚洲av高清一级| 欧美日韩av久久| 国产成人影院久久av| 国产精品久久久久久精品古装| 一区在线观看完整版| 欧美+亚洲+日韩+国产| 国产97色在线日韩免费| 午夜福利乱码中文字幕| 最新美女视频免费是黄的| 国产伦理片在线播放av一区| 亚洲色图综合在线观看| 欧美日韩亚洲综合一区二区三区_| 欧美国产精品va在线观看不卡| 69av精品久久久久久 | 热99re8久久精品国产| 亚洲美女黄片视频| 成人18禁在线播放| 久久久水蜜桃国产精品网| 色综合婷婷激情| 久久人妻熟女aⅴ| 欧美精品亚洲一区二区| 国产精品成人在线| 老熟女久久久| 欧美另类亚洲清纯唯美| 18禁国产床啪视频网站| 少妇猛男粗大的猛烈进出视频| 老司机影院毛片| 欧美人与性动交α欧美软件| 一边摸一边抽搐一进一小说 | 99国产精品免费福利视频| 99国产精品99久久久久| 中文字幕人妻丝袜制服| 精品少妇一区二区三区视频日本电影| 制服诱惑二区| 精品亚洲成a人片在线观看| 99国产精品一区二区三区| 电影成人av| 建设人人有责人人尽责人人享有的| 美女主播在线视频| 一个人免费看片子| 国产高清国产精品国产三级| 国产一区二区三区综合在线观看| 黑人操中国人逼视频| 99国产精品一区二区蜜桃av | 国产91精品成人一区二区三区 | 少妇裸体淫交视频免费看高清 | 丰满人妻熟妇乱又伦精品不卡| 日本精品一区二区三区蜜桃| 丝袜美足系列| av网站免费在线观看视频| 热re99久久精品国产66热6| 精品国产一区二区久久| 巨乳人妻的诱惑在线观看| 欧美大码av| 在线天堂中文资源库| 午夜福利影视在线免费观看| 亚洲av国产av综合av卡| 亚洲 国产 在线| 777米奇影视久久| 欧美激情 高清一区二区三区| 少妇猛男粗大的猛烈进出视频| 18禁美女被吸乳视频| 一边摸一边做爽爽视频免费| 精品国产一区二区久久| 亚洲精品久久成人aⅴ小说| 色视频在线一区二区三区| 99精国产麻豆久久婷婷| 精品福利观看| 操出白浆在线播放| 亚洲国产看品久久| 青青草视频在线视频观看| 国产精品自产拍在线观看55亚洲 | 日韩欧美一区二区三区在线观看 | 久久精品国产综合久久久| 亚洲黑人精品在线| 国产高清视频在线播放一区| 成人国语在线视频| 亚洲,欧美精品.| 亚洲午夜理论影院| 国产伦人伦偷精品视频| 热99国产精品久久久久久7| 男女之事视频高清在线观看| a级毛片黄视频| 亚洲精华国产精华精| 热re99久久精品国产66热6| 国产精品美女特级片免费视频播放器 | 亚洲av成人一区二区三| 女性生殖器流出的白浆| 亚洲专区字幕在线| cao死你这个sao货| 变态另类成人亚洲欧美熟女 | 99香蕉大伊视频| 动漫黄色视频在线观看| 后天国语完整版免费观看| 在线观看www视频免费| 欧美成狂野欧美在线观看| 成人av一区二区三区在线看| 高清毛片免费观看视频网站 | 精品第一国产精品| 男女床上黄色一级片免费看| 高潮久久久久久久久久久不卡| 啦啦啦 在线观看视频| 窝窝影院91人妻| 女性生殖器流出的白浆| 热re99久久国产66热| 国内毛片毛片毛片毛片毛片| 老司机深夜福利视频在线观看| 国产成人精品久久二区二区91| 久久精品国产综合久久久| 啦啦啦 在线观看视频| 久久国产亚洲av麻豆专区| 久久精品人人爽人人爽视色| 一区二区三区国产精品乱码| 91成年电影在线观看| 欧美亚洲日本最大视频资源| 18在线观看网站| 不卡一级毛片| 免费观看a级毛片全部| 婷婷成人精品国产| 成人手机av| 国产精品久久久久久精品古装| 桃红色精品国产亚洲av| 飞空精品影院首页| 国产极品粉嫩免费观看在线| 亚洲熟女毛片儿| 亚洲天堂av无毛| 自线自在国产av| 一本—道久久a久久精品蜜桃钙片| 精品一区二区三卡| 免费不卡黄色视频| 波多野结衣av一区二区av| 婷婷丁香在线五月| 日韩欧美一区二区三区在线观看 | 久久精品国产99精品国产亚洲性色 | 一级片免费观看大全| 久久久久久人人人人人| 免费久久久久久久精品成人欧美视频| 欧美激情极品国产一区二区三区| 老熟妇乱子伦视频在线观看| 99精品久久久久人妻精品| 成年动漫av网址| 午夜福利视频精品| 久久精品成人免费网站| 女人高潮潮喷娇喘18禁视频| 男女高潮啪啪啪动态图| 亚洲av日韩精品久久久久久密| 亚洲午夜精品一区,二区,三区| 后天国语完整版免费观看| 欧美av亚洲av综合av国产av| 高潮久久久久久久久久久不卡| 亚洲国产精品一区二区三区在线| 成人av一区二区三区在线看| 亚洲精品国产一区二区精华液| 在线天堂中文资源库| 成年版毛片免费区| 国产麻豆69| 国产精品免费一区二区三区在线 | 在线观看免费视频日本深夜| 日本wwww免费看| 蜜桃在线观看..| 国产亚洲精品一区二区www | 91老司机精品| 精品免费久久久久久久清纯 | 五月开心婷婷网| 一二三四在线观看免费中文在| 性少妇av在线| 欧美黑人欧美精品刺激| 在线 av 中文字幕| 人人妻人人澡人人看| 久久中文看片网| 国产精品久久久久成人av| 国产av又大| 人妻 亚洲 视频| 成人永久免费在线观看视频 | 精品亚洲成国产av| 国产成人av激情在线播放| 国产精品久久久久久人妻精品电影 | 欧美亚洲 丝袜 人妻 在线| av天堂在线播放| 日本av手机在线免费观看| 日韩人妻精品一区2区三区| 丁香六月天网| 窝窝影院91人妻| 又大又爽又粗| 欧美在线黄色| 亚洲国产欧美网| 亚洲熟女精品中文字幕| 两个人免费观看高清视频| 欧美日韩国产mv在线观看视频| 波多野结衣一区麻豆| 国产精品国产av在线观看|