張永紅,陳軍
1. 中國測繪科學(xué)研究院,北京 100036;2. 國家基礎(chǔ)地理信息中心,北京 100830
地面沉降是一種可由多種因素引起的地面標(biāo)高緩慢降低的環(huán)境地質(zhì)現(xiàn)象,嚴(yán)重時會成為一種地質(zhì)災(zāi)害(薛禹群等,2003)。引起地面沉降的原因分為自然因素與人為因素。自然因素主要包括地震活動、地質(zhì)構(gòu)造作用以及地表松散層或半松散層在重力作用下的致密硬化過程。人為因素包括地下水(含地下熱水)、石油、天然氣及地下固體礦產(chǎn)資源開采,建筑物荷載與工程建設(shè)等。對于目前我國地面沉降而言,人為因素的作用已多于自然因素(何慶成等,2006b)。我國的主要沉降區(qū)是由過量開采地下水造成的(薛禹群等,2003)。
地面沉降是一種累進(jìn)性的緩變地質(zhì)災(zāi)害。其發(fā)展過程是不可逆的,一旦形成便難以恢復(fù)(殷躍平等,2005)。地面沉降對資源環(huán)境安全的影響是多方面的,主要有:①地面沉降可導(dǎo)致建筑物墻面開裂、墻面傾斜、地面開裂、地基下沉、建筑物傾斜甚至倒塌等;②地面沉降可對通過沉降區(qū)的管網(wǎng)設(shè)施產(chǎn)生剪切作用,使其斷裂或破損;③影響鐵路、公路、地鐵、機場等交通設(shè)施及河道航運;④加劇洪澇災(zāi)害,包括加劇內(nèi)澇、導(dǎo)致防洪堤壩的抗洪能力下降、引起沿海地區(qū)的海水入侵等;⑤破壞生態(tài)環(huán)境,造成地表塌陷、地裂縫、土地荒漠化等生態(tài)環(huán)境傷害。
由于地面沉降廣泛性和造成的巨大損失,我國政府對地面沉降的監(jiān)測和防治高度重視,先后出臺了《全國地面沉降防治規(guī)劃(2011—2020 年)》《京津冀平原地面沉降綜合防治總體規(guī)劃(2019—2035年)》等重要規(guī)劃。地面沉降的準(zhǔn)確監(jiān)測是進(jìn)行地面沉降預(yù)防和治理的基礎(chǔ)。經(jīng)過多年的科技攻關(guān)和發(fā)展,我國的地面沉降監(jiān)測已從以往的以水準(zhǔn)測量、分層標(biāo)測量等地面手段為主,發(fā)展到以衛(wèi)星合成孔徑雷達(dá)干涉測量(interferometric synthetic aperture radar,InSAR)遙感技術(shù)為主,能夠快速獲取大區(qū)域、高精度、長時序的地面沉降信息。例如,京津冀平原區(qū)、北京市、西安市等地區(qū)數(shù)十年地面沉降監(jiān)測,以及江蘇全省地面沉降監(jiān)測等(Zhang等,2016,2021;Dong 等,2023;Li 等,2023)。
空間規(guī)劃是社會經(jīng)濟(jì)發(fā)展到一定階段后,為解決空間矛盾采取的政策措施與工具(蔡玉梅和高平,2013)。國土空間規(guī)劃是優(yōu)化生態(tài)環(huán)境保護(hù)、保障糧食安全、促進(jìn)國土資源集約利用的國家意志導(dǎo)向,是推進(jìn)生態(tài)文明建設(shè)的關(guān)鍵舉措(陳磊和姜海,2021)。地面沉降作為一種重要的自然資源時空信息,對于科學(xué)編制國土空間規(guī)劃、開展山水林田湖草沙海的整體保護(hù)、系統(tǒng)修復(fù)和綜合治理、保障國家生態(tài)安全等具有重要意義(陳軍等,2022)。在地面沉降信息主要依靠水準(zhǔn)、分層標(biāo)等地面手段獲取的時代,由于地面測點的稀疏,得到的地面沉降信息無論是在空間覆蓋還是沉降值域方面都是不完整的,很難真正被利用。2016 年以后,InSAR技術(shù)被大規(guī)模應(yīng)用于全國幾大主要平原區(qū)的地面沉降定期監(jiān)測,獲取了區(qū)域全覆蓋、高空間分辨率(幾米至幾十米)、高精度(毫米級)的地面沉降信息。這樣的地面沉降信息才真正成為一種重要的自然資源時空信息,具備服務(wù)國土空間規(guī)劃編制等行業(yè)應(yīng)用的潛力。從已有研究來看,將地質(zhì)災(zāi)害信息用于國土空間規(guī)劃編制或者資源環(huán)境承載力評價和國土空間開發(fā)適宜性評價的研究僅有少部分。鄭續(xù)等(2020)納入地質(zhì)災(zāi)害因子開展了延川縣國土空間開發(fā)適宜性評價;趙增玉等(2021)將地質(zhì)環(huán)境安全作為限制性因子完成了泰州市國土空間開發(fā)適宜性評價;武漢市基于城市地質(zhì)調(diào)查成果編制了建設(shè)用地地質(zhì)環(huán)境(彭漢發(fā)等,2018);上海市將地面沉降、地下水環(huán)境等生態(tài)地質(zhì)環(huán)境信息用于《上海市國土空間生態(tài)修復(fù)規(guī)劃》編制(黎兵等,2023)。這些成果為不同層級(省、市、縣等)的局地案例,而且地質(zhì)災(zāi)害信息參與應(yīng)用的深淺不一。因此,如何在全國范圍內(nèi)不同層級的國土空間規(guī)劃編制的理論和實踐中充分利用地面沉降信息,是一個亟待深入研究的復(fù)雜課題。
本文首先對我國的地面沉降狀況及主要的沉降區(qū)信息進(jìn)行概述,其次從理論上指出地面沉降與國土空間規(guī)劃的關(guān)聯(lián),最后以北京市為例,分析國土空間規(guī)劃的變更及其與地面沉降空間分布之間的關(guān)系,以期為推動地面沉降信息在國土空間規(guī)劃編制中大規(guī)模深入應(yīng)用提供一點參考。
據(jù)中國地質(zhì)調(diào)查局自然資源航空物探遙感中心參與編寫的《全國地面沉降防治工作報告(2018)》,我國地面沉降易發(fā)區(qū)面積達(dá)117.7 萬km2,其中,高易發(fā)區(qū)21.5 萬km2,中易發(fā)區(qū)12.2 萬km2,主要分布于華北平原、長江三角洲、汾渭盆地、珠江三角洲等平原地區(qū)。全國共有21 個?。▍^(qū)、市)的103 個地級行政區(qū)發(fā)生地面沉降。
京津冀地區(qū)是當(dāng)前我國地面沉降發(fā)展速度最快、發(fā)生范圍最廣、直接影響最嚴(yán)重的地區(qū)(何慶成等,2006a)。自20世紀(jì)90年代出現(xiàn)快速地面下沉以來,北京市、天津市周邊區(qū)縣及河北平原出現(xiàn)大面積連片沉降,最大速率均超過100 mm/a,部分沉降中心達(dá)到180 mm/a(Zhang 等,2016;Zhou 等,2020)。地面沉降中心與地下水水位降落漏斗分布一致,均呈現(xiàn)出不斷擴展的趨勢。
2016 年以來的監(jiān)測表明,隨著南水北調(diào)和地下水控采措施的實施,北京市、天津市地面沉降防控取得了顯著成效,地面沉降狀況逐年減緩,京津兩市周邊地面沉降嚴(yán)重區(qū)(沉降速率超過50 mm/a 的地區(qū))面積大幅減少(Du 等,2021;Dong 等,2023)。河北省地面沉降主要由農(nóng)業(yè)灌溉開采地下水引起,因此,地面沉降狀況隨著氣候條件的變化呈現(xiàn)較大的年際波動。豐水年降水多,地下水開采減少,地面沉降減緩;反之,枯水年份則開采加大,地面沉降加劇。2016~2019 年河北省地面沉降嚴(yán)重區(qū)面積依次為8784 km2、21905 km2、19556 km2、23224 km2(中國測繪科學(xué)研究院,2020)。
長江三角洲地區(qū)地面沉降發(fā)育較早,20 世紀(jì)20 年代上海市出現(xiàn)地面沉降,70 年代起,長江三角洲的蘇錫常、杭嘉湖等主要城市也相繼出現(xiàn)地面沉降,其中,至2007 年上海中心城區(qū)最大累計沉降量超過3 m,蘇錫常等地區(qū)因嚴(yán)重的差異性沉降作用出現(xiàn)地裂縫(龔士良等,2008;王光亞等,2019)。2011 年以來,長江三角洲經(jīng)濟(jì)區(qū)內(nèi)通過區(qū)域聯(lián)防聯(lián)控措施,地面沉降得到了有效控制,以往大于30 mm/a 的江陰—武進(jìn)、吳江市盛澤、浙江慈溪等沉降區(qū)均明顯減緩。
當(dāng)前,長江三角洲地區(qū)地面沉降呈零星分布。江蘇徐州、連云港填海造陸區(qū)、鹽城化工區(qū)、南通沿海地區(qū),上海浦東機場和臨港新城區(qū),浙江溫黃平原、溫瑞平原等均為局部沉降中心。2016~2018年,徐州礦區(qū)因地下工作面開采,出現(xiàn)了年度沉降超過1 m 的情況,此外,連云港徐圩港區(qū)填海造地形成的短期最大地面沉降速率超過 200 mm/a(Zhang 等,2021)。長江三角洲其他沉降中心的沉降速率絕大多數(shù)為15~80 mm/a。2016~2018 年長三角地區(qū)地面沉降嚴(yán)重區(qū)面積為55 km2。
汾渭盆地地面沉降總體較為嚴(yán)重。2019 年度,汾渭盆地地面沉降嚴(yán)重區(qū)面積達(dá)到550 km2,分布于太原盆地和運城盆地,大型沉降區(qū)帶明顯,因地下水開采與汾渭盆地特殊的地質(zhì)條件,伴隨多條順構(gòu)造線展布的地裂縫。
山西省境內(nèi)主要沉降區(qū)位于太原市區(qū)南部、榆次—清徐—太谷—祁縣、文水、孝義—介休等地區(qū);運城市與夏縣交界地區(qū)沉降范圍擴大,臨汾市稷山—新絳沉降區(qū)發(fā)展為山區(qū)平原沉降區(qū)帶。
西安市受地面沉降、地裂縫影響顯著,市區(qū)南部、西南部沉降區(qū)面積近年來逐年減緩。西安市魚化寨沉降中心一直以來屬西安市最嚴(yán)重的沉降區(qū),平均速率超過100 mm/a,2019 年因城中村拆遷改造完成,出現(xiàn)地面回彈,最大回彈量超過20 mm/a(Li 等,2023)。
珠江三角洲地區(qū)受軟土自然固結(jié)影響,存在多個零星分布小幅沉降區(qū)。該地區(qū)為珠江河口堆積形成的復(fù)合型三角洲,平原區(qū)淤泥類軟弱土層較厚,在城市建設(shè)、地下水開采、地下開發(fā)等作用下,土體易于變形,地面出現(xiàn)緩慢下沉(謝先明等,2021)。
2019 年度,珠江三角洲地區(qū)年沉降量大于30 mm 的區(qū)域總面積為8.5 km2,沉降中心數(shù)量達(dá)20 個以上,主要分布于江門市、中山市、揭陽市、汕頭市與珠海市等地。部分沉降中心出現(xiàn)了建筑物開裂、地基與路面交接開裂等現(xiàn)象。
2019 年 5 月發(fā)布的 《中共中央 國務(wù)院關(guān)于建立國土空間規(guī)劃體系并監(jiān)督實施的若干意見》明確了建立國土空間規(guī)劃體系并監(jiān)督實施的“頂層設(shè)計”,指出“在資源環(huán)境承載能力和國土空間開發(fā)適宜性評價的基礎(chǔ)上,科學(xué)有序統(tǒng)籌布局生態(tài)、農(nóng)業(yè)、城鎮(zhèn)等功能空間,劃定生態(tài)保護(hù)紅線、永久基本農(nóng)田、城鎮(zhèn)開發(fā)邊界等空間管控邊界以及各類海域保護(hù)線”??梢?,資源環(huán)境承載力評價和國土空間開發(fā)適宜性評價(即“雙評價”)是編制實施國土空間規(guī)劃的科學(xué)基礎(chǔ)。
資源環(huán)境承載力是從分類到綜合的承載力概念統(tǒng)稱,是包括各類資源要素承載力、環(huán)境要素承載力、生態(tài)承載力以及各要素承載力綜合與集成評價的集合概念。開展資源環(huán)境承載力評價的主要目的是通過開展土地、水、生態(tài)、環(huán)境、災(zāi)害等要素的單項評價和集成評價,了解區(qū)域的資源稟賦和環(huán)境本底,分析區(qū)域發(fā)展的優(yōu)勢與制約因素(郝慶等,2019)。
根據(jù)自然資源部2020 年發(fā)布的《資源環(huán)境承載能力與國土空間開發(fā)適宜性評價指南(試行)》,國土空間開發(fā)適宜性定義為“在維系生態(tài)系統(tǒng)健康和國土安全的前提下,綜合考慮資源環(huán)境等要素條件,特定國土空間進(jìn)行農(nóng)業(yè)生產(chǎn)、城鎮(zhèn)建設(shè)等人類活動的適宜程度”。指標(biāo)的確定是國土空間開發(fā)適宜性評價的關(guān)鍵,雖然不同評價框架側(cè)重點不同,但主要評價指標(biāo)都集中在生態(tài)環(huán)境狀態(tài)、災(zāi)害風(fēng)險、水土資源存量、社會經(jīng)濟(jì)發(fā)展基礎(chǔ)與潛力、交通區(qū)位等方面(喻忠磊等,2015)。其中的災(zāi)害風(fēng)險就主要包括地面沉降、滑坡、坍塌、地裂縫等地質(zhì)災(zāi)害。
可以看出,地面沉降作為一種地質(zhì)災(zāi)害,是“雙評價”都應(yīng)該考慮的因素。但是,地質(zhì)成果(包括地面沉降調(diào)查成果)在支撐服務(wù)國土空間規(guī)劃上仍然面臨內(nèi)容不明確、程序缺失、產(chǎn)品難以得到有效應(yīng)用的窘境(黎兵等,2023)。造成這種局面的原因是多方面的。在體制機制層面,目前尚沒有強制性法律法規(guī)來明確要求國土空間規(guī)劃必須開展相關(guān)地質(zhì)工作(程光華等,2018)。在技術(shù)層面,缺乏針對國土空間規(guī)劃的地質(zhì)調(diào)查工作的要求和標(biāo)準(zhǔn),如針對不同層級、不同類型的規(guī)劃,地質(zhì)調(diào)查的尺度和精度分別是多少并不明確。
北京市有記載的地面沉降最早出現(xiàn)在1935 年的東單至西單一帶,后逐步擴展至平原地區(qū)(賈三滿等,2006)。結(jié)合賈三滿等(2006)、張永紅等(2016)已有成果,本文將北京市地面沉降分為五個階段,即形成階段(1955~1973 年)、發(fā)展階段(1973~1983 年)、擴展階段(1987~1999 年)、快速發(fā)展階段(1999~2016 年)、減緩階段(2016 年至今)。圖1 展示了利用時序InSAR 技術(shù)分別獲取的2012~2019 年六個時段的北京市地面沉降分級區(qū)劃,按地面沉降速率從小到大,依次為穩(wěn)定區(qū)(沉降速率小于10 mm/a)、一般區(qū)(沉降速率大于10 mm/a且小于30 mm/a)、較嚴(yán)重區(qū)(沉降速率大于30 mm/a且小于50 mm/a)、嚴(yán)重區(qū)(沉降速率大于50 mm/a 且小于100 mm/a)、非常嚴(yán)重區(qū)(沉降速率大于100 mm/a)。在前三個時段,即2012 年1 月~2016 年11 月,北京市地面沉降處于不斷發(fā)展階段。沉降漏斗主要分布在朝陽區(qū)東部、通州區(qū)西北部,以及海淀區(qū)北部、昌平區(qū)南部、順義區(qū)南部。此外大興區(qū)南部也有零星分布,三環(huán)以內(nèi)的主城區(qū)長期較穩(wěn)定。自2017年以后,北京市地面沉降呈現(xiàn)穩(wěn)步減緩的趨勢,地面沉降嚴(yán)重區(qū)和非常嚴(yán)重區(qū)的面積大幅減少。
圖1 2012~2019 年北京市地面沉降分級區(qū)劃圖Fig.1 Land subsidence zoning map in Beijing from 2012 to 2019
為了適應(yīng)首都現(xiàn)代化建設(shè)的需要,2002 年5 月北京市第九次黨代會提出了修編北京城市總體規(guī)劃的工作任務(wù),2004 年《北京城市總體規(guī)劃(2004 年—2020 年)》(簡稱北京04 規(guī)劃)被批復(fù)。根據(jù)該規(guī)劃,北京市將構(gòu)建“兩軸–兩帶–多中心”的城市空間結(jié)構(gòu)。其中,兩軸為沿長安街的東西軸和傳統(tǒng)中軸線的南北軸;兩帶是包括通州、順義、亦莊、懷柔、密云、平谷的東部發(fā)展帶,以及包括大興、房山、昌平、延慶、門頭溝的西部發(fā)展帶,如圖2 所示。
圖2 《北京城市總體規(guī)劃(2004 年—2020 年)》中“兩軸–兩帶–多中心”的空間結(jié)構(gòu)Fig.2 Spatial structure of “two axes – two belts – multiple centers” in the Beijing Overall Spatial Planning (2004-2020)
東部發(fā)展帶中的平原地區(qū)稱為東部次區(qū)域。對比圖2、圖1,直到2018 年,北京市地面沉降嚴(yán)重區(qū)主要分布在東部次區(qū)域。為了防止東部次區(qū)域地面沉降進(jìn)一步擴展,勢必要采取地下水控采、限制高層建筑密集建設(shè)等措施,這與東部次區(qū)域作為“未來重點發(fā)展的地區(qū)”“承接中心城人口、職能疏解和新的產(chǎn)業(yè)集聚的主要地區(qū)”等發(fā)展定位存在一定的矛盾。
北京04 規(guī)劃批復(fù)后,北京的經(jīng)濟(jì)社會發(fā)展和城鄉(xiāng)面貌發(fā)生了巨大變化,部分之前設(shè)定的2020年發(fā)展指標(biāo)已經(jīng)被突破,一些深層次矛盾和問題也逐步顯現(xiàn)。2014 年,北京市政府啟動了新的城市總體規(guī)劃的編制。2017 年9 月《北京城市總體規(guī)劃(2016 年—2035 年)》(簡稱北京新規(guī)劃)被批復(fù)。北京新規(guī)劃中提出構(gòu)建“一核一主一副、兩軸多點一區(qū)”的城市空間結(jié)構(gòu),將北京04 規(guī)劃的“兩帶–多中心”改成了“多點一區(qū)”,東部發(fā)展帶不再被提及,并明確指出“綜合評估地面沉降、活動斷裂、巖溶塌陷、砂土液化、地下有害氣體、地下采空、地下水位變化等災(zāi)害因素,消除災(zāi)害隱患,確保地上地下空間安全”。同時,為保護(hù)和修復(fù)自然生態(tài)系統(tǒng),優(yōu)化生態(tài)空間格局,北京新規(guī)劃將北京市域按照生態(tài)控制區(qū)、集中建設(shè)區(qū)、限制建設(shè)區(qū)(“三區(qū)”)進(jìn)行了詳細(xì)的空間區(qū)分,如圖3 所示。
圖3 北京新規(guī)劃中的“三區(qū)”分布圖Fig.3 Distribution map of ecological control zone, construction restriction zone, and construction zone in Beijing Overall Spatial Planning (2016-2035)
將圖3、圖1(f)進(jìn)行對比,結(jié)果如圖4 所示,可以看到,北京市全部地面沉降嚴(yán)重區(qū)和非常嚴(yán)重區(qū)中很大部分在北京新規(guī)劃中被劃為生態(tài)控制區(qū)、限制建設(shè)區(qū),這在北京市東部沉降帶(圖4(a)的黑框區(qū)域)更加明顯。將圖3 和圖1(f)進(jìn)行疊加分析,統(tǒng)計了北京市域沉降速率超過30 mm/a 的區(qū)域(含沉降較嚴(yán)重區(qū)、嚴(yán)重區(qū)、非常嚴(yán)重區(qū))、沉降速率超過50 mm/a 的區(qū)域(含沉降嚴(yán)重區(qū)、非常嚴(yán)重區(qū))中三區(qū)的面積占比,結(jié)果如表1 所示。
表1 北京市2019 年地面沉降區(qū)中生態(tài)控制區(qū)、限制建設(shè)區(qū)、集中建設(shè)區(qū)面積占比Tab. 1 Percentage of ecological control areas,restricted construction areas,and concentrated construction zone in the subsiding regions in Beijing in 2019%
圖4 2019 年北京市主要沉降區(qū)(a)與對應(yīng)的北京新規(guī)劃中的“三區(qū)”分布(部分;b)Fig.4 Main subsidence areas in Beijing in 2019 (a) and the corresponding distribution of the “three zones” in Beijing Overall Spatial Planning (2016-2035) (Part; b)
沉降速率超過30 mm/a 的區(qū)域中,36.6%的區(qū)域在北京新規(guī)劃中被劃定為生態(tài)控制區(qū),22.6%的區(qū)域被劃定為限制建設(shè)區(qū),二者合計占59.2%,只有40.8%的面積為集中建設(shè)區(qū)。沉降速率超過50 mm/a的區(qū)域中,24.1%的區(qū)域在北京新規(guī)劃中被劃定為生態(tài)控制區(qū),44.3%的區(qū)域被劃定為限制建設(shè)區(qū),二者合計占比高達(dá)68.4%,只有31.6%的面積為集中建設(shè)區(qū)??梢钥闯?,北京新規(guī)劃的編制過程中充分考慮了地面沉降這一災(zāi)害風(fēng)險對北京“建設(shè)國際一流的和諧宜居之都”的城市發(fā)展目標(biāo)的限制性影響,將存在較嚴(yán)重區(qū)及以上地面沉降的地區(qū)大部分劃為了生態(tài)控制區(qū)和限制建設(shè)區(qū),切實降低了平原地區(qū)開發(fā)強度。這些舉措對地面沉降的減緩和治理及地面沉降風(fēng)險的控制都會產(chǎn)生直接效果,對推動北京市可持續(xù)發(fā)展、高質(zhì)量發(fā)展、生態(tài)文明建設(shè)等具有重要意義。
地面沉降作為影響我國經(jīng)濟(jì)社會可持續(xù)發(fā)展的主要地質(zhì)災(zāi)害之一,可能會在未來一段時間內(nèi)持續(xù)存在,這本質(zhì)上是由經(jīng)濟(jì)社會發(fā)展對各類資源(包括地下水、煤炭、石油、天然氣等礦產(chǎn)和能源資源,以及土地資源)的持續(xù)需求所決定的,如河北平原作為我國主要的糧食生產(chǎn)區(qū),農(nóng)業(yè)灌溉對地下水的需求是長期持續(xù)的。因此,地面沉降作為一種地質(zhì)災(zāi)害因素,在國土空間規(guī)劃編制中必須加以考慮。
從北京市國土空間規(guī)劃與地面沉降狀況看,北京04 規(guī)劃中的東部發(fā)展帶是北京市地面沉降嚴(yán)重區(qū)的主要分布地。北京新規(guī)劃中,地面沉降速率超過30 mm/a和50 mm/a的地區(qū),分別有59.2%、68.4%的國土被劃為生態(tài)控制區(qū)和限制建設(shè)區(qū)。這充分表明地面沉降信息在北京市國土空間規(guī)劃編制中發(fā)揮了重要作用。
自《中共中央 國務(wù)院關(guān)于建立國土空間規(guī)劃體系并監(jiān)督實施的若干意見》于2019 年發(fā)布后,《全國國土空間規(guī)劃綱要(2021—2035 年)》于2023 年印發(fā)。目前,多個省級國土空間規(guī)劃相繼獲批,國土空間詳細(xì)規(guī)劃編制工作已經(jīng)啟動。在各級國土空間總體規(guī)劃和詳細(xì)規(guī)劃編制過程中,如何充分利用InSAR 技術(shù)獲取的高精度、區(qū)域全覆蓋、高分辨率的地面沉降信息,進(jìn)一步提升規(guī)劃的科學(xué)性,是下一步亟待研究和解決的關(guān)鍵問題。
致謝本文部分資料來源于陳軍院士牽頭編制的《基于地理空間數(shù)據(jù)的國土空間規(guī)劃分析報告之三——全國地面沉降狀況分析》。中國自然資源航空物探航遙中心葛大慶研究員、長安大學(xué)趙超英教授、中南大學(xué)胡駿教授分別提供了有關(guān)全國、汾渭盆地和珠江三角洲的地面沉降資料。在此一并表示感謝。