劉銳金，劉 東，何長(zhǎng)輝，伍 薇，陳幫乾

（中國(guó)熱帶農(nóng)業(yè)科學(xué)院橡膠研究所，?？?571101）

一個(gè)地區(qū)內(nèi)城鎮(zhèn)和農(nóng)村的人口結(jié)構(gòu)是反映經(jīng)濟(jì)社會(huì)發(fā)展?fàn)顩r的一項(xiàng)重要指標(biāo)（厲以寧，2012），協(xié)調(diào)推進(jìn)城鎮(zhèn)化是實(shí)現(xiàn)現(xiàn)代化的重大戰(zhàn)略選擇（李克強(qiáng)，2012）。中國(guó)城市化過(guò)程中存在“化地不化人”的不合理傾向（文貫中 等，2012），相比土地城鎮(zhèn)化，人口城鎮(zhèn)化是一個(gè)更為復(fù)雜的過(guò)程，受到諸如戶(hù)口、市民化成本、產(chǎn)業(yè)結(jié)構(gòu)、農(nóng)村產(chǎn)權(quán)制度改革等眾多因素影響，因此滯后于土地城鎮(zhèn)化的發(fā)展（呂添貴 等，2016；吳一凡 等，2018）。新型城鎮(zhèn)化背景下，地區(qū)發(fā)展更應(yīng)重視以人為本?！丁笆奈濉毙滦统擎?zhèn)化實(shí)施方案》（發(fā)改規(guī)劃〔2022〕960 號(hào)）明確，要深入推進(jìn)以人為核心的新型城鎮(zhèn)化戰(zhàn)略，計(jì)劃到2025年全國(guó)常住人口城鎮(zhèn)化率穩(wěn)步提高，戶(hù)籍人口城鎮(zhèn)化率要明顯提高，戶(hù)籍人口城鎮(zhèn)化率與常住人口城鎮(zhèn)化率差距明顯縮小，提高城鎮(zhèn)化的質(zhì)量（國(guó)家發(fā)展改革委，2022）。因此，推動(dòng)高質(zhì)量城鄉(xiāng)融合發(fā)展是未來(lái)較長(zhǎng)一段時(shí)期的重要任務(wù)，尤其是農(nóng)業(yè)轉(zhuǎn)移人口就業(yè)、農(nóng)民工留守家屬隨遷進(jìn)城落戶(hù)。

資源豐裕程度和經(jīng)濟(jì)發(fā)展水平之間的關(guān)系是發(fā)展經(jīng)濟(jì)學(xué)所涉及的重要研究命題之一，資源豐富并不總是促進(jìn)經(jīng)濟(jì)發(fā)展，還可能是阻礙作用，即經(jīng)濟(jì)學(xué)研究中的“資源詛咒”，在某些特殊情境下，資源相對(duì)豐裕但經(jīng)濟(jì)發(fā)展反而相對(duì)緩慢（Ross,1999）。資源與發(fā)展的關(guān)系呈現(xiàn)明顯的空間異質(zhì)性，挪威、加拿大等國(guó)很好地實(shí)現(xiàn)資源促進(jìn)發(fā)展，從石油等資源獲得額外收入用于長(zhǎng)期性投資，實(shí)現(xiàn)長(zhǎng)期經(jīng)濟(jì)社會(huì)發(fā)展，但非洲等地區(qū)部分資源豐富國(guó)家卻長(zhǎng)期陷入資源“詛咒”（Rosa et al., 2012），政府腐敗、產(chǎn)業(yè)結(jié)構(gòu)失衡，陷入長(zhǎng)期停滯。從中國(guó)省市區(qū)層面看，經(jīng)濟(jì)發(fā)達(dá)的省份大多集中在東部沿海地區(qū)，相較于西部省份，這些地區(qū)煤礦、石油等自然資源并不豐裕。具體到海南島，熱帶氣候資源豐富，但經(jīng)濟(jì)發(fā)展相對(duì)滯后，人口城鎮(zhèn)化任務(wù)艱巨，不僅受到制度層面的隱性影響，還有產(chǎn)業(yè)規(guī)模及其結(jié)構(gòu)的作用，地區(qū)就業(yè)空間有限（徐翠枚，2014）。天然橡膠是重要工業(yè)原料，橡膠種植屬于勞動(dòng)密集型農(nóng)業(yè)產(chǎn)業(yè)，也是海南島種植面積最大的農(nóng)業(yè)產(chǎn)業(yè)，形成了較為完整的產(chǎn)業(yè)和技術(shù)體系，同時(shí)還是海南島230萬(wàn)鄉(xiāng)村人口的重要收入來(lái)源之一。2020年末，海南島橡膠種植面積為51.92萬(wàn)hm2，占海南島陸地面積的比例14.8% （海南省統(tǒng)計(jì)局 等，2021），是中國(guó)天然橡膠的兩大產(chǎn)區(qū)之一。按照《關(guān)于建立糧食生產(chǎn)功能區(qū)和重要農(nóng)產(chǎn)品生產(chǎn)保護(hù)區(qū)的指導(dǎo)意見(jiàn)》（國(guó)發(fā)〔2017〕24 號(hào)）要求，中國(guó)劃定120 萬(wàn)hm2天然橡膠生產(chǎn)保護(hù)區(qū)（國(guó)務(wù)院，2017），其中海南56 萬(wàn)hm2，主要分布在海南島中西部山區(qū)。海南正在全島推進(jìn)建設(shè)自由貿(mào)易港，建設(shè)和管護(hù)天然橡膠生產(chǎn)保護(hù)區(qū)，明晰天然橡膠種植生產(chǎn)與人口城鎮(zhèn)化之間的關(guān)系，對(duì)實(shí)現(xiàn)二者協(xié)調(diào)發(fā)展具有重要現(xiàn)實(shí)意義。

橡膠樹(shù)經(jīng)歷定植及8 a左右的非生產(chǎn)期撫管后，進(jìn)入長(zhǎng)達(dá)30 a的割膠生產(chǎn)期，屬于較為典型的長(zhǎng)期性生物資產(chǎn)，且產(chǎn)品無(wú)滯銷(xiāo)情況。割膠是日常性生產(chǎn)活動(dòng)，類(lèi)似于工業(yè)生產(chǎn)，一定規(guī)模的橡膠園意味著一個(gè)就業(yè)崗位，橡膠種植對(duì)人口城鎮(zhèn)化的可能作用渠道主要有：1）土地利用。橡膠種植可能擠占其他農(nóng)業(yè)生產(chǎn)用地，加之橡膠園受到較為嚴(yán)格的管制，抑制了其他需要大量土地的產(chǎn)業(yè)發(fā)展，通過(guò)本地區(qū)財(cái)政收入等方式影響城鎮(zhèn)化。2）勞動(dòng)力。橡膠種植屬于勞動(dòng)密集型產(chǎn)業(yè)，橡膠面積占比越大，意味著農(nóng)村吸附越多的勞動(dòng)力，向城市轉(zhuǎn)移的勞動(dòng)力較少，農(nóng)民收入可能更依賴(lài)農(nóng)業(yè)經(jīng)營(yíng)。3）溢出效應(yīng)。橡膠種植強(qiáng)度變化可能影響周邊地區(qū)，勞動(dòng)力等生產(chǎn)要素可能隨天然橡膠價(jià)格變動(dòng)等因素，在橡膠種植面積占比高和低的地區(qū)流動(dòng)。為了探尋橡膠種植是否對(duì)地區(qū)城鄉(xiāng)就業(yè)與人口結(jié)構(gòu)造成顯著影響，本研究擬運(yùn)用海南島市縣級(jí)面板數(shù)據(jù)，借助中介效應(yīng)模型、空間計(jì)量經(jīng)濟(jì)模型、耦合協(xié)調(diào)度模型，闡釋橡膠種植對(duì)人口城鎮(zhèn)化的影響機(jī)制，找出空間差異。以期在穩(wěn)定橡膠種植面積的情況下，推動(dòng)天然橡膠產(chǎn)業(yè)科學(xué)合理布局，提出促進(jìn)非農(nóng)就業(yè)、加快人口城鎮(zhèn)化協(xié)調(diào)發(fā)展的路徑與措施，使橡膠產(chǎn)業(yè)更好地服務(wù)于地方經(jīng)濟(jì)發(fā)展。

1 文獻(xiàn)綜述

1.1 人口城鎮(zhèn)化與土地城鎮(zhèn)化等指標(biāo)的耦合協(xié)調(diào)關(guān)系

中國(guó)經(jīng)歷了“工業(yè)城鎮(zhèn)化—土地城鎮(zhèn)化—人口城鎮(zhèn)化”3 個(gè)階段（周飛舟 等，2018），人口城鎮(zhèn)化對(duì)經(jīng)濟(jì)增長(zhǎng)的綜合效應(yīng)和投資效應(yīng)大于土地城鎮(zhèn)化（朱紀(jì)廣 等，2020）。按照常住人口統(tǒng)計(jì)的人口城鎮(zhèn)化率并不科學(xué)，進(jìn)城務(wù)工人口雖然大部分時(shí)間居住在城市，但由于戶(hù)籍仍保留在農(nóng)村，并不能很好地享受城市的公共服務(wù)福利，無(wú)法融入城市社會(huì)（寧越敏 等，2019）。因此，常住人口城鎮(zhèn)化率很難準(zhǔn)確反映實(shí)際情況，需要更多關(guān)注質(zhì)量指標(biāo)，如城市生活空間要素、人文體現(xiàn)（薛德升 等，2016）。非農(nóng)就業(yè)是推進(jìn)人口城鎮(zhèn)化的重要環(huán)節(jié)，兩者之間具有長(zhǎng)期穩(wěn)定的關(guān)系（曾令華 等，2007），在經(jīng)濟(jì)發(fā)達(dá)地區(qū)，人口、土地、產(chǎn)業(yè)非農(nóng)化轉(zhuǎn)型的耦合度和協(xié)調(diào)度處于相對(duì)高的水平（楊忍 等，2015），但這是個(gè)長(zhǎng)期的過(guò)程，即使在長(zhǎng)三角地區(qū)，“人口-土地-經(jīng)濟(jì)-社會(huì)”四維協(xié)調(diào)耦合度仍有一定的提升空間（王鋒 等，2018）。一方面，產(chǎn)業(yè)發(fā)展提供就業(yè)機(jī)會(huì)，促進(jìn)人口城鎮(zhèn)化（文樂(lè) 等，2017），另一方面，人口城鎮(zhèn)化能推動(dòng)產(chǎn)業(yè)結(jié)構(gòu)升級(jí)，有正向的溢出效應(yīng)（宋麗敏，2017）。具體到海南島，“人口-土地”城鎮(zhèn)化發(fā)展處于協(xié)調(diào)狀態(tài)的市縣多位于沿海地帶（楊新 等，2022）。

1.2 影響人口城鎮(zhèn)化的主要因素

在勞動(dòng)收益、作業(yè)環(huán)境、保障體系和獲得感方面，外來(lái)務(wù)工人口與本地勞動(dòng)人口之間差距較大，且呈逐年擴(kuò)大趨勢(shì)（楊超 等，2022），其關(guān)鍵原因在于人力資本積累的差異化，但隨著教育發(fā)展更加均衡（謝童偉 等，2021），這種差距會(huì)逐漸縮小。從有效市場(chǎng)和有為政府的關(guān)系看，市場(chǎng)自發(fā)的產(chǎn)業(yè)結(jié)構(gòu)調(diào)整，更有利于人口城鎮(zhèn)化（郭樹(shù)華 等，2021），對(duì)于顯著影響人口跨省遷移的城鄉(xiāng)社保體系，則有賴(lài)于政府的強(qiáng)有力推進(jìn)（城鄉(xiāng)社保一體化與人口遷移課題組，2021）。當(dāng)前，亟待削弱對(duì)戶(hù)籍的歧視甚至完全放開(kāi)對(duì)戶(hù)籍的限制，推動(dòng)農(nóng)業(yè)轉(zhuǎn)移人口落戶(hù)，促進(jìn)區(qū)域人口流動(dòng)（王永進(jìn) 等，2022），改善城市勞動(dòng)者福利（劉瑋瑋，2021）。

1.3 資源豐裕度對(duì)經(jīng)濟(jì)社會(huì)發(fā)展的效應(yīng)

當(dāng)?shù)貎?yōu)勢(shì)資源性產(chǎn)品處于市場(chǎng)繁榮期，所獲得額外收益會(huì)給不同地區(qū)帶來(lái)不同的結(jié)果，部分地區(qū)仍受制于“資源詛咒”，需要避免此類(lèi)資源性行業(yè)對(duì)其他行業(yè)的擠出效應(yīng)，構(gòu)建相互促進(jìn)、互利共生的機(jī)制（張智光，2022）。當(dāng)經(jīng)濟(jì)社會(huì)發(fā)展水平很低時(shí)，農(nóng)產(chǎn)品市場(chǎng)繁榮可大幅度提高家庭和國(guó)家收入，農(nóng)業(yè)發(fā)展有助于降低瘧疾發(fā)生率（Pan et al,2019），但往往難以形成市場(chǎng)以維持社會(huì)長(zhǎng)期可持續(xù)發(fā)展。例如，冬蟲(chóng)夏草的市場(chǎng)不斷擴(kuò)大，顯著地提高了西藏采收地家庭的收入和生活水平，但并沒(méi)有促進(jìn)支撐長(zhǎng)期增長(zhǎng)的有效投資和人力資源積累（Wang et al., 2018）。中美洲和非洲地區(qū)，以咖啡和可可為主要出口的國(guó)家，如埃塞俄比亞、科特迪瓦、洪都拉斯等，其2011—2013年出口額中可可或咖啡占比超過(guò)30%，在咖啡和可可國(guó)際市場(chǎng)價(jià)格大幅度上漲期間，出口收入增量大，但這些國(guó)家的經(jīng)濟(jì)發(fā)展水平卻一直處于較低水平；其他例子包括貝寧、多哥、巴基斯坦的棉花，馬拉維、津巴布韋的煙草等。事實(shí)上，成功實(shí)現(xiàn)經(jīng)濟(jì)轉(zhuǎn)型的案例較少，其取決于不同國(guó)家采取的不同策略，推進(jìn)市場(chǎng)化、優(yōu)化要素配置效率是緩解資源詛咒的一個(gè)有效途徑（邵帥 等，2011）?？铺氐贤邔?guó)內(nèi)可可豆收購(gòu)價(jià)格維持在低水平，種植戶(hù)不但無(wú)法得到國(guó)際市場(chǎng)繁榮所帶來(lái)的紅利，還要承受價(jià)格低谷時(shí)的收入損失（Benjamin et al, 1993），可可豆的大部分收益被精英集團(tuán)在國(guó)外揮霍。新疆作為油氣資源豐富的地區(qū)之一，其資源開(kāi)發(fā)過(guò)程中受“詛咒效應(yīng)”影響，綠色全要素生產(chǎn)率反而下降（許志明 等，2022）。相反，19世紀(jì)中期美國(guó)南北戰(zhàn)爭(zhēng)推高了棉花價(jià)格，埃及棉花出口的額外收益被用于建設(shè)交通設(shè)施和港口，這些基礎(chǔ)設(shè)施很好地促進(jìn)了埃及20世紀(jì)初的工業(yè)化進(jìn)程（Dodwell, 1931）。過(guò)往文獻(xiàn)對(duì)森林資源“詛咒現(xiàn)象”的實(shí)證研究結(jié)論不一，其因研究對(duì)象、時(shí)間、樣本、模型和指標(biāo)選取的不同而變化，對(duì)森林資源詛咒的傳導(dǎo)機(jī)制缺乏實(shí)證檢驗(yàn)（謝煜 等，2019）。

橡膠樹(shù)屬高大喬木，是資源性資產(chǎn)，種植面積占國(guó)土面積的比重在較大程度上代表區(qū)域氣候等資源的豐裕度。人口城鎮(zhèn)化是經(jīng)濟(jì)社會(huì)發(fā)展的重要結(jié)果，適宜用橡膠種植探討資源豐裕度對(duì)海南島人口城鎮(zhèn)化的影響。通過(guò)文獻(xiàn)梳理，尚未發(fā)現(xiàn)有成果探討橡膠種植與人口城鎮(zhèn)化之間的關(guān)系，故本研究借助多種計(jì)量方法識(shí)別上述二者的關(guān)系，同時(shí)將其空間關(guān)聯(lián)因素納入考量。

2 數(shù)據(jù)與研究區(qū)域概況

2.1 數(shù)據(jù)來(lái)源與預(yù)處理

本研究所使用的數(shù)據(jù)主要來(lái)自1989—2021年的《海南省統(tǒng)計(jì)年鑒》（海南省統(tǒng)計(jì)局 等，1989—2021），由于時(shí)期跨度大，部分指標(biāo)統(tǒng)計(jì)口徑有所調(diào)整。為比較不同年度用貨幣表示的指標(biāo)，借鑒文獻(xiàn)（陳詩(shī)一，2011）的做法，利用海南省居民消費(fèi)價(jià)格指數(shù)（1978年=100）對(duì)貨幣表示的指標(biāo)進(jìn)行平減，以1991—2001年平均值為基數(shù)，將貨幣表示的指標(biāo)折算為可比價(jià)格。2015 年《海南省統(tǒng)計(jì)年鑒》公布農(nóng)村常住居民人均可支配收入，取代之前的農(nóng)村居民家庭人均純收入，二者均可反映農(nóng)村居民的收入狀況，且沒(méi)有跳躍的情況，無(wú)需調(diào)整。農(nóng)墾在海南島占有重要地位，曾屬于相對(duì)獨(dú)立的經(jīng)濟(jì)社會(huì)單元，經(jīng)過(guò)多輪改革，逐漸融入地方。2011 年開(kāi)始，海南農(nóng)墾橡膠種植面積不再單獨(dú)統(tǒng)計(jì)，納入各市縣，為綜合整體反映各地區(qū)橡膠種植面積占國(guó)土面積的比例，將2011年之前逐年的農(nóng)墾橡膠種植面積按農(nóng)場(chǎng)場(chǎng)部所在位置，并入各市縣。農(nóng)場(chǎng)橡膠種植面積等數(shù)據(jù)來(lái)自海南農(nóng)墾歷年經(jīng)濟(jì)社會(huì)統(tǒng)計(jì)資料①海南農(nóng)墾經(jīng)濟(jì)社與社會(huì)統(tǒng)計(jì)資料統(tǒng)計(jì)了各國(guó)有農(nóng)場(chǎng)的發(fā)展情況，不對(duì)外發(fā)行，從原海南農(nóng)墾總局獲取1993—2010年的資料。。由于部分國(guó)有農(nóng)場(chǎng)土地跨市縣分布，按農(nóng)場(chǎng)場(chǎng)部所在地并入市縣，可能存在誤差，但由于跨市縣分布的面積較小，不會(huì)對(duì)結(jié)果造成實(shí)質(zhì)性影響。

2.2 海南島人口城鎮(zhèn)化進(jìn)展

圖1為海南自建省辦特區(qū)以來(lái)的逐年人口城鎮(zhèn)化率（1988—2020 年），2004 年戶(hù)籍人口中非農(nóng)戶(hù)口比例較上年增加了10.54 個(gè)百分點(diǎn)，造成突變的主要原因是海南農(nóng)墾國(guó)有農(nóng)場(chǎng)職工及家屬統(tǒng)一辦理“農(nóng)轉(zhuǎn)非”②其他地方也有存在這種情況，據(jù)《盤(pán)錦日?qǐng)?bào)》（2006-01-25，第1版）報(bào)道，2005年盤(pán)錦市墾區(qū)近4萬(wàn)戶(hù)、12萬(wàn)人實(shí)現(xiàn)“農(nóng)轉(zhuǎn)非”；據(jù)《青海日?qǐng)?bào)》（2005-02-04，第1版）報(bào)道，2005年青海墾區(qū)3萬(wàn)人實(shí)現(xiàn)“農(nóng)轉(zhuǎn)非”，納入城鎮(zhèn)低保覆蓋范圍。，全省非農(nóng)戶(hù)口的人數(shù)較2005 年增加89.17 萬(wàn)人；2005—2015 年，戶(hù)籍人口城鎮(zhèn)化率沒(méi)有增長(zhǎng)，反而下降了1.21 個(gè)百分點(diǎn)，此后逐年增長(zhǎng)，2020 年達(dá)44.58%，與全國(guó)水平相當(dāng)。2010 年海南省常住人口城鎮(zhèn)化率與全國(guó)常住人口平均數(shù)基本一致，此后發(fā)展速度相對(duì)緩慢，2020年與全國(guó)常住人口相差3.61個(gè)百分點(diǎn)。

圖1 海南省及全國(guó)人口城鎮(zhèn)化率的變化（1988—2020年）Fig.1 Changes of population urbanization ratios in Hainan and China during 1988-2020

2004年農(nóng)墾國(guó)有農(nóng)場(chǎng)職工及其家屬的戶(hù)口轉(zhuǎn)變對(duì)多個(gè)市縣戶(hù)籍人口城鎮(zhèn)化率產(chǎn)生顯著的影響，保亭縣、瓊中縣、白沙縣的非農(nóng)戶(hù)人口占比分別提高了29.4、27.5、26.3 個(gè)百分點(diǎn)，屯昌縣、萬(wàn)寧市、澄邁縣分別提高了19.6、15.4、15.2 個(gè)百分點(diǎn)（圖2）。農(nóng)墾國(guó)有農(nóng)場(chǎng)人口轉(zhuǎn)入非農(nóng)戶(hù)口后，農(nóng)墾橡膠種植面積、從事種植業(yè)的勞動(dòng)人口越多，該地區(qū)的戶(hù)籍人口城鎮(zhèn)化率可能相對(duì)越高，這實(shí)際由政策調(diào)整造成，應(yīng)予以處理。為了避免2003 與2004 年之間戶(hù)籍人口城鎮(zhèn)化率的跳躍，又考慮到每年均有部分農(nóng)業(yè)戶(hù)籍人口轉(zhuǎn)為城鎮(zhèn)人口，將2004年較上年減少的農(nóng)業(yè)戶(hù)籍人口69.08 萬(wàn)人中的95%，以各市縣減少的農(nóng)業(yè)人口為權(quán)重，調(diào)減2003年各個(gè)市縣的農(nóng)業(yè)戶(hù)籍人口數(shù)，即2003年全省戶(hù)籍人口從570.43萬(wàn)人調(diào)整為504.80 萬(wàn)人，減少65.63 萬(wàn)人。調(diào)整后，2004 年各市縣的戶(hù)籍人口城鎮(zhèn)化率相較于2003 年的變動(dòng)幅度均在4%以?xún)?nèi)。

圖2 海南農(nóng)墾國(guó)有農(nóng)場(chǎng)職工及其家屬“農(nóng)轉(zhuǎn)非”對(duì)海南島市縣人口城鎮(zhèn)化率的影響Fig.2 Impacts of workers and their family members in state-owned farms from rural "Hukou" to urban on population urbanization of counties and cities in Hainan Island

2.3 海南島天然橡膠產(chǎn)業(yè)發(fā)展歷程

1990年以來(lái)海南島橡膠種植業(yè)經(jīng)歷了較為明顯的3 個(gè)階段（圖3）。1）第一階段（1992—2002年）：中國(guó)天然橡膠市場(chǎng)管制逐漸放開(kāi)，1995 年前后完成市場(chǎng)化改革后，迅速融入國(guó)際市場(chǎng)，在亞洲金融危機(jī)沖擊下，天然橡膠價(jià)格大幅度下跌，1998年較1995年下跌了43.21%，2001年較1998年下跌11.29%；農(nóng)業(yè)稅、農(nóng)林特產(chǎn)稅仍在征收，且是地方財(cái)政稅收的主要來(lái)源之一，尤其是在以農(nóng)業(yè)為主的山區(qū)，如海南省白沙縣，其2002年農(nóng)稅收入占地方財(cái)政收入的比例達(dá)56%（陳瓊列 等，2006）；橡膠種植規(guī)模保持穩(wěn)定，2002 年末種植面積36.91 萬(wàn)hm2，與1992 年持平。2）第二階段（2003—2013年）：2003 年海南省瓊山市正式撤市設(shè)區(qū)，并入?？谑校?002年農(nóng)業(yè)稅占全國(guó)稅收總額的比例下降至2.6%（季程遠(yuǎn)，2021），涉及天然橡膠產(chǎn)業(yè)的稅收逐漸減少或取消，減輕了農(nóng)戶(hù)負(fù)擔(dān)。但同時(shí)依賴(lài)于農(nóng)業(yè)稅的主產(chǎn)區(qū)市縣，地方財(cái)政收入也相應(yīng)減少。該階段天然橡膠價(jià)格整體處于上漲周期，2011年價(jià)格 是2003 年 的2.73 倍，2013 年 較2011 年 下 跌41.17%，但仍是2003 年的1.61 倍。持續(xù)上漲的價(jià)格刺激生產(chǎn)，海南島橡膠種植面積從2003 年的37.98 萬(wàn)hm2增加到2013 年的54.02 萬(wàn)hm2，增長(zhǎng)主要來(lái)自民營(yíng)橡膠種植業(yè)。2003年海南農(nóng)墾橡膠種植面積24.61 萬(wàn)hm2，2008 年增加至26.58 萬(wàn)hm2，此后保持相對(duì)穩(wěn)定；2003年海南島民營(yíng)膠園種植面積13.37 萬(wàn)hm2，根據(jù)原農(nóng)業(yè)部農(nóng)墾局統(tǒng)計(jì)③來(lái)自原農(nóng)業(yè)部農(nóng)墾局提供的《全國(guó)熱帶、南亞熱帶作物生產(chǎn)情況（2013年）》，未公開(kāi)發(fā)行。2012年開(kāi)始，《海南省統(tǒng)計(jì)年鑒》不再公布海南農(nóng)墾橡膠種植面積。，2013 年增加至28.92萬(wàn)hm2，較2003年增長(zhǎng)了116.31%。3）第三階段（2014—2020 年）：天然橡膠價(jià)格持續(xù)低位運(yùn)行，2014 年價(jià)格較2011 年下跌56.42%，2018年較2014 年又下跌29.30%，2020 年仍未恢復(fù)到2014年的水平。此階段勞動(dòng)力價(jià)格上漲，橡膠種植業(yè)整體進(jìn)入穩(wěn)定調(diào)整期，種植面積呈減少趨勢(shì)，2020 年較2014 年減少了4.3%。由于3 個(gè)時(shí)間階段具有顯著分異特征，因此在實(shí)證分析部分?jǐn)M采用分階段回歸策略。由散點(diǎn)圖（圖4）可知，2020年各市縣橡膠種植占國(guó)土面積比例和戶(hù)籍人口城鎮(zhèn)化率兩者總體呈負(fù)相關(guān)。

圖3 天然橡膠價(jià)格變化與海南島橡膠種植面積（1990—2020年）Fig.3 Prices of natural rubber and rubber-planted areas in Hainan (1990—2020)

圖4 橡膠種植占海南國(guó)土面積比例和戶(hù)籍人口城鎮(zhèn)化率（2020年）Fig.4 Shares of rubber-planted area and the registered population urbanization ratios (2020) in Hainan

3 變量選擇與研究方法

3.1 變量選擇

1）因變量。最直接反映人口城鎮(zhèn)化的指標(biāo)是城鎮(zhèn)戶(hù)籍或常住人口占比，但部分研究認(rèn)為就業(yè)情況也可以反映人口城鎮(zhèn)化程度（中國(guó)人口與發(fā)展研究中心課題組，2012），變量設(shè)計(jì)具體如下：

①戶(hù)籍人口城鎮(zhèn)化率。城鄉(xiāng)融合發(fā)展是個(gè)復(fù)雜的過(guò)程（陳錫文，2021），農(nóng)民市民化成本高、遷入地政府負(fù)擔(dān)重、產(chǎn)業(yè)人口落戶(hù)難等問(wèn)題，導(dǎo)致農(nóng)民工市民化緩慢（劉銳 等，2014），戶(hù)籍人口城鎮(zhèn)化率顯著滯后于常住人口城鎮(zhèn)化率（戚偉 等，2017）。大部分農(nóng)民工未能實(shí)現(xiàn)市民化，無(wú)法很好地融入城市，使用戶(hù)籍人口城鎮(zhèn)化率更能說(shuō)明城鄉(xiāng)融合的質(zhì)量，但也有文獻(xiàn)采用城區(qū)常住人口占年末總?cè)丝诒嚷时硎救丝诔擎?zhèn)化（劉瓊 等，2018；吳一凡 等，2018），其考慮的是城鎮(zhèn)常住人口也能享受到部分基本公共服務(wù)，因此2個(gè)方面的指標(biāo)都需考慮。2013年起，海南省統(tǒng)計(jì)局和國(guó)家統(tǒng)計(jì)局海南調(diào)查總隊(duì)開(kāi)始公布各市縣的常住人口，這之前僅有戶(hù)籍人口，為了保證數(shù)據(jù)一致性，先使用戶(hù)籍非農(nóng)業(yè)戶(hù)口作為人口城鎮(zhèn)化率的指標(biāo)，再使用2014—2020年常住人口城鎮(zhèn)化率進(jìn)行穩(wěn)健性檢驗(yàn)。

②就業(yè)指標(biāo)。人口流動(dòng)與就業(yè)緊密相關(guān)，引入非農(nóng)就業(yè)勞動(dòng)力占比、職工人數(shù)占比2 個(gè)指標(biāo)。1992—1995年非農(nóng)就業(yè)勞動(dòng)力占比通過(guò)農(nóng)村勞動(dòng)者人數(shù)、社會(huì)勞動(dòng)者人數(shù)計(jì)算得到；1996年社會(huì)勞動(dòng)者人數(shù)變更為從業(yè)人員數(shù)，2009年農(nóng)村勞動(dòng)者人數(shù)改為農(nóng)村家庭從業(yè)人員數(shù)，2013年開(kāi)始采用的是鄉(xiāng)村從業(yè)人員數(shù)。此外，2005年之前，職工人數(shù)指年末全部職工數(shù)量，2005—2013年職工人數(shù)為在崗職工年末人數(shù)，2013年以后指非私營(yíng)單位在崗職工年末人數(shù)，計(jì)算占比的除數(shù)與非農(nóng)就業(yè)勞動(dòng)力占比相同。

2012年起，海南農(nóng)墾從業(yè)人員數(shù)納入所在市縣統(tǒng)一計(jì)算，這使得大部分市縣職工人數(shù)占從業(yè)人員數(shù)的比例有所提高，三亞市2012年市縣職工人數(shù)占比較2011年增加5.2個(gè)百分點(diǎn)，澄邁縣增加4.4個(gè)百分點(diǎn)、瓊中縣增加3.1 個(gè)百分點(diǎn)。由于城鄉(xiāng)從業(yè)人員數(shù)量提高，而鄉(xiāng)村從業(yè)人員數(shù)基本穩(wěn)定，使得鄉(xiāng)村從業(yè)人員占比降低，保亭縣2013年鄉(xiāng)村從業(yè)人員占比較2012年降低9.0個(gè)百分點(diǎn)，白沙縣降低6.4個(gè)百分點(diǎn)、瓊海市降低5.2 個(gè)百分點(diǎn)、三亞市整體降低3.8個(gè)百分點(diǎn)。

2）關(guān)鍵解釋變量。橡膠種植強(qiáng)度包含多維度的含義，如從業(yè)人員占比、產(chǎn)值占比，根據(jù)數(shù)據(jù)可獲得性和指標(biāo)直觀(guān)性，借鑒已有研究（封志明 等，2013），選擇橡膠種植年末種植面積占本地區(qū)國(guó)土面積比例作為橡膠種植強(qiáng)度的代理變量。

3）中介變量。自民營(yíng)膠園快速發(fā)展以來(lái)，天然橡膠是主產(chǎn)區(qū)農(nóng)民的收入來(lái)源之一，其種植面積占比會(huì)對(duì)地方財(cái)政造成影響，特別是存在農(nóng)業(yè)稅的時(shí)期，種植面積占比的影響會(huì)更直接。農(nóng)民收入與地方財(cái)政分別是家庭和政府2個(gè)主體進(jìn)行決策的基礎(chǔ)，擬運(yùn)用中介模型，識(shí)別通過(guò)農(nóng)民收入、地方財(cái)政收入產(chǎn)生影響的中介效應(yīng)（圖5）。

圖5 橡膠種植強(qiáng)度對(duì)人口城鎮(zhèn)化的影響機(jī)制Fig.5 Influencing mechanism of rubber-planted intensity on population urbanization

4）控制變量。已有研究認(rèn)為人口城鎮(zhèn)化與經(jīng)濟(jì)發(fā)展水平（趙金華 等，2009）、產(chǎn)業(yè)結(jié)構(gòu)等緊密相關(guān)，如人均GDP、二三產(chǎn)業(yè)產(chǎn)值水平（秦佳 等，2013）等。結(jié)合數(shù)據(jù)可獲得性，選擇人均GDP、農(nóng)業(yè)和工業(yè)總產(chǎn)值的自然對(duì)數(shù)作為控制變量。

3.2 研究方法

中介效應(yīng)分析參考溫忠麟等（2022）的研究，通過(guò)方程（1）（2）（3）進(jìn)行檢驗(yàn)。地區(qū)的溢出效應(yīng)是指橡膠種植對(duì)周邊地區(qū)產(chǎn)生正向或負(fù)向的作用，采用方程（4） 的空間杜賓模型（Spatial Durbin Model, SDM）檢驗(yàn)。

式中：α、β、γ、ρ、φ為估計(jì)參數(shù)；ε是誤差項(xiàng)；x、y、z如圖5所示，H為控制變量。

以面板數(shù)據(jù)模型基礎(chǔ)，判斷橡膠種植強(qiáng)度是否會(huì)對(duì)人口城鎮(zhèn)化造成“詛咒”效應(yīng)，將不同距離矩陣的空間計(jì)量模型都納入考量。面板數(shù)據(jù)模型能在一定程度上控制無(wú)法觀(guān)測(cè)的異質(zhì)性，有助于解決因遺漏變量而導(dǎo)致的內(nèi)生性問(wèn)題（陳強(qiáng)，2014），其他導(dǎo)致內(nèi)生性的原因有度量誤差、聯(lián)立性。人口城鎮(zhèn)化是經(jīng)濟(jì)社會(huì)多變量綜合作用的結(jié)果，自然條件對(duì)橡膠種植的影響總體要大于經(jīng)濟(jì)社會(huì)，尤其在民營(yíng)膠園快速發(fā)展的第二階段，各市縣農(nóng)業(yè)勞動(dòng)力相對(duì)充足，割膠生產(chǎn)的預(yù)期收入較高，風(fēng)害、降雨等氣候因素的作用更突出，反向因果雖有可能存在，但估計(jì)不足以導(dǎo)致內(nèi)生性問(wèn)題。依據(jù)文獻(xiàn)（王淑佳等，2021）關(guān)于耦合協(xié)調(diào)度的計(jì)算方法，將橡膠種植面積占比、職工人數(shù)占比、非農(nóng)就業(yè)勞動(dòng)力占比、戶(hù)籍人口城鎮(zhèn)化率進(jìn)行標(biāo)準(zhǔn)化處理，仍用原符號(hào)表示，x為橡膠種植面積占國(guó)土面積的比例，y代表人口城鎮(zhèn)化的3項(xiàng)指標(biāo)。假設(shè)指標(biāo)的重要性是一樣的，包含2個(gè)子系統(tǒng)的耦合協(xié)調(diào)度D計(jì)算公式為：

4 實(shí)證分析

借助計(jì)量經(jīng)濟(jì)學(xué)方法，確定橡膠種植與人口城鎮(zhèn)化之間關(guān)系屬于正向還是負(fù)向，再探討兩者之間的耦合協(xié)調(diào)關(guān)系。

4.1 基于面板數(shù)據(jù)模型的分析

使用Stata軟件識(shí)別面板數(shù)據(jù)模型，運(yùn)用豪斯曼檢驗(yàn)來(lái)確定使用固定效應(yīng)還是隨機(jī)效應(yīng)模型。將3個(gè)階段的虛擬變量加入方程（1），以戶(hù)籍人口城鎮(zhèn)化率、非農(nóng)勞動(dòng)力占比、職工人數(shù)占比為因變量，虛擬變量的系數(shù)均顯著不為零。3 個(gè)階段的海南島各市縣橡膠種植面積占比平均分別為10.3%、13.3%、15.5%，以12%為斷點(diǎn)對(duì)前2個(gè)階段進(jìn)行回歸檢驗(yàn)，使用Stata 指令rd 估計(jì)，戶(hù)籍人口城鎮(zhèn)化率、非農(nóng)勞動(dòng)力占比的斷點(diǎn)效應(yīng)分別在5%和10%水平下通過(guò)檢驗(yàn)；以14%為斷點(diǎn)對(duì)后2個(gè)階段檢驗(yàn)，戶(hù)籍人口城鎮(zhèn)化率、非農(nóng)勞動(dòng)力占比的斷點(diǎn)效應(yīng)分別在10%和1%水平下顯著；職工人數(shù)占比均通過(guò)檢驗(yàn)。綜合看，階段劃分是必要且合理的。

表1顯示，1992—2002年橡膠種植強(qiáng)度對(duì)農(nóng)民收入的回歸系數(shù)接近零，且統(tǒng)計(jì)上不顯著，這一階段大部分橡膠園屬于國(guó)有農(nóng)場(chǎng)，民營(yíng)膠園占比小，天然橡膠價(jià)格沒(méi)有形成上漲趨勢(shì)（見(jiàn)圖3），對(duì)增加農(nóng)民收入的作用非常有限；2003—2013年橡膠種植面積擴(kuò)大可有效地增加農(nóng)民收入，從整體上看，橡膠種植占國(guó)土面積的比例提高了1個(gè)百分點(diǎn)，農(nóng)民收入約可增長(zhǎng)0.8%；2014 年后，情況發(fā)生根本性改變，擁有更多膠園的地區(qū)，當(dāng)?shù)剞r(nóng)民可支配收入反而更低，主產(chǎn)區(qū)種植戶(hù)平均橡膠收入從2012年的2.1 萬(wàn)元，下降至2018 年的1.2 萬(wàn)元（劉銳金 等，2022）。橡膠種植強(qiáng)度較高的地區(qū)，人均地方財(cái)政收入更低，即使在有農(nóng)業(yè)稅、農(nóng)林特產(chǎn)稅的第一階段和價(jià)格持續(xù)上漲的第二階段也是如此，第三階段的橡膠種植強(qiáng)度與人均地方財(cái)政收入更是顯著負(fù)相關(guān)，且影響幅度明顯增大。從回歸系數(shù)看，第二階段的系數(shù)絕對(duì)值最小，主要原因可能是持續(xù)上漲的天然橡膠價(jià)格給種植戶(hù)帶來(lái)更高的現(xiàn)金收入。根據(jù)本研究設(shè)置的橡膠種植戶(hù)固定觀(guān)察點(diǎn)數(shù)據(jù)，2011年海南主產(chǎn)區(qū)割膠日收入約是當(dāng)年城鎮(zhèn)非私營(yíng)單位在崗職工平均工資的2倍，較高的農(nóng)民收入活躍了農(nóng)村經(jīng)濟(jì)，促進(jìn)區(qū)域經(jīng)濟(jì)發(fā)展，在一定程度上帶動(dòng)地方財(cái)政收入增加。

表1 海南橡膠種植強(qiáng)度、農(nóng)民收入、財(cái)政收入與人口城鎮(zhèn)化的回歸結(jié)果Table 1 Regressions of rubber-planted intensity, farmer income, and fiscal revenue on population urbanization in Hainan

橡膠種植占國(guó)土面積的比例對(duì)人口城鎮(zhèn)化3項(xiàng)指標(biāo)的回歸系數(shù)全部是負(fù)數(shù)，但回歸系數(shù)隨時(shí)間變化呈現(xiàn)不同的特征。對(duì)非私營(yíng)單位職工占比的作用在不同階段沒(méi)有明顯的變化；對(duì)非農(nóng)就業(yè)人員占比的影響逐漸減弱，橡膠種植面積占比增加10個(gè)百分點(diǎn)，第一階段非農(nóng)就業(yè)勞動(dòng)力占比約減少10.3個(gè)百分點(diǎn)，由于青年人就業(yè)更加多元化，割膠勞動(dòng)力偏老齡化，第三階段下降至3.3 個(gè)百分點(diǎn)，且回歸系數(shù)在統(tǒng)計(jì)上不顯著；相反，對(duì)戶(hù)籍人口城鎮(zhèn)化率的回歸系數(shù)絕對(duì)變大，統(tǒng)計(jì)上也更顯著，尤其是在第三階段。從整體上看，橡膠種植面積占比較大的白沙、瓊中、屯昌、儋州等市縣的經(jīng)濟(jì)發(fā)展相對(duì)滯后，部分曾屬于國(guó)家或省級(jí)貧困縣，2020年戶(hù)籍人口城鎮(zhèn)化率僅35%左右（如圖4）。從內(nèi)在激勵(lì)上看，由于國(guó)家高度重視農(nóng)業(yè)農(nóng)村發(fā)展，村集體經(jīng)濟(jì)發(fā)展，土地制度改革深入，農(nóng)村資源資產(chǎn)價(jià)值實(shí)現(xiàn)渠道多元化，農(nóng)民“農(nóng)轉(zhuǎn)非”的意愿減弱，不少農(nóng)民外出工作但不轉(zhuǎn)戶(hù)口。這些因素可能導(dǎo)致橡膠種植面積占比對(duì)戶(hù)籍人口城鎮(zhèn)化率的估計(jì)系數(shù)中，第三階段的絕對(duì)值遠(yuǎn)大于第二階段，負(fù)向影響加強(qiáng)，但其對(duì)本區(qū)域就業(yè)結(jié)構(gòu)的作用卻明顯變?nèi)酢?/p>

農(nóng)民可支配（純）收入水平、地方人均財(cái)政收入對(duì)人口城鎮(zhèn)化相關(guān)指標(biāo)的影響在不同階段呈現(xiàn)不同特征。第一階段，橡膠種植對(duì)非私營(yíng)單位職工占比、戶(hù)籍人口城鎮(zhèn)化率的回歸系數(shù)統(tǒng)計(jì)上均不顯著，農(nóng)民收入、地方財(cái)政收入對(duì)橡膠種植與非鄉(xiāng)村從業(yè)人員占比關(guān)系的貢獻(xiàn)也接近零。第二階段，農(nóng)民收入、財(cái)政收入對(duì)橡膠種植與就業(yè)結(jié)構(gòu)的負(fù)向關(guān)系起“遮掩”效應(yīng)，主要原因可能是處于上漲通道的天然橡膠價(jià)格，主產(chǎn)區(qū)吸引更多勞動(dòng)力留在農(nóng)村割膠生產(chǎn)，在一定程度上緩解了橡膠種植對(duì)非農(nóng)就業(yè)的負(fù)相關(guān)系；該階段農(nóng)民收入、財(cái)政收入對(duì)戶(hù)籍人口城鎮(zhèn)化率的系數(shù)均顯著，其作為中介的效應(yīng)比較弱。第三階段，橡膠種植對(duì)本地區(qū)城鄉(xiāng)就業(yè)結(jié)構(gòu)的作用不明顯，從中介效應(yīng)估計(jì)結(jié)果看，橡膠種植強(qiáng)度對(duì)職工人數(shù)占比的負(fù)向關(guān)系中，農(nóng)民收入變量貢獻(xiàn)了11.64%，財(cái)政收入變量貢獻(xiàn)了12.48%；對(duì)戶(hù)籍人口城鎮(zhèn)化率的負(fù)相關(guān)系中，農(nóng)民收入和財(cái)政收入分別貢獻(xiàn)了6.71%和17.04%。使用常住人口城鎮(zhèn)化率替代2014—2020年戶(hù)籍人口城鎮(zhèn)化率，回歸系數(shù)為-0.542 3，在1%水平下顯著，絕對(duì)值相較于戶(hù)籍人口城鎮(zhèn)化率明顯變小，再次表明橡膠種植對(duì)區(qū)域內(nèi)就業(yè)多元化的影響減弱。

橡膠種植面積占比對(duì)人口城鎮(zhèn)化的影響具有很強(qiáng)的延續(xù)性，故將橡膠種植面積占比滯后一期作為工具變量，使用兩階段最小二乘法（2SLS）重新估計(jì)。以戶(hù)籍人口城鎮(zhèn)化率為例，3 個(gè)階段的豪斯曼檢驗(yàn)結(jié)果表明，沒(méi)有工具變量的估計(jì)方程的統(tǒng)計(jì)量是顯著的；忽略時(shí)間維度，將面板數(shù)據(jù)轉(zhuǎn)化為截面數(shù)據(jù)，DWH 檢驗(yàn)結(jié)果均不顯著，即無(wú)法拒絕變量外生性的原假設(shè)，可以排除內(nèi)生性問(wèn)題的影響。

4.2 基于空間計(jì)量模型的進(jìn)一步檢驗(yàn)

橡膠種植分布存在空間關(guān)聯(lián)性，采用空間計(jì)量模型進(jìn)一步檢驗(yàn)。橡膠種植分布具有較強(qiáng)的空間限制特征，如自然災(zāi)害頻發(fā)區(qū)域的農(nóng)戶(hù)和農(nóng)場(chǎng)往往會(huì)減少橡膠種植，而人口城鎮(zhèn)化與經(jīng)濟(jì)社會(huì)發(fā)展緊密相關(guān)，因此設(shè)計(jì)了相鄰矩陣、地理距離矩陣以及經(jīng)濟(jì)矩陣分析橡膠種植空間分布特征。相鄰矩陣的對(duì)角線(xiàn)元素為零，若兩個(gè)市縣行政區(qū)域邊界上相鄰設(shè)為1，否則為0，與瓊中縣相鄰的市縣達(dá)9個(gè)，即海南島內(nèi)一半市縣均與其相鄰。地理距離矩陣使用百度地圖測(cè)距工具獲得，采用直接距離數(shù)據(jù)表示，對(duì)角線(xiàn)為0，其他元素是距離的倒數(shù)。經(jīng)濟(jì)矩陣設(shè)置參照林光平等（2005）的做法，使用相鄰矩陣與人均GDP 差距絕對(duì)值倒數(shù)矩陣的乘積得到，人均GDP 是可比價(jià)格下1992—2002、2003—2020 年的平均值，經(jīng)濟(jì)矩陣不是對(duì)稱(chēng)矩陣，其中1992—2002年的矩陣包含瓊山市。運(yùn)用莫蘭指數(shù)I（Moran'sI）、吉爾里指數(shù)c（Geary'sc），對(duì)2000、2010和2020年的數(shù)據(jù)進(jìn)行空間相關(guān)性檢驗(yàn)，從表2可知，人均地方財(cái)政收入、戶(hù)籍人口城鎮(zhèn)化率、非鄉(xiāng)村從業(yè)人員占比在不同矩陣及方法下，均統(tǒng)計(jì)不顯著，其余變量都存在一定程度的空間相關(guān)性。橡膠種植面積占比具有明顯的空間相關(guān)性，橡膠種植從海南島中間由西往東的儋州、白沙、瓊中、瓊海等市縣，向南、向北逐步減少。2000年及以前，橡膠種植以國(guó)有農(nóng)場(chǎng)為主，分布具有一定的局限性；而后逐漸轉(zhuǎn)為以民營(yíng)膠園逐漸增加，相鄰矩陣下的空間相關(guān)性有所增強(qiáng)。

表2 個(gè)別年份空間相關(guān)性的檢驗(yàn)結(jié)果Table2 Tested results of spatial correlations in some years

采用空間杜賓模型，即方程（4），重新檢驗(yàn)橡膠種植面積占比對(duì)地區(qū)就業(yè)結(jié)構(gòu)、人口城鎮(zhèn)化率的影響，通過(guò)豪斯曼檢驗(yàn)確定固定和隨機(jī)效應(yīng)模型。從ρ的估計(jì)值及其P值可知，人口城鎮(zhèn)化的指標(biāo)在大部分情形下存在較為明顯的正向空間自相關(guān)（表3）。

表3 橡膠種植強(qiáng)度對(duì)海南人口城鎮(zhèn)化的空間計(jì)量模型估計(jì)結(jié)果Table3 Estimated spatial econometric models of rubber-planted intensity on urbanization in Hainan

橡膠種植強(qiáng)度對(duì)人口城鎮(zhèn)化3個(gè)指標(biāo)的直接效應(yīng)均是負(fù)向的，全省天然橡膠產(chǎn)業(yè)集中在種植和初加工環(huán)節(jié)，對(duì)非農(nóng)就業(yè)的貢獻(xiàn)有限，與人口城鎮(zhèn)化指標(biāo)呈負(fù)向關(guān)系，與表1的估計(jì)一致，但橡膠種植強(qiáng)度對(duì)人口城鎮(zhèn)化指標(biāo)具有一定的正向空間溢出效應(yīng)。在第三階段的相鄰矩陣、經(jīng)濟(jì)矩陣下，橡膠種植強(qiáng)度對(duì)戶(hù)籍、常住人口城鎮(zhèn)化率的間接效用和總效用估計(jì)值均為正，且對(duì)常住人口城鎮(zhèn)化率的方程還通過(guò)了統(tǒng)計(jì)顯著性檢驗(yàn)。在第二階段的相鄰矩陣、地理距離矩陣下，橡膠種植強(qiáng)度對(duì)戶(hù)籍人口城鎮(zhèn)化率有正向的間接效應(yīng)、總效應(yīng)，并在5%顯著性水平通過(guò)檢。第三階段的3個(gè)距離矩陣下，橡膠種植強(qiáng)度對(duì)非農(nóng)就業(yè)人員占比均有正向溢出效應(yīng)，第二階段的相鄰矩陣下，間接效應(yīng)同樣是正向的。這表明若周邊市縣橡膠種植面積占比增大，可能有利于吸納更多勞動(dòng)力就業(yè)，在一定程度上使得本地區(qū)人口城鎮(zhèn)化水平相對(duì)高。橡膠種植總體上從海南島中部地區(qū)向四周延伸，人口城鎮(zhèn)化率較高的地區(qū)幾乎都分布在沿海市縣，2020年戶(hù)籍人口城鎮(zhèn)化率＞65%、常住人口城鎮(zhèn)化率＞70%的?？谑泻腿齺喪校约皯?hù)籍人口城鎮(zhèn)化率＞40%（圖4）的文昌市、萬(wàn)寧市、澄邁縣、東方市、昌江縣皆屬于沿海地區(qū)，而上述地區(qū)周邊市縣的橡膠種植強(qiáng)度總體相對(duì)高。其中一個(gè)原因或?yàn)橄鹉z種植強(qiáng)度較大的地區(qū)，部分割膠勞動(dòng)力放棄割膠，人口向周邊流出后，成為相鄰市縣的城鎮(zhèn)常住人口、城鎮(zhèn)從業(yè)人員，進(jìn)而提升流入地區(qū)的人口城鎮(zhèn)化水平。

4.3 耦合協(xié)調(diào)關(guān)系分析

計(jì)量結(jié)果顯示，橡膠種植占本地國(guó)土面積的比例對(duì)人口城鎮(zhèn)化的作用是逆向的，因此將橡膠種植強(qiáng)度按照負(fù)向指標(biāo)進(jìn)行標(biāo)準(zhǔn)化處理。通過(guò)式（5）分別計(jì)算橡膠種植面積占比與非私營(yíng)單位職工占比、非鄉(xiāng)村從業(yè)人員占比、戶(hù)籍人口城鎮(zhèn)化率的耦合協(xié)調(diào)度，各市縣3個(gè)指標(biāo)之間的耦合協(xié)調(diào)度沒(méi)有明顯差異，選擇橡膠種植面積占比與戶(hù)籍人口城鎮(zhèn)化率進(jìn)行分析（圖6）。海口市和三亞市的耦合協(xié)調(diào)度分別達(dá)到0.96和0.87，橡膠種植強(qiáng)度低，經(jīng)濟(jì)發(fā)展水平較高；南部地區(qū)五指山市、昌江縣、東方市、保亭縣的耦合協(xié)調(diào)度＞0.6。相較于1992—2002年，2003—2020年平均耦合協(xié)調(diào)度均有不同程度的提高，提高幅度從東往西遞減；儋州市是海南島最大的產(chǎn)區(qū)，耦合協(xié)調(diào)度幾乎沒(méi)有變化，橡膠種植面積占比和戶(hù)籍人口城鎮(zhèn)化率分別提高7.5 個(gè)百分點(diǎn)和13.5個(gè)百分點(diǎn)；而萬(wàn)寧市、瓊海市在天然橡膠種植規(guī)模穩(wěn)定的同時(shí)，人口城鎮(zhèn)化率分別增加了16.4個(gè)百分點(diǎn)和12.5個(gè)百分點(diǎn)?？傮w而言，人口城鎮(zhèn)化率提高幅度較大、橡膠種植占比相對(duì)穩(wěn)定或有所降低的市縣，2 個(gè)變量間的耦合協(xié)調(diào)度提升相對(duì)較快（見(jiàn)圖6-a、b）。海南島橡膠種植業(yè)已進(jìn)入穩(wěn)定發(fā)展階段，可重點(diǎn)依托該產(chǎn)業(yè)延伸產(chǎn)業(yè)鏈，為農(nóng)村勞動(dòng)力提供更多非農(nóng)就業(yè)機(jī)會(huì)，進(jìn)而提高人口城鎮(zhèn)化率。

圖6 海南島各市縣橡膠種植強(qiáng)度（a）及其與人口城鎮(zhèn)化率的耦合協(xié)調(diào)度（b）分布Fig.6 Distribution of rubber-planted intensity(a) and coupling coordination between it and population urbanization ratios(b) in the Hainan Island

5 結(jié)論與討論

通過(guò)構(gòu)建中介效應(yīng)、空間計(jì)量模型、耦合協(xié)調(diào)度模型，探討海南橡膠種植強(qiáng)度對(duì)城鄉(xiāng)就業(yè)與人口結(jié)構(gòu)的影響，主要結(jié)論如下：

1）橡膠種植業(yè)曾對(duì)農(nóng)民增收起積極作用，2003—2013 年，較好地促進(jìn)主產(chǎn)區(qū)農(nóng)村經(jīng)濟(jì)繁榮，但與人均地方財(cái)政收入的關(guān)系卻始終是負(fù)向的，橡膠主產(chǎn)區(qū)財(cái)政收入偏低，即使有農(nóng)業(yè)稅存在，也是如此。

2）橡膠種植對(duì)非私營(yíng)單位職工占比、非鄉(xiāng)村從業(yè)人員占比、戶(hù)籍和常住人口城鎮(zhèn)化率均為負(fù)向關(guān)聯(lián)，面板數(shù)據(jù)模型和空間杜賓模型的直接效應(yīng)估計(jì)結(jié)果都證明了這點(diǎn)，既有各市縣區(qū)域經(jīng)濟(jì)發(fā)展戰(zhàn)略不同的原因，也有主產(chǎn)市縣自身發(fā)展動(dòng)力不足的影響。

3）農(nóng)民收入、地方財(cái)政收入作為中介變量在不同時(shí)期發(fā)揮的作用不同，在天然橡膠價(jià)格持續(xù)上漲期間，對(duì)人口城鎮(zhèn)化率是正向貢獻(xiàn)。

4）橡膠種植對(duì)周邊地區(qū)具有一定的正向溢出效應(yīng)，橡膠種植強(qiáng)度大的市縣對(duì)周邊地區(qū)促進(jìn)非農(nóng)就業(yè)、提高人口城鎮(zhèn)化率都具有一定的促進(jìn)作用，橡膠面積占比與人口城鎮(zhèn)化的耦合協(xié)調(diào)度總體在提高。

促進(jìn)天然橡膠更好地服務(wù)地區(qū)經(jīng)濟(jì)發(fā)展，尤其是通過(guò)促進(jìn)地區(qū)非農(nóng)就業(yè)帶動(dòng)人口城鎮(zhèn)化，存在2個(gè)關(guān)鍵點(diǎn)。1）提高勞動(dòng)生產(chǎn)率，采用低頻割膠技術(shù)或更先進(jìn)的割膠工具，大幅度提高每個(gè)割膠工人負(fù)責(zé)的株數(shù)，將更多勞動(dòng)力從割膠生產(chǎn)中解放出來(lái)，同時(shí)提供更多的非農(nóng)就業(yè)機(jī)會(huì)。2）延長(zhǎng)產(chǎn)業(yè)鏈，天然橡膠作為重要工業(yè)原料，后端產(chǎn)業(yè)鏈長(zhǎng)，基于原料的優(yōu)勢(shì)發(fā)展下游產(chǎn)業(yè)，發(fā)展天然橡膠新材料和制品業(yè)，不僅能創(chuàng)造非農(nóng)就業(yè)機(jī)會(huì)，還能增加地方財(cái)政收入，將資源優(yōu)勢(shì)轉(zhuǎn)化為經(jīng)濟(jì)優(yōu)勢(shì)，促使更多勞動(dòng)力在城鎮(zhèn)安家立業(yè)，提高人口城鎮(zhèn)化率。以馬來(lái)西亞為例，早期是一個(gè)天然橡膠原料出口國(guó)，種植端占絕對(duì)主體，但目前轉(zhuǎn)型為天然橡膠原料、橡膠制品、橡膠木三足鼎立的發(fā)展格局，附加值提升，非農(nóng)就業(yè)崗位增加，這對(duì)海南島天然橡膠產(chǎn)業(yè)轉(zhuǎn)型具有明顯借鑒意義。但需要推動(dòng)差異化發(fā)展，不同市縣發(fā)展定位不同，并不是都要致力于延伸產(chǎn)業(yè)鏈，白沙、瓊中、五指山等市縣處于海南島生態(tài)核心區(qū)內(nèi)，制造業(yè)不是其發(fā)展重點(diǎn)，需從全島同城化框架下協(xié)調(diào)推進(jìn)發(fā)展橡膠新材料和制品制造業(yè)發(fā)展。本文解釋了橡膠種植強(qiáng)度與人口城鎮(zhèn)化水平的多維度關(guān)系，發(fā)現(xiàn)橡膠種植面積占比對(duì)本區(qū)域推動(dòng)人口城鎮(zhèn)化具有抑制作用，對(duì)平衡上級(jí)政府穩(wěn)定面積和本區(qū)域經(jīng)濟(jì)社會(huì)發(fā)展的目標(biāo)沖突提出見(jiàn)解，為完善天然橡膠產(chǎn)業(yè)扶持政策提供參考，如天然橡膠主產(chǎn)區(qū)利益補(bǔ)償機(jī)制。本研究重點(diǎn)關(guān)注了人口城鎮(zhèn)化，但實(shí)際上處于海南島國(guó)家熱帶雨林公園區(qū)域內(nèi)和周邊的市縣，發(fā)展優(yōu)先目標(biāo)重點(diǎn)在于生態(tài)保護(hù)，其經(jīng)濟(jì)社會(huì)發(fā)展和橡膠產(chǎn)業(yè)也應(yīng)服務(wù)于此，本文未對(duì)差異化目標(biāo)進(jìn)行分析，未來(lái)將考慮區(qū)域定位差異，設(shè)定不同發(fā)展目標(biāo)，深入探討橡膠種植業(yè)與發(fā)展目標(biāo)的協(xié)調(diào)性，有助于提出更適合特定地區(qū)的對(duì)策措施。