余慧容，杜鵬飛

（1.清華大學(xué)環(huán)境學(xué)院， 北京100084；2.環(huán)境模擬與污染控制國家重點(diǎn)實(shí)驗(yàn)室，北京 100084）

1 引言

耕地資源是保障食物安全的核心要素，直接關(guān)系到植物性食物的生產(chǎn)及畜禽水產(chǎn)飼養(yǎng)所需飼料的供給[1-2]。隨著社會經(jīng)濟(jì)的快速發(fā)展，耕地資源的可持續(xù)利用面臨著前所未有的壓力：一方面來自生產(chǎn)端耕地的被占用與功能消退，另一方面源于消費(fèi)端的人口食物需求總量增加及居民食物消費(fèi)結(jié)構(gòu)變化?，F(xiàn)有關(guān)耕地資源可持續(xù)利用方面的研究大多聚焦在生產(chǎn)端[3-5]，而針對消費(fèi)端的研究則相對較少。隨著中國城鄉(xiāng)居民對食物的需求日益增長且日趨多樣化，消費(fèi)端對耕地資源的影響日益凸顯。因此，非常有必要從消費(fèi)端探究膳食對耕地資源的影響。在中國，城鄉(xiāng)發(fā)展不平衡不協(xié)調(diào)是制約當(dāng)前經(jīng)濟(jì)社會發(fā)展的最突出矛盾[6]，在消費(fèi)端尤為明顯[7]。中共十九大報(bào)告明確提出“建立健全城鄉(xiāng)融合發(fā)展體制機(jī)制和政策體系”。故此，在中國重塑城鄉(xiāng)關(guān)系的大背景下，亟需厘清基于消費(fèi)端的城鄉(xiāng)居民膳食需求差異及其與耕地資源的相互關(guān)系，這對于中國耕地保護(hù)、食物安全、城鄉(xiāng)融合、國民健康乃至國家健康可持續(xù)發(fā)展均具有重要意義。

現(xiàn)階段，已有不少學(xué)者對膳食與耕地資源之間的關(guān)系進(jìn)行研究。國際上對膳食與耕地資源之間相互關(guān)系的研究興起于21世紀(jì)初，GERBENS-LEENES等于2002年利用統(tǒng)計(jì)數(shù)據(jù)和家庭調(diào)查數(shù)據(jù)，詳細(xì)核算荷蘭1990年100多種食物在其生產(chǎn)過程中所需的土地面積[8]。該方法對數(shù)據(jù)要求高、數(shù)據(jù)收集耗時多、運(yùn)算量大，難以運(yùn)用于長時間尺度的研究。為了克服此缺陷，KASTNER和NONHEBEL在GERBENSLEENES等研究的基礎(chǔ)上，引入等效農(nóng)作物的概念，進(jìn)一步簡化計(jì)算方法與計(jì)算流程，并借助FAO的食物平衡表計(jì)算分析了菲律賓1919—2003年長時間尺度的耕地需求變化[9]。此后，該方法在全世界范圍內(nèi)得到廣泛的應(yīng)用[10-13]。國內(nèi)對膳食與耕地資源的關(guān)系研究起步相對較晚，2010年ZHEN等首次將GERBENS-LEENES等的方法運(yùn)用于中國西部的寧夏固原地區(qū)，借助實(shí)地調(diào)查數(shù)據(jù)計(jì)算2006年當(dāng)?shù)剞r(nóng)村居民基于食物消費(fèi)的耕地需求量[14]。由于該方法對數(shù)據(jù)的高要求，已有文獻(xiàn)僅將其運(yùn)用于局域尺度特定人群在特定年限的膳食特征與耕地資源關(guān)系研究，如高利偉將其運(yùn)用于2010年西藏“一江兩河”流域的農(nóng)村居民[15]，LI等則針對旅游人員這一特殊群體，用GERBENS-LEENES等人的方法研究2013—2015年拉薩旅游人員的膳食耕地需求量[16]。而KASTNER和NONHEBEL的改進(jìn)方法則往往運(yùn)用于長時間尺度和較大空間尺度的研究：在省級空間尺度上，王佳友等運(yùn)用該改進(jìn)方法借助統(tǒng)計(jì)年鑒數(shù)據(jù)對廣東省居民自20世紀(jì)80年代以來食物消費(fèi)結(jié)構(gòu)對耕地需求的影響進(jìn)行分析[17]。在全國尺度上，楊旺明等[18]、劉春霞與王芳[19]、張翠玲等[20]均運(yùn)用KASTNER和NONHEBEL的改進(jìn)方法，假設(shè)市場上的供給量和需求量相等，借助FAO中國食物均衡表（Food Balance Sheet）中的食物總供給量數(shù)據(jù)，分別從飲食影響因素、居民食物消費(fèi)模式、食物消費(fèi)結(jié)構(gòu)優(yōu)化等角度研究近三四十年中國居民總體膳食特征與耕地資源之間的相互關(guān)系。FAO中國食物均衡表來源于聯(lián)合國糧農(nóng)組織網(wǎng)站①參見https://www.fao.org/faostat/en/#data。,記錄著中國自20世紀(jì)60年代以來的各類食物總供給量、總生產(chǎn)量等信息，將其用于膳食與耕地資源關(guān)系分析具有數(shù)據(jù)完整、易獲取等優(yōu)點(diǎn)，缺點(diǎn)在于：（1）其將食物總供給量作為消費(fèi)量進(jìn)行分析，供給量與消費(fèi)量相等的前提條件是市場完全競爭，這與不完全競爭市場的現(xiàn)實(shí)情形有所偏差；（2）食物總供給量無法體現(xiàn)城鄉(xiāng)居民的食物消費(fèi)差異。而長時間尺度的中國城鄉(xiāng)居民的食物消費(fèi)數(shù)據(jù)需通過收集整理匯總國內(nèi)不同年代的統(tǒng)計(jì)資料而得。故此，本文通過中國統(tǒng)計(jì)資料館、圖書館、中國經(jīng)濟(jì)社會大數(shù)據(jù)研究平臺、EPS全球統(tǒng)計(jì)數(shù)據(jù)/分析平臺、布瑞克農(nóng)業(yè)數(shù)據(jù)平臺等多個渠道收集整理20世紀(jì)80年代以來的中國城鄉(xiāng)居民食物消費(fèi)數(shù)據(jù)和農(nóng)產(chǎn)品生產(chǎn)數(shù)據(jù)，并在此基礎(chǔ)上對城鄉(xiāng)居民膳食結(jié)構(gòu)與營養(yǎng)變遷對耕地資源的影響進(jìn)行長時間尺度的分析，并結(jié)合最新《中國居民膳食指南》提出的不同能量水平基本理想膳食模型分析未來城鄉(xiāng)膳食變遷方向及其對耕地資源的影響，旨在為城鄉(xiāng)融合、國民健康及耕地保護(hù)決策提供相應(yīng)的科學(xué)支撐。

2 研究方法

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

城鄉(xiāng)居民膳食數(shù)據(jù)來源于中國統(tǒng)計(jì)局，是由歷年的《中國統(tǒng)計(jì)年鑒》（1981—2020年）、《中國住戶調(diào)查年鑒》（2011—2020年）、《中國城市（鎮(zhèn)）生活與價格年鑒》（2006—2012年）、《中國價格及城鎮(zhèn)居民家庭收支調(diào)查統(tǒng)計(jì)年鑒》（1996—2005年）、《中國城鎮(zhèn)居民家庭收支調(diào)查資料》（1987—1991年、1994年）、《“六五”期間我國城鎮(zhèn)居民家庭收支調(diào)查資料》、《農(nóng)村住戶調(diào)查資料提要1954—1988年》、《農(nóng)村住戶調(diào)查資料》（1990—1991年）、《中國農(nóng)村住戶調(diào)查年鑒》（1992年、2000—2010年）及各省的統(tǒng)計(jì)年鑒綜合整理而得，時間尺度為1981—2019年。其中，2013年以前城鎮(zhèn)的糧食消費(fèi)量是按加工糧統(tǒng)計(jì)，本文中依據(jù)相關(guān)文獻(xiàn)研究成果[21-22]，按0.75的轉(zhuǎn)換系數(shù)將其折算為原糧，糧食消費(fèi)量統(tǒng)一用原糧計(jì)算。由于國家統(tǒng)計(jì)局的食物消費(fèi)數(shù)據(jù)主要是指城鄉(xiāng)居民的居家消費(fèi)量，并未考慮外出就餐及其他食物來源的消費(fèi)量[18]，本文綜合《2021中國與全球食物政策報(bào)告》[23]和蔡祖聰[24]等研究成果中的城鄉(xiāng)居民在外消費(fèi)比例系數(shù)對城鄉(xiāng)居民的人均食物消費(fèi)量進(jìn)行修正，獲得最終的城鄉(xiāng)居民膳食數(shù)據(jù)。

此外，油料作物中菜籽油與大豆油的總產(chǎn)量與總壓榨量、糖料作物中的甘蔗與甜菜總產(chǎn)量與總含糖量均來源于布瑞克數(shù)據(jù)庫②參見http://www.agdata.cn/。，動物性食物的耗糧量來源于歷年的《全國農(nóng)產(chǎn)品成本收益年鑒》，肉蛋奶產(chǎn)量、各農(nóng)作物播種面積、耕地面積來自EPS數(shù)據(jù)庫③參見https://www.epsnet.com.cn。及相關(guān)文獻(xiàn)[25-26]。其中，部分年份缺失的數(shù)據(jù)按其最近的兩個年份均值進(jìn)行估算。

2.2 研究方法

2.2.1 膳食營養(yǎng)分析

本文選取能量以及蛋白質(zhì)、脂肪、膽固醇、膳食纖維及碳水化合物5種重要的宏量營養(yǎng)素來對城鄉(xiāng)居民的膳食營養(yǎng)進(jìn)行分析。城鄉(xiāng)居民的每日膳食營養(yǎng)計(jì)算公式如下：

式（1）中：Nutrienturban or rural為城鎮(zhèn)居民（或農(nóng)村居民）的每日膳食營養(yǎng)（kcal或g）；Fi為人均膳食量，即人均第i類食物的消費(fèi)量（kg/人）；ni，j為第i類食物的第j種營養(yǎng)素代表值，其取值基于《中國食物成分表》[27]。j為能量與主要的宏觀營養(yǎng)素（包括蛋白質(zhì)、脂肪、膽固醇、膳食纖維及碳水化合物）。i為膳食結(jié)構(gòu)中的食物種類，文中涉及的主要食物種類有11類，分別為糧食、蔬菜、水果、植物油、食糖、豬肉、牛羊肉、禽肉、水產(chǎn)品、蛋類和奶類。

2.2.2 膳食耕地足跡

本文基于KASTNER和NONHEBEL的耕地需求核算方法[9]，借助膳食耕地足跡來分析膳食變遷對耕地資源的影響。其中，膳食耕地足跡是指膳食結(jié)構(gòu)中各類食物在生產(chǎn)過程中所占用的耕地面積[9,28]，其計(jì)算流程如圖1所示。

圖1 居民膳食與耕地資源的轉(zhuǎn)換關(guān)系Fig.1 Conversion from diet to cultivated land footprint

具體的計(jì)算公式為：

式（2）—式（5）中：L為人均膳食耕地足跡（m2/人）；LTotal為總膳食足跡（m2），即人均膳食足跡（L）與人口（Pop）的乘積；Ai為人均所需食物等效農(nóng)作物播種面積（m2/人）；I為復(fù)種指數(shù)，為食物等效農(nóng)作物總播種面積與對應(yīng)的耕地面積（總耕地面積扣除棉麻煙的種植面積）的比率[29]；Pi為人均食物等效農(nóng)作物需求量（kg/人）；Yi為食物種類i對應(yīng)等效農(nóng)作物的單位面積產(chǎn)量（kg/ m2）；Fi為人均膳食量，即人均第i類食物的消費(fèi)量（kg/人）。在本文所涉及的11類食物種類中，水果作物包括瓜果和園林水果，按照土地利用分類的最新官方文件《國土空間調(diào)查、規(guī)劃、用途管制用地用海分類指南（試行）》 中耕地的定義，園林水果屬于園地出產(chǎn)，故在計(jì)算膳食耕地足跡時僅考慮瓜果部分。瓜果的消費(fèi)量缺乏相應(yīng)統(tǒng)計(jì)數(shù)據(jù)，其數(shù)值由瓜果產(chǎn)量占水果總產(chǎn)量的比例進(jìn)行估算。Ci為食物種類與對應(yīng)等效農(nóng)作物類型之間的轉(zhuǎn)換系數(shù)。糧食、蔬菜、瓜果、植物油、食糖及動物性食物所對應(yīng)的等效農(nóng)作物分別為糧食作物、蔬菜、瓜果、油料作物、糖料作物和用于飼料的糧食作物。這里的轉(zhuǎn)換系數(shù)主要反映食物的加工情況。其中，糧食（以“原糧”計(jì)算）、蔬菜和瓜果由其等效農(nóng)作物直接產(chǎn)出，故糧食、蔬菜類及水果類的轉(zhuǎn)換系數(shù)為1。植物油與油料作物的轉(zhuǎn)換系數(shù)取決于油料作物的出油率，數(shù)值上等于出油率的倒數(shù)。本文以菜籽油和大豆油為代表（因菜籽油和大豆油是我國最主要的兩種食用植物油，二者消費(fèi)份額占植物油總消費(fèi)份額的一半以上），根據(jù)每年菜籽油與大豆油的總產(chǎn)量與總壓榨量之比計(jì)算出每年相應(yīng)的出油率。食糖與糖料作物的轉(zhuǎn)換系數(shù)為糖料作物出糖率的倒數(shù)，數(shù)值上等于甘蔗與甜菜的總產(chǎn)量與總含糖量的比值。動物性食物與飼料糧作物的轉(zhuǎn)換系數(shù)取決于豬肉、牛羊肉、禽肉、水產(chǎn)品、蛋類和奶類各自的耗糧系數(shù)。耗糧系數(shù)為各類動物性食物的耗糧量（耗用的各種飼料折成糧食的數(shù)量）與凈肉（或蛋或奶）產(chǎn)量之比。由于統(tǒng)計(jì)年鑒上的肉類產(chǎn)量以酮體（禽畜、水產(chǎn)等屠宰后的含骨頭的軀干部分）重量為準(zhǔn)，需將其乘上凈肉系數(shù)以獲得凈肉產(chǎn)量（酮體去骨頭后的肉量）。其中，豬肉、牛肉、羊肉、禽肉、水產(chǎn)品、蛋類和奶類的凈肉系數(shù)分別為0.83、0.80、0.70、0.75、0.55、0.95和1[30]。本文中涉及糧食部分（包括飼料用糧）均統(tǒng)一折算為“原糧”。

3 結(jié)果與分析

3.1 城鄉(xiāng)居民膳食變遷特征

3.1.1 膳食結(jié)構(gòu)變遷

由圖2可知，1981—2019年中國城鄉(xiāng)居民的人均膳食結(jié)構(gòu)呈現(xiàn)出如下特征：（1）城鎮(zhèn)居民的動物性食物比例均高于農(nóng)村居民；（2）食物日趨均衡化，城鄉(xiāng)居民膳食結(jié)構(gòu)差異在縮小。城鄉(xiāng)居民的植物性食物比例均日趨減少、動物性食物比例均日趨上升。其中，城鎮(zhèn)居民的動物性食物比例由20世紀(jì)80年代初的10.1%上升至2019年的23.9%，農(nóng)村居民的動物性食物比例由20世紀(jì)80年代初的3.2%上升至2019年的19.3%。這與我國城鄉(xiāng)居民收入增長、膳食消費(fèi)水平提升密切相關(guān)。

圖2 城鄉(xiāng)居民人均膳食結(jié)構(gòu)變遷Fig.2 Changes in the per capita dietary structure of urban and rural residents

3.1.2 膳食能量與營養(yǎng)變遷

1981—2019 年，城鄉(xiāng)居民的人均每日膳食營養(yǎng)整體呈現(xiàn)出如下特征：在膳食能量方面（圖3（a）），農(nóng)村居民高于城鎮(zhèn)居民，城鎮(zhèn)居民膳食能量有所上升，農(nóng)村居民所有波動且近幾年呈上升趨勢；在膳食蛋白質(zhì)方面（圖3（b）），城鄉(xiāng)居民均有所上升，城鎮(zhèn)居民的膳食蛋白質(zhì)增速高于農(nóng)村居民；在膳食脂肪與膽固醇方面（圖3（c）和圖3（d）），城鄉(xiāng)居民均呈現(xiàn)出較大幅度的增長；在膳食纖維方面（圖3（e）），城鄉(xiāng)居民均呈現(xiàn)出先減少后增加的趨勢；在碳水化合物方面（圖3（f）），城鄉(xiāng)居民均有所減少，農(nóng)村居民的減少幅度高于城鎮(zhèn)居民。這些主要受城鄉(xiāng)居民膳食結(jié)構(gòu)變動的影響：依據(jù)《中國食物成分表》，動物性食物的蛋白含量低于豆類但高于糧食、水果、蔬菜，脂類主要存在于油脂及動物性食物中，而膳食纖維和碳水化合物主要存在于植物性食物中；故此，隨著城鄉(xiāng)居民動物性食物食用比例上升、植物性食物食用比例下降，其動物性蛋白攝入量增加，脂肪與膽固醇攝入量則相應(yīng)增高，而膳食纖維和碳水化合物的攝入量則相應(yīng)的減少。這些膳食營養(yǎng)在2013年有波動，這是由于受城鄉(xiāng)住戶調(diào)查一體化改革影響，國家統(tǒng)計(jì)局的統(tǒng)計(jì)指標(biāo)與統(tǒng)計(jì)范圍有所變化而致[31]，但該波動對整體營養(yǎng)變遷趨勢研究影響不大。

圖3 城鄉(xiāng)居民人均每日膳食營養(yǎng)變遷Fig.3 Changes in the per capita dietary nutrition of urban and rural residents

3.2 城鄉(xiāng)膳食變遷對耕地資源的影響

從城鄉(xiāng)人均膳食耕地足跡來看（圖（4a）—圖（4b）），城鎮(zhèn)居民的整體人均膳食耕地足跡大于農(nóng)村居民。1981—2019年，城鄉(xiāng)居民的人均膳食耕地足跡均呈現(xiàn)出先下降后上升的趨勢。前期城鄉(xiāng)居民人均膳食耕地足跡的減少主要受中國食物等效農(nóng)作物生產(chǎn)水平提高的影響；而后期城鄉(xiāng)居民人均膳食耕地足跡的增加主要受膳食結(jié)構(gòu)變化的影響，城鄉(xiāng)居民人均動物性食物攝入量較大幅度增加，人均動物性食物耕地足跡增加的幅度大于其植物性食物足跡減少的幅度。其中，2000—2003年、2009年出現(xiàn)斷層式增長，前者是由于國家農(nóng)業(yè)政策調(diào)整、農(nóng)民生產(chǎn)積極性下降及自然災(zāi)害比較嚴(yán)重等綜合因素導(dǎo)致此階段國家糧食單產(chǎn)減少、耕地足跡增加，后者是源于第二次全國土地調(diào)查的耕地面積大于詳查數(shù)據(jù)而致。

圖4 城鄉(xiāng)膳食耕地足跡演變特征Fig.4 Changes in the dietary cultivated land footprints of urban and rural residents

從城鄉(xiāng)總膳食耕地足跡來看（圖（4c）—圖（4e）），1981—2019年城鄉(xiāng)居民的總膳食耕地足跡波動中有所上升，由1981年的12 700億m2增加到2019年的12 493億m2。其中，城鎮(zhèn)居民的總膳食耕地足跡由1981年的2 821億m2增加到2019年的7 889億m2，而農(nóng)村居民的總膳食耕地足跡則從1981年的9 879億m2減少到2019年的4 605億m2。這主要受到城鄉(xiāng)人口變動的影響，隨著我國城鎮(zhèn)化的推進(jìn)，特別是1995年之后城鎮(zhèn)人口急劇上升、農(nóng)村人口急劇減少，城鎮(zhèn)人口從1981年的17 245萬人增加至2019年的84 843萬人，農(nóng)村人口由79 901萬人減少至55 162萬人。

超聲造影是通過向體內(nèi)注入造影劑來提高超聲顯像效果的技術(shù)，可以較好地反映組織的血流灌注情況，是彩色多普勒成像之后超聲醫(yī)學(xué)的再一次革命，具有里程碑式的意義。與常規(guī)超聲相比，超聲造影可以檢出常規(guī)超聲顯示不清的微弱血流信息和微小病灶，葉桂宏等[4]研究發(fā)現(xiàn)，超聲造影與增強(qiáng)CT掃描在肝癌病灶檢出率方面無差異，且超聲造影對＜1cm的肝癌病灶的診斷比增強(qiáng)CT更有優(yōu)勢。

3.3 城鄉(xiāng)膳食變遷趨勢與耕地資源保護(hù)

本文基于最新《中國居民膳食指南》所提出的理想膳食模型及城鄉(xiāng)融合的發(fā)展方向，選取其低、中、高三種能量膳食模型作為未來的城鄉(xiāng)健康營養(yǎng)膳食趨勢目標(biāo)，并結(jié)合《國家人口發(fā)展規(guī)劃（2016—2030年）》中預(yù)測2030年人口14.5億人、城鎮(zhèn)化70%等指標(biāo)測算未來三種理想膳食模型的耕地足跡，并將其與2019年城鄉(xiāng)居民膳食現(xiàn)狀進(jìn)行對比，所得的結(jié)果如圖5、表1所示。

圖5 城鄉(xiāng)居民膳食現(xiàn)狀及未來理想膳食模型Fig.5 Current situation and future ideal dietary patterns of urban and rural residents

表1 城鄉(xiāng)居民膳食變遷趨勢及其對耕地影響Tab.1 Dietary changing trends of urban and rural residents and the impacts on cultivated land

在膳食結(jié)構(gòu)方面（圖5、表1），城鄉(xiāng)居民的動物性食物比例現(xiàn)狀與理想膳食模型相比仍然偏低，現(xiàn)階段城鎮(zhèn)居民與農(nóng)村居民的動物性食物比例分別為24%和19%，而低、中、高能量三種理想膳食模型的動物性食物比例均大于25%，分別為34%、30%和27%。按照膳食指南，可以預(yù)見未來城鄉(xiāng)居民的動物性食物膳食比例將進(jìn)一步增加。此外，城鄉(xiāng)居民的食物消費(fèi)品種搭配有待進(jìn)一步優(yōu)化：在植物性食物中，現(xiàn)狀城鄉(xiāng)居民的水果食物比例為15%和11%，低于理想膳食模型16%～20%；在動物性食物中，現(xiàn)狀城鄉(xiāng)居民的畜肉食物消費(fèi)比例偏高而奶類則偏低，其中畜肉食物7%～8%的比例高于理想膳食模型的2%～3%，奶類食物2%～4%的比例則遠(yuǎn)低于理想膳食模型的15%～24%。

在能量與膳食營養(yǎng)方面（表1），城鄉(xiāng)居民膳食能量現(xiàn)狀趨近于中能量理想膳食模型（約2 200 kcal/d）。與中能量理想膳食模型相比，城鎮(zhèn)居民的脂肪攝入偏高、碳水化合物與膳食纖維攝入偏低，這與其畜肉食物比例偏高、果蔬食物比例偏低有直接的關(guān)系；農(nóng)村居民較高的糧食消費(fèi)比例一定程度上彌補(bǔ)碳水化合物攝入量，但在膳食纖維方面仍有待進(jìn)一步提高，增加果蔬食物比例及優(yōu)化動物性食物消費(fèi)品種是其未來膳食完善方向。

在膳食耕地足跡方面（表1），城鄉(xiāng)居民現(xiàn)有的人均膳食耕地足跡分別為929.81 m2、834.74 m2，高于低能量理想膳食模型（623.76 m2），趨近于中能量理想膳食模型（865.66 m2），低于高能量理想膳食模型（1 060.21 m2）。在未來保持現(xiàn)有農(nóng)作物生產(chǎn)水平情況下，隨著總?cè)丝诘脑黾蛹俺擎?zhèn)化水平的提升，預(yù)計(jì)到2030年中高能量理想膳食模型的耕地需求量將超過現(xiàn)在的12 493.40億m2（18.74億畝），分別達(dá)到12 552.09億m2（18.83億畝）、15 373.07億 m2（23.06億畝），其中高能量理想膳食模型的耕地需求量將超過中國現(xiàn)有耕地面積（19.179億畝）?；诮鼛啄晟攀衬芰康脑鲩L趨勢，預(yù)計(jì)未來膳食能量將進(jìn)一步提高，其所需的耕地需求量亦將進(jìn)一步增加，高能量膳食模型將對中國耕地保護(hù)提出更高的要求。為緩解耕地資源供需壓力，一方面立足國內(nèi)大循環(huán)，保護(hù)現(xiàn)有耕地資源、提高耕地生產(chǎn)力，另一方面充分利用國際大循環(huán)，借助“農(nóng)業(yè)走出去”、國際貿(mào)易等戰(zhàn)略提高食物供給能力?？梢?，耕地保護(hù)對于保障國民健康營養(yǎng)膳食具有關(guān)鍵性的作用。

4 討論

城鄉(xiāng)居民人均膳食耕地足跡先減少后增加的趨勢（圖4（a）—圖4（b）、圖6（a））主要受膳食結(jié)構(gòu)和農(nóng)作物生產(chǎn)水平變化的影響。由圖6中1981年與1999年的數(shù)據(jù)對比可以看出，城鄉(xiāng)居民在人均膳食量變化不大的情況下（圖6（b）），人均膳食耕地足跡有了較大幅度的減少，主要受其膳食等效農(nóng)作物生產(chǎn)水平大幅度提高的影響（圖6（c））。而從1999年與2019年的數(shù)據(jù)對比可以看出（圖6（a）—圖6（b）），城鄉(xiāng)居民人均動物性食物膳食消費(fèi)量與結(jié)構(gòu)比重的大幅度增加是促使其人均膳食耕地足跡增加的主要緣由。

圖6 城鄉(xiāng)居民膳食耕地足跡與食物等效農(nóng)作物的單產(chǎn)變化Fig.6 Dietary cultivated land footprints of urban and rural residents and yield changes of crop equivalents

需要說明的是，本文計(jì)算過程中存在一定的不確定性。該不確定性主要受數(shù)據(jù)有限性影響，主要體現(xiàn)在復(fù)種指數(shù)的估算和食物與等效農(nóng)作物轉(zhuǎn)換系數(shù)的確定兩個方面。前者由于統(tǒng)計(jì)口徑的不統(tǒng)一，耕地面積數(shù)據(jù)存在序列不連續(xù)的問題，本文按照文獻(xiàn)[29]的修訂方法對耕地面積進(jìn)行訂正，訂正后的數(shù)據(jù)更趨近于現(xiàn)實(shí)，但同時也帶來了不確定性；后者由于等效農(nóng)作物種類繁多，本文選取了代表性農(nóng)作物作為等效農(nóng)作物，這里存在著一定的不確定性。相信隨著調(diào)查統(tǒng)計(jì)數(shù)據(jù)覆蓋面與精確度的提高，這些不確定性將會降低。

5 結(jié)論與展望

可持續(xù)的耕地資源及健康營養(yǎng)的城鄉(xiāng)居民膳食是中國社會經(jīng)濟(jì)穩(wěn)定發(fā)展的重要基礎(chǔ)。本文基于1981—2019年長時間尺度中國城鄉(xiāng)居民食物消費(fèi)數(shù)據(jù)，運(yùn)用膳食營養(yǎng)分析及耕地足跡方法分析中國城鄉(xiāng)居民的膳食結(jié)構(gòu)與能量營養(yǎng)演變特征及其對耕地資源的影響。主要結(jié)論如下：

（1）城鄉(xiāng)居民膳食結(jié)構(gòu)與營養(yǎng)差距在1981—2019年間逐漸縮小。雖然城鎮(zhèn)居民的動物性食物比例均高于農(nóng)村居民，但城鄉(xiāng)居民膳食結(jié)構(gòu)差異在縮小，食物日趨均衡化，植物性食物比例日趨減少、動物性食物比例日趨上升。膳食結(jié)構(gòu)變動帶來膳食能量與營養(yǎng)的變遷，城鄉(xiāng)膳食能量、蛋白質(zhì)、脂肪與膽固醇近幾年均呈上升趨勢，膳食纖維城鄉(xiāng)居民均呈現(xiàn)出先減少后增加的趨勢，城鄉(xiāng)居民在碳水化合物方面均有所減少，農(nóng)村居民的碳水化合物減少幅度高于城鎮(zhèn)居民。

（2）城鄉(xiāng)居民膳食變遷直接影響到耕地足跡。城鄉(xiāng)居民的人均膳食耕地足跡較于20世紀(jì)80年代均有所減少，主要源于農(nóng)業(yè)生產(chǎn)水平的提高。但隨著城鎮(zhèn)化的推進(jìn)、動物性食物比例上升、總?cè)丝谠黾?、城?zhèn)人口上升、農(nóng)村人口減少，1981—2019年城鄉(xiāng)居民的總膳食耕地足跡波動中有所上升，城鎮(zhèn)居民的總膳食耕地足跡增加，農(nóng)村居民的總膳食耕地足跡有所減少。

（3）未來城鄉(xiāng)居民膳食仍需進(jìn)一步優(yōu)化，耕地保護(hù)對于保障國民健康營養(yǎng)膳食具有關(guān)鍵性的作用?，F(xiàn)城鄉(xiāng)居民的水果及奶類食物比例較理想膳食模型仍偏低，農(nóng)村居民偏低更為明顯，而城鄉(xiāng)居民的畜肉食物則偏高，食物消費(fèi)品種搭配有待進(jìn)一步優(yōu)化。與中能量理想膳食模型相比，城鎮(zhèn)居民的脂肪攝入偏高、碳水化合物與膳食纖維攝入偏低，農(nóng)村居民在膳食纖維方面仍有待進(jìn)一步提高，增加果蔬食物比例及優(yōu)化動物性食物消費(fèi)品種是其未來城鄉(xiāng)膳食均衡的完善方向?；诮鼛啄晟攀衬芰康脑鲩L趨勢，未來膳食能量將進(jìn)一步提高，預(yù)計(jì)到2030年中高能量理想膳食模型的耕地足跡將進(jìn)一步增加，其中高能量理想膳食模型的耕地足跡將超過中國現(xiàn)有耕地面積，這意味著高能量膳食模型的耕地需求量將超過實(shí)際耕地供給量，中國耕地保護(hù)的壓力將進(jìn)一步增加。保護(hù)現(xiàn)有耕地資源、提高耕地生產(chǎn)力，同時借助“農(nóng)業(yè)走出去”、國際貿(mào)易等戰(zhàn)略充分利用國際資源是緩解中國耕地資源供需壓力的重要舉措。