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

    京津冀城市群熱島定量評估

    2017-10-13 12:33:05劉勇洪房小怡欒慶祖權(quán)維俊
    生態(tài)學(xué)報(bào) 2017年17期

    劉勇洪,房小怡,張 碩,欒慶祖,權(quán)維俊

    1 北京市氣候中心,北京 100089 2 中國氣象局北京城市氣象研究所, 北京 100089 3 中國氣象局京津冀環(huán)境氣象預(yù)報(bào)預(yù)警中心, 北京 100089

    京津冀城市群熱島定量評估

    劉勇洪1,2,*,房小怡1,張 碩1,欒慶祖1,權(quán)維俊3

    1 北京市氣候中心,北京 100089 2 中國氣象局北京城市氣象研究所, 北京 100089 3 中國氣象局京津冀環(huán)境氣象預(yù)報(bào)預(yù)警中心, 北京 100089

    在采用城鄉(xiāng)二分法估算區(qū)域范圍內(nèi)多個城市地表熱島強(qiáng)度(Surface Urban Heat Island,SUHI)時(shí),如何選擇城鎮(zhèn)化影響最小的周邊鄉(xiāng)村背景是一個技術(shù)難點(diǎn),提出了一種基于地形、土地利用、植被覆蓋和城市夜間燈光指數(shù)來確定鄉(xiāng)村背景的SUHI估算方法,并建立了基于SUHI和熱島比例指數(shù)(Urban Heat Island Proportion Index,UHPI)的城市熱島強(qiáng)度定量評估方法。利用上述方法,基于長時(shí)間序列MODIS和NOAA衛(wèi)星資料,開展了京津冀城市群11個平原城市熱島時(shí)空變化分析與評估,并開展了社會經(jīng)濟(jì)驅(qū)動因子對城市熱島大小的影響評估研究,同時(shí)結(jié)合未來京津冀一體化發(fā)展提出相應(yīng)參考建議。研究結(jié)果表明:(1)建立的SUHI估算方法能有效監(jiān)測城市群熱島年/季和晝/夜變化,近5年(2010—2014)年均SUHI≥3℃的熱島總面積1926 km2,但在熱島最強(qiáng)的夏季白天可達(dá)7386 km2(占行政區(qū)域面積的5.8%);排名前四的分別是北京(2351 km2)、天津(1883 km2)、唐山(889 km2)和石家莊(611 km2),顯示出超大、特大城市及資源性城市貢獻(xiàn)了大部分城市群熱島面積;各中心城區(qū)平均SUHI和UHPI分別為3.0℃和0.61,熱島評估達(dá)到較嚴(yán)重等級以上的城市占到73%,表明當(dāng)前城市群整體熱島處于嚴(yán)峻現(xiàn)狀;(2)1994、2004年和2014年夏季白天城市群強(qiáng)熱島面積分別為190、1975、4539 km2,各中心城區(qū)平均SUHI分別為1.2、2.6、3.2℃,UHPI分別為0.29、0.58和0.69,熱島評估等級分別為“一般”、“較嚴(yán)重”和“嚴(yán)重”等級,反映了20年來京津冀城市群熱島迅速增強(qiáng)增大事實(shí);(3)各城市年均熱島面積增加2—86 km2/a,強(qiáng)熱島面積增加主要發(fā)生在超大城市,北京、天津強(qiáng)熱島區(qū)之間的最短空間距離從1994年的94 km逐步縮減到2014年的52 km,未來存在形成“京津區(qū)域熱島群”的可能,建議在京津之間建立“綠色生態(tài)屏障”來消除這種可能性;(4)城鎮(zhèn)人口數(shù)、國內(nèi)生產(chǎn)總值和用電量都極大地影響著京津冀城市熱島大小,擬合模型決定系數(shù)R2分別為0.9097、0.912和0.9661,意味著在未來京津冀一體化城市發(fā)展中可采取控制城市人口規(guī)模、減少能源消耗等措施減緩熱島效應(yīng)。

    鄉(xiāng)村背景;地表熱島強(qiáng)度;熱島比例指數(shù);定量評估;經(jīng)濟(jì)驅(qū)動因子;京津冀城市群

    Abstract: When estimating surface urban heat island intensity (SUHI) of different cities on a regional scale using the urban-rural difference method, it is technically difficult to select the rural districts or reference areas that exhibit minimal impact by urbanization. In this study, we propose a method of estimating SUHI based on land topography, land use, vegetation coverage, and urban nighttime light index to define the rural districts. Then, a quantitative evaluation method can be established to determine heat island intensity based on SUHI and the urban heat island proportion index (UHPI). Using this method, the temporal and spatial variation of urban heat island for the Beijing-Tianjin-Hebei urban agglomeration region was quantitatively evaluated based on long time series MODIS and NOAA satellite data. An impact assessment for socio-economic factors driving urban heat islands was also conducted, and some recommendations for Beijing-Tianjin-Hebei economic unity development are proposed. The results show that the SUHI estimation method can effectively reflect annual, seasonal, and diurnal variations of heat island in urban agglomerations. The annual mean heat island area with SUHI greater than or equal to 3.0℃ from 2010—2014 was 1,926 km2. It is highest on summer days at 7386 km2(5.8% of total urban group area). The four largest area cities are Beijing at 2351 km2, Tianjin at 1883 km2, Tangshan at 889 km2, and Shijiazhang at 611 km2, which represent megacities, extra-large cities, and resource-dependent cities contributing to most of the heat island areas. The mean SUHI and UHPI of all central urban areas are 3.0℃ and 0.61, respectively. Seventy-three percent of cities have heat island assessment levels that are more than relatively severe with SUHI greater than or equal to 3.0℃. This indicates that the phenomenon of heat islands in urban agglomeration is currently severe. During summer days in 1994, 2004, and 2014, the areas with SUHI greater than or equal to 5.0℃ were 190 km2, 1975 km2, and 4539 km2; the mean SUHI were 1.2℃, 2.6℃, and 3.2℃; and the UHPI were 0.29, 0.58, and 0.69. The heat island assessment levels were normal, relatively severe, and severe, respectively, which reflects the fact that the heat island effect for urban agglomerations has strengthened and expanded in the last 20 years. The annual growth rate of heat island areas is 2—86 km2/a, and the increasing areas of high-level heat islands in urban group regions occur mainly in megacity behemoths. The shortest distance between the Beijing heat island region and the Tianjin heat island region was reduced from 94 km in 1994 to 52 km in 2014. Thus, in future, the Beijing and Tianjin heat island regions will probably connect to form a new Beijing Tianjin heat island group, which may be an ecological disaster for urban agglomeration development. Therefore, an ecological green barrier between Beijing and Tianjin should be set up to halt the trend. Heat island area with SUHI greater than or equal to 3.0℃ for various cities correlate notably with the urban population, gross domestic product, and electricity consumption, and theR2are 0.9097, 0.912, and 0.9661, which indicates the important role of population control and energy consumption reduction in the development of the Beijing-Tianjin-Hebei economic unit.

    KeyWords: rural background;surface urban heat island intensity;heat island proportion index;quantitative evaluation;socio-economic driving factors;Beijing-Tianjin-Hebei urban agglomeration

    近年來,在全球增溫和高速城市化的背景下,城市氣候被認(rèn)為是主導(dǎo)城市生態(tài)環(huán)境的重要因素之一,城市氣候最明顯的特征就是城市熱島效應(yīng)(urban heat island effect,簡稱UHI)。UHI已成為21世紀(jì)面臨的重要生態(tài)環(huán)境問題[1- 2],UHI降低了人們生活的舒適度并加劇了大氣污染,嚴(yán)重影響了居民生活質(zhì)量,如何定量地監(jiān)測UHI的動態(tài)變化并開展UHI成因分析已成為當(dāng)前UHI研究的重要內(nèi)容。目前,常用的UHI研究手段有3種:氣象觀測、數(shù)值模擬和衛(wèi)星遙感。氣象觀測由于具有觀測資料時(shí)效長、定點(diǎn)、準(zhǔn)確、定量的優(yōu)勢,一直被作為基礎(chǔ)手段用于研究UHI的時(shí)空演變規(guī)律[3-7]。遙感則因具有時(shí)間同步性好、覆蓋范圍廣、空間結(jié)構(gòu)直觀定量等特點(diǎn),近20年已成為UHI研究的普遍手段[8-14]。此外,數(shù)值模擬也由于能深人研究UHI現(xiàn)象和成因之間的物理本質(zhì),已成為UHI機(jī)理形成、模擬及預(yù)測研究的一種重要手段[15-18]。

    UHI最初定義為城市中的氣象臺站觀測的氣溫明顯高于外圍郊區(qū)或鄉(xiāng)村的現(xiàn)象,近年來已擴(kuò)展至遙感觀測的城市陸表溫度高于郊區(qū)或鄉(xiāng)村的現(xiàn)象[19- 20],通常把這種主要通過衛(wèi)星觀測的城市熱島稱之為地表城市熱島(Surface Urban Heat Island,簡稱SUHI)[21-22]??紤]到UHI是一個城市高于鄉(xiāng)村的相對溫度距平值,由于涉及到城市區(qū)(或城市氣象代表站點(diǎn))與鄉(xiāng)村背景區(qū)(或鄉(xiāng)村氣象代表站點(diǎn))選擇的復(fù)雜性,不斷有新的UHI或SUHI計(jì)算方法被提出。如利用城市氣象站距離遠(yuǎn)近結(jié)合土地利用類型選擇選擇鄉(xiāng)村氣象站點(diǎn)來估算UHI[23],在城市內(nèi)部選擇大面積植被覆蓋的公園作為城市氣象代表站來反映城市化引起的熱島效應(yīng)[24],選擇城市中最低氣溫的站點(diǎn)作為背景參考站估算UHI等[25]?;谛l(wèi)星資料估算SUHI最常用的方法有城鄉(xiāng)二分法[26- 28],均值-標(biāo)準(zhǔn)差法[29- 31],正規(guī)化法或歸一化[32- 33]。在采用城鄉(xiāng)二分法計(jì)算SUHI時(shí),難點(diǎn)在于郊區(qū)或農(nóng)村等參考或背景地區(qū)的確定[21]。近年來,中國城鎮(zhèn)化迅速,分布于大城市周邊的郊區(qū)衛(wèi)星城鎮(zhèn)擴(kuò)展迅速,如何提取城鎮(zhèn)化影響最小的鄉(xiāng)村背景是一個復(fù)雜的技術(shù)難點(diǎn)。王建凱等[26]在利用MODIS影像研究北京城市熱島(冷島)時(shí),根據(jù)高分辨率影像結(jié)合土地覆蓋類型人為選定遠(yuǎn)離北京五環(huán)外的一塊平原農(nóng)田區(qū)域作為郊區(qū)背景,Imhoff等[34]在利用MODIS影像估算美國38個城市熱島強(qiáng)度時(shí),基于高分辨率美國國家土地覆蓋數(shù)據(jù)庫(NLCD)中的不透水蓋度數(shù)據(jù)集、中心城市距離和地形因子來確定不同城市的鄉(xiāng)村背景,Schwarz等[35]則采用了離城市周圍20 km區(qū)域作為鄉(xiāng)村背景,葉彩華等[36]人為選取了遠(yuǎn)離城區(qū)的兩塊5×5像元大小的長期種植結(jié)構(gòu)無變化的農(nóng)田作為鄉(xiāng)村背景,林奕桐等[28]則基于遙感反演的植被覆蓋度高值區(qū)來選定郊區(qū)農(nóng)田參照點(diǎn)來估算南寧城市熱島強(qiáng)度,劉勇洪等[25]在計(jì)算京津唐城市群熱島時(shí),主要根據(jù)植被覆蓋度和不透水蓋度結(jié)合土地利用類型來確定各城市的鄉(xiāng)村背景。上述方法在選擇鄉(xiāng)村背景時(shí),在選擇農(nóng)田的基礎(chǔ)上更多地考慮了植被覆蓋與不透水蓋度,由于這兩個指標(biāo)計(jì)算的不確定性以及國內(nèi)尚沒有建立全國性的標(biāo)準(zhǔn)植被覆蓋度和不透水蓋度數(shù)據(jù)集,使得國內(nèi)依據(jù)這些指標(biāo)選擇鄉(xiāng)村背景具有諸多不確定性,本文將探尋一種考慮簡便性與易獲取性來選擇鄉(xiāng)村背景的SUHI估算方法。

    SUHI具有明顯的時(shí)空變化特征[37-38],其大小和形態(tài)不但與城市核心區(qū)(或城市氣象代表站)與鄉(xiāng)村背景區(qū)(或郊區(qū)氣象代表站)選擇有關(guān),還與研究的時(shí)空尺度有關(guān)。在城市宏觀尺度上氣候類型[34]、天氣條件(包括溫度、云、風(fēng)、逆溫等)[39-41]、物候(晝/夜、季節(jié)、生長期)[31,39- 42]和土地利用類型[23,31]等因子對SUHI影響較明顯,在城市微觀尺度上地形條件[43-44]、植被豐度與綠地景觀結(jié)構(gòu)[16,45-46]、城市不透水蓋度[22,34]、城市建筑結(jié)構(gòu)[47-48]等因子對SUHI有較大影響;而人口、國內(nèi)生產(chǎn)總值(GDP)、能耗、人為熱等等作為社會經(jīng)濟(jì)驅(qū)動因子,對SUHI的形成和擴(kuò)大具有重要作用[22,42,48-49]。當(dāng)前研究大多關(guān)注單一城市或者是處于不同氣候帶的多個大城市之間局地SUHI之間的相互比較,而對處于同一氣候類型較大范圍城市連片發(fā)展的城市群SUHI效應(yīng)研究相對較少,由于城市群下墊面土地利用的變化對區(qū)域氣候的增溫效應(yīng)明顯[50],而且由于城市群下墊面的不均勻性,局地增溫效應(yīng)在城市群內(nèi)部以及城市群之間都存在明顯差異[51],因此有必要開展城市群內(nèi)各城市之間的SUHI效應(yīng)研究,為城市群的整體發(fā)展規(guī)劃提供參考。

    京津冀城市群作為中國未來要發(fā)展成一個世界級城市群的重要經(jīng)濟(jì)地區(qū),近30年城鎮(zhèn)化進(jìn)入了空前快速發(fā)展階段,城市規(guī)模不斷擴(kuò)大,生態(tài)環(huán)境惡化、能源消耗高以及SUHI帶來的不利影響和問題日益突顯。在京津冀一體化發(fā)展規(guī)劃和北京城市功能布局調(diào)整中,減緩、降低SUHI及其影響,成為改善整個京津冀地區(qū)城市生態(tài)環(huán)境亟待解決的問題,因此對京津冀城市群的SUHI監(jiān)測及評估將為改善城市生態(tài)環(huán)境提供基礎(chǔ)。本文將探討如何在大尺度區(qū)域范圍內(nèi)建立較為簡單、數(shù)據(jù)易獲取的SUHI估算和評估方法,并開展城市驅(qū)動經(jīng)濟(jì)因子包括經(jīng)濟(jì)規(guī)模(國內(nèi)生產(chǎn)總值,簡稱GDP)、人口規(guī)模、能源消耗(以電力為代表)、人為熱排放(以機(jī)動車保有量為代表)對SUHI大小的影響評估,為未來京津冀一體化發(fā)展及改善城市生態(tài)環(huán)境提供建議。

    1 研究區(qū)概況

    京津冀地區(qū)位于36.07—42.65°N,113.46—119.79°E(圖1),包括北京、天津直轄市與河北省。整個地區(qū)地勢西北高,東南低,燕山山脈和太行山脈自東向西又向南形成弧形山脈,山脈東部及南部為河北平原區(qū),境內(nèi)海拔相差懸殊,最高海拔為小五臺山(約2840 m),平原地區(qū)多數(shù)海拔為50 m以下。該地區(qū)屬溫帶季風(fēng)氣候區(qū),大部分地區(qū)四季分明,其特點(diǎn)是夏季炎熱多雨,春季干燥少雨多風(fēng)沙,秋季天氣晴爽,冬季寒冷少雪。年日照時(shí)數(shù)2400—3100 h,年無霜期120—200 d;年均降水量300—800 mm;年平均氣溫在4—13℃之間。本文主要研究京津冀城市群地級以上中等平原城市,包括北京、天津、石家莊、唐山、廊坊、保定、滄州、衡水、邢臺、邯鄲和承德等11個平原城市,秦皇島、張家口由于位于山區(qū)及山地,不做熱島研究,其中北京、天津和石家莊為省會城市,其余為地級城市;按照2013年12月中國城市規(guī)模新劃分標(biāo)準(zhǔn)(http://politics.people.com.cn/n/2014/1121/c70731- 26064695.html),結(jié)合2013年城區(qū)常住人口統(tǒng)計(jì)資料,北京、天津?yàn)槌蟪鞘?城區(qū)常住人口在1000萬以上),石家莊為特大城市(城區(qū)常住人口在500—1000萬之間),唐山、保定、滄州、形態(tài)、邯鄲為Ⅰ型大城市(城區(qū)常住人口在300—500萬之間),廊坊、衡水、秦皇島為Ⅱ型大城市(城區(qū)常住人口在100—300萬之間)。另外選擇各城市市轄區(qū)作為城市中心城區(qū)熱島狀況進(jìn)行評估。

    2 資料和方法

    2.1 資料

    (1)MODIS衛(wèi)星資料 MODIS衛(wèi)星陸表溫度產(chǎn)品(Land Surface Temperature, LST),由于其較好的質(zhì)量、誤差在1 K之內(nèi)[52],已在SUHI監(jiān)測中得到廣泛應(yīng)用[34-35,38-39,53]。本研究選用京津冀地區(qū)2003—2014年全年MODIS/Terra和MODIS/Aqua白天和夜晚8天合成平均LST產(chǎn)品(MOD11A2,MYD11A2),來源于美國國家航空與航天局陸表數(shù)據(jù)中心網(wǎng)站(https://lpdaac.usgs.gov/data_access/data_pool),可免費(fèi)下載。利用LST有效像元均值法可計(jì)算得到2003—2014京津冀地區(qū)逐年平均及春季(3—5月)、夏季(6—8月)、秋季(9—11月)、冬季(12—2月)LST影像。

    (2)NOAA衛(wèi)星資料 NOAA衛(wèi)星具有與MODIS衛(wèi)星類似的熱紅外通道(第4通道(10.3—11.3μm)和第5通道(11.5—12.5μm)),空間分辨率為1 km左右,利用其反演的陸表溫度數(shù)據(jù)可以有效監(jiān)測城市熱島[54-56]。本文選用京津冀地區(qū)1994年NOAA11/AVHRR LST逐日產(chǎn)品數(shù)據(jù)集[57],LST的反演采用改進(jìn)型的Becker分裂窗方法,反演結(jié)果與國際上廣泛應(yīng)用的MODIS陸表溫度產(chǎn)品具有較好的一致性[58]。

    (3)DMSP/OLS衛(wèi)星資料 美國軍事氣象衛(wèi)星Defense Meteorological Satellite Program(DMSP)搭載的Operational Linescan System(OLS)傳感器能夠探測到城市燈光甚至小規(guī)模居民地、車流等發(fā)出的低強(qiáng)度燈光,這一特征使其系列數(shù)據(jù)產(chǎn)品已被廣泛而成功地應(yīng)用于城市化相關(guān)研究[59-61]。其中DMSP/OLS數(shù)據(jù)灰度值(Digital Number,DN值有效范圍為1—63)高值區(qū)多用來提取城市建成區(qū)反映人類活動影響范圍,使之明顯區(qū)別于低值區(qū)的黑暗鄉(xiāng)村背景[62-63],因此在本研究中可利用DMSP/OLS數(shù)據(jù)的低值區(qū)提取幾乎不受人類活動影響的黑暗鄉(xiāng)村農(nóng)田背景區(qū)域。本研究選用京津冀地區(qū)2013年DMSP/ OLS夜間燈光數(shù)據(jù),來源于美國國家航空與航天局國家環(huán)境信息數(shù)據(jù)中心網(wǎng)站(http://ngdc.noaa.gov/eog/dmsp/downloadV4composites.html#AVSLCFC),可免費(fèi)下載。

    (4)植被指數(shù)資料 遙感估算的歸一化植被指數(shù)(NDVI)定義為近紅外波段與可見光波段數(shù)值之差和這兩個波段數(shù)值之和的比值,是植被生長狀態(tài)及植被覆蓋的良好指示因子。本文選用京津冀地區(qū)2013年1 km空間分辨率MODIS月歸一化植被指數(shù)NDVI產(chǎn)品(MOD13A3),來源于美國國家地質(zhì)調(diào)查局?jǐn)?shù)據(jù)陸表數(shù)據(jù)中心網(wǎng)站(https://lpdaac.usgs.gov/data_access/data_pool),可免費(fèi)下載。采用最大值法可獲取2013年最大NDVI影像。

    (5)社會經(jīng)濟(jì)資料 京津冀城市群11個城市2013年社會經(jīng)濟(jì)數(shù)據(jù),包括國內(nèi)生產(chǎn)總值、城鎮(zhèn)人口數(shù)、電力用量、機(jī)動車保有量,來源于當(dāng)?shù)亟y(tǒng)計(jì)局。

    (6)其他資料 京津冀地區(qū)90 m空間分辨率的數(shù)字高程影像(DEM),來源于國際科學(xué)數(shù)據(jù)服務(wù)平臺((http://datamirror.csdb.cn/dem/search.jsp),可免費(fèi)下載。

    2.2 研究方法

    2.2.1 SUHI估算

    本文選用城鄉(xiāng)二分法,由于UHI是一個城市高于鄉(xiāng)村的相對溫度距平值,而這個距平值無論是地表城市熱島,還是氣象熱島,至今沒有統(tǒng)一的閾值標(biāo)準(zhǔn)[1,34]。在這里,采用葉彩華等[36]提出的SUHI指標(biāo)及等級劃分方法來估算遙感地表熱島強(qiáng)度,即:

    (1)

    式中,SUHIi為圖象上第i個象元所對應(yīng)的熱島強(qiáng)度(℃),Ti是陸表溫度(℃),n為鄉(xiāng)村背景內(nèi)的有效象元數(shù),Tcrop為鄉(xiāng)村背景內(nèi)的陸表溫度(℃)。并按熱島強(qiáng)度值的大小,劃分為7級熱島強(qiáng)度:強(qiáng)冷島、較強(qiáng)冷島、弱冷島、無熱島、弱熱島、較強(qiáng)熱島和強(qiáng)熱島等,分別賦值為1、2、3、4、5、6和7。SUHIi的具體等級劃分參見表1,本文主要為季、年SUHI等級劃分。

    表1 地表熱島強(qiáng)度(SUHI)等級劃分及含義

    2.2.2 鄉(xiāng)村背景的確定

    按照城市熱島的定義,SUHI是城市與周邊鄉(xiāng)村或郊區(qū)陸表溫度的差異,因此對大尺度空間區(qū)域內(nèi)城市群的多個城市不能用同一個鄉(xiāng)村背景值來計(jì)算各城市熱島強(qiáng)度。考慮估算指標(biāo)的簡便性、易獲取性和通用性,本文確定4個指標(biāo)來選取鄉(xiāng)村背景區(qū)域。

    (1)平原(平壩) 與城市所在區(qū)域高程差別≤50 m,通過該閾值的設(shè)置可以忽略由于地形高度造成的兩地溫度差異,減少由于地形引起的熱島強(qiáng)度誤差,該指標(biāo)可通過當(dāng)?shù)財(cái)?shù)字高程圖像(DEM)計(jì)算得到。

    (2)農(nóng)田類型 大多研究采用城鄉(xiāng)二分法計(jì)算城市熱島強(qiáng)度時(shí)選擇的鄉(xiāng)村或郊區(qū)背景參考區(qū)域一般都是農(nóng)田地區(qū),本文繼續(xù)遵循這一原則。該指標(biāo)一般可通過當(dāng)?shù)赝恋乩妙愋瞳@取。本文采用2014年京津冀地區(qū)MODIS土地覆蓋分類產(chǎn)品(MCD12Q1),產(chǎn)品來源于美國國家地質(zhì)調(diào)查局?jǐn)?shù)據(jù)陸表數(shù)據(jù)中心網(wǎng)站(https://lpdaac.usgs.gov/data_access/data_pool),可免費(fèi)下載?;谠撏恋馗采w類型產(chǎn)品,結(jié)合2007年華北地區(qū)的植被覆蓋分類圖[64]與2011年京津唐地區(qū)土地利用類型[22]以及部分Landsat-TM高分辨率遙感資料,進(jìn)行了部分類別如城鎮(zhèn)用地、農(nóng)田類型的目視解譯修正,總體分類精度在90%以上,可得到京津冀地區(qū)2013年土地利用類型(圖2)。

    圖2 京津冀地區(qū)2013年土地利用類型、夜間燈光強(qiáng)度指數(shù)值、年最大植被指數(shù)及提取的郊區(qū)農(nóng)田背景Fig.2 Landuse types,DMSP/OLS DN values,annual maximal NDVI and rural crop background areas in Beijing-Tianjin-Hebei region in 2013

    (3)夜間燈光強(qiáng)度指數(shù)值≤15,通過該指標(biāo)可以確定基本不受人類活動影響(可以由夜間燈光反映)的區(qū)域,包含了一般不存在熱島效應(yīng)的鄉(xiāng)村,該指標(biāo)可以從美國衛(wèi)星DMSP/OLS資料中的灰度值(Digital Number,DN值有效范圍為1—63)計(jì)算得到,如圖2為提取的京津冀地區(qū)2013年夜間燈光指數(shù)圖,與土地利用圖2中的城鎮(zhèn)用地區(qū)域相比,夜間燈光指數(shù)高值區(qū)(指數(shù)值在50—63)面積明顯大于城鎮(zhèn)用地區(qū)域,顯示城市與衛(wèi)星小城鎮(zhèn)周邊郊區(qū)雖然也是農(nóng)田植被,但受城鎮(zhèn)化影響較大,已不適合作為郊區(qū)農(nóng)田背景,因此需要在夜間燈光指數(shù)低值區(qū)來選擇鄉(xiāng)村背景農(nóng)田。

    (4)年最大植被指數(shù)NDVI≥0.7,通過該指標(biāo)可以確定不包含水體且為高植被覆蓋的純植被,也排除了城市周邊郊區(qū)的大量鄉(xiāng)鎮(zhèn)及衛(wèi)星城鎮(zhèn)等低植被覆蓋人類活動區(qū),可以根據(jù)MODIS月歸一化植被指數(shù)產(chǎn)品(MOD13A3)求取年最大值獲得,如圖2所示:城鎮(zhèn)地區(qū)和沿海水體NDVI值較低,一般在0.50以下,而山區(qū)林地和平原農(nóng)田NDVI值較高,普遍在0.60以上,值得注意的是西部太行山山前的植被NDVI普遍不高,一般在0.40—0.60之間。

    根據(jù)上述4個指標(biāo),可提取得到1 km空間分辨率尺度上的京津冀地區(qū)的鄉(xiāng)村農(nóng)田背景區(qū)域,如圖2d所示:京津冀地區(qū)北部張家口市、承德市大部分地區(qū)處于山區(qū)或山地,幾乎沒有合適的郊區(qū)農(nóng)田背景,而北京由于城鎮(zhèn)化程度高,僅在南部與東部地區(qū)有少量的局部農(nóng)田能作為北京城市的郊區(qū)農(nóng)田背景,而其他城市地區(qū)由于地處平原,農(nóng)田面積大,城鎮(zhèn)化程度不高,有大量的農(nóng)田可入選為郊區(qū)農(nóng)田背景。

    2.2.3 城市熱島比例指數(shù)

    在這里,采用葉彩華等[36]提出的地表熱島比例指數(shù)(Urban Heat Island Proportion Index, UHPI)的計(jì)算方法來定量估算各城市空間單元內(nèi)的熱島強(qiáng)度,可以比較各城市不同時(shí)空尺度下的熱島強(qiáng)度大小。

    (2)

    式中,UHPI為城市熱島比例指數(shù);m為熱島強(qiáng)度等級數(shù);i為城區(qū)溫度高于郊區(qū)溫度等級序號;n為城區(qū)溫度高于郊區(qū)溫度的等級數(shù);wi為第i級的權(quán)重,取等級值,pi為第i級所占的面積百分比。UHPI值在0—1.0之間,該值越大,熱島現(xiàn)象越嚴(yán)重。其值為0時(shí),表明此地沒有熱島現(xiàn)象,值為1時(shí),表明此地均處于強(qiáng)熱島范圍。由前面定義的熱島強(qiáng)度等級可知,m=7,n=3。

    2.2.4 空間單元熱島強(qiáng)度評估

    (1)SUHI均值評估

    對各城市采用中心城區(qū)空間單元的SUHI平均值(單位為℃)定量評估該城市的熱島強(qiáng)度大小。

    (2)熱島面積評估

    對各城市區(qū)域范圍內(nèi)統(tǒng)計(jì)弱熱島、較強(qiáng)熱島及強(qiáng)熱島面積(單位為km2)及占行政區(qū)域內(nèi)面積百分比,可以有效評估城市化影響的熱島范圍大小。

    (3)UHPI評估

    由式(2)可知,UHPI通過具有熱島效應(yīng)的多個等級熱島強(qiáng)度在空間單元所占的面積百分比加權(quán)和綜合反映該城市空間單元熱島強(qiáng)度大小,較SUHI均值和熱島面積評估指標(biāo)更具客觀性與科學(xué)性,通過對UHPI分級可對各城市空間單元總體熱島狀況進(jìn)行定量評估(表2)。

    表2 城市空間單元熱島評估等級劃分

    可以看出,SUHI均值、熱島面積、UHPI估算的區(qū)域并不一致,SUHI均值針對所有等級熱島范圍的均值狀況進(jìn)行評估,而熱島面積、UHPI則主要從該區(qū)域內(nèi)產(chǎn)生的熱島范圍進(jìn)行總體評估。由于UHPI兼具了熱島等級和熱島范圍特征,較SUHI均值和熱島面積更具科學(xué)性與合理性,本文主要評估中心城區(qū)(Central Urban Area,簡稱CA)內(nèi)的UHPI。評估熱島強(qiáng)度大小的目的是分析城市建筑等人類活動影響所引起的城市氣候環(huán)境變化,因此對于海拔較高溫度明顯低于附近平原受地形影響的山區(qū)和山地,不符合做熱島評估要求,因此在估算各城市空間單元UHPI時(shí),需要排除山區(qū)、山地等區(qū)域,只對平原地區(qū)進(jìn)行評估。

    在這里統(tǒng)一用字母簡化文中出現(xiàn)的各種評估指標(biāo):熱島面積(指弱熱島以上等級面積,Heat Island Area),簡稱HA;較強(qiáng)熱島以上面積(Relative High Above Heat Island Area),簡稱RHA;強(qiáng)熱島面積(High Above Heat Island Area),簡稱HHA;中心城區(qū)平均SUHI,簡稱SUHI_CA;中心城區(qū)UHPI,簡稱UHPI_CA;熱島評估等級(Assessment Grades),簡稱AG。

    3 結(jié)果與分析

    3.1 年均熱島分析

    利用2010—2014年MODIS衛(wèi)星逐日4個時(shí)次的LST產(chǎn)品可計(jì)算得到近5年京津冀地區(qū)全年、白天和夜晚平均SUHI圖,如圖3所示。結(jié)合土地利用圖(圖2)分析可知:就全年平均SUHI(圖3)來看,京津冀城市群大部分平原農(nóng)田地區(qū)處于無熱島狀態(tài),大部分山區(qū)林地則呈現(xiàn)冷島現(xiàn)象,局部山區(qū)例如河北西部的太行山山前地區(qū)受地形和較低植被覆蓋影響,出現(xiàn)弱熱島現(xiàn)象;各大中城市地區(qū)均出現(xiàn)熱島效應(yīng),部分特大以上城市如北京、天津和石家莊等城市出現(xiàn)強(qiáng)熱島。就全年白天平均SUHI(圖3)來看,熱島空間分布趨勢與全年一致,但出現(xiàn)熱島的區(qū)域以弱熱島為主,部分大中城市中心局部地區(qū)會出現(xiàn)較強(qiáng)熱島。就全年夜晚平均SUHI(圖3)來看,熱島空間分布趨勢與全年基本一致,但各城市熱島強(qiáng)度等級和范圍明顯大于白天和全年平均。其中水體例如北京的密云水庫和天津沿海等區(qū)域具有明顯熱島效應(yīng),這與未結(jié)冰的水體在夜晚具有明顯熱島效應(yīng)這一結(jié)論[55]具有一致性。對比圖3可知,京津冀地區(qū)各城市年平均夜晚熱島強(qiáng)度明顯大于白天,這與眾多學(xué)者利用氣象資料對北京、天津和石家莊等城市熱島研究的“夜晚城市熱島強(qiáng)度大于白天”具有一致性[65-68],表明利用本文建立的SUHI方法可有效監(jiān)測京津冀城市群年均熱島狀況。

    圖3 京津冀地區(qū)2010—2014年均全年、白天、夜晚SUHIFig.3 Annual mean SUHI for all year, day, night in Beijing-Tianjin-Hebei region in 2010—2014

    熱島評估表明(表3):就年均熱島面積來看,京津冀城市群總的RHA為1926 km2,占行政區(qū)域總面積的1.5%,其中排名前3位的分別是天津(640 km2)、北京(433 km2)和唐山(243 km2),石家莊與邯鄲分列第四(223 km2)和第五(134 km2)。就SUHI_CA來看,排名前3位的分別是石家莊(3.6℃)、保定(3.3℃)和邢臺(2.6℃),北京和天津并列第4(2.3℃),廊坊、衡水、秦皇島都比較低(0.10℃以下)。就UHPI_CA來看,排名前3位的與SUHI_CA排名相同,分別是石家莊(0.82)、保定(0.74)和邢臺(0.72),天津和北京分別為0.64和0.62,此外邯鄲、唐山、滄州也在0.50以上??梢钥闯?SUHI_CA排名與UHPI_CA排名基本相同,但北京和天津的SUHI_CA一樣,而天津的UHPI_CA卻略大于北京,同時(shí)參考RHA值,表明天津的年均熱島很可能較北京嚴(yán)重。城市群各城市UHPI_CA平均值為0.50,總體AG為3級(較嚴(yán)重),其中5級(非常嚴(yán)重)和4級(嚴(yán)重)的城市有6個,占到54%,表明京津冀城市群中心城區(qū)年均熱島總體都比較嚴(yán)重。

    表3 京津冀各城市年平均RHA、SUHI_CA、UHPI_CA和AG對比(2010—2014年)

    3.2 季節(jié)熱島分析

    利用2010—2014年MODIS衛(wèi)星逐日LST產(chǎn)品可計(jì)算得到近5年京津冀地區(qū)春、夏、秋、冬四季白天和夜晚平均SUHI圖,如圖4和圖5所示。

    圖4 京津冀地區(qū)2010—2014春季、夏季、秋季、冬季平均白天SUHIFig.4 Seasonal day SUHI for spring,summer,autumn,winter in Beijing-Tianjin-Hebei area in 2010—2014

    對白天而言(圖4),京津冀地區(qū)各城市四季SUHI分布差異明顯。春季白天熱島空間差異最為明顯(圖4),大部分城市城鎮(zhèn)地區(qū)存在熱島效應(yīng),北京以弱熱島為主,天津、石家莊出現(xiàn)較強(qiáng)以上熱島,但許多平原農(nóng)田地區(qū)以及西部太行山山前地帶存在明顯熱島區(qū)域,這主要是由于京津冀地區(qū)春季降水少,平原既分布著大量的較高植被覆蓋的冬小麥(低溫冷島區(qū)),又分布著大量的白天受熱升溫較快的干燥裸土區(qū)(高溫?zé)釐u區(qū)),而太行山山前受地形遮擋春季西北冷空氣影響加之植被覆蓋較低存在較多干燥裸土區(qū)也存在一高溫?zé)釐u區(qū)。夏季白天(圖4),植被覆蓋普遍較高,大部分平原農(nóng)田為無熱島區(qū)域,城鎮(zhèn)地區(qū)顯現(xiàn)出明顯熱島效應(yīng),北京、天津、唐山、石家莊、邯鄲、滄州、邢臺等城市均有明顯強(qiáng)熱島發(fā)生,其中強(qiáng)熱島面積排名前3的分別是北京(1501 km2)、天津(1198 km2)和唐山(362 km2),石家莊排第4(307 km2)。秋季白天(圖4),大部分城市城鎮(zhèn)所在區(qū)域以弱熱島為主,局部有較強(qiáng)熱島。冬季白天(圖4),與春季、夏季、秋季不同,大部分城市城鎮(zhèn)所在區(qū)域出現(xiàn)冷島效應(yīng),這是由于冬季白天,郊區(qū)地面裸露土壤較為干燥,蒸發(fā)量較小,其升溫速度高于城區(qū)地表;另一方面原因是冬季較為嚴(yán)重的城區(qū)污染物對太陽輻射有散射和吸收作用,很大程度削弱了到達(dá)地表的太陽輻射[26]。因此對白天而言,京津冀大部分城市熱島在城鎮(zhèn)區(qū)域內(nèi)均有夏季最強(qiáng)、冬季最弱特點(diǎn),這與眾多學(xué)者利用遙感監(jiān)測城市白天熱島結(jié)果具有一致性[26,34,49]。但這與氣象觀測熱島不一致,主要原因是白天受太陽輻射影響,衛(wèi)星觀測與氣象觀測的城市熱島原理不同,兩者相關(guān)關(guān)系普遍較差[25,38]。

    對夜晚而言(圖5),京津冀地區(qū)各城市四季熱島分布較為相似,各城鎮(zhèn)地區(qū)均存在明顯熱島效應(yīng),水體具有“暖島效應(yīng)”。春季夜晚SUHI(圖5),多數(shù)城市均出現(xiàn)較強(qiáng)以上熱島,局部強(qiáng)熱島,RHA排名前3位的分別是天津(1555 km2)、北京(779 km2)和唐山(730 km2),石家莊為232 km2。夏季夜晚SUHI(圖5),大部分城市以弱熱島為主,熱島范圍也明顯小于其他季節(jié),RHA排名前3位的分別是天津(770 km2)、北京(585 km2)和唐山(272 km2),石家莊為160 km2。秋季夜晚SUHI(圖5)與春季分布相似,大部分城鎮(zhèn)地區(qū)出現(xiàn)熱島,RHA排名前3位的分別是天津(1729 km2)、北京(958 km2)和唐山(632 km2),石家莊為276 km2。冬季夜晚SUHI(圖5),大部分城市出現(xiàn)強(qiáng)熱島,熱島強(qiáng)度及范圍也明顯大于其他季節(jié),RHA排名前3位的分別是天津(2033 km2)、北京(1186 km2)和唐山(773 km2),石家莊為310 km2。值得一提的是四季夜晚天津的較強(qiáng)熱島面積均為各城市中最大,這與夜晚天津沿海分布著大面積的水體具有“暖島”效應(yīng)有關(guān)。對比四季SUHI(圖5)可知,對夜晚SUHI而言,京津冀大部分城市在城鎮(zhèn)區(qū)域內(nèi)均有冬季最強(qiáng)、夏季最弱特點(diǎn)。這與北京市氣象局氣候資料室[69]及謝莊等[65]利用氣溫資料對北京地區(qū)的研究結(jié)論相似:即冬季熱島最強(qiáng),夏季最弱。這是由于衛(wèi)星遙感監(jiān)測的地溫與氣溫在夜間具有良好的正相關(guān)性,而這種高度相關(guān)可以由夜間地表可以看作一個均質(zhì)的單一熱源來解釋[70]。即對城市熱島監(jiān)測而言,遙感監(jiān)測只有在夜晚才與氣溫監(jiān)測具有較好的一致性[22]。

    圖5 京津冀地區(qū)2010—2014春季、夏季、秋季、冬季平均夜晚SUHIFig.5 Seasonal night SUHI for spring,summer,autumn,winter in Beijing-Tianjin-Hebei region in 2010—2014

    由圖4和圖5對比可知,不考慮京津冀地區(qū)地形對熱島影響和沿海水體夜晚的“暖島效應(yīng)”,夏季白天各城市SUHI強(qiáng)度和范圍最大,而且夏季白天也是高溫易發(fā)對人居環(huán)境及人體健康影響最大的時(shí)段,在這里對夏季白天各城市RHA、較強(qiáng)以上熱島面積比例(占該行政區(qū)域面積的百分比,簡稱P_RHA)、SUHI_CA、UHPI_CA及AG進(jìn)行分析,如表4所示:RHA排名前三的分別是北京(2351 km2)、天津(1883 km2)和唐山(889 km2),石家莊排第四(611 km2),分別占該行政區(qū)域面積的14%、16%、7%和4%,其余各城市RHA均在400 km2以下,面積百分比為4%以下,整個京津冀城市群RHA為7386 km2,占行政區(qū)域面積的5.8%;11個平原城市SUHI_CA和UHPI_CA平均值分別為3.0℃和0.61,但排名前三的SUHI_CA和UHPI_CA城市并不一致,前三的SUHI_CA是保定(4.5℃)、北京和石家莊(并列第二,4.3℃)、天津(3.6℃),前三的UHPI_CA為石家莊(0.87)、保定(0.83)、北京和邢臺(并列,0.81),UHPI評估為5級(非常嚴(yán)重);唐山、天津、邯鄲也在0.60以上。就AG而言,5級(非常嚴(yán)重)的城市有4個,4級(嚴(yán)重)的有3個,3級(較嚴(yán)重)的有1個,2級(一般)的有3個,3級(較嚴(yán)重)以上城市占到73%,表明夏季白天京津冀城市群各中心城區(qū)熱島總體都很嚴(yán)重。

    3.3 熱島時(shí)空變化現(xiàn)狀分析

    由前面分析可知,夏季白天是京津冀城市群SUHI最強(qiáng)的時(shí)段,利用1994年NOAA11、2004年和2014 MODIS衛(wèi)星資料可估算京津冀地區(qū)不同年代夏季白天SUHI_CA、UHPI_CA和HHA(圖6)。

    圖6 京津冀地區(qū)不同年代1994年、2004年、2014年夏季白天SUHI變化Fig.6 SUHI temporal change in summer day in 1994,2004,2014 in Beijing-Tianjin-Hebei region

    1994年 各城市熱島強(qiáng)度和范圍比較弱(圖6),以弱熱島和較強(qiáng)熱島為主,大部分UHPI_CA在0.60以下,排名前三位的分別是石家莊(0.67)、唐山(0.46)、北京(0.37),大部分城市HHA很小,一般在100 km2以下。

    2004年 各城市熱島強(qiáng)度和范圍明顯增加,出現(xiàn)明顯強(qiáng)熱島區(qū)域(圖6),各UHPI_CA明顯增加,UHPI_CA在0.60以上的城市占到45%,各城市HHA均增加2倍以上,其中北京、天津、石家莊HHA分別增加到866、373、189 km2。

    2014年 各城市熱島強(qiáng)度和范圍較2004年繼續(xù)擴(kuò)大(圖6),各城市UHPI_CA呈持續(xù)增加趨勢,UHPI_CA在0.60以上的城市上升到73%,其中排名前三的分別是石家莊(0.93)、保定(0.88)、北京(0.85)與邢臺(0.85),AG均達(dá)5級(非常嚴(yán)重);天津、唐山、滄州、邯鄲UHPI_CA分別為0.79、0.78、0.73和0.66,AG均達(dá)4級(嚴(yán)重),表明2014年京津冀城市群大部分中心城區(qū)夏季白天熱島已經(jīng)很嚴(yán)重;HHA也迅速增加,北京、天津分別已達(dá)到1346km2和1342km2,明顯高于其他城市,石家莊排第三,達(dá)562 km2,唐山也達(dá)486 km2,保定、廊坊、邢臺、邯鄲的HHA為124—254 km2,其他城市在100 km2以下,該結(jié)果也顯示出超大城市(北京、天津)、特大城市(石家莊)、資源型城市(唐山)貢獻(xiàn)了京津冀城市群大部分熱島面積。

    表4 京津冀各城市夏季白天RHA、P_RHA、SUHI_CA、UHPI_CA和AG對比(2010—2014)

    對比3個時(shí)期的熱島變化分析看出,經(jīng)過20a的城市發(fā)展,11個平原城市夏季白天SUHI_CA平均值為1.2℃,UHPI_CA平均值在1994年僅為0.29,AG為2級(一般),總的HHA為190 km2;到2004年SUHI_CA和UHPI_CA平均值分別增加到2.6℃和0.58,AG為3級(較嚴(yán)重),總的HHA增加到1975 km2;到2014年SUHI_CA和UHPI_CA平均值分別增加到3.2℃和0.69,AG為4級(嚴(yán)重),總的HHA增加到4539 km2,占行政區(qū)域總面積的3.6%,這種變化反映了京津冀城市群20年來城市擴(kuò)展導(dǎo)致的熱島效應(yīng)增強(qiáng)增大事實(shí)。

    3.4 熱島時(shí)空變化趨勢

    利用MODIS衛(wèi)星資料估算2003—2014逐年城市年均熱島面積(HA)與夏季白天強(qiáng)熱島面積(HHA),對熱島面積時(shí)間變化進(jìn)行趨勢線性擬合(模型達(dá)到0.05顯著水平),可以得到京津冀各城市近12年熱島面積變化量:京津冀各城市HA呈增加趨勢,平均增加45 km2/a,除了秦皇島增加很少(2 km2/a),其余城市增加23—86 km2/a,天津HA增加最大(86 km2/a),其次是北京(61 km2/a),石家莊為43 km2/a。對夏季白天HHA而言,京津冀各城市平均增加15 km2/a,但主要增加在天津和北京,分別為69 km2/a和63 km2/a,唐山增加12km2/a,其余城市在6 km2/a以下,而衡水和秦皇島幾乎沒有增加。

    從前面分析可知,不考慮受河北太行山山前地區(qū)影響的熱島區(qū)域,無論時(shí)間上還是空間上,京津冀城市群熱島強(qiáng)度與范圍呈持續(xù)增加趨勢,其中北京城市熱島范圍呈現(xiàn)向東、南、背三面擴(kuò)展的趨勢,而天津熱島呈向東北、東擴(kuò)展趨勢。根據(jù)主要強(qiáng)熱島區(qū)域最外圍位置估算,北京與天津強(qiáng)熱島最近空間距離1994年為94 km(圖6),到2004年縮減到79 km(圖6),到2014年則已縮減到52 km(圖6),雖然兩個超大城市中間隔著的廊坊城市熱島范圍擴(kuò)展并不十分明顯,但由京津冀夜間燈光指數(shù)圖(圖2)可知廊坊的城鎮(zhèn)化也非常明顯,隨著未來北京城市重心東移以及東南部的北京第二國際機(jī)場建設(shè),將會大大促進(jìn)周邊廊坊市的城鎮(zhèn)發(fā)展,北京、廊坊、天津大范圍熱島連片的可能性大大增加,這很可能將給京津冀生態(tài)圈環(huán)境帶來惡劣后果。許多研究指出增加植被覆蓋是有效防止城市熱島的重要措施[16,46,58,63],因此,在未來的京津冀一體化發(fā)展城市規(guī)劃中,需要在北京與天津之間建立有效的綠色生態(tài)屏障,防止大面積熱島連片成為超大的“京津區(qū)域熱島群”。

    3.5 經(jīng)濟(jì)驅(qū)動因子對SUHI的影響

    利用京津冀地區(qū)各城市近5年平均(2010—2014)夏季白天RHA(km2)與2013各城市國內(nèi)生產(chǎn)總值(GDP(萬元))、城鎮(zhèn)人口數(shù)(UP(萬人))、用電量(EC(kW·h))、機(jī)動車數(shù)量(VP(萬輛))進(jìn)行線性擬合分析(圖7)。京津冀城市群RHA與GDP、UP、EC、VP存在顯著正相關(guān)關(guān)系,模型決定系數(shù)R2分別為0.9097、0.912、0.9661和0.7098(均達(dá)到0.01顯著水平),顯示京津冀城市群經(jīng)濟(jì)規(guī)模、人口規(guī)模、用電量和機(jī)動車數(shù)量與城市熱島有明顯的正相關(guān)關(guān)系,顯示了這些社會經(jīng)濟(jì)驅(qū)動因子對城市熱島大小的重要影響,意味著控制人口規(guī)模和減少能耗可有效緩解熱島效應(yīng)。

    一般來說,經(jīng)濟(jì)規(guī)模、電力消耗、機(jī)動車數(shù)量均是以一定的城鎮(zhèn)人口規(guī)模為基礎(chǔ),根據(jù)城鎮(zhèn)人口數(shù)與熱島面積之間模型(圖7)擬合結(jié)果,可估算出京津冀Ⅱ型大城市、Ⅰ型大城市、特大城市、超大城市的RHA分別為38—280、280—524、524—1133、≥1133 km2,這為京津冀城市群內(nèi)各城市發(fā)展時(shí),如何通過控制城市規(guī)模來避免當(dāng)?shù)爻鞘袩釐u大小提供了一個參考。

    圖7 2013年京津冀各城市GDP、UP、EC、VP與2010—2014年夏季白天RHA_HHA關(guān)系Fig.7 Relation between GDP,UP,EC,VP of cities in 2013 and the RHA_HHA in summer day in 2010—2014P<0.01表示線性擬合模型達(dá)到0.01顯著水平

    同時(shí)為避免未來“京津區(qū)域熱島群”出現(xiàn)的可能,一方面需要控制或減少超大城市的城鎮(zhèn)人口規(guī)模來控制熱島面積,在2015年4月我國發(fā)布的《京津冀協(xié)同發(fā)展規(guī)劃綱要》中已提到解決“大城市病”的一條重要措施就是人口疏解,這有利于緩解超大城市的熱島效應(yīng);另一方面則需要做好京津冀一體化的城市規(guī)劃,尤其是城市綠地和土地利用規(guī)劃,已有研究表明通過城市的土地利用有效規(guī)劃能明顯緩解熱島效應(yīng)[56]。

    4 結(jié)論與討論

    利用遙感定量表達(dá)城市熱島強(qiáng)弱是個復(fù)雜的科學(xué)問題,其中關(guān)于郊區(qū)或鄉(xiāng)村背景以及城市區(qū)選擇來估算熱島強(qiáng)度是一個重要的技術(shù)問題。本文建立的SUHI估算方法能有效監(jiān)測京津冀城市群熱島,各大中城市熱島時(shí)空分布與前人的研究結(jié)果相似,城市群SUHI不但存在年/季、晝/夜變化,還受土地利用類型、地形影響,例如平原農(nóng)田一般無熱島,大部分山區(qū)林地呈冷島,水體在夜間呈“暖湖效應(yīng)”,河北太行山山前地區(qū)易出現(xiàn)熱島。另外,雖然在夜晚由衛(wèi)星遙感估算的SUHI與氣溫?zé)釐u具有較好一致性,但由于白天陸表溫度與氣溫形成原理上的差異,使得兩種方式估算的熱島強(qiáng)度并不具有一致性,許多研究也提到這種明顯差異[25,39- 40],同時(shí)SUHI(或UHI)的分級也沒有統(tǒng)一的標(biāo)準(zhǔn),這給城市熱島強(qiáng)度的合理、準(zhǔn)確描述及城市熱島比較分析帶來困惑。因此,如何統(tǒng)一、協(xié)調(diào)好二者來更加合理、準(zhǔn)確表達(dá)一個城市的熱島強(qiáng)度是今后值得研究的內(nèi)容。另外,本文提出的熱島強(qiáng)度估算方法在屬于一個氣候帶里的京津冀城市群進(jìn)行了有效應(yīng)用,由于熱島強(qiáng)度是一個相對溫度距平指標(biāo)值,它考慮了一定范圍內(nèi)(即氣候條件相似)地形、土地利用類型、植被覆蓋、人類活動因素,且各因素指標(biāo)容易獲取,因此該方法可在其他氣候類型區(qū)進(jìn)行推廣應(yīng)用。但受氣候與地形差異,不同氣候帶城市熱島很可能呈現(xiàn)與京津冀城市群不一致的時(shí)空分布,如楊英寶利用衛(wèi)星資料[71]對處于亞熱帶季風(fēng)氣候區(qū)的南京城市熱島研究發(fā)現(xiàn)熱島強(qiáng)度最高值出現(xiàn)在秋季,李曉敏[72]對同處于亞熱帶季風(fēng)氣候區(qū)但又受盆地影響的成都城市熱島研究發(fā)現(xiàn)熱島強(qiáng)度最高值出現(xiàn)在春季,而處于溫帶海洋性氣候的荷蘭鹿特丹城市熱島強(qiáng)度最高值出現(xiàn)在冬季[38]。因此,針對不同氣候類型地區(qū)的熱島時(shí)空格局還需要考慮多種因素進(jìn)一步的研究。

    對京津冀城市群熱島定量監(jiān)測評估表明:城市群近5年(2010—2014)年均較強(qiáng)以上熱島面積RHA(SUHI≥3℃)為1926 km2,但在熱島最強(qiáng)的夏季白天,RHA達(dá)到7386 km2(占行政區(qū)域面積的5.8%),排名前三的分別是北京(2351 km2)、天津(1087 km2)和唐山(889 km2),石家莊排第四(611 km2);11個平原城市中心城區(qū)平均SUHI和UHPI分別為3.0℃和0.61,熱島評估較嚴(yán)重(UHPI大于0.6)以上城市占到73%,評估結(jié)果顯示城市群各中心城區(qū)熱島總體都很嚴(yán)重,同時(shí)顯示超大城市(北京、天津)、特大城市(石家莊)和以重化工業(yè)代表的資源型城市(唐山)占據(jù)了城市群熱島面積的大部分,這對未來京津冀城市群一體化發(fā)展,如何控制城市規(guī)模、把握城市定位、保護(hù)生態(tài)環(huán)境提供一個參考。

    對比3個不同時(shí)期發(fā)現(xiàn),京津冀城市群11個平原城市夏季白天SUHI和UHPI平均值,在1994年僅為1.2℃和0.29(熱島評估為“一般”等級),強(qiáng)熱島面積HHA(SUHI≥5℃)為190 km2;到2004年分別增加到2.6℃和0.58(熱島評估為“較嚴(yán)重”等級),HHA增加到1975 km2;到2014年則增加到3.2℃和0.69(熱島評估為“嚴(yán)重”等級),HHA增加到4539 km2,占行政區(qū)域總面積的3.6%,這種變化反映了京津冀城市迅速擴(kuò)展導(dǎo)致的熱島增強(qiáng)增大事實(shí)。同時(shí)分析表明北京、天津強(qiáng)熱島區(qū)的空間距離日趨縮小,未來隨著北京城市重心東移和北京第二機(jī)場的建設(shè),北京、天津兩個強(qiáng)熱島區(qū)存在集中連片為“京津區(qū)域熱島群”的可能,這是未來京津冀城市群“一體化”發(fā)展需要極力避免的重大生態(tài)環(huán)境問題,可在京津之間建立“綠色生態(tài)屏障”有效措施防止熱島蔓延趨勢。

    社會經(jīng)濟(jì)驅(qū)動因子分析表明,京津冀城市經(jīng)濟(jì)規(guī)模(GDP)、人口數(shù)量和用電量都極大地影響著京津冀各城市熱島大小(擬合模型決定系數(shù)R2均在0.90以上),意味著控制城市規(guī)模、減少能源消耗可有效減緩熱島效應(yīng),這為《京津冀協(xié)同發(fā)展規(guī)劃綱要》中提到“人口疏解”和“促進(jìn)綠色循環(huán)低碳發(fā)展”以及合理地規(guī)劃城市規(guī)模大小以減緩熱島并避免“京津熱島群”出現(xiàn)的可能提供了參考。同時(shí),城市熱島形成的機(jī)制很復(fù)雜,本文只是從社會經(jīng)濟(jì)驅(qū)動因子來描述其對城市熱島的可能影響,但并沒有從機(jī)理上也難以從衛(wèi)星遙感的監(jiān)測現(xiàn)狀去解釋各種因素對熱島的具體影響過程,而數(shù)值模擬則為熱島形成機(jī)制提供了一種較好解決手段[15-16,50-51,53]。因此通過衛(wèi)星遙感獲取精細(xì)化地表參數(shù),融合數(shù)值模式和能量平衡模型,結(jié)合實(shí)地試驗(yàn)觀測數(shù)據(jù)來研究城市熱島的形成機(jī)制是未來的重要研究方向,這必將對城市的熱島效應(yīng)有深入理解,促進(jìn)各城市采取有力措施以減緩城市熱島效應(yīng)給自然環(huán)境、區(qū)域氣候造成的負(fù)面影響。

    [1] Voogt J A. Urban heat island//Douglas I, Ed. Encyclopedia of Global Environmental Change Volume 3: Causes and Consequences of Global Environmental Change. New York: John Wiley & Sons, 2002: 660- 666.

    [2] Rizwan A M, Dennis L Y C, Liu C H. A review on the generation, determination and mitigation of urban heat island. Journal of Environmental Sciences, 2008, 20(1): 120- 128.

    [3] 張光智, 徐祥德, 王繼志, 楊元琴. 北京及周邊地區(qū)城市尺度熱島特征及其演變. 應(yīng)用氣象學(xué)報(bào), 2002, 13(Z1): 43- 50.

    [4] 季崇萍, 劉偉東, 軒春怡. 北京城市化進(jìn)程對城市熱島的影響研究. 地球物理學(xué)報(bào), 2006, 49(1): 69- 77.

    [5] 郭勇, 龍步菊, 劉偉東, 李明. 北京城市熱島效應(yīng)的流動觀測和初步研究. 氣象科技, 2006, 34(6): 656- 661.

    [6] Liu W D, Ji C P, Zhong X, Jiang X, Zheng Z. Temporal characteristics of the Beijing urban heat island. Theoretical and Applied Climatology, 2007, 87(1/4): 213- 221.

    [7] 馬玉霞, 王式功, 魏海茹. 蘭州市近50年城市熱島強(qiáng)度變化特征. 氣象科技, 2009, 37(6): 660- 664.

    [8] 陳云浩, 史培軍, 李曉兵. 基于遙感和GIS的上海城市空間熱環(huán)境研究. 測繪學(xué)報(bào), 2002, 31(2): 139- 144.

    [9] Voogt J A, Oke T R. Thermal remote sensing of urban climates. Remote Sensing of Environment, 2003, 86(3): 370- 384.

    [10] 蘇偉忠, 楊英寶, 楊桂山. 南京市熱場分布特征及其與土地利用/覆被關(guān)系研究. 地理科學(xué), 2005, 25(6): 697- 703.

    [11] 宋艷暾, 余世孝, 李楠, 龔友夫. 深圳特區(qū)表面溫度空間分異特征. 生態(tài)學(xué)報(bào), 2007, 27(4): 1489- 1498.

    [12] 江學(xué)頂, 夏北成. 珠江三角洲城市群熱環(huán)境空間格局動態(tài). 生態(tài)學(xué)報(bào), 2007, 27(4): 1461- 1470.

    [13] 胡姝婧, 胡德勇, 李小娟, 王榮華. 京津冀都市圈熱環(huán)境空間格局遙感分析. 國土資源遙感, 2009, (3): 94- 99.

    [14] Li C F, Yin J Y. Study on urban thermal field of Shanghai using multi-source remote sensing data. Journal of the Indian Society Remote Sensing, 2013, 41(4): 1009- 1019.

    [15] 楊玉華, 徐祥德, 翁永輝. 北京城市邊界層熱島的日變化周期模擬. 應(yīng)用氣象學(xué)報(bào), 2003, 14(1): 61- 68.

    [16] 佟華, 劉輝志, 李延明, 桑建國, 胡非. 北京夏季城市熱島現(xiàn)狀及楔形綠地規(guī)劃對緩解城市熱島的作用. 應(yīng)用氣象學(xué)報(bào), 2005, 6(3): 357- 366.

    [17] 江曉燕, 張朝林, 高華, 苗世光. 城市下墊面反照率變化對北京市熱島過程的影響——個例分析. 氣象學(xué)報(bào), 2007, 65(2): 301- 307.

    [18] Sarkar A, De Ridder K. The urban heat island intensity of Paris: A case study based on a simple urban surface parametrization. Boundary-Layer Meteorology, 2011, 138(3): 511- 520.

    [19] Stewart I D. A systematic review and scientific critique of methodology in modern urban heat island literature. International Journal of Climatology, 2011, 31(2), 200- 217.

    [20] De Lucena A J, Filho O C R, de Almeida Fran?a J R, de Faria Peres L. Urban climate and clues of heat island events in the metropolitan area of Rio de Janeiro. Theoretial and Applied Climatology, 2013, 111(3/4): 497- 511.

    [21] Martin P, Baudouin Y, Gachon P. An alternative method to characterize the surface urban heat island. International Journal of Biometeorology, 2015, 59(7): 849- 861.

    [22] 劉勇洪, 欒慶祖, 權(quán)維俊, 張碩. 基于多源衛(wèi)星資料的京津唐城市群熱環(huán)境研究. 生態(tài)環(huán)境學(xué)報(bào), 2015, 24(7): 1150- 1158.

    [23] He J F, Liu J Y, Zhuang D F, Zhang W, Liu M L. Assessing the effect of land use/land cover change on the change of urban heat island intensity. Theoretical and Applied Climatology, 2007, 90(3/4): 217- 226.

    [24] Gaffin S R, Rosenzweig C, Khanbilvardi R, Parshall L, Mahani S, Glickman H, Goldberg R, Blake R, Slosberg R B, Hillel D. Variations in New York city′s urban heat island strength over time and space. Theoretical and Applied Climatology, 2008, 94(1/2): 1- 11.

    [25] Mohan M, Kikegawa Y, Gurjar B R, Bhati S, Kolli N R. Assessment of urban heat island effect for different land use-land cover from micrometeorological measurements and remote sensing data for megacity Delhi. Theoretical and Applied Climatology, 2013, 112(3/4): 647- 658.

    [26] 王建凱, 王開存, 王普才. 基于MODIS地表溫度產(chǎn)品的北京城市熱島(冷島)強(qiáng)度分析. 遙感學(xué)報(bào), 2007, 11(3): 330- 339.

    [27] 鄧玉嬌, 匡耀求, 黃鋒. 基于Landsat/TM資料研究廣州城市熱島現(xiàn)象. 氣象, 2010, 36(1): 26- 30.

    [28] 林奕桐, 葉駿菲, 林開平, 張家安. 南寧市熱島效應(yīng)的遙感研究. 災(zāi)害學(xué), 2014, 29(4): 192- 197.

    [29] 王天星, 陳松林, 閻廣建. 地表參數(shù)反演及城市熱島時(shí)空演變分析. 地理科學(xué), 2009, 29(5): 697- 702.

    [30] 徐賽萍, 沙晉明, 劉玉琴, 祝坤子. 城市熱島效應(yīng)對植被生態(tài)環(huán)境質(zhì)量的響應(yīng)研究——以福州市馬尾區(qū)為例. 福建師范大學(xué)學(xué)報(bào): 自然科學(xué)版, 2014, 30(3): 137- 146.

    [31] Walawender J P, Szymanowski M, Hajto M J, Bokwa A. Land surface temperature patterns in the urban agglomeration of Krakow (Poland) derived from Landsat- 7/ETM+Data. Pure and Applied Geophysics, 2014, 171(6): 913- 940.

    [32] 徐涵秋, 陳本清. 不同時(shí)相的遙感熱紅外圖像在研究城市熱島變化中的處理方法. 遙感技術(shù)與應(yīng)用, 2003, 18(3): 129- 133.

    [33] 潘衛(wèi)華, 張春桂. 泉州市城市化進(jìn)程中的熱島效應(yīng)遙感研究. 國土資源遙感, 2006, (4): 50- 54.

    [34] Imhoff M L, Zhang P, Wolfe R E, Bounoua L. Remote sensing of the urban heat island effect across biomes in the continental USA. Remote Sensing of Environment, 2010, 114(3): 504- 513.

    [35] Schwarz N, Lautenbach S, Seppelt R. Exploring indicators for quantifying surface urban heat islands of European cities with MODIS land surface temperatures. Remote Sensing of Environment, 2011, 115(12): 3175- 3186.

    [36] 葉彩華, 劉勇洪, 劉偉東, 劉誠, 權(quán)維俊. 城市地表熱環(huán)境遙感監(jiān)測指標(biāo)研究及應(yīng)用. 氣象科技, 2011, 39(1): 95- 101.

    [37] Wang J, Huang B, Fu D, Atkinson P M. Spatiotemporal variation in surface urban heat Island intensity and associated determinants across major Chinese cities. Remote Sensing, 2015, 7(4): 3670- 3689.

    [38] Van Hove L W A, Jacobs C M J, Heusinkveld B G, Elbers J A, Van Driel B L, Holtslag A A M. Temporal and spatial variability of urban heat island and thermal comfort within the Rotterdam agglomeration. Building and Environment, 2015, 83: 91- 103.

    [39] Cheval S, Dumitrescu A, Bell A. The urban heat island of Bucharest during the extreme high temperatures of July 2007. Theoretical and Applied Climatology, 2009, 97(3/4): 391- 401.

    [40] Morris C J G, Simmonds I, Plummer N. Quantification of the influences of wind and cloud on the nocturnal urban heat island of a large city. Journal of Applied Meteorology, 2001, 40(2): 169- 182.

    [41] Hu L Q, Brunsell N A. The impact of temporal aggregation of land surface temperature data for surface urban heat island (SUHI) monitoring. Remote Sensing of Environment, 2013, 134: 162- 174.

    [42] Cheval S, Dumitrescu A. The summer surface urban heat island of Bucharest (Romania) retrieved from MODIS images. Theoretical and Applied Climatology, 2015, 121(3/4): 631- 640.

    [43] Mohsin T, Gough W A. Characterization and estimation of urban heat island at Toronto: impact of the choice of rural sites. Theoretical and Applied Climatology, 2012, 108(1/2): 105- 117.

    [44] Bokwa A, Hajto M J, Walawender J P, Szymanowski M. Influence of diversified relief on the urban heat island in the city of Kraków, Poland. Theoretical and Applied Climatology, 2015, 122(1/2): 365- 382.

    [45] Ren Z B, Zheng H F, He X Y, Zhang D, Yu X Y. Estimation of the relationship between urban vegetation configuration and land surface temperature with remote sensing. Journal of the Indian Society of Remote Sensing, 2015, 43(1): 89- 100.

    [46] Connors J P, Galletti C S, Chow W T L. Landscape configuration and urban heat island effects: assessing the relationship between landscape characteristics and land surface temperature in Phoenix, Arizona. Landscape Ecology, 2013, 28(2): 271- 283.

    [47] Yamashita S, Sekine K, Shoda M, Yamashita K, Hara Y. On relationships between heat island and sky view factor in the cities of Tama River basin, Japan. Atmospheric Environment, 1986, 20(4): 681- 686.

    [48] Roth M. Review of urban climate research in (sub)tropical regions. International Journal of Climatology, 2007, 27(14): 1859- 1873.

    [49] Oke T R. City size and the urban heat island. Atmospheric Environment, 1973, 7(8): 769- 779.

    [50] 崔耀平, 劉紀(jì)遠(yuǎn), 張學(xué)珍, 秦耀辰, 董金瑋. 京津唐城市群土地利用變化的區(qū)域增溫效應(yīng)模擬. 生態(tài)學(xué)報(bào), 2015, 35(4): 993- 1003.

    [51] 周莉, 江志紅, 李肇新, 楊修群. 中國東部不同區(qū)域城市群下墊面變化氣候效應(yīng)的模擬研究. 大氣科學(xué), 2015, 39(3): 596- 610.

    [52] Wan Z M. New refinements and validation of the MODIS land-surface temperature/emissivity products. Remote Sensing of Environment, 2008, 112(1): 59- 74.

    [53] Rajasekar U, Weng Q H. Urban heat island monitoring and analysis using a non-parametric model: A case study of Indianapolis. International Journal of Photogrammetry and Remote Sensing, 2009, 64(1): 86- 96.

    [54] Prata A J. Land surface temperatures derived from the Advanced Very High Resolution Radiometer and the along-track scanning radiometer: 1. Theory. Journal of Geophysical Research, 1993, 98(D9): 16689- 16702.

    [55] 劉勇洪, 軒春怡, 權(quán)維俊. 基于衛(wèi)星資料的北京陸表水體的熱環(huán)境效應(yīng)分析. 湖泊科學(xué), 2013, 25(1):73- 81.

    [56] 劉勇洪, 徐永明, 馬京津, 權(quán)維俊. 北京城市熱島的定量監(jiān)測及規(guī)劃模擬研究. 生態(tài)環(huán)境學(xué)報(bào), 2014, 23(7): 1156- 1163.

    [57] 權(quán)維俊. 長序列AVHRR的地表參數(shù)反演及其在青藏高原地區(qū)時(shí)空變化的分析[D]. 北京: 中國科學(xué)院大學(xué), 2014.

    [58] Quan W J, Chen H B, Han X Z, Liu Y H, Ye C H. A modified Becker′s split-window approach for retrieving land surface Temperature from AVHRR and VIRR. Acta Meteorologica Sinica, 2012, 26(2): 229- 240.

    [59] 卓莉, 李強(qiáng), 史培軍, 陳晉, 鄭璟, 黎夏. 基于夜間燈光數(shù)據(jù)的中國城市用地?cái)U(kuò)展類型. 地理學(xué)報(bào), 2006, 61(2): 169- 178.

    [60] 楊眉, 王世新, 周藝, 王麗濤. DMSP/OLS夜間燈光數(shù)據(jù)應(yīng)用研究綜述. 遙感技術(shù)與應(yīng)用, 2011, 26(1): 45- 51.

    [61] Liu Q P, Yang Y C, Tian H Z, Zhang B, Gu L. Assessment of human impacts on vegetation in built-up areas in China based on AVHRR, MODIS and DMSP_OLS nighttime light data, 1992- 2010. Chinese Geographical Science, 2014, 24(2): 231- 244.

    [62] 楊洋, 何春陽, 趙媛媛, 李通, 喬云偉. 利用DMSP/OLS穩(wěn)定夜間燈光數(shù)據(jù)提取城鎮(zhèn)用地信息的分層閾值法研究. 中國圖象圖形學(xué)報(bào), 2011, 16(4): 666- 673.

    [63] 王翠平, 王豪偉, 李春明, 董仁才. 基于DMSP/OLS影像的我國主要城市群空間擴(kuò)張?zhí)卣鞣治? 生態(tài)學(xué)報(bào), 2012, 32(3): 942- 954.

    [64] 劉勇洪, 葉彩華, 權(quán)維俊, 高燕虎.華北植被的光能利用率及其時(shí)空變化分析. 高原氣象, 2010, 29(3): 696- 702.

    [65] 謝莊, 崔繼良, 陳大剛, 胡保昆. 北京城市熱島效應(yīng)的晝夜變化特征分析. 氣候與環(huán)境研究, 2006, 11(1): 69- 75.

    [66] 鄭祚芳, 劉偉東, 王迎春. 北京地區(qū)城市熱島的時(shí)空分布特征. 南京氣象學(xué)院學(xué)報(bào), 2006, 29(5): 694- 699.

    [67] 黃利萍, 苗峻峰, 劉月琨. 天津城市熱島效應(yīng)的時(shí)空變化特征. 大氣科學(xué)學(xué)報(bào), 2012, 35(5): 620- 632.

    [68] 卞韜, 任國玉, 張翠華, 閻訪. 石家莊氣象站記錄的城市熱島效應(yīng)及其趨勢變化.南京信息工程大學(xué)學(xué)報(bào):自然科學(xué)版, 2012, 4(5): 402- 408.

    [69] 北京市氣象局氣候資料室. 北京城市氣候. 北京: 氣象出版社, 1992: 2- 7.

    [70] Wang W H, Liang S L, Meyers T. Validating MODIS land surface temperature products using long-term nighttime ground measurements. Remote Sensing of Environment, 2008, 112(3): 623- 635.

    [71] 楊英寶, 蘇偉忠, 江南. 南京市熱島效應(yīng)時(shí)空特征的遙感分析. 遙感技術(shù)與應(yīng)用, 2006, 21(6): 488- 492.

    [72] 李曉敏, 曾勝蘭. 成都、重慶城市熱島效應(yīng)特征對比. 氣象科技, 2015, 43(5): 888- 897.

    ResearchonquantitativeevaluationsofheatislandsfortheBeijing-Tianjin-HebeiUrbanAgglomeration

    LIU Yonghong1,2,*,FANG Xiaoyi1, ZHANG Shuo1, LUAN Qingzu1, QUAN Weijun3

    1BeijingMunicipalClimateCenter,Beijing100089,China2InstituteofBeijingUrbanMeteorology,ChinaMeteorologicalAdiministration,Beijing100089,China3EnvironmentalMeteorologyForecastCenterofBeijing-Tianjin-Hebei,ChinaMeteorologicalAdiministration,Beijing100089,China

    高分氣象遙感應(yīng)用示范系統(tǒng)(一期);中國氣象局氣候變化專項(xiàng)(CCSF201618);北京市科技計(jì)劃課題(Z161100001216011);北京市氣象局城市氣候評估創(chuàng)新團(tuán)隊(duì);中國氣象局氣候變化創(chuàng)新團(tuán)隊(duì)

    2016- 06- 03; < class="emphasis_bold">網(wǎng)絡(luò)出版日期

    日期:2017- 04- 24

    10.5846/stxb201606031074

    *通訊作者Corresponding author.E-mail: lyh7414@163.com

    劉勇洪,房小怡,張碩,欒慶祖,權(quán)維俊.京津冀城市群熱島定量評估.生態(tài)學(xué)報(bào),2017,37(17):5818- 5835.

    Liu Y H, Fang X Y, Zhang S, Luan Q Z, Quan W J.Research on quantitative evaluations of heat islands for the Beijing-Tianjin-Hebei Urban Agglomeration.Acta Ecologica Sinica,2017,37(17):5818- 5835.

    深夜精品福利| 99久久精品国产亚洲精品| 欧美日韩乱码在线| 欧美在线一区亚洲| 欧美丝袜亚洲另类 | 黄色毛片三级朝国网站| 丝袜在线中文字幕| 国产精品二区激情视频| 亚洲人成电影观看| 欧美日韩亚洲国产一区二区在线观看| 咕卡用的链子| 欧美日本亚洲视频在线播放| 女性生殖器流出的白浆| 欧美日韩福利视频一区二区| xxxhd国产人妻xxx| 欧美亚洲日本最大视频资源| 成年人免费黄色播放视频| a在线观看视频网站| 午夜久久久在线观看| 日日爽夜夜爽网站| 99久久久亚洲精品蜜臀av| 亚洲男人天堂网一区| 久久精品成人免费网站| 国内久久婷婷六月综合欲色啪| 又黄又爽又免费观看的视频| 99re在线观看精品视频| 欧美在线一区亚洲| 多毛熟女@视频| 国产精品国产高清国产av| 无遮挡黄片免费观看| 日本精品一区二区三区蜜桃| 久久九九热精品免费| 高清av免费在线| 韩国av一区二区三区四区| 一级黄色大片毛片| 欧美激情久久久久久爽电影 | 久久热在线av| 黄色怎么调成土黄色| 热re99久久国产66热| 在线国产一区二区在线| 香蕉久久夜色| 老司机在亚洲福利影院| 一区二区三区激情视频| 三上悠亚av全集在线观看| 女人高潮潮喷娇喘18禁视频| 久久国产乱子伦精品免费另类| 亚洲成人精品中文字幕电影 | 日本撒尿小便嘘嘘汇集6| 一级毛片女人18水好多| 女人被躁到高潮嗷嗷叫费观| 精品乱码久久久久久99久播| 国产精华一区二区三区| 精品久久久久久久久久免费视频 | 欧美成人性av电影在线观看| 国产亚洲欧美98| 亚洲av成人一区二区三| 国产av在哪里看| 成在线人永久免费视频| 首页视频小说图片口味搜索| 日韩免费高清中文字幕av| 亚洲成人精品中文字幕电影 | 成人精品一区二区免费| 久久久久亚洲av毛片大全| 精品久久久久久久毛片微露脸| 久99久视频精品免费| 免费观看精品视频网站| 97碰自拍视频| 日日干狠狠操夜夜爽| 老司机福利观看| 精品卡一卡二卡四卡免费| 乱人伦中国视频| 日韩三级视频一区二区三区| 啪啪无遮挡十八禁网站| 久久精品成人免费网站| 国产成人精品在线电影| 88av欧美| 久久久精品国产亚洲av高清涩受| 超碰成人久久| 欧美午夜高清在线| 在线观看免费日韩欧美大片| 色综合婷婷激情| 波多野结衣高清无吗| 午夜福利免费观看在线| 国产精品亚洲av一区麻豆| 欧美国产精品va在线观看不卡| 国产又爽黄色视频| 视频在线观看一区二区三区| 久久 成人 亚洲| 波多野结衣高清无吗| 香蕉丝袜av| 中出人妻视频一区二区| 变态另类成人亚洲欧美熟女 | 久久久久久久午夜电影 | 色在线成人网| 18禁裸乳无遮挡免费网站照片 | 久久香蕉激情| 久久精品亚洲精品国产色婷小说| avwww免费| 国产精品 国内视频| 夫妻午夜视频| 精品国产超薄肉色丝袜足j| 国产欧美日韩精品亚洲av| 亚洲成a人片在线一区二区| 亚洲自偷自拍图片 自拍| 精品国产一区二区三区四区第35| 一边摸一边抽搐一进一小说| 99精品在免费线老司机午夜| 婷婷精品国产亚洲av在线| 欧美激情高清一区二区三区| 级片在线观看| 91字幕亚洲| 99精国产麻豆久久婷婷| 狂野欧美激情性xxxx| 亚洲国产精品合色在线| 宅男免费午夜| 99久久国产精品久久久| 不卡av一区二区三区| av在线播放免费不卡| 国产黄色免费在线视频| 99精品在免费线老司机午夜| 国产成人系列免费观看| 在线观看午夜福利视频| 在线观看免费高清a一片| 精品人妻1区二区| 最近最新中文字幕大全免费视频| 人妻丰满熟妇av一区二区三区| 国产精品综合久久久久久久免费 | 免费搜索国产男女视频| 亚洲欧美日韩高清在线视频| 国产亚洲精品第一综合不卡| www.www免费av| 亚洲色图综合在线观看| 亚洲自拍偷在线| 99国产精品一区二区蜜桃av| 看免费av毛片| 国产不卡一卡二| 久久久久九九精品影院| 日日夜夜操网爽| 国产一区二区三区视频了| 午夜福利影视在线免费观看| 搡老乐熟女国产| 高潮久久久久久久久久久不卡| 在线观看66精品国产| 女人高潮潮喷娇喘18禁视频| 精品一品国产午夜福利视频| 91成年电影在线观看| 成人黄色视频免费在线看| www.精华液| 精品卡一卡二卡四卡免费| 色哟哟哟哟哟哟| 91麻豆精品激情在线观看国产 | 天堂影院成人在线观看| 午夜福利影视在线免费观看| 国产精品综合久久久久久久免费 | 国产男靠女视频免费网站| 日韩三级视频一区二区三区| 成人三级做爰电影| 999精品在线视频| 成在线人永久免费视频| 欧美大码av| 亚洲全国av大片| 亚洲av熟女| 久久国产精品人妻蜜桃| 国产男靠女视频免费网站| 电影成人av| 亚洲在线自拍视频| 无遮挡黄片免费观看| 18美女黄网站色大片免费观看| 99精品在免费线老司机午夜| 国产蜜桃级精品一区二区三区| 久久午夜亚洲精品久久| 无限看片的www在线观看| 亚洲视频免费观看视频| 久久久久久亚洲精品国产蜜桃av| 不卡一级毛片| cao死你这个sao货| 91精品国产国语对白视频| 丁香欧美五月| 亚洲精华国产精华精| 国产99久久九九免费精品| 无遮挡黄片免费观看| a在线观看视频网站| 久久久久久久久中文| 男人舔女人下体高潮全视频| 精品国产一区二区久久| 激情在线观看视频在线高清| 午夜福利在线免费观看网站| 亚洲在线自拍视频| 国产又色又爽无遮挡免费看| 精品国产一区二区久久| 亚洲熟女毛片儿| 国产免费现黄频在线看| 亚洲欧洲精品一区二区精品久久久| 亚洲成人久久性| 久久精品国产亚洲av高清一级| 国产精品永久免费网站| 女人高潮潮喷娇喘18禁视频| a在线观看视频网站| 久久精品国产综合久久久| 欧美午夜高清在线| 亚洲免费av在线视频| 国产精品免费视频内射| 91av网站免费观看| 国产精品电影一区二区三区| 久久天堂一区二区三区四区| 亚洲一区高清亚洲精品| 亚洲三区欧美一区| 亚洲精品中文字幕一二三四区| 国产人伦9x9x在线观看| 91成人精品电影| 嫁个100分男人电影在线观看| 久久久久久久久中文| 日韩欧美一区二区三区在线观看| 两人在一起打扑克的视频| 亚洲专区国产一区二区| 日韩欧美免费精品| 国产99久久九九免费精品| av在线播放免费不卡| 日本一区二区免费在线视频| 精品一区二区三区视频在线观看免费 | 国产一区二区三区综合在线观看| √禁漫天堂资源中文www| 999久久久国产精品视频| 久久久久久久精品吃奶| 色播在线永久视频| 长腿黑丝高跟| 正在播放国产对白刺激| 欧美在线黄色| av片东京热男人的天堂| 亚洲 欧美 日韩 在线 免费| 欧美在线一区亚洲| 一本综合久久免费| 亚洲精品国产精品久久久不卡| 久久久国产一区二区| 超碰97精品在线观看| 满18在线观看网站| 色综合婷婷激情| 黄片小视频在线播放| 国产区一区二久久| 精品电影一区二区在线| 久久久久久久久久久久大奶| 亚洲va日本ⅴa欧美va伊人久久| 男女床上黄色一级片免费看| 日本三级黄在线观看| 99久久人妻综合| 一级,二级,三级黄色视频| 看黄色毛片网站| 欧美午夜高清在线| 国产激情久久老熟女| 国产高清国产精品国产三级| 免费看a级黄色片| 99riav亚洲国产免费| 午夜福利影视在线免费观看| 国产成人精品在线电影| 欧美最黄视频在线播放免费 | 久久草成人影院| 久久婷婷成人综合色麻豆| 久久中文字幕人妻熟女| 中文字幕av电影在线播放| 国产精品久久久人人做人人爽| 极品人妻少妇av视频| 亚洲欧美激情在线| 亚洲一码二码三码区别大吗| 一本综合久久免费| 亚洲人成77777在线视频| 天堂影院成人在线观看| 国产亚洲精品第一综合不卡| 无人区码免费观看不卡| 每晚都被弄得嗷嗷叫到高潮| 青草久久国产| 精品少妇一区二区三区视频日本电影| 大型av网站在线播放| 91精品国产国语对白视频| 国产精品免费一区二区三区在线| 亚洲男人的天堂狠狠| 国产精品影院久久| 热99re8久久精品国产| 成人特级黄色片久久久久久久| 99国产精品一区二区蜜桃av| 欧美国产精品va在线观看不卡| 淫秽高清视频在线观看| 中文字幕av电影在线播放| 18禁国产床啪视频网站| 国内久久婷婷六月综合欲色啪| 天天躁狠狠躁夜夜躁狠狠躁| 精品久久久久久电影网| 国产欧美日韩一区二区三区在线| 日本免费a在线| 亚洲精品国产一区二区精华液| 免费日韩欧美在线观看| 欧美成人午夜精品| 他把我摸到了高潮在线观看| 99在线视频只有这里精品首页| 久久国产亚洲av麻豆专区| 99国产极品粉嫩在线观看| 天堂中文最新版在线下载| 国产激情久久老熟女| 亚洲成人免费电影在线观看| 五月开心婷婷网| 成人国语在线视频| 91精品国产国语对白视频| 亚洲三区欧美一区| 日本黄色日本黄色录像| 亚洲精华国产精华精| 最近最新中文字幕大全电影3 | 欧美乱妇无乱码| 香蕉久久夜色| 国产精品av久久久久免费| 亚洲 欧美一区二区三区| 国产91精品成人一区二区三区| 亚洲avbb在线观看| 午夜福利欧美成人| 国产亚洲欧美精品永久| 一a级毛片在线观看| 最新在线观看一区二区三区| www.自偷自拍.com| 国产极品粉嫩免费观看在线| 少妇粗大呻吟视频| 国产熟女午夜一区二区三区| 超碰97精品在线观看| 日韩大码丰满熟妇| 亚洲熟女毛片儿| 欧美在线一区亚洲| 国产精品乱码一区二三区的特点 | 中文字幕色久视频| 巨乳人妻的诱惑在线观看| 精品一区二区三区四区五区乱码| 12—13女人毛片做爰片一| 一区二区三区精品91| 欧美精品亚洲一区二区| 天堂中文最新版在线下载| 国产伦一二天堂av在线观看| 日本免费a在线| 咕卡用的链子| 欧美日韩乱码在线| 岛国在线观看网站| 亚洲国产欧美一区二区综合| 少妇 在线观看| 久久亚洲真实| 欧美精品亚洲一区二区| 中国美女看黄片| 精品福利观看| 国产不卡一卡二| 欧美日韩福利视频一区二区| 欧美日韩亚洲高清精品| av福利片在线| 亚洲全国av大片| 日韩有码中文字幕| 伊人久久大香线蕉亚洲五| 男人操女人黄网站| 亚洲国产欧美日韩在线播放| 精品熟女少妇八av免费久了| 麻豆久久精品国产亚洲av | 国产高清视频在线播放一区| 别揉我奶头~嗯~啊~动态视频| 欧美亚洲日本最大视频资源| 操出白浆在线播放| 日韩精品青青久久久久久| 操出白浆在线播放| 人人妻,人人澡人人爽秒播| 很黄的视频免费| 天天躁夜夜躁狠狠躁躁| 久久久精品国产亚洲av高清涩受| 亚洲精品在线观看二区| 亚洲一区高清亚洲精品| 99re在线观看精品视频| 国产精品永久免费网站| 免费在线观看完整版高清| 最新美女视频免费是黄的| 久久中文看片网| 黑人欧美特级aaaaaa片| www.999成人在线观看| xxxhd国产人妻xxx| 90打野战视频偷拍视频| 91九色精品人成在线观看| 俄罗斯特黄特色一大片| 午夜福利影视在线免费观看| 91av网站免费观看| 欧美不卡视频在线免费观看 | ponron亚洲| 欧美色视频一区免费| 每晚都被弄得嗷嗷叫到高潮| 日本wwww免费看| 国产深夜福利视频在线观看| 亚洲成国产人片在线观看| 成人三级黄色视频| 久久久久国产一级毛片高清牌| 国产又爽黄色视频| 高清在线国产一区| 国产精品一区二区三区四区久久 | 成年人黄色毛片网站| 999精品在线视频| 亚洲人成电影免费在线| 免费看a级黄色片| 久久久国产成人精品二区 | 欧美久久黑人一区二区| 啦啦啦 在线观看视频| 可以在线观看毛片的网站| 亚洲av日韩精品久久久久久密| 色播在线永久视频| 亚洲七黄色美女视频| 手机成人av网站| 中亚洲国语对白在线视频| 一二三四在线观看免费中文在| 夜夜爽天天搞| 国产亚洲av高清不卡| 久久九九热精品免费| 成人av一区二区三区在线看| 色婷婷av一区二区三区视频| 美女午夜性视频免费| 成人三级做爰电影| 18禁国产床啪视频网站| 亚洲精品av麻豆狂野| 国产黄a三级三级三级人| 国产精品野战在线观看 | 精品福利观看| 91国产中文字幕| 免费在线观看影片大全网站| 99久久精品国产亚洲精品| 成人三级黄色视频| 国产片内射在线| 男女床上黄色一级片免费看| 免费高清视频大片| 免费在线观看日本一区| 黄色视频,在线免费观看| 99国产精品一区二区三区| 免费搜索国产男女视频| 999久久久精品免费观看国产| 日韩有码中文字幕| 天天躁狠狠躁夜夜躁狠狠躁| 高清在线国产一区| 黄色视频不卡| 高清欧美精品videossex| 亚洲aⅴ乱码一区二区在线播放 | 亚洲免费av在线视频| 美女午夜性视频免费| 国产片内射在线| 一二三四社区在线视频社区8| 一区二区三区国产精品乱码| 中文字幕另类日韩欧美亚洲嫩草| 国产精品免费视频内射| 亚洲一卡2卡3卡4卡5卡精品中文| 国产男靠女视频免费网站| 成人三级做爰电影| 亚洲男人天堂网一区| 久久 成人 亚洲| 国产有黄有色有爽视频| 97人妻天天添夜夜摸| 青草久久国产| 婷婷精品国产亚洲av在线| 精品午夜福利视频在线观看一区| 宅男免费午夜| 久久久久久久久中文| 亚洲精品中文字幕一二三四区| 亚洲专区国产一区二区| 亚洲人成电影免费在线| 天天影视国产精品| 国产人伦9x9x在线观看| 视频区欧美日本亚洲| 两性午夜刺激爽爽歪歪视频在线观看 | xxx96com| 亚洲视频免费观看视频| 人妻久久中文字幕网| 国产精品综合久久久久久久免费 | 欧美日韩视频精品一区| 欧美最黄视频在线播放免费 | 欧美人与性动交α欧美精品济南到| 一二三四在线观看免费中文在| 国产精品 国内视频| 亚洲一区二区三区色噜噜 | 中文字幕高清在线视频| 免费高清在线观看日韩| 曰老女人黄片| 亚洲国产欧美一区二区综合| 欧美精品亚洲一区二区| 亚洲人成电影观看| 亚洲精品在线观看二区| 久久伊人香网站| 日韩欧美一区二区三区在线观看| 久久 成人 亚洲| 精品免费久久久久久久清纯| 国产精品亚洲一级av第二区| 久久久精品国产亚洲av高清涩受| 精品一品国产午夜福利视频| 99精品在免费线老司机午夜| 亚洲欧美一区二区三区黑人| 黄网站色视频无遮挡免费观看| 久99久视频精品免费| 美女大奶头视频| 国产人伦9x9x在线观看| 午夜福利在线免费观看网站| 一二三四在线观看免费中文在| 精品熟女少妇八av免费久了| 日韩 欧美 亚洲 中文字幕| 国产在线精品亚洲第一网站| 午夜福利在线免费观看网站| 九色亚洲精品在线播放| 女警被强在线播放| av欧美777| 精品一区二区三区av网在线观看| 黑丝袜美女国产一区| 两性午夜刺激爽爽歪歪视频在线观看 | 精品人妻在线不人妻| 夜夜夜夜夜久久久久| 欧美+亚洲+日韩+国产| 高清毛片免费观看视频网站 | 免费高清视频大片| 亚洲精品久久午夜乱码| 免费高清视频大片| 免费一级毛片在线播放高清视频 | 久久久久国内视频| 国产视频一区二区在线看| 18美女黄网站色大片免费观看| 黑人巨大精品欧美一区二区mp4| 免费看十八禁软件| 国产精品亚洲一级av第二区| 这个男人来自地球电影免费观看| 在线观看免费视频日本深夜| 午夜免费激情av| 黄网站色视频无遮挡免费观看| 久久精品亚洲av国产电影网| 在线观看www视频免费| 欧美色视频一区免费| 女人被躁到高潮嗷嗷叫费观| 女警被强在线播放| 神马国产精品三级电影在线观看 | 亚洲精品在线观看二区| 国产高清videossex| 午夜福利免费观看在线| 一级黄色大片毛片| 97超级碰碰碰精品色视频在线观看| 男女做爰动态图高潮gif福利片 | 黄色成人免费大全| 18禁黄网站禁片午夜丰满| 欧美黑人精品巨大| 亚洲国产欧美网| 又大又爽又粗| 精品人妻1区二区| 国产aⅴ精品一区二区三区波| 国产精品爽爽va在线观看网站 | 伦理电影免费视频| 极品教师在线免费播放| 激情视频va一区二区三区| 国产精品亚洲av一区麻豆| 欧美黑人欧美精品刺激| 无限看片的www在线观看| 香蕉国产在线看| 国产又爽黄色视频| 亚洲精品av麻豆狂野| 母亲3免费完整高清在线观看| 男人舔女人的私密视频| 法律面前人人平等表现在哪些方面| 成人免费观看视频高清| 亚洲五月色婷婷综合| 夜夜夜夜夜久久久久| 人人妻人人爽人人添夜夜欢视频| 热99国产精品久久久久久7| 麻豆国产av国片精品| 国产精品香港三级国产av潘金莲| 成人精品一区二区免费| 18禁裸乳无遮挡免费网站照片 | 欧美日韩瑟瑟在线播放| 操出白浆在线播放| 国产精品国产av在线观看| 亚洲午夜理论影院| a在线观看视频网站| 亚洲欧美日韩无卡精品| 亚洲 欧美一区二区三区| 一个人观看的视频www高清免费观看 | 亚洲专区国产一区二区| 两人在一起打扑克的视频| 国产成年人精品一区二区 | 成人黄色视频免费在线看| 男人舔女人下体高潮全视频| 看黄色毛片网站| 桃红色精品国产亚洲av| 日本欧美视频一区| 国产精品电影一区二区三区| 国产精品久久视频播放| 日本五十路高清| 国产成人av激情在线播放| 亚洲精品中文字幕一二三四区| 久久久久久人人人人人| 成人亚洲精品av一区二区 | 黄片播放在线免费| 久久人妻福利社区极品人妻图片| 在线观看66精品国产| av福利片在线| 一级a爱片免费观看的视频| 我的亚洲天堂| 欧美日韩亚洲综合一区二区三区_| 精品少妇一区二区三区视频日本电影| 国产麻豆69| 在线免费观看的www视频| 精品高清国产在线一区| 国产成人精品久久二区二区91| 成人特级黄色片久久久久久久| 女人精品久久久久毛片| 欧美激情 高清一区二区三区| 色播在线永久视频| 在线观看午夜福利视频| 国产av在哪里看| 久久国产精品人妻蜜桃| 香蕉丝袜av| 999久久久精品免费观看国产| 999久久久国产精品视频| 法律面前人人平等表现在哪些方面| 琪琪午夜伦伦电影理论片6080| 亚洲欧美一区二区三区久久| 午夜免费激情av| 乱人伦中国视频| 制服人妻中文乱码| 99国产综合亚洲精品| 日本黄色视频三级网站网址| 久久久国产成人精品二区 |