去城市化作用前后黑龍江省積溫對比分析*

李亞濱，陳 莉**，魏 磊，王一達，劉 棟，孫 爽，王麗巖，黃英偉

去城市化作用前后黑龍江省積溫對比分析*

李亞濱1，陳 莉1**，魏 磊1，王一達2，劉 棟1，孫 爽1，王麗巖3，黃英偉4

（1.黑龍江省氣候中心，哈爾濱 150030；2.哈爾濱市氣象臺,哈爾濱 150030；3.中國氣象局華風集團，北京 100081；4.黑龍江省公共氣象服務中心，哈爾濱 150030）

肩負糧食安全重要職責的黑龍江省，地處中高緯度，溫度資源是限制該省農業生產的重要因素。≥10℃積溫及積溫帶的變化，對農業品種布局、農業種植結構調整、產量形成都起著決定性作用。本研究采用1961?2020年黑龍江省31個國家氣象站實測日平均氣溫、月平均氣溫和“基于城市化偏差訂正的中國大陸月平均地面氣溫網格數據集”資料，進行了去城市化作用前后黑龍江省≥10℃積溫對比分析。研究表明，與去城市化作用前相比，去城市化作用后近60a黑龍江省積溫線性變化趨勢和變化傾向率有所放緩；1991?2020年近30a積溫平均值普遍降低，城市化作用導致的積溫增量孫吳站超過200℃·d，呼瑪、黑河、尚志、北安站超過100℃·d，其他大部地區在20～100℃·d；黑龍江省6條積溫帶中有5條因去城市化作用而調整，第一積溫帶范圍大幅縮小，進而導致第二積溫帶范圍顯著增大，第三、四、五積溫帶有所南移。

黑龍江省；積溫；積溫帶；城市化作用；對比分析

IPCC《氣候變化與土地特別報告》[1]指出，在全球變暖背景下，氣候變化和極端氣候事件加重了土地壓力，并嚴重影響全球糧食安全。氣候變化對低緯地區糧食產量產生負面影響，而對較高緯地區一些作物（如玉米、小麥等）產量起到了積極作用[2]。在全球變暖大背景下，地處中國最高緯度的黑龍江省變暖顯著，1961?2018年黑龍江省年平均氣溫呈顯著上升趨勢，上升速率為0.34℃·10a?1，高于全國平均水平[3]。由于其所處緯度較高，溫度資源一直是制約該糧食大省農業生產的重要因素，氣候變暖造成的溫度資源增加，會對農業結構、種植制度和農作物產量產生重要影響。

熱量是農業生物生存、生長發育必需的外界環境因子和能量來源[4]。溫度資源在傳統上常用積溫、無霜期等指標來衡量。積溫，又稱度·日，為某一時期內大于或小于某一界限溫度的日平均溫度的總和[5]。它在農作物育種、選種、生產和生長發育評估模型中被廣泛使用。前人從積溫的概念及物理意義[6?9]、變化特征[10?16]、模擬[17?19]、空間格點化方法[20]等方面進行了大量研究。尤其對地處中高緯度、農業生產受溫度資源制約的黑龍江省積溫變化，開展了不少研究。曹萌萌等[21]研究表明，1981?2012年黑龍江省≥10℃積溫年代際變化呈顯著增加趨勢，各積溫帶較20世紀90年代劃分的積溫帶大致向北移、東擴了一個積溫帶。季生太等[22?23]指出，1961?2005年黑龍江省積溫呈增加趨勢，90年代增加顯著，1976?2005年積溫帶分布較1961?1990年北移東擴明顯，尤其是第一積溫帶北移平均達0.5個緯距，第三積溫帶東擴達 2個經距。袁湘玲等[11]研究指出，1961?2017年黑龍江省穩定通過10℃積溫以27～101℃·d·10a?1的趨勢顯著增加。朱紅蕊等[15]研究得出，1961?2014年黑龍江省≥10℃積溫呈顯著上升趨勢，傾向率達67.7℃·d·10a?1，1981?2014年黑龍江省6個積溫帶區域均有變化，第一積溫帶所占全省面積由1961?1990年的1%增加到8%，第六積溫帶范圍縮小8個百分點。這些研究揭示了黑龍江省積溫變化特征，但基于城市化偏差的積溫變化研究還較為少見。同時，東北平原廣袤的農田大多距離城市較遠，近幾十年受城市化影響程度遠小于氣象站所受到的影響。更為關鍵的是，研究表明城市化偏差是中國地面氣溫觀測記錄中最大的系統性偏差[24?28]，由此必然帶來因城市化因素引起的氣溫資料的不確定性。本研究旨在分離出溫度資源增加中的城市化作用，對黑龍江省溫度資源增加量進行科學評估，以便為農業結構、種植制度和農作物產量決策提供更科學依據。研究采用黑龍江省31個國家氣象站1961?2020年實測逐日平均氣溫和“基于城市化偏差訂正的中國大陸月平均地面氣溫網格數據集”中的月平均氣溫資料，進行去城市化作用前后積溫對比分析，對于厘清溫度資源的真實變化，指導當地農業生產有明顯的科學意義，同時對生長發育進程測算模型中積溫的訂正使用提供參考。

1 資料與方法

1.1 資料及其來源

黑龍江省31個國家基本/基準站逐日平均氣溫觀測資料和月平均氣溫資料時間范圍為1961?2020年，來自黑龍江省氣象數據中心。

經城市化偏差訂正后各站月平均氣溫資料取自“基于城市化偏差訂正的中國大陸月平均地面氣溫網格數據集”[26]，時間范圍為1961?2015年。該數據集綜合考慮了觀測連續性、遷站次數、觀測點附近聚落人口數量、測站周圍12km2內人工建筑比率、歷史資料的起止年份、與附近城鎮中心的距離等因素，確定了143個中國大陸地面氣溫參考站[25?27]，而后采取自東往西迭代訂正的方法，以各參考站與其目標站年均氣溫相關系數的平方為權重，得到目標站年和月平均氣溫參考序列，利用目標站及其參考氣溫序列趨勢之差作為該目標站城市化偏差訂正值。該數據集以2015年為基準年，向前逐漸遞加平均每年的城市化影響，訂正過的序列把站點當前及未來幾年的氣溫作為固定值，使新的氣溫資料在未來數年內具有更好的可延展性。因此，把2016?2020年實測月平均氣溫數據續接到該序列中，形成黑龍省31個氣象站基于城市化偏差訂正的1961?2020年月平均氣溫序列。

地理信息數據采用1:25萬的DEM數據，包括省、市、縣各級行政邊界和高程數據，資料來自中國氣象局。

1.2 基于日平均氣溫數據計算≥10℃活動積溫的方法

日平均氣溫穩定通過10℃為黑龍江省農作物正常生長的適宜溫度，通常用穩定通過10℃期間的積溫作為衡量作物生長的熱量指標[29]。首先采用5日滑動平均法確定≥10℃的初/終日，初日為滑動平均溫度中穩定通過10℃的第一個日平均氣溫≥10℃的日期，終日為滑動平均溫度中第一個小于10℃的前一個5d中最后一個日平均氣溫大于或等于10℃的日期，持續天數為初日至終日之間的天數[29]。活動積溫即≥10℃初日和終日之間的日平均氣溫之和。

1.3 基于月平均氣溫數據計算≥10℃活動積溫的方法

1.3.1 ≥10℃初/終日計算方法

由表1可知，與各站基于日平均氣溫得到的1961?2020年≥10℃初日和終日平均值相比較，由實測月平均氣溫計算得到的各站初日普遍偏早、終日普遍偏晚。

表1 黑龍江省各站基于實測月平均氣溫得到的≥10℃初日、終日誤差天數(d)

1.3.2 ≥10℃積溫計算方法

初、終日所在月≥10℃積溫值R為

式中，T為≥10℃初日所在月的最后一天日平均氣溫或≥10℃終日所在月的第一天日平均氣溫（℃）；h為初日或終日所在月≥10℃的天數（d）。

T由式（6）求得，即

由初日、終日所在月≥10℃積溫之和，再加上初日、終日其間各完整月的積溫值，則可求得全年≥10℃積溫值。

與基于各站日平均氣溫得到的各站活動積溫相比較，由實測月平均氣溫計算得到的活動積溫普遍偏小，偏小數值在35～96℃·d，占全年積溫值的2.1%～3.5%，二者相關系數為0.9996，為高度相關。

1.4 基于城市化偏差訂正的近30a積溫平均值（1991?2020年）計算流程

根據以上方法，可以得到基于城市化偏差訂正的近30a積溫平均值（1991?2020年）。其計算流程為從實測日、月平均氣溫出發，采用五日滑動平均方法分別計算≥10℃初、終日，隨后完成對康德拉公式的修正；然后把基于城市化偏差訂正的月平均氣溫帶入修正的康德拉公式中，分別計算相應的≥10℃初、終日和積溫值；最后結合實測日平均氣溫得到的前30a（1961?1990年）積溫平均值，得到基于城市化偏差訂正的近30a積溫平均值（1991?2020年），具體流程如圖1所示。

1.5 80%保證率積溫計算方法

在實際指導農業生產過程中，排除生物學零度以下的無效積溫，使用80%保證率數據來表征作物發育所需要的熱量更為貼切和符合實際[29]。在積溫帶劃分中采用經驗頻率法計算各站80%保證率下穩定通過10℃活動積溫數據[30]，保證率（P）計算公式為

圖1 基于城市化偏差訂正的近30a（1991?2020年）積溫均值計算流程

式中，n為樣本序列數（年），m為任一序號，m=1，2，…，n。

1.6 積溫帶劃分方法和積溫均值計算方法

農業、氣象部門通常以30a內80%保證率下≥10℃積溫進行積溫帶劃分。20世紀90年代中期，黑龍江省劃分了6條積溫帶[31]。本研究基于1961?1990年實測逐日平均氣溫數據，確定了第一至第六積溫帶的空間位置。31個氣象站所屬積溫帶情況如表2所示，各積溫帶積溫均值，即各積溫帶所包含的氣象站積溫平均值。

1.7 數據處理

利用趨勢分析方法[32?33]分析時間序列的趨勢和線性變化傾向率。數據分析和處理應用Fortran、SPSS、EXCEL軟件，繪圖采用了GIS軟件。

2結果與分析

2.1 去城市化作用前后黑龍江省≥10℃積溫近60a變化趨勢對比

1961?2020年黑龍江省≥10℃積溫序列用31個氣象站≥10℃積溫逐年平均值計算得到。基于實測日平均氣溫計算得到的黑龍江省≥10℃積溫近60a呈顯著增加趨勢，趨勢系數為0.6654，通過了0.01水平的顯著性檢驗，線性變化傾向率為66.0℃·d·10a?1。基于城市化偏差訂正的月平均氣溫計算得到的積溫近60a來趨勢系數為0.6236，線性變化傾向率為48.5℃·d·10a?1（圖2）。可見，去城市化作用后，黑龍江省積溫線性增加趨勢和增加速率均有所放緩。

表2 黑龍江≥10℃積溫帶中氣象站分布情況（1961?1990年）

圖2 1961?2020年去城市化作用前后黑龍江省≥10℃積溫變化趨勢

2.2 去城市化作用前后黑龍江省近30a≥10℃活動積溫空間分布對比

圖3a為去城市化作用前后近30a（1991?2020年）黑龍江省≥10℃積溫均值空間分布圖。由圖可見，去城市化作用后，黑龍江省各積溫等值線普遍南移，其中大興安嶺北部2100℃·d以下積溫值范圍變化不大，但2100～2400℃·d和2400～2700℃·d積溫區域范圍明顯增大，黑河部分市縣的≥10℃積溫均值由原來的2400～2700℃·d減至2100～2400℃·d；三江平原東部一帶在去城市化作用后2700℃·d以上的積溫范圍明顯收縮，由原來的大于2700℃·d調整為2400～2700℃·d，積溫差在100～200℃·d。由此可見，近30a黑河部分市縣和三江平原東部區域城市化作用影響較大。從去城市化作用前后的積溫均值變化量來看（圖3b），城市化作用導致的積溫增加量大部地區在20～100℃·d，其中除漠河以外的大興安嶺地區、黑河、齊齊哈爾西部、哈爾濱中西部及三江平原大部市縣增加量在50℃·d以上，尤其是呼瑪、黑河、尚志、北安的增量超過100℃·d，孫吳超過200℃·d。

2.3 去城市化作用前后黑龍江省積溫帶分布對比

積溫帶是農業部門進行農作物種植結構調整和品種選擇的重要依據。以1961?1990年積溫帶空間分布為標準（表2），分析近30a（1991?2020年）去城市化作用前后黑龍江省六條積溫帶積溫均值情況。結果見表3。由表中可見，各積溫帶積溫均值因去城市化作用有所減小，第一積溫帶均值減少23.8℃·d，第六積溫帶均值減少72.3℃·d。

對比去城市化作用前后的近30a積溫帶分布圖（圖4a和圖4b）可見，第一積溫帶范圍大幅縮小，原處于第一積溫帶的三江平原中部、牡丹江西部及哈爾濱東南部市縣，在去城市化作用后均處于第二積溫帶，由此導致第二積溫帶范圍顯著增大；第三、四、五積溫帶有所南移，尤其是黑河北部在去城市化之前為第四積溫帶，去城市化之后調整為第五積溫帶。盡管第六積溫帶積溫均值在去城市化作用后減少較明顯，但去城市化前后積溫均值都在1900℃·d以下，所以該積溫帶空間位置基本無變化。

圖3 去城市化作用前后黑龍江省近30a積溫均值（a）和城市化作用引起的積溫增量（b）的空間分布

注：圖3a 中，不同灰度的分界線為去城市化作用之前積溫均值，黑色實線為去城市化作用之后。

Note: The contour map filled with different grayscale is before de-urbanization and the solid line is after de-urbanization in Fig.2a.

表3 去城市化作用前后黑龍江≥10℃積溫帶近30a平均值(℃·d)

圖4 去城市化作用之前（a）和之后（b）黑龍江省近30a積溫帶分布圖

3 結論與討論

3.1 討論

氣候變暖給處于中高緯度、農業生產受溫度資源制約的地區帶來了難得的溫度資源，科學評估增加的溫度資源是正確決策的依據。最近幾十年的變暖，使黑龍江省積溫帶北移東擴[15,21?23]，本研究基于城市化偏差訂正月平均氣溫數據集，揭示了城市化作用在積溫帶北移東擴中有不可忽視的貢獻，尤其像三江平原中部、牡丹江西部、北部的黑河地區，其貢獻更為明顯。如果忽視了城市化偏差造成農作物品種越區種植，會造成低溫冷害風險加大，可能會給當地農業穩產、高產和農民增收帶來巨大挑戰。

目前城市化偏差訂正的月平均氣溫數據集資料年份為2015年，采用的是向前訂正方法，因此，本研究把2016?2020年數據直接續接到該數據集中是科學而合理的。但隨著時間的推移，距離2015年越來越遠，未來幾年的數據直接續接會產生較大的誤差，因此，有必要建立新的城市化偏差訂正月平均氣溫數據集，才能保證在未來幾年或幾十年相關的城市化作用研究得以開展。同時該數據集時間分辨率為月，如果在當地布設自動氣象站開展城市化作用監測，把城市化偏差數據訂正到日尺度，由此開展積溫研究會得到更精確的結果。因此，開展氣候變化背景下黑龍江省乃至東北地區城市化偏差精細化訂正和積溫帶變化研究，對于調整農業品種布局、種植結構和產量形成都有重要意義。

3.2 結論

（1）1961?2020年黑龍江省≥10℃積溫線性增加趨勢和增加速率有所放緩。在去城市化作用前，趨勢系數為0.6654，通過了0.01水平的顯著性檢驗，線性變化傾向率為66.0℃·d·10a?1，去城市化作用后，趨勢系數為0.6236，線性變化傾向率為48.5℃·d·10a?1。

（2）與去城市化作用前相比，黑龍江省各站近30a（1991?2020年）積溫均值在去城市化作用后普遍降低。城市化作用導致的積溫增加量孫吳站超過200℃·d，呼瑪、黑河、尚志、北安站增量超過100℃·d，其他大部地區在100℃·d以下。

（3）與去城市化作用前相比，黑龍江省6條積溫帶中有5條因去城市化作用而調整。各條積溫帶積溫均值因去城市化作用而有所減小。第一積溫帶范圍大幅縮小，原處于第一積溫帶的三江平原中部、牡丹江西部及哈爾濱東南部市縣，在去城市化作用后均處于第二積溫帶，由此也導致第二積溫帶范圍顯著增大；第三、四、五積溫帶有所南移。

[1] Valérie M D, Zhai P M, P?rtner H O, et al. Summary for policymakers[M/OL]//IPCC. Climate change and land: an IPCC special report on climate change, desertification, land degradation, sustainable land management, food security, and greenhouse gas fluxes in terrestrial ecosystems. https://www.ipcc.ch/srccl/chapter/summary-for-policymakers/, 2021-08-25.

[2] Cheikh M, Cynthia R. Food security [M/OL]://IPCC. Climate change and land: an IPCC special report on climate change, desertification, land degradation, sustainable land management, food security, and greenhouse gas fluxes in terrestrial ecosystems. https://www.ipcc.ch/site/assets/uploads/ sites/4/2021/02/08_Chapter-5_3.pdf,2021-08-25.

[3] 陳莉,王波,劉玉蓮,等.黑龍江省氣候變化監測公報(2019)[M].北京:氣象出版社,2020:6.

Chen L, Wang B, Liu Y L, et al. Heilongjiang climate change monitoring bulletin (2019) [M]. Beijing: China Meteorological Press,2020:6.(in Chinese)

[4] 《中國農業百科全書》委員會.中國農業百科全書·農業氣象卷[M].北京:農業出版社,1986:237.

Committee of Agricultural Encyclopedia of China. Agriculture encyclopedia of China. Volume of agricultural meteorology[M]. Beijing: China Agricultural Press,1986: 237. (in Chinese)

[5] 顧鈞禧.大氣科學辭典[M].北京:氣象出版社,1994.

Gu J X. Dictionary of atmospheric sciences[M].Beijing: China Meteorological Press,1994.(in Chinese)

[6] 張子源,鄭大瑋,潘宇鷹,等.積溫及熱量資源概念的科學性問題與改進[J].中國農業氣象,2021,42(8):686-692.

Zhang Z Y, Zheng D W, Pan Y Y, et al. The scientific problem and improvement of the concepts of accumulated temperature and heat resource[J].Chinese Journal of Agrometeorology,2010,42(8):686-692.(in Chinese)

[7] 高洪蛟,李奇.淺談積溫[J].黑龍江氣象,2017,34(2):15.

Gao H J, Li Q. Discussion on thermal time[J]. Heilongjiang Meteorology,2017,34(2):15.(in Chinese)

[8] 馮樹常,汪擴軍.論濕積溫[J].氣象科學,1992,12(1):58-63.

Feng S C, Wang K J. Discussion on wet thermal time[J]. Journal of the Meteorological Sciences,1992,12(1):58-63. (in Chinese)

[9] 鄭大瑋,孫忠富.關于積溫一詞及其度量單位科學性問題的討論[J].中國農業氣象,2010,31(2):165-169.

Zheng D W, Sun Z F. Discussion on scientificalness problem of accumulated temperature and its unit[J]. Chinese Journal of Agrometeorology,2010,31(2):165-169. (in Chinese)

[10] 邱新法,王喆,曾燕,等.1960-2013年中國≥10℃積溫時空變化特征及其主導因素分析[J].江蘇農業科學,2017,45 (2):220-225.

Qiu X F, Wang Z, Zeng Y, et al. Change of ≥10℃ accumulated temperature's spatiotemporal characteristics and its changing factors during 1960-2013[J]. Jiangsu Agricultural Sciences,2017,45(2):220-225.(in Chinese)

[11] 袁湘玲,劉春生,魏磊.黑龍江熱量資源時空變化特征[J].中國農學通報,2019,35(10):121-127.

Yuan X L, Liu C S, Wei L. Temporal and spatial variation of heat resources in Heilongjiang[J]. Chinese Agricultural Science Bulletin,2019,35(10):121-127.(in Chinese)

[12] 李帥,張勃,馬彬,等.基于格點數據的中國1961?2016年≥5℃、≥10℃有效積溫時空演變[J].自然資源學報,2020, 35(5):1216-1227.

Li S, Zhang B, Ma B, et al. Spatiotemporal evolution of effective accumulated temperatures of ≥5℃ and ≥10℃ based on grid data in China from 1961 to 2016[J]. Journal of Natural Resources,2020,35(5):1216-1227.(in Chinese)

[13] 白磊,張帆,尚明,等.基于格點數據的1961-2018年中國多種積溫時空變化研究[J].地球信息科學學報,2021,23(8): 1446-1460.

Bai L, Zhang F, Shang M, et al. Evolution of the multiple accumulated temperature across mainland China in 1961-2018 with the gridded meteorological dataset[J]. Journal of Geo-Information Science,2021,23(8):1446-1460. (in Chinese)

[14] 蘆艷珍,韓小英,張蕾,等.基于GIS的山西有效積溫時空變化分析[J].山西農業科學,2020,48(11):1829-1834.

Lu Y Z, Han X Y, Zhang L, et al. Analysis of temporal and spatial variation of effective accumulated temperature in Shanxi province based on GIS[J]. Journal of Shanxi Agricultural Sciences,2020,48(11):1829-1834.(in Chinese)

[15] 朱紅蕊,尹嫦姣,張洪玲,等.黑龍江省≥10℃積溫時空變化特征[J].冰川凍土,2015,37(6):1473-1479.

Zhu H R, Yin C J, Zhang H L, et al. The spatial-temporal change of active cumulated temperature of ≥10℃ in Heilongjiang Province[J].Journal of Glaciology and Geocryology,2015,37(6):1473-1479.(in Chinese)

[16] 郝成元,李明秋.1960-2009年間河南省≥0℃和≥10℃活動積溫時空變化[J].中國農業資源與區劃,2019,40(6): 89-95.

Hao C Y, Li M Q. Spatio-temporal characters of ≥0℃ and ≥10℃ accumulated temperatures From 1960 to 2009 in Henan province[J].Chinese Journal of Agricultural and Regional Planning,2019,40(6):89-95.(in Chinese)

[17] 邱新法,卞宗雅,曾燕,等.重慶山地界限溫度起止日期和持續日數的分布式模擬[J].自然科學通報,2009,19(7): 746-753.

Qiu X F, Bian Z Y, Zeng Y, et al. Distributed simulation of starting and ending date and duration of Chongqing mountain boundary temperature[J]. Prog Nat Sci,2009, 19(7): 746-753.(in Chinese)

[18] 姜會飛,郭勇,張玉瑩.不同下限基點溫度對積溫模型模擬效果的影響[J].中國農業大學學報,2018,23(5):131-141.

Jiang H F, Guo Y, Zhang Y Y, et al. Impact of base temperature on the growing degree-day and simulation effect of GDD model[J]. Journal of China Agricultural University,2018,23(5):131-141.(in Chinese)

[19] 朱海霞,李秀芬,王萍,等.黑龍江省水稻生長季積溫計算方法[J].應用氣象學報,2017,28(2):247-256.

Zhu H X, Li X F, Wang P, et al. Methods of accumulated temperature during rice growing stage in Heilongjiang province[J]. Quarterly Journal of Applied Meteorology, 2017,28(2):247-256.(in Chinese)

[20] 李飛,王春,趙軍,等.中國陸地多年平均積溫空間化研究[J].自然資源學報,2010,25(5):778-784.

Li F, Wang C, Zhao J, et al. The spatialization of multi-year average accumulated temperature in China[J]. Journal of Natural Resources,2010,25(5):778-784.(in Chinese)

[21] 曹萌萌,李俏,張立友,等.黑龍江省積溫時空變化及積溫帶的重新劃分[J].中國農業氣象,2014,35(5):492-496.

Cao M M, Li Q, Zhang L Y, et al. Accumulated temperature variation and accumulated temperature rezone in Heilongjiang Province[J]. Chinese Journal of Agrometeorology,2014,35(5):492-496.(in Chinese)

[22] 季生太,楊明,紀仰慧,等.黑龍江省近45年積溫變化及積溫帶的演變趨勢[J].中國農業氣象,2009,30(2):133-137.

Ji S T, Yang M, Ji Y H, et al. Change of accumulated temperature and evolution trends of accumulated temperature zone over last 45 years in Heilongjiang Province[J]. Chinese Journal of Agrometeorology,2009, 30(2):133-137.(in Chinese)

[23] 閆平,楊明,王萍,等.基于GIS的黑龍江省積溫帶精細化區劃[J].黑龍江氣象,2009,26(1):26-27.

Yan P, Yang M, Wang P, et al. Fine zoning of thermal time zone in Heilongjiang Province on GIS[J]. Heilongjiang Meteorology,2009,26(1):26-27.(in Chinese)

[24] 司鵬,李慶祥,李偉.城市化進程對中國東北部氣溫增暖的貢獻檢測[J].氣象,2010,36(2):13-21.

Si P, Li Q X, Li W. Impact of urbanization on surface air temperature warming in Northeast China[J]. Meteorological Monthly,2010,36(2):13-21.(in Chinese)

[25] 任玉玉,任國玉,張愛英,等.城市化對地面氣溫變化趨勢影響研究綜述[J].地理科學進展,2010,29(11):1301-1310.

Ren Y Y, Ren G Y, Zhang A Y. An overview of researches of urbanization effect on land surface air temperature trends[J].Progress in Geography,2010,29(11):1301-1310. (in Chinese)

[26] 溫康民,任國玉,李嬌,等.國家基本/基準站地面氣溫資料城市化偏差訂正[J].地理科學進展,2019,38(4):600-611.

Wen K M, Ren G Y, Li J, et al. Adjustment of urbanization bias in surface air temperature over the mainland of China[J].Prog Geogr,2019,38(4):600-611.(in Chinese)

[27] 張愛英,任國玉,周江興,等.中國地面氣溫變化趨勢中的城市化影響偏差[J].氣象學報,2010,68(6):957-966.

Zhang A Y, Ren G Y, Zhou J X, et al. On the urbanization effect on surface air temperature records over China[J]. Acta Meteorologica Sinica,2010,68(6):957-966.(in Chinese)

[28] 周雅清,任國玉.華北地區地表氣溫觀測中城鎮化影響的檢測和訂正[J].氣候與環境研究,2005,10(4):743-753.

Zhou Y Q, Ren G Y. Identifying and correcting urban bias for regional surface air temperature series of North China over time period of 1961-2000[J].Climatic and Environmental Research,2005,10(4):743-753.(in Chinese)

[29] 胡毅,李萍.應用氣象學[M].北京:氣象出版社,2007.

Hu Y, Li P. Applied meteorology[M]. Beijing: China Meteorological Press,2007.(in Chinese)

[30] 任國玉.美國中東部溫帶穩定≥10℃積溫[J].地理研究,1987(4):70-76.

Ren G Y. Thermal time above 10℃ in the central and eastern in the United States[J]. Geographical Research, 1987(4):70-76.(in Chinese)

[30] 北京農業大學農業氣象專業農業氣候教學組.農業氣候學[M].北京:農業出版社,1987:80.

Agroclimatology Teaching Group of Agricultural Meteorology of Beijing Agricultural University. Agroclimatology [M]. Beijing: Agriculture Press,1987:80.(in Chinese)

[31] 黑龍江省農業科學院.黑龍江省農作物品種積溫區劃[M/OL].http://www.hljagri.org.cn/nygk/nzwpzjwqh/200807/P020080702301995009250.jpg,2021-09-12.

Heilongjiang Academy of Agricultural Sciences. Regionalization of thermal time of crop varieties in Heilongjiang Province [M/OL]. http://www.hljagri.org.cn/nygk/nzwpzjwqh/200807/ P020080702301995009250.jpg,2021-09-12.(in Chinese)

[32] 馬開玉,丁裕國,屠其璞,等.氣候統計原理與方法[M].氣象出版社,1993.

Ma K Y, Ding Y G, Tu Q P, et al. Principles and methods of climate statistics[M].Beijing: China Meteorological Press,1993.(in Chinese)

[33] 魏鳳英.現代氣候統計診斷與預測技術.2版[M].北京:氣象出版社,2007.

Wei F Y. Modern climate statistics diagnosis and prediction technology.2nd ed[M].Beijing: China Meteorological Press,2007.(in Chinese)

Comparative Analysis of Thermal Time above 10℃ in Heilongjiang Province before and after De-urbanization

LI Ya-bin1, CHEN Li1, WEI Lei1, WANG Yi-da2, LIU Dong1, SUN Shuang1, WANG Li-yan3, HUANG Ying-wei4

(1. Heilongjiang Climate Center, Harbin 150030, China; 2. Harbin Meteorological Observatory, Harbin 150030; 3. Huafeng Group of China Meteorological Administration, Beijing 100081; 4. Public Meteorological Service Center of Heilongjiang Province, Harbin 150030)

Heilongjiang province, with an important responsibility for food security, is in the middle-high latitudes. Thermal resource is an important factor restricting agricultural production in the province. Meanwhile its climate warming has been significantly higher than the national average in the past 60 years. Changes in thermal time and thermal time zones play decisive roles in the distribution of agricultural varieties, the adjustment of agricultural planting structure, and the formation of output. Using daily and monthly average temperature data of 31 national meteorological stations in Heilongjiang province from 1961 to 2020 and “Grid dataset of monthly average surface temperature in mainland China revised based on the urbanization deviation”, changes of thermal time above 10℃ in Heilongjiang province before and after de-urbanization were analyzed. The results show that after removing the effects of urbanization, the linear change trend and rate of change tendency of thermal time from 1961 to 2020 have slowed down in Heilongjiang province. In the past 30 years during 1991 to 2020, the average thermal time generally decreased after de-urbanization. The increase thermal time caused by urbanization was below 100℃·d in most areas. The amount of Huma, Heihe, Shangzhi and Bei'an exceeds 100℃·d, and Sunwu exceeds 200℃·d. Five of six thermal time zones in Heilongjiang province were adjusted due to the removal of urbanization. The range of First Zone was greatly reduced, which led to a significant increase in Second Zone, and Third, Fourth, and Fifth Zones moved south.

Heilongjiang province; Thermal time above 10℃; Zone of thermal time; Urbanization; Comparative analysis

10.3969/j.issn.1000-6362.2022.06.002

李亞濱,陳莉,魏磊,等.去城市化作用前后黑龍江省積溫對比分析[J].中國農業氣象,2022,43(6):440-449

收稿日期：2021?09?28

黑龍江省自然科學基金聯合引導項目（LH2020D015）

通訊作者：陳莉，正高級工程師，博士，主要從事氣候與氣候變化研究，E-mail:climate_chenli@163.com

李亞濱，E-mail:113931077@qq.com