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

    基于MaxEnt模型的白鶴越冬棲息地潛在適生區(qū)分布研究

    2024-10-09 00:00:00劉小燕張增信李軍陳娟華軍彭冶晏欣邱健

    摘要:【目的】近年來白鶴(Grus leucogeranus)的適生區(qū)分布因氣候變化和人類活動等因素受到了極大影響,故運用最大熵模型(MaxEnt)模擬白鶴在中國不同時期的適生區(qū)分布,以預測不同氣候變化情景下白鶴越冬棲息地的潛在適生區(qū)變化,為瀕臨滅絕的水鳥保護提供數(shù)據(jù)支撐和理論依據(jù)?!痉椒ā扛鶕?jù)白鶴在中國越冬棲息地適生區(qū)的分布變化,基于白鶴在中國區(qū)域的分布點位、植被、地形、實測氣溫、降水以及全球氣候模式(GCMs)等數(shù)據(jù),分析其適宜生境與各環(huán)境因子的關系,利用MaxEnt開展白鶴越冬棲息地潛在適生區(qū)預測研究。【結果】①MaxEnt模型在研究白鶴越冬棲息地適生區(qū)分布和氣候適宜性方面具有較高的適用性,其中海拔、最濕月降水量、坡度和歸一化植被指數(shù)(NDVI)是影響白鶴適生區(qū)分布的主要環(huán)境因素。②白鶴的適生區(qū)主要分布在長江中下游地區(qū),其中鄱陽湖流域是最主要的高適生區(qū),而中、低適生區(qū)也主要圍繞鄱陽湖周邊分布。此外,長江中下游干流區(qū)、渤海灣沿岸、海河、遼河以及嫩江和松花江流域等地也有低適生區(qū)分布。③全球氣候模式預測結果顯示,在未來20年,白鶴越冬區(qū)低適生區(qū)面積呈上升趨勢,但中、高適生區(qū)面積普遍下降。其中,中適生區(qū)面積將收縮2 500~25 700 km2,高適生區(qū)面積將收縮3 800~12 200 km2?!窘Y論】不同溫室氣體排放情景對珍稀物種白鶴分布產(chǎn)生的影響有所差異。在全球氣候變暖背景下,白鶴越冬棲息地中、高適生區(qū)面積可能會出現(xiàn)明顯萎縮,對白鶴的保護帶來嚴峻挑戰(zhàn),應加強氣候變化背景下白鶴棲息地的研究,為白鶴保護和棲息地管理提供數(shù)據(jù)支撐。

    關鍵詞:白鶴;適生區(qū)分布;氣候情景;MaxEnt模型;棲息地;中國

    中圖分類號:S718"""""" 文獻標志碼:A開放科學(資源服務)標識碼(OSID):

    文章編號:1000-2006(2024)05-0181-08

    An prediction of" the potential distribution of suitable habitat" for Grus leucogeranus " using the MaxEnt model

    LIU Xiaoyan1, ZHANG Zengxin1*, LI Jun2, CHEN Juan2, HUA Jun2, PENG Ye3, YAN Xin1, QIU Jian1

    (1. Co-Innovation Center for Sustainable Forestry in Southern China,College of Water and Soil Conservation, College of Forestry and Grassland, Nanjing Forestry University, Nanjing 210037, China; 2. Administrative Committee of Yunlong Lake Scenic Area of Xuzhou City, Xuzhou 221018, China; 3. College of Life Sciences, Nanjing Forestry University, Nanjing 210037, China)

    Abstract: 【Objective】The middle and lower reaches of the Yangtze River basin are the most important wintering habitats for the Grus leucogeranus. However, in recent years, the distribution of suitable habitats for the G. leucogeranus in this area has been greatly altered by the combined effects of climate change and various other factors. Therefore, conducting research on the habitat suitability of G. leucogeranus in this area is of great significance for the protection and scientific management of this endangered species.【Method】In this study, based on the changes in the distribution of suitable areas of the G. leucogeranus wintering habitat in China, we used the Maximum Entropy Model (MaxEnt) to predict the potential suitable areas of the G. leucogeranus wintering habitat. The prediction was based on the G. leucogeranus distribution points, vegetation, topography, observed air temperature, precipitation, Global Climate Models" (GCMs) data, and other various environmental factors.【Result】(1) The MaxEnt model proved" highly effective in predicting the distribution of suitable habitats and assessing climate suitability for G. leucogeranus in their wintering grounds, with an AUC value of 0.978. The most influential environmental factors affecting the distribution of suitable habitat for" G. leucogeranus " were elevation, precipitation during the wettest month, slope, and NDVI. (2) The contemporary suitable habitat of the G. leucogeranus " mainly distributed in the middle and lower reaches of the Yangtze River, with the the Poyang Lake being the most important high suitable" area for the species. The" medium and low suitable" areas"" also mainly distributed around the the Poyang Lake. In addition,"" low suitable" areas were" found in the middle and lower reaches of the Yangtze River, along the coast of the Bohai Bay, and in the Haihe River, the" Liaohe River," the" Nenjiang, and" the Songhua River basins. ③The global(3) Global climate model prediction results showed" predictions indicate that"" the next 20 years, although the low suitable"" areas for" G. leucogeranus during wintering is on the rise," the medium and high suitable"" areas"" generally decreasing. Among them, compared" to the contemporary period, the medium suitable living area"" are projected to shrink by 2 500"" to 25 700 km2," and the high suitable living area will shrink by" 3" 800" to 12" 200 km2.【Conclusion】Different greenhouse gas emission scenarios will have different varying impacts on the distribution of G. leucogeranus. In the context of global warming, the wintering habitat of"" G. leucogeranus may shrink significantly, which will pose a serious challenge to the conservation of this rare species, and we should strengthen the" research on G. leucogeranus habitat in the context of climate change," to provide data support for the conservation of G. leucogeranus and the management of their habitats.

    Keywords:Grus leucogeranus; suitable habitat distribution; climate scenario; MaxEnt model; habitat; China

    白鶴(Grus leucogeranus)是國家一級重點保護野生動物,被世界自然保護聯(lián)盟(IUCN)紅色名錄列入“極危(CR)”等級。據(jù)統(tǒng)計,全球白鶴種群數(shù)量約為4 000只,它們的越冬棲息地主要分布在中國,少部分在印度[1]。白鶴在我國有東、中、西等3個遷徙種群,其中,東部地區(qū)為主的東部遷徙種群數(shù)量最多,超過99%,而西部和中部遷徙種群已接近滅絕[2]。白鶴對棲息地生境質(zhì)量要求極高,常依賴于開闊的濕地水域生存,作為濕地生態(tài)系統(tǒng)環(huán)境變化的重要指示物種之一,對整個棲息地生態(tài)系統(tǒng)及其他水鳥的保護具有重要意義[3]。棲息地的適宜性則代表了環(huán)境為野生動物提供適宜條件的能力,同時也直接反映了野生動物對棲息地的適宜性[4]。隨著氣候變化和人類活動的加劇,棲息地的喪失已成為一個全球性問題,特別是位于東亞—澳大利亞遷徙路線上的沿江沿海濕地已發(fā)生巨大變化[5],這對依賴于濕地生存的白鶴來說無疑是巨大的威脅。我國在1978—2008年濕地面積減少了約33 000" km2,與此同時,超過一半的鳥類種群數(shù)量明顯減少[6]。進入21世紀以來,我國濕地面積仍在不同程度上受到氣候變化和人類活動的影響而發(fā)生退化喪失[7],嚴重影響了遷徙鳥類的生存和繁衍,導致水鳥棲息地生境質(zhì)量下降和水鳥種群數(shù)量減少。因此,開展白鶴越冬棲息地潛在適生區(qū)分布的預測研究,是保護白鶴等珍稀瀕危物種和維護生物多樣性的基礎。

    近年來,物種分布模型(SDMs)已成為一種重要的物種適生區(qū)分布預測工具,它可以通過降雨量、相對濕度、最高溫度、最低溫度、地形和其他環(huán)境因子來預測物種的潛在分布范圍,從而建立物種發(fā)生位點與所選環(huán)境變量之間的聯(lián)系[8]。最大熵(MaxEnt)模型結合了機器學習和最大熵原理,被許多學者認為是目前預測物種潛在分布比較理想的模型之一,具有預測精度高、輸入數(shù)據(jù)規(guī)范等特點,已被廣泛用于生物地理學、入侵生物學、保護生物學和全球氣候變化對物種分布影響等研究領域[9]。如Singh[10]基于MaxEnt模型,評估了未來氣候和森林覆蓋變化對東南亞保護區(qū)受威脅鳥類的棲息地分布影響;Maleki等[11]采用MaxEnt模型和多準則評價方法評估了上游人為活動干擾加劇了氣候變化對下游水鳥棲息地的影響;白雪紅等[12]運用MaxEnt模型預測了不同情景下京津冀地區(qū)多種瀕危水鳥的潛在適宜生境分布,并對水鳥多樣性保護空缺進行了深入分析。鑒于此,本研究以白鶴的地理分布數(shù)據(jù)及當代(1970—2000年)和未來(2021—2040年)氣候數(shù)據(jù)為基礎,運用MaxEnt模型模擬白鶴在中國不同時期的適生區(qū)分布,分析當代白鶴越冬棲息地適生區(qū)分布規(guī)律,并預測不同氣候變化情景下白鶴越冬棲息地的潛在適生區(qū)變化,為瀕臨滅絕的水鳥保護提供數(shù)據(jù)支撐和理論依據(jù)。

    1 材料與方法

    1.1 數(shù)據(jù)來源與處理

    數(shù)據(jù)來源分為2部分,一部分來源于網(wǎng)站查找,另一部分來自資料收集。從全球生物多樣性信息設施(GBIF, http://www.gbif.org)、國際農(nóng)業(yè)與生物科學中心(CABl, http://www.cabi.org)和國家動物標本資源中心(http://museum.ioz.ac.cn/specimens.aspx)下載白鶴的相關位點數(shù)據(jù)(包括經(jīng)度和緯度),結合文獻資料,共收集了全球722個白鶴位點分布記錄。通過剔除缺乏詳細地理位置數(shù)據(jù)的重復位點和無效分布位點等,最終獲取了105個有效位點數(shù)據(jù)。實地調(diào)查區(qū)域位于鄱陽湖、洞庭湖、上海崇明東灘鳥類國家級自然保護區(qū)和徐州市云龍湖風景名勝區(qū)等地,獲得了209個分布位點。最終得到314個白鶴分布位點。

    1.2 環(huán)境變量獲取與選擇

    環(huán)境變量的獲取,當代(1970—2000年)的氣候數(shù)據(jù)來源于世界氣候數(shù)據(jù)庫(https://www.worldclim.org/),主要有19個用于模擬物種潛在分布最重要的生物氣候變量(bio_1—bio_19),包括氣溫和降水,空間分辨率約1 km[13]。以第2代國家(北京)氣候中心氣候系統(tǒng)模式(BCC-CSM2-MR)、加拿大氣候中心氣候系統(tǒng)模式(CanESM5)和法國氣候中心氣候系統(tǒng)模式(CNRM-CM6-1)作為未來(2021—2040年)氣候系統(tǒng)模式,分別采用3種不同溫室氣體排放濃度和社會發(fā)展水平下的SSP1-2.6、SSP2-4.5和SSP5-8.5氣候情景,用于預測未來白鶴的適生區(qū)分布。這3種氣候情景分別代表全球氣候變暖由低到高3個等級的變化趨勢,可以對未來氣候變化進行較為準確的預估[14]。其中,SSP1-2.6表示低強迫情景,代表2100年輻射強迫穩(wěn)定在約2.6 W/m2;SSP2-4.5表示中等強迫情景,代表2100年輻射強迫穩(wěn)定在約4.5 W/m2;SSP5-8.5表示高強迫情景,代表2100年輻射強迫穩(wěn)定在約8.5 W/m2。

    地形數(shù)據(jù)選取數(shù)字高程模型(DEM)來自中國科學院地理科學與資源研究所環(huán)境科學與數(shù)據(jù)中心(http://www.resdc.cn/),作為其中1個地形因子,空間分辨率為1 km×1 km,并基于海拔層在ArcGIS 10.2中提取2個地形因子(坡度和坡向)。植被數(shù)據(jù)選用從GIMMS NDVIl 3g數(shù)據(jù)集中獲取的歸一化植被指數(shù)(NDVI)變量1個。該數(shù)據(jù)集涵蓋了1982—2000年,空間分辨率為8 km×8 km,再通過最大值合成法,得到多年植被指數(shù)的平均值。最后,將23個因子(其中包括海拔、坡度、坡面、氣溫、降水量以及歸一化植被指數(shù)等)作為定義白鶴分布點的重要變量,用于構建模型。

    1.3 分布模型構建與精度評價

    使用ArcGIS 10.2軟件對所有環(huán)境變量進行坐標定義,地理坐標統(tǒng)一設置為:CCS_WCS_1984。投影坐標統(tǒng)一設置為:WGS_1984_UTM_Zone_48N。批量處理提取出中國范圍,將所有數(shù)據(jù)重采樣至8 km×8 km,最終將圖層轉(zhuǎn)換為ASCII格式。在建模過程中,將白鶴的314個分布數(shù)據(jù)樣本隨機選取75%作為訓練數(shù)據(jù)(training data),剩余的25%樣本用于測試數(shù)據(jù)(test data)。迭代運算500次,重復運行次數(shù)為10,其余參數(shù)設置默認值。為避免變量的多重共線性和MaxEnt的過度擬合[15],使用SPSS 20.0對23個環(huán)境變量進行Pearson相關性分析,考慮到兩個變量的相關系數(shù)大于0.8時生態(tài)學意義較低[16],最終選擇10個較為顯著的環(huán)境變量用于MaxEnt建模。

    MaxEnt模型利用受試者工作特征曲線(receive operating characteristic curve,ROC)分析法來檢驗模型精度。ROC曲線和橫坐標之間的面積是平均訓練AUC值,用來判斷預測結果的準確性[17]。當MaxEnt模型重復運行10次時,自動計算AUC值,AUC值范圍從0到1,值越高,模型預測效果越好[18]。通常來說,AUC值劃分為5類,分別為無效[0.5,0.6)、效果差[0.6,0.7)、效果一般[0.7,0.8)、效果良好[0.8,0.9)、效果優(yōu)秀[0.9,1.0)[19]。該模型能夠從數(shù)據(jù)中推斷或預測氣候、海拔、植被和其他環(huán)境因素與物種分布點及研究區(qū)域之間的非隨機關系,在一定的約束條件下,找到熵最大的概率區(qū)作為預測物種的最佳適生區(qū)[20]。

    1.4 潛在適生區(qū)劃分

    氣候適宜性等級劃分參考政府間氣候變化專門委員會(IPCC)第5次評估報告中關于評估可能性的劃分標準[21]和相關學者的劃分方法[22],再依據(jù)物種豐富度從多到少,使用ArcGIS 10.2軟件中的空間分析工具和自然斷點法對區(qū)域進行重分類劃分,將當代和未來白鶴適生區(qū)存在值的標準分為4類:非適生區(qū)[0,0.1)、低適生區(qū)[0.1,0.3)、中適生區(qū)[0.3,0.5)和高適生區(qū)[0.5,1.0)。繪制出不同時期白鶴在中國的適生區(qū)分布預測圖,并利用ArcGIS的空間統(tǒng)計功能,計算出不同適生區(qū)的面積以及在不同的未來氣候情景下當代與未來各適生區(qū)面積的擴張和收縮情況。

    2 結果與分析

    2.1 模型預測結果檢驗

    MaxEnt模型預測結果檢驗如圖1所示,AUC值達0.978±0.009,預測結果達到優(yōu)秀水平,表明運用MaxEnt模型能夠很好地預測白鶴在中國的適生區(qū)分布。

    2.2 環(huán)境變量重要性及響應曲線分析

    運用Jackknife測試法計算每個變量的相對重要性和貢獻率,結果如表1所示:海拔、最濕月份降水量、坡度和歸一化植被指數(shù)對模型貢獻率最大,分別占43.9%、19.0%、10.1%和9.2%,累計貢獻率達82.2%。

    由表1可知,這4個環(huán)境因子是影響白鶴分布的主要環(huán)境因子。根據(jù)環(huán)境變量與對應的s物種存在概率,可獲得相應的關鍵變量響應曲線,以存在概率gt;0.5為適宜區(qū)域。進一步分析(圖2)表明:白鶴適生區(qū)的海拔在0~500 m,最適宜的降水量為270 mm;坡度在0°左右最為適宜;歸一化植被指數(shù)在0.4~0.7,其中0.52為最佳。

    2.3 當代白鶴在中國的適生區(qū)分布

    白鶴僅分布在我國經(jīng)濟較為發(fā)達的東部地區(qū),且越冬種群過于集中,其中,主要以長江中下游地區(qū)的鄱陽湖流域為聚集區(qū)(圖3),而鄱陽湖濕地生境質(zhì)量的好壞也將直接關系到這一物種的生存。白鶴當代的適生區(qū)主要分布在長江口沿岸和以鄱陽湖為中心的區(qū)域,其中,低適生區(qū)主要分布在江蘇沿江地區(qū)、東北部渤海灣沿岸、海河、遼河以及嫩江和松花江匯合流域。江蘇沿江地區(qū)又主要以東部太湖平原為低適生區(qū)分布地,包括太湖、陽澄湖和滆湖等大型湖泊,這些停歇點濕地是遷飛白鶴重要的能量補給地。另外,天津濱海地以及嫩江和松花江匯合處也有極少量的中適生區(qū)分布。在整個中國區(qū)域白鶴的非適生區(qū)面積為9 133 000 km2,約占國土面積的95%,而低、中、高適生區(qū)面積分別為354 700、77 200和34 300 km2。由此可見,鶴的適生區(qū)范圍極其狹窄,如果再加上人為因素的影響,其適生區(qū)將更容易遭到破壞進而導致白鶴種群數(shù)量下降。

    2.4 不同氣候情景下白鶴的適生區(qū)分布范圍和空間格局的變化

    從白鶴在中國的適生區(qū)分布面積變化來看(圖4),在未來3種氣候模式情景下,中、高適生區(qū)面積均有所下降,下降規(guī)模為2 500~25 700 km2,極少數(shù)情況下中適生區(qū)面積的增長可忽略不計。從低適生區(qū)和非適生區(qū)面積的變化來看,白鶴的適生區(qū)情況則有所好轉(zhuǎn)。在BCC-CSM2-MR和CanESM5氣候模式下,白鶴的適生區(qū)面積變化差異較大,除SSPS-8.5情景外,非適生區(qū)面積減少,低適生區(qū)面積增加趨勢明顯。其中,BCC-CSM2-MR模式SSP5-8.5情景下總適生區(qū)面積減少最多,減少了78 800 km2;非適生區(qū)面積增加最多,增加了92 800 km2。CanESM5模式SSP1-2.6情景下非適生區(qū)面積減少最多,減少了98 700 km2;低適生區(qū)面積增加最多,增加了136 900 km2,且中、高適生區(qū)面積在3種情景下減少最多。CNRM-CM6-1模式各情景適生區(qū)面積變化趨勢較為平緩,中、高適生區(qū)面積均呈下降趨勢。

    與當代白鶴適生區(qū)面積相比(圖5)可知,2021—2040年, BCC-CSM2-MR模式顯示,在SSP1-2.6情景下,適生區(qū)面積減少了261 700 km2。經(jīng)分析,高適生區(qū)面積損失了9 400 km2,損失地主要集中在近鄱陽湖流域的贛江支流地區(qū),該區(qū)是白鶴越冬的主要棲息地。中適生區(qū)面積損失了30 100 km2,損失地主要集中在湖南洞庭湖流域。低適生區(qū)面積損失最多為222 200 km2,損失地主要包括重慶西南部、山東西北角、河北東南部以及吉林西北角;在SSP2-4.5情景下,適生區(qū)面積減少了288 700 km2。經(jīng)分析,高適生區(qū)面積損失最多,達11 500 km2,損失地主要集中在江西中北部鄱陽湖流域周邊。中適生區(qū)面積損失最大,達45 900 km2,主要集中湖南洞庭湖流域以南城市較為密集的地區(qū),包括岳陽、長沙、湘潭、衡陽以及永州等地。低適生區(qū)面積損失了231 300 km2,主要集中在重慶西南角、渤海灣的海河流域、遼河流域;在SSP5-8.5情景下,適生區(qū)面積減少了143 300 km2,是所有情景中最低的。經(jīng)分析,高適生區(qū)面積損失了6 200 km2,中適生區(qū)面積損失了20 700 km2,低適生區(qū)面積損失了116 400 km2,損失地主要集中在遼河流域。而CanESM5模式和CNRM-CM6-1模式也顯示未來不同情景下適生區(qū)面積雖有差異,但總體來看,未來白鶴的非適生區(qū)面積與當代白鶴適生區(qū)面積相比主要以下降為主,其中,非適生區(qū)與低適生區(qū)面積之間的轉(zhuǎn)變較為明顯,而中、高適生區(qū)面積則有所下降??傮w來說,在SSP2-4.5情景下白鶴中、高適生區(qū)下降的趨勢較為緩和,而在SSP5-8.5情景下白鶴的非適生區(qū)則有所擴張。

    白鶴在中國的適生區(qū)分布范圍區(qū)劃面積如表2所示。2021—2040年,在BCC-CSM2-MR氣候模式的3種排放情景下,白鶴在SSP2-4.5情景下總適生區(qū)面積最高,達到了578 400 km2,低、中、高適生區(qū)面積均最高;在SSP5-8.5情景下總適生區(qū)面積最低,為387 300 km2,低、中、高適生區(qū)面積均最低。在CanESM5氣候模式的3種排放情景下,白鶴在SSP1-2.6情景下總適生區(qū)面積最高,達到了578 800 km2;而高適生區(qū)面積在SSP2-4.5情景下達到最高,為29 600 km2。在CNRM-CM6-1氣候模式的3種排放情景下,白鶴在SSP1-2.6情景下總適生區(qū)面積最低,為408 100 km2,低、中、高適生區(qū)面積均最低;在SSP5-8.5情景下總適生區(qū)面積最高,為486 500 km2,除高適生區(qū)外,低、中適生區(qū)面積均最高。

    3 討 論

    本研究重點探討了氣候、地形和歸一化植被指數(shù)等非生物因素對白鶴分布的影響,其中,海拔、最濕月份降水量、坡度和歸一化植被指數(shù)是影響白鶴適生區(qū)分布的最主要影響因素,而地形因子在物種分布中發(fā)揮著不可替代的作用,濕地生境質(zhì)量的高低直接影響著水鳥對濕地棲息地的利用和行為適應[23]。吳慶明等[24]認為長江中下游地區(qū),雖然海拔高差較小,但海拔與水位分布、植被分布及白鶴隱蔽性等具有極大相關性。最濕月份降水量是白鶴適生區(qū)分布的第二大影響因子,而降水對于營造良好的棲息環(huán)境同樣至關重要。Li等[25]研究表明,雨季水位變化對水鳥適生區(qū)分布影響顯著。當汛期來臨時,洪澇災害的長期作用會引起清水下泄,加劇河床侵蝕和下游水位降低,導致鄱陽湖和洞庭湖的枯水期提前、旱季延長,導致越冬水鳥獲取食物資源的難度增加。此外,作為遷徙水鳥之一的白鶴,其集中停歇的濕地分布在我國東北地區(qū)(黑龍江流域、遼河中下游流域)和黃河濕地,而資料顯示,水壩的大量建設加速了自然濕地的喪失,使白鶴在東北和長江濕地之間的停歇區(qū)域大為減少[26]。Hou等[27]認為,自然濕地中食物資源的短缺會直接影響到水鳥種群對人工濕地的利用。特別是在2020年夏季極端枯水位的影響下,鄱陽湖濕地食物資源豐富度急劇下降,導致鶴類和雁類種群轉(zhuǎn)向湖區(qū)周邊稻田和藕塘等人工濕地覓食[28]。除海拔和降水外,坡度和歸一化植被指數(shù)也是影響白鶴適生區(qū)分布的重要因子。歸一化植被指數(shù)能夠反映區(qū)域植被資源的豐富程度,當NDVI>0時,表示地表具有一定數(shù)量的植被覆蓋,且隨NDVI值而增加[29]。植被覆蓋度一方面代表食物資源的可利用性,另一方面也反映植被覆蓋能否為水鳥提供適宜的隱蔽條件,這對白鶴的適生區(qū)分布具有重要作用。

    未來氣候情景下,多模式多情景均顯示白鶴越冬棲息地中、高適生區(qū)面積普遍呈下降趨勢,不利于白鶴的越冬棲息。其中,白鶴棲息地中、高適生區(qū)主要分布在長江中下游流域,特別是鄱陽湖、洞庭湖流域,而中、低適生區(qū)主要圍繞高適生區(qū)呈環(huán)形向外擴散。相對而言,在各模式中SSP2-4.5情景下,未來白鶴適生區(qū)面積變化趨于穩(wěn)定,并且局部地區(qū)呈增加趨勢,對適生區(qū)有增加的區(qū)域也應得到重視,適當提前做好白鶴越冬棲息地的保護工作;而在SSP5-8.5高排放情景下,中、高適生區(qū)面積下降趨勢明顯,說明不同溫室氣體排放情景對水鳥分布也將產(chǎn)生不同影響[30]。已有研究表明,高緯度地區(qū)氣溫升高比低緯度地區(qū)表現(xiàn)的更加明顯[31],白鶴夏季繁殖區(qū)也將受到氣候變化影響。

    綜上,本研究通過對白鶴在中國的適生區(qū)分布與環(huán)境因素的關系分析,以及MaxEnt生態(tài)位模型的應用,提高了對白鶴分布情況的準確性和預測能力。這對于越冬水鳥的分布預測和棲息地管理具有重要的科學價值。未來可以進一步加強氣候變化對白鶴適生區(qū)影響研究,特別是加強對中、高適生區(qū)的保護和研究,并探究白鶴等其他瀕危物種的適生區(qū)變化情況,加強環(huán)境保護措施,保護和維護生物多樣性。

    參考文獻(reference):

    [1]MEINE C,ARCHIBALD G.The cranes:status survey and conservation action plan[M].Gland,Switzerland:IUCN,1996.

    [2]KANAI Y,UETA M,GERMOGENOV N,et al.Migration routes and important resting areas of Siberian cranes (Grus leucogeranus) between northeastern Siberia and China as revealed by satellite tracking[J].Biol Conserv,2002,106(3):339-346.DOI: 10.1016/S0006-3207(01)00259-2.

    [3]曹壘,孟凡娟,趙青山.基于前沿監(jiān)測技術探討“大開發(fā)”對鳥類遷徙及其棲息地的影響[J].中國科學院院刊,2021,36(4):436-447.CAO L,MENG F J,ZHAO Q S.Understanding effects of large-scale development on bird migration and habitats through cutting edge avian monitoring techniques[J].Bull Chin Acad Sci,2021,36(4):436-447.DOI: 10.16418/j.issn.1000-3045.20210309002.

    [4]TANG X G,LI H P,XU X B,et al.Changing land use and its impact on the habitat suitability for wintering Anseriformes in Chinas Poyang Lake region[J].Sci Total Environ,2016,557/558:296-306.DOI: 10.1016/j.scitotenv.2016.03.108.

    [5]MACKINNON J,VERKUIL Y I,MURRAY N J.IUCN situation analysis on east and southeast Asian intertidal habitats, with particular reference to the Yellow Sea (including the Bohai Sea)[M]. Gland, Switzeyland: IUCN, 2012.

    [6]NIU Z G,ZHANG H Y,WANG X W,et al.Mapping wetland changes in China between 1978 and 2008[J].Chin Sci Bull,2012,57(22):2813-2823.DOI: 10.1007/s11434-012-5093-3.

    [7]耿國彪.我國濕地保護形勢不容樂觀:第二次全國濕地資源調(diào)查結果公布[J].綠色中國,2014(3):8-11.GENG G B.The situation of wetland protection in China is not optimistic[J].Green China,2014(3):8-11.

    [8]PHILLIPS S J,ANDERSON R P,SCHAPIRE R E.Maximum entropy modeling of species geographic distributions[J].Ecol Model,2006,190(3/4):231-259.DOI: 10.1016/j.ecolmodel.2005.03.026.

    [9]MEROW C,SILANDER J A.A comparison of Maxlike and Maxent for modelling species distributions[J].Methods Ecol Evol,2014,5(3):215-225.DOI: 10.1111/2041-210x.12152.

    [10]SINGH M.Evaluating the impact of future climate and forest cover change on the ability of Southeast (SE) Asias protected areas to provide coverage to the habitats of threatened avian species[J].Ecol Indic,2020,114:106307.DOI: 10.1016/j.ecolind.2020.106307.

    [11]MALEKI S,RAHDARI V,BAGHDADI N.Humans in the upstream can exacerbate climate change impacts on water birds habitat in the downstream[J].Sci Rep,2021,11(1):20203.DOI: 10.1038/s41598-021-99822-8.

    [12]白雪紅,王文杰,蔣衛(wèi)國,等.氣候變化背景下京津冀地區(qū)瀕危水鳥潛在適宜區(qū)模擬及保護空缺分析[J].環(huán)境科學研究,2019,32(6):1001-1011.BAI X H,WANG W J,JIANG W G, et al.Simulation of potential suitable distribution of endangered waterfowl and its gap analysis of national nature reserves in Beijing-Tianjin-Hebei region under climate change scenarios[J].Res Environ Sci,2019,32(6):1001-1011.DOI: 10.13198/j.issn.1001-6929.2018.08.15.

    [13]MOLLOY S W,DAVIS R A,VAN ETTEN E J B.Species distribution modelling using bioclimatic variables to determine the impacts of a changing climate on the western ringtail possum (Pseudocheirus occidentals;Pseudocheiridae)[J].Envir Conserv,2014,41(2):176-186.DOI: 10.1017/s0376892913000337.

    [14]翁宇威,蔡聞佳,王燦.共享社會經(jīng)濟路徑(SSPs)的應用與展望[J].氣候變化研究進展,2020,16(2):215-222.WENG Y W,CAI W J,WANG C.The application and future directions of the shared socioeconomic pathways(SSPs)[J].Clim Change Res,2020,16(2):215-222.DOI: 10.12006/j.issn.1673-1719.2019.078.

    [15]FOTHERINGHAM A S,OSHAN T M.Geographically weighted regression and multicollinearity:dispelling the myth[J].J Geogr Syst,2016,18(4):303-329.DOI: 10.1007/s10109-016-0239-5.

    [16]翟新宇,申宇芳,朱圣華,等.未來氣候變化對孑遺植物鵝掌楸地理分布的影響[J].熱帶亞熱帶植物學報,2021,29(2):151-161.ZHAI X Y,SHEN Y F,ZHU S H,et al.Potential impacts of climate change in future on the geographical distributions of relic Liriodendron chinense[J].J Trop Subtrop Bot,2021,29(2):151-161.DOI: 10.11926/jtsb.4322.

    [17]李麗鶴,劉會玉,林振山,等.基于MAXENT和ZONATION的加拿大一枝黃花入侵重點監(jiān)控區(qū)確定[J].生態(tài)學報,2017,37(9):3124-3132.LI L H,LIU H Y,LIN Z S,et al.Identifying priority areas for monitoring the invasion of Solidago canadensis based on MAXENT and ZONATION[J].Acta Ecol Sin,2017,37(9):3124-3132.DOI: 10.5846/stxb201601260182.

    [18]SWETS J A.Measuring the accuracy of diagnostic systems[J]. Science,1988,240(4857):1285-1293.DOI: 10.1126/science.3287615.

    [19]ZHAN P,WANG F Y,XIA P G,et al.Assessment of suitable cultivation region for Panax notoginseng under different climatic conditions using MaxEnt model and high-performance liquid chromatography in China[J].Ind Crops Prod,2022,176:114416.DOI: 10.1016/j.indcrop.2021.114416.

    [20]LI Z J,LIU Y M,ZENG H.Application of the MaxEnt model in improving the accuracy of ecological red line identification:a case study of Zhanjiang,China[J].Ecol Indic,2022,137:108767.DOI: 10.1016/j.ecolind.2022.108767.

    [21]孫穎,秦大河,劉洪濱.IPCC第五次評估報告不確定性處理方法的介紹[J].氣候變化研究進展,2012,8(2):150-153.SUN Y,QIN D H,LIU H B.Introduction to treatment of uncertainties for IPCC fifth assessment report[J].Progressus Inquisitiones DE Mutat Clim,2012,8(2):150-153.DOI: 10.3969/j.issn.1673-1719.2012.02.012.

    [22]ZHENG H F,SHEN G Q,SHANG L Y,et al.Efficacy of conservation strategies for endangered oriental white storks (Ciconia boyciana) under climate change in northeast China[J].Biol Conserv,2016,204:367-377.DOI: 10.1016/j.biocon.2016.11.004.

    [23]王新建,周立志,陳錦云,等.長江下游沿江濕地升金湖越冬水鳥覓食集團結構及生態(tài)位特征[J].湖泊科學,2021,33(2):518-528,651.WANG X J,ZHOU L Z,CHEN J Y,et al.Foraging guild structure and niche characteristics of wintering waterbirds at Lake Shengjin,lower reaches of the Yangtze River[J].J Lake Sci,2021,33(2):518-528,651.DOI: 10.18307/2021.0216.

    [24]吳慶明,王磊,朱瑞萍,等.基于MAXENT模型的丹頂鶴營巢生境適宜性分析:以扎龍保護區(qū)為例[J].生態(tài)學報,2016,36(12):3758-3764.WU Q M,WANG L,ZHU R P,et al.Nesting habitat suitability analysis of red-crowned crane in Zhalong Nature Reserve based on MAXENT modeling[J].Acta Ecol Sin,2016,36(12):3758-3764.DOI: 10.5846/stxb201410101997.

    [25]LI Y K,QIAN F W,SILBERNAGEL J,et al.Community structure,abundance variation and population trends of waterbirds in relation to water level fluctuation in Poyang Lake[J].J Gt Lakes Res,2019,45(5):976-985.DOI: 10.1016/j.jglr.2019.08.002.

    [26]馬柱國,符淙斌,周天軍,等.黃河流域氣候與水文變化的現(xiàn)狀及思考[J].中國科學院院刊,2020,35(1):52-60.MA Z G,F(xiàn)U C B,ZHOU T J,et al.Status and ponder of climate and hydrology changes in the Yellow River basin[J].Bull Chin Acad Sci,2020,35(1):52-60.DOI: 10.16418/j.issn.1000-3045.20191223002.

    [27]HOU J J,LIU Y F,F(xiàn)RASER J D,et al.Drivers of a habitat shift by critically endangered Siberian cranes:evidence from long-term data[J].Ecol Evol,2020,10(20):11055-11068.DOI: 10.1002/ece3.6720.

    [28]HOU J J,LI L,WANG Y F,et al.Influences of submerged plant collapse on diet composition,breadth,and overlap among four crane species at Poyang Lake,China[J].Front Zool,2021,18(1):24.DOI: 10.1186/s12983-021-00411-2.

    [29]高浩翔,申立泉,劉瑞,等.基于最大熵模型的野生馬麝夏季生境適宜性研究[J].生態(tài)學報,2023,43(1):441-448.GAO H X,SHEN L Q,LIU R,et al.Summer habitat suitability of wild alpine musk deer based on MaxEnt model[J].Acta Ecol Sin,2023,43(1):441-448.DOI: 10.5846/stxb202112043432.

    [30]LEHIKOINEN A,JAATINEN K,VHTALO A V,et al.Rapid climate driven shifts in wintering distributions of three common waterbird species[J].Glob Chang Biol,2013,19(7):2071-2081.DOI: 10.1111/gcb.12200.

    [31]王會軍,唐國利,陳海山,等.“一帶一路”區(qū)域氣候變化事實、影響及可能風險[J].大氣科學學報,2020,43(1):1-9.WANG H J,TANG G L,CHEN H S,et al.The Belt and Road region climate change:facts,impacts and possible risks[J].Trans Atmos Sci,2020,43(1):1-9.DOI: 10.13878/j.cnki.dqkxxb.20191110003.

    (責任編輯 王國棟)

    基金項目:國家重點研發(fā)計劃(2019YFC0409004);國家自然科學基金項目(41971025);江蘇高校優(yōu)勢學科建設工程資助項目(PAPD)。

    第一作者:劉小燕(15070701057@163.com)。

    *通信作者:張增信(nfuzhang@163.com),教授。

    引文格式:劉小燕, 張增信, 李軍,等.基于MaxEnt模型的白鶴越冬棲息地潛在適生區(qū)分布研究[J]. 南京林業(yè)大學學報(自然科學版),2024,48(5):181-188.

    LIU X Y, ZHANG Z X, LI J, et al.An prediction of the potential distribution of suitable habitat" for Grus leucogeranus " using the MaxEnt model[J]. Journal of Nanjing Forestry University (Natural Sciences Edition),2024,48(5):181-188.

    亚洲熟女精品中文字幕| 男人添女人高潮全过程视频| 女人久久www免费人成看片| 丝袜人妻中文字幕| 亚洲av电影在线观看一区二区三区| 又黄又粗又硬又大视频| 午夜免费成人在线视频| 国产精品.久久久| 如日韩欧美国产精品一区二区三区| 国产男人的电影天堂91| 王馨瑶露胸无遮挡在线观看| 久久精品国产a三级三级三级| 日本黄色日本黄色录像| 中文欧美无线码| 男女边吃奶边做爰视频| 一边摸一边抽搐一进一出视频| 97在线人人人人妻| 精品国产一区二区三区久久久樱花| 精品国产乱码久久久久久小说| 欧美久久黑人一区二区| 亚洲精品中文字幕在线视频| 十八禁网站网址无遮挡| 后天国语完整版免费观看| 成人影院久久| 国产视频一区二区在线看| 日韩欧美一区视频在线观看| 大型av网站在线播放| 狂野欧美激情性xxxx| netflix在线观看网站| 久久 成人 亚洲| 亚洲中文字幕日韩| 中文字幕人妻丝袜一区二区| 成人黄色视频免费在线看| 国产精品一区二区在线不卡| 又紧又爽又黄一区二区| 熟女av电影| 亚洲精品自拍成人| 久久午夜综合久久蜜桃| 亚洲国产欧美日韩在线播放| 亚洲成av片中文字幕在线观看| 国产精品偷伦视频观看了| 久久ye,这里只有精品| 免费日韩欧美在线观看| 丁香六月天网| 伊人亚洲综合成人网| 大片电影免费在线观看免费| 久热爱精品视频在线9| 日韩欧美一区视频在线观看| 老鸭窝网址在线观看| 成人三级做爰电影| 王馨瑶露胸无遮挡在线观看| 国产精品麻豆人妻色哟哟久久| 国产av精品麻豆| 精品国产一区二区久久| 夜夜骑夜夜射夜夜干| 国产野战对白在线观看| 99国产精品一区二区蜜桃av | 亚洲中文字幕日韩| 啦啦啦 在线观看视频| 国产精品一国产av| 啦啦啦视频在线资源免费观看| 每晚都被弄得嗷嗷叫到高潮| 免费高清在线观看日韩| 一级片'在线观看视频| 老司机亚洲免费影院| 国产精品久久久av美女十八| 纵有疾风起免费观看全集完整版| 亚洲精品国产色婷婷电影| 国产黄色视频一区二区在线观看| 蜜桃在线观看..| 亚洲人成电影观看| av在线播放精品| 晚上一个人看的免费电影| 一二三四社区在线视频社区8| 91国产中文字幕| 亚洲欧洲精品一区二区精品久久久| 九色亚洲精品在线播放| 少妇被粗大的猛进出69影院| 婷婷色麻豆天堂久久| svipshipincom国产片| 久久久久网色| 宅男免费午夜| 波多野结衣一区麻豆| 人人妻人人添人人爽欧美一区卜| 国产精品国产三级国产专区5o| 午夜福利乱码中文字幕| 亚洲av国产av综合av卡| 国产精品免费大片| 一级a爱视频在线免费观看| 在线观看免费视频网站a站| 天天操日日干夜夜撸| 亚洲国产欧美在线一区| 97精品久久久久久久久久精品| 日韩精品免费视频一区二区三区| 高清黄色对白视频在线免费看| 一本色道久久久久久精品综合| 丰满少妇做爰视频| 男女午夜视频在线观看| 在线看a的网站| 啦啦啦在线观看免费高清www| 国产精品熟女久久久久浪| 免费少妇av软件| 亚洲七黄色美女视频| 可以免费在线观看a视频的电影网站| 国产在视频线精品| 国产日韩欧美在线精品| 日韩人妻精品一区2区三区| 午夜久久久在线观看| 亚洲精品日本国产第一区| 日本五十路高清| 一级毛片黄色毛片免费观看视频| 各种免费的搞黄视频| 亚洲自偷自拍图片 自拍| 国产精品麻豆人妻色哟哟久久| 男女之事视频高清在线观看 | 妹子高潮喷水视频| 精品久久久久久久毛片微露脸 | 男女国产视频网站| 悠悠久久av| 波野结衣二区三区在线| 亚洲国产欧美在线一区| 欧美亚洲日本最大视频资源| 国产三级黄色录像| 2018国产大陆天天弄谢| 欧美亚洲日本最大视频资源| 欧美xxⅹ黑人| 校园人妻丝袜中文字幕| 女人高潮潮喷娇喘18禁视频| 麻豆av在线久日| 黄色a级毛片大全视频| 男人舔女人的私密视频| 国产成人精品久久二区二区91| svipshipincom国产片| 久久久久久免费高清国产稀缺| 19禁男女啪啪无遮挡网站| 少妇人妻久久综合中文| 欧美精品啪啪一区二区三区 | 国产精品九九99| 国产成人免费观看mmmm| 欧美成人午夜精品| 性高湖久久久久久久久免费观看| 国产精品免费大片| 久久人人97超碰香蕉20202| 中文字幕另类日韩欧美亚洲嫩草| 美女扒开内裤让男人捅视频| 热re99久久精品国产66热6| 五月开心婷婷网| 久久久久久亚洲精品国产蜜桃av| 亚洲五月色婷婷综合| 99re6热这里在线精品视频| 精品国产一区二区久久| 久久人人爽人人片av| 国产欧美亚洲国产| 欧美日韩亚洲国产一区二区在线观看 | 18禁观看日本| 欧美在线一区亚洲| 男男h啪啪无遮挡| 韩国精品一区二区三区| 在线天堂中文资源库| 操美女的视频在线观看| 亚洲国产成人一精品久久久| 成年美女黄网站色视频大全免费| 色视频在线一区二区三区| 欧美+亚洲+日韩+国产| 日本a在线网址| 国产高清videossex| 999精品在线视频| 香蕉丝袜av| 国产精品 欧美亚洲| 国产无遮挡羞羞视频在线观看| 岛国毛片在线播放| 国产精品免费大片| 亚洲精品一区蜜桃| 久久99精品国语久久久| 超碰成人久久| 捣出白浆h1v1| 国产精品三级大全| 男人舔女人的私密视频| 超色免费av| 国产av国产精品国产| 少妇精品久久久久久久| 亚洲av国产av综合av卡| 99国产精品一区二区三区| 久久av网站| 亚洲国产av影院在线观看| 美女国产高潮福利片在线看| 成人免费观看视频高清| 久久ye,这里只有精品| 亚洲av成人不卡在线观看播放网 | 狂野欧美激情性xxxx| 亚洲熟女毛片儿| videosex国产| 日韩av免费高清视频| 国产无遮挡羞羞视频在线观看| 国产熟女欧美一区二区| 亚洲精品国产色婷婷电影| 国产深夜福利视频在线观看| 亚洲人成电影观看| 国产片特级美女逼逼视频| 尾随美女入室| www.精华液| 午夜精品国产一区二区电影| 久久天躁狠狠躁夜夜2o2o | 最近最新中文字幕大全免费视频 | 中文字幕另类日韩欧美亚洲嫩草| 国产av国产精品国产| 男女无遮挡免费网站观看| 欧美成狂野欧美在线观看| 国语对白做爰xxxⅹ性视频网站| 国产精品久久久人人做人人爽| 最黄视频免费看| 免费日韩欧美在线观看| 天天添夜夜摸| 成人影院久久| av欧美777| bbb黄色大片| 午夜激情av网站| 看免费成人av毛片| 国产精品三级大全| 国产在线观看jvid| 人人妻人人澡人人看| 国产精品一区二区精品视频观看| 国产在线免费精品| 精品第一国产精品| 国产成人系列免费观看| 18禁黄网站禁片午夜丰满| 亚洲欧美一区二区三区国产| 午夜福利免费观看在线| 9色porny在线观看| 一区福利在线观看| 一本—道久久a久久精品蜜桃钙片| 在线观看www视频免费| 国产有黄有色有爽视频| 亚洲中文日韩欧美视频| 久久天堂一区二区三区四区| 免费不卡黄色视频| 少妇粗大呻吟视频| av网站免费在线观看视频| 久久99热这里只频精品6学生| 国产日韩欧美在线精品| av福利片在线| 国产无遮挡羞羞视频在线观看| 在线看a的网站| 中文字幕av电影在线播放| 午夜福利乱码中文字幕| 精品欧美一区二区三区在线| 男女午夜视频在线观看| 18禁黄网站禁片午夜丰满| 亚洲五月色婷婷综合| 欧美激情高清一区二区三区| 真人做人爱边吃奶动态| 一本一本久久a久久精品综合妖精| xxx大片免费视频| 久久久久久久大尺度免费视频| cao死你这个sao货| 人妻人人澡人人爽人人| 五月开心婷婷网| 亚洲精品国产色婷婷电影| 日本91视频免费播放| 香蕉国产在线看| a级毛片黄视频| 午夜福利免费观看在线| 高清不卡的av网站| 久久久久精品人妻al黑| 久久国产精品人妻蜜桃| 午夜福利一区二区在线看| av片东京热男人的天堂| 三上悠亚av全集在线观看| a级毛片在线看网站| 99国产精品99久久久久| www.999成人在线观看| 国产精品欧美亚洲77777| 国产亚洲av片在线观看秒播厂| 亚洲美女黄色视频免费看| a级毛片黄视频| 最近中文字幕2019免费版| 看免费av毛片| 在线精品无人区一区二区三| 丰满迷人的少妇在线观看| 美女中出高潮动态图| 一区二区av电影网| 成年av动漫网址| 久久精品亚洲av国产电影网| 亚洲欧美日韩另类电影网站| cao死你这个sao货| 王馨瑶露胸无遮挡在线观看| 成人亚洲精品一区在线观看| 啦啦啦中文免费视频观看日本| 日韩欧美一区视频在线观看| 免费久久久久久久精品成人欧美视频| 美女午夜性视频免费| 亚洲国产av新网站| 91字幕亚洲| 亚洲五月婷婷丁香| 天堂俺去俺来也www色官网| 国产成人精品久久二区二区91| 亚洲欧美一区二区三区久久| 看十八女毛片水多多多| av国产精品久久久久影院| 成人午夜精彩视频在线观看| 国产亚洲欧美精品永久| 国产又爽黄色视频| 亚洲成国产人片在线观看| 亚洲精品国产av蜜桃| 成人手机av| 欧美精品亚洲一区二区| 国产女主播在线喷水免费视频网站| 一区二区三区乱码不卡18| 看免费成人av毛片| 高清不卡的av网站| 高清视频免费观看一区二区| 黑人欧美特级aaaaaa片| 宅男免费午夜| 久久热在线av| 尾随美女入室| 黄片播放在线免费| www.999成人在线观看| 欧美大码av| 午夜福利乱码中文字幕| 多毛熟女@视频| 欧美精品高潮呻吟av久久| 视频区图区小说| 蜜桃国产av成人99| 欧美精品一区二区免费开放| 777久久人妻少妇嫩草av网站| 国产91精品成人一区二区三区 | 国产精品麻豆人妻色哟哟久久| 99re6热这里在线精品视频| av又黄又爽大尺度在线免费看| 水蜜桃什么品种好| 午夜免费男女啪啪视频观看| 大香蕉久久成人网| 中文字幕av电影在线播放| 波多野结衣av一区二区av| 人成视频在线观看免费观看| 成人影院久久| www.自偷自拍.com| 亚洲av日韩在线播放| 国产人伦9x9x在线观看| 亚洲 国产 在线| 久久九九热精品免费| 亚洲专区国产一区二区| 亚洲精品日韩在线中文字幕| 99热全是精品| 亚洲国产精品一区三区| 赤兔流量卡办理| 国产精品免费大片| 午夜免费男女啪啪视频观看| 丰满迷人的少妇在线观看| 国产又色又爽无遮挡免| 国产精品一区二区免费欧美 | 欧美亚洲 丝袜 人妻 在线| a 毛片基地| 国产熟女欧美一区二区| 欧美久久黑人一区二区| 亚洲国产精品国产精品| 久久天躁狠狠躁夜夜2o2o | 免费在线观看完整版高清| www.熟女人妻精品国产| 热99久久久久精品小说推荐| 成人影院久久| 欧美人与善性xxx| 精品一区二区三区av网在线观看 | 超碰97精品在线观看| 这个男人来自地球电影免费观看| 天堂俺去俺来也www色官网| 一区二区三区精品91| 一区二区三区激情视频| 亚洲天堂av无毛| 老司机靠b影院| 69精品国产乱码久久久| 嫩草影视91久久| 精品一区二区三区四区五区乱码 | 欧美人与性动交α欧美精品济南到| 在线观看www视频免费| 两性夫妻黄色片| 性色av一级| 免费看十八禁软件| 欧美日韩黄片免| 国产精品一国产av| 伊人久久大香线蕉亚洲五| 久久这里只有精品19| 精品人妻1区二区| 久久久久久亚洲精品国产蜜桃av| 美女午夜性视频免费| 日韩av免费高清视频| 一区二区日韩欧美中文字幕| 搡老乐熟女国产| 免费在线观看完整版高清| 色视频在线一区二区三区| 国产一区二区三区综合在线观看| 91精品伊人久久大香线蕉| 午夜福利在线免费观看网站| 免费在线观看视频国产中文字幕亚洲 | 男女下面插进去视频免费观看| 成年人午夜在线观看视频| 欧美精品人与动牲交sv欧美| 一边摸一边做爽爽视频免费| 91国产中文字幕| 国产一级毛片在线| 亚洲国产av新网站| 免费少妇av软件| 日韩av免费高清视频| av视频免费观看在线观看| 亚洲精品国产色婷婷电影| 黄片小视频在线播放| 青青草视频在线视频观看| 国产一区二区三区综合在线观看| 99久久99久久久精品蜜桃| 人体艺术视频欧美日本| 亚洲午夜精品一区,二区,三区| 大型av网站在线播放| 亚洲专区国产一区二区| 成人黄色视频免费在线看| 天天躁夜夜躁狠狠躁躁| 首页视频小说图片口味搜索 | 午夜免费鲁丝| 久久鲁丝午夜福利片| 国产成人免费无遮挡视频| 大陆偷拍与自拍| 热99国产精品久久久久久7| 精品福利永久在线观看| 一边摸一边做爽爽视频免费| 精品熟女少妇八av免费久了| 极品少妇高潮喷水抽搐| 可以免费在线观看a视频的电影网站| 美女中出高潮动态图| 国产高清videossex| 一本一本久久a久久精品综合妖精| 脱女人内裤的视频| 精品久久久久久久毛片微露脸 | 丰满饥渴人妻一区二区三| 麻豆国产av国片精品| 自线自在国产av| 国产精品亚洲av一区麻豆| 国产又爽黄色视频| 国产日韩一区二区三区精品不卡| 黄色视频不卡| 女人爽到高潮嗷嗷叫在线视频| 中文字幕另类日韩欧美亚洲嫩草| 国产黄色视频一区二区在线观看| 后天国语完整版免费观看| 丝瓜视频免费看黄片| 两人在一起打扑克的视频| 丝袜喷水一区| 一级黄色大片毛片| 最黄视频免费看| 国产日韩欧美亚洲二区| 国产精品 欧美亚洲| 高清视频免费观看一区二区| 国产97色在线日韩免费| 十八禁高潮呻吟视频| 国产一区亚洲一区在线观看| 亚洲精品久久成人aⅴ小说| 中文字幕最新亚洲高清| 黄色一级大片看看| 国产精品熟女久久久久浪| 麻豆av在线久日| 美女视频免费永久观看网站| 性少妇av在线| 人妻 亚洲 视频| 久久热在线av| 国产男女超爽视频在线观看| 亚洲国产精品成人久久小说| 最新的欧美精品一区二区| 国产视频一区二区在线看| 久久九九热精品免费| 亚洲av成人不卡在线观看播放网 | 人妻人人澡人人爽人人| 日韩av不卡免费在线播放| 国产一级毛片在线| 国产精品亚洲av一区麻豆| 精品一品国产午夜福利视频| 午夜免费鲁丝| 午夜福利免费观看在线| 欧美日韩精品网址| 在线观看免费日韩欧美大片| 亚洲av成人不卡在线观看播放网 | 久久人人爽av亚洲精品天堂| 色网站视频免费| 777久久人妻少妇嫩草av网站| 如日韩欧美国产精品一区二区三区| 老熟女久久久| 国产亚洲午夜精品一区二区久久| www.999成人在线观看| 人妻人人澡人人爽人人| 国产精品 国内视频| 欧美成人精品欧美一级黄| 18禁裸乳无遮挡动漫免费视频| 另类亚洲欧美激情| 丁香六月欧美| 纵有疾风起免费观看全集完整版| 黄色 视频免费看| 亚洲免费av在线视频| 久久毛片免费看一区二区三区| 丁香六月欧美| 人成视频在线观看免费观看| 欧美 日韩 精品 国产| 18禁裸乳无遮挡动漫免费视频| 久久天堂一区二区三区四区| 好男人视频免费观看在线| 亚洲国产欧美一区二区综合| 晚上一个人看的免费电影| 成年女人毛片免费观看观看9 | 后天国语完整版免费观看| 久久久久久久国产电影| 又粗又硬又长又爽又黄的视频| 高清视频免费观看一区二区| 一级毛片我不卡| 久久久亚洲精品成人影院| 精品卡一卡二卡四卡免费| 51午夜福利影视在线观看| 亚洲精品一卡2卡三卡4卡5卡 | 亚洲成国产人片在线观看| 欧美老熟妇乱子伦牲交| 亚洲国产欧美网| 色精品久久人妻99蜜桃| 国产高清视频在线播放一区 | 老司机亚洲免费影院| 香蕉国产在线看| 久久久国产欧美日韩av| 中文字幕最新亚洲高清| 美女脱内裤让男人舔精品视频| 一级a爱视频在线免费观看| 视频区图区小说| 国产主播在线观看一区二区 | 满18在线观看网站| 欧美日韩黄片免| 少妇的丰满在线观看| 宅男免费午夜| av国产久精品久网站免费入址| 1024视频免费在线观看| 伊人久久大香线蕉亚洲五| 国产亚洲午夜精品一区二区久久| 成人免费观看视频高清| 九色亚洲精品在线播放| 精品久久久久久电影网| 亚洲国产看品久久| 我的亚洲天堂| 大型av网站在线播放| 国产日韩欧美视频二区| 国产成人91sexporn| 宅男免费午夜| 一本大道久久a久久精品| 欧美成狂野欧美在线观看| 国产精品久久久av美女十八| av天堂久久9| 欧美日韩综合久久久久久| 久久国产精品影院| 少妇人妻久久综合中文| 黄片小视频在线播放| 欧美人与性动交α欧美软件| 亚洲色图 男人天堂 中文字幕| 国产免费又黄又爽又色| xxxhd国产人妻xxx| 男的添女的下面高潮视频| 成人国产一区最新在线观看 | 久久精品人人爽人人爽视色| 国产精品99久久99久久久不卡| 操出白浆在线播放| 亚洲人成网站在线观看播放| 王馨瑶露胸无遮挡在线观看| 亚洲五月婷婷丁香| 搡老岳熟女国产| 久久青草综合色| 国产精品久久久久久人妻精品电影 | 黑人欧美特级aaaaaa片| 黄色一级大片看看| 一级黄色大片毛片| 成人国语在线视频| 亚洲国产日韩一区二区| 亚洲国产精品成人久久小说| 国产熟女午夜一区二区三区| 精品视频人人做人人爽| 99久久综合免费| 国产黄色免费在线视频| 丁香六月天网| 高潮久久久久久久久久久不卡| 欧美xxⅹ黑人| 男人舔女人的私密视频| 国产精品一区二区在线不卡| 男女之事视频高清在线观看 | 只有这里有精品99| 国产精品久久久av美女十八| 高清不卡的av网站| 18在线观看网站| 欧美精品高潮呻吟av久久| 国产av国产精品国产| 中国美女看黄片| 波多野结衣av一区二区av| av国产精品久久久久影院| 日韩一本色道免费dvd| 91成人精品电影| 亚洲色图综合在线观看| 中国美女看黄片| 亚洲国产毛片av蜜桃av| 日日摸夜夜添夜夜爱| 1024视频免费在线观看| 亚洲国产毛片av蜜桃av| 少妇猛男粗大的猛烈进出视频| 免费高清在线观看视频在线观看| 在线观看一区二区三区激情| 国产99久久九九免费精品| 麻豆av在线久日| 国产精品久久久av美女十八| 在线观看免费高清a一片| 只有这里有精品99| 男女无遮挡免费网站观看| av国产久精品久网站免费入址| 如日韩欧美国产精品一区二区三区| 在线观看免费视频网站a站| 精品久久久精品久久久| 99香蕉大伊视频| 欧美变态另类bdsm刘玥| 2021少妇久久久久久久久久久|