下载:
-
生态系统服务是指人类从生态系统获得的各种惠益，如供给服务、调节服务、支持服务和文化服务等[1]。生态系统提供了人类几乎所有的福祉要素[2]，人类的生存和发展依赖于生态系统所提供的服务。目前，生态系统的退化已经对人类福祉造成了日益严重的损伤，定量研究生态系统服务价值，有助于提高对生态系统服务价值的认识，同时为政府政策权衡提供依据，有利于对自然生态系统进行可持续开发和利用，从而对自然生态系统的进行有效的管理与保护[3]。近些年来，国内外对生态系统服务价值的计算研究逐渐日趋增多[4]，相关理论与方法也不断涌现。国外学者在全球和区域尺度、流域尺度、单个生态系统尺度、单项服务价值方面开展了大量的研究土作[5]。如Sutton等[6]的全球尺度生态系统服务评价，Grêt-Regamey等[7]的关于区域尺度生态系统服务评价，Pattanayak[8]的关于流域尺度服务价值，Lal[9]的关于单个生态系统尺度生态服务价值核算，以及de Mendonca[10]关于单项生态系统服务价值的研究。国内关于生态系统价值核算的研究起步较晚，欧阳志云等[11]在1999年首次使用生态系统服务的概念并对中国陆地生态系统的6种服务价值进行了初步评估，也是中国较早研究生态系统服务价值的论著之一。随后，赵景柱等[12]对13个国家生态系统服务价值的测算进行比较分析，赵同谦等[13-15]对中国陆地地表水、中国草地、中国森林、海南岛等大尺度区域的生态系统服务进行了价值评估，谢高地等[16]在Costanza等提出的评价模型的基础上，对国内200多位生态学学者进行了问卷调查，提出了“中国生态系统服务价值当量因子表”。21世纪以来，国内外的研究范围更加广阔，但对昆虫生态系统服务的研究依然很少涉及。昆虫生态系统服务是指昆虫类群在生态系统服务中发生的作用，以及为人类提供的各种收益，一般包括有形收益的产品和无形收益的服务[17]。国内对昆虫生态服务的研究主要停留在定性的服务作用分析上，例如传粉、病虫害控制、物质分解等[18]。国外对蜜蜂Apis spp.等昆虫的生态服务价值（如传粉等）已经开始进行定性的研究和定量核算，但尚未形成共识的理论体系和方法。本研究以河南省蜣螂Catharsius molossus为例，初步对其生态系统服务价值的进行测算，以期推动国内昆虫服务价值评价研究与发展，同时对生态服务价值与政府政策权衡进行简单评述，为政府政策的生态意识提供参考。
-
河南省首位于中国中东部，地处沿海与中西部结合部，31°23′~36°22′N，110°21′~116°39′E，属于暖温带-亚热带、湿润-半湿润季风气候。冬季较为寒冷雨雪少，春季干旱风沙多，而夏季炎热雨水丰沛，秋季晴和日照充足。全省的年平均气温为12~16 ℃，1月-3.0~3.0 ℃，7月24.0~29.0 ℃，大体东高西低，南高北低，山地与平原间差异相对比较明显。在中国动物地理区划上，则处于华北区和华中区交接地带。因此，河南省动物区系构成呈现复杂化、多样化的局面。
-
河南省由于特殊的地质历史、交错的地理位置和过渡性的气候环境，生物资源非常丰富。河南省昆虫物种有7 387种，是目前全国记录昆虫种类最多的省份，为河南省生物产业的发展提供了宝贵的资源储备，也为世界生物多样性和中国生物资源宝库增添了科学财富[19]，创造了巨大的生态服务价值。蜣螂在河南省分布极广，各个地市均有分布，具有较强的地域代表性。
-
对生态服务价值评定一般采用的是Costanza和Chee的评价方法[20]，即物质量评价法（PAM）和价值量评价法（VAM）。物质量评价法主要是从物质量的角度对生态系统提供的服务进行整体评价，而价值量评价法主要是价值量的角度对生态系统提供的服务进行整体评价[21]。本研究主要运用价值量评价法中的价值替代法进行评价。本研究选取蜣螂作为评价对象，对其生态服务价值进行核算。蜣螂的生态系统服务功能主要有生粪便分解、牧草污染处理、氮循环和有害生物（蚊蝇等）控制等[22]，其生态系统服务范围主要集中在粪便中。我们主要以蜣螂对牛粪的生粪便处理、牧草污染控制和氮循环这3种生态系统服务价值进行核算。对具体生态服务价值指标的选取并不是以其在生态系统服务中所起的作用，而是以数据的可获取性和方法的可操作性为原则。
-
蜣螂可以将固体牛粪破碎化并分解入地下，但仅限于草地和土地（简称自然地块），蜣螂活动能力在水泥地和铺装路面上受限。农村牛类养殖基本上都是自然地块个体放养，而养殖场的混合地块较多。另外，一些养殖场采用药物对粪便进行杀菌消毒，蜣螂排斥杀菌消毒处理的牛粪。据国外研究，在自然情况下牛粪完全分解的平均时间为（28.14±0.71）个月。在有蜣螂作用的情况下，分解平均时间为（22.74±0.64）个月[24]。因此，可以得出蜣螂将粪便分解速率提升了19%。采用利用价值替代法，将蜣螂处理粪便的服务价值转化为人工处理粪便等固体垃圾所需要的价值量。
$${V_{\rm{d}}} = T \times {R_1} \times {P_1} \times D \times I \times {M_{\rm{d}}}.$$ (1) 式（1）中：Vd为蜣螂处理粪便的价值，元；T为河南省牛的数量，头；R1为自然地块的牛的比例；P1为粪便没有被药物处理的牛的比例；D为牛所产生的粪便量，m3·头-1；I=19%；Md为当地人工处理土方量的单位价值，元·m-3。
-
自然地块上的牛粪堆积能提高牛粪堆下水肥含量，改善土壤理化性质，利于植物生长。但是新鲜牛粪质地黏稠，分解速度缓慢，造成牛粪下植物缺氧、少光而死亡。牲畜拒绝取食牛粪周围的草地面积高达牛粪面积的17.15倍。在缺少蜣螂的处理分解下，牛取食范围缩小，造成牛肉减产达7.63 kg·头-1[24]。采用价值替代法计算如下：
$${V_{\rm{f}}} = T \times C \times {M_{\rm{b}}}$$ (2) 式（2）中：Vf为蜣螂处理牧草（饲料）污染的生态系统服务价值，元；T为河南省牛的数量，头；C=7.63 kg·头-1；Mb为当地牛肉价格，元·kg-1。
-
以牛粪中的氮循环为例。格拉德研究，氮在牛粪中的含量约为2%，自然状态下没有蜣螂处理的粪便大约80%的氮直接释放到空中。然而蜣螂可以将牛粪输入地下，为植物提供养分。同时牛大约能产生氮27 kg·头-1·a-1[25]，转化为人工施氮肥的成本进行替代如下：
$${V_{\rm{n}}} = T \times {M_{\rm{n}}} \times {P_2} \times {R_2}.$$ (3) 式（3）中：Vn为氮循环价值，T为河南省牛的数量，头；Mn为当地人工施氮肥的成本，元·kg-1；R2取2%，P2取80%。
-
D参考的数值来自于文献[26]，约为21 m3。T参考的数值来自于河南省畜牧业统计数据[22]，为9.55×106头。P1和R1来自河南省统计年鉴结合抽样调查数据P1=72%，R1=80%，I=19%[24]，Md为当地人工处理牛粪土方量单位价值25.00元·m-3。代入式（1）中得到5.49×108元。即河南省牛粪处理的生态价值约为5.49×108元。这里仅以蜣螂处理牛粪的生态服务价值进行核算，实际中蜣螂关于粪便的处理还包括马、猪、鸡等其他家庭驯养类动物和兔子、狼等野生动物。家庭驯养类的可以根据此法结合相关数据类推，野生动物的数据尚需要长期观测，暂未能深入研究。因此，蜣螂处理粪便的生态服务价值远远大于5.49×108元。
-
式（2）中Mb的参考值来自于当地牛肉市场价格50.00元·kg-1；T取自于河南省畜牧业统计数据为9.55×106头；C取7.63 kg·头-1。代入式（2）得出蜣螂对牧草污染处理的价值Vf为3.64×109元。
本研究在计算处理牧草污染的价值时，将其转化为牛肉减产的产量价值，也可以采用牛类拒绝使用的牧草面积转化为牧草价值来替代计算。
-
式（3）人工施氮肥的成本取值：以河南省的氮肥以及复合肥的销量为基础，计算出它们的相互比例，同时结合氮元素在其中的比例以及销售价格，计算出氮元素的价格（元·kg-1），结合河南省平均人工费用，得出Mn为1.50元·kg-1。R2取2%，P2取80%。代入式（3）中得出蜣螂对碳循环的处理价值Vn为9.09×108元。在牛粪的氮释放过程中，除了蜣螂作用外，还有其他微生物作用，计算结果比实际偏大。
-
通过对蜣螂的3项生态系统服务（生粪便分解、牧草污染处理、氮循环）价值的评估，得出其在河南省的价值量约为5.10×109元。结果仅是蜣螂对生牛粪处理的生态服务价值，同时也是蜣螂若干生态服务价值中的一部分。而蜣螂只是昆虫种类之一，昆虫的种类和数量巨多，因此，昆虫生态服务价值远大于本研究数值。
长期以来，人类对其依赖和生存的环境所提供的生态系统服务定量不准，尤其是昆虫的生态服务价值的核算几乎没有涉及，在经济发展和环境保护的困惑中，无法精确衡量其中的利弊关系。本研究通过价值替代法对蜣螂部分的生态系统服务进行了粗略的价值统计。价值替代法对于大范围的区域尺度的生态系统价值核算提供了便利性，同时避免了相对复杂的生态过程分析，对缺乏难获得的生态检测数据的生态系统服务价值核算问题的处理具有一定的积极意义；但由于评价中缺乏对生态系统结构、生态过程与服务功能的深入分析，导致结果往往缺乏直接的生态学依据。同样不同地区的物价、劳动力价值水平存在差异，影响到评价结果的一致性。更为准确的昆虫生态服务价值的核算有以下几个方面可以讨论：
生态系统是一个开放的系统，昆虫的迁徙和活动不受行政区划的控制，因此，对数据的取得和计算应进行基于昆虫活动过程和范围的全局考虑。本研究采用的一些数据和研究方法以国外为主，国外的一些昆虫习性和地理特征等相关因素与国内存在差异，因此，计算结果有一定偏差，精确的测算需要加强对本土昆虫的生态系统服务进行长期的观测和深入的过程研究。
关于昆虫对野生动植物以及微生物，即对生物多样性的影响所带来直接或间接的对人类的生态系统服务价值，因限于数据的取得，本研究没有进行具体核算。这将是昆虫生态服务影响中的重点和难点。
-
目前，生态服务价值评估在全球范围都是一个热点，之前的研究侧重于全球和区域范围，最近的研究偏向于基于生态过程的单项生态系统。昆虫生态系统服务价值核算研究也出于此背景下。定量生态系统服务价值的意义在于有助于提高公众的环境意识，同时也是政府政策决策权衡的重要依据：①促进各级政府树立科学的发展观，正确权衡和评价每项措施的影响。政府在制定任何项目前期，必须以科学、准确的生态系统价值评估为依据，将项目带来的经济收益与项目带来的生态系统服务价值丧失进行对比，或者是在进行生态系统保护工程规划时权衡项目自身花费成本与带来的生态系统保护价值的增加量进行对比，来决定政策或项目的实施与否。②研究生态系统服务的边际价值，确定精确的生态系统补偿标准和方案。尽管生态系统服务的总体价值巨大，但对于具体的问题分析显得过于笼统，因此，研究优化或恶化导致生态服务功能的增加或减少的边际价值更具有实际意义。如生态系统购买者的活动导致某种生境丧失引起昆虫生态系统服务价值的降低，必须由生态系统购买者对生态系统服务提供者进行生态补偿，因此，边际价值的确定不仅需要生态系统提供服务的真实值，还需要量化生态服务功能的变化对人们偏好的满足程度来确定使之成为一个长期和谐的持续机制。③推进环境与国民经济的科学评价：当前环境评价往往以某类指标的高低作为某种环境因素的评价标准，如细颗粒物（pm 2.5）浓度的高低反映雾霾的严重程度，不能表达环境与生态系统服务变化所产生的具体价值影响，即雾霾天气对人类健康带来的具体经济损失等。将生态价值与环境评价联系起来，有助于具体的环境科学评价，促进环境成本纳入国民经济核算体系，现行的国民经济核算体系以国民生产总值或国内生产总值作为主要指标，只重视经济产值及其增长速度的核算，而忽视国民经济赖以发展的生态资源基础和环境成本的计算。
Research and discussion on the valuation of insects' ecosystem services
-
摘要: 以蜣螂Catharsius molossus为例, 采用Costanza和Chee的价值替代法对河南省蜣螂的生态服务如生粪便分解、牧草污染处理和氮循环的价值统计, 得出其生态服务价值为5.10×109元, 阐述昆虫提供的巨大生态服务价值。初步提出了昆虫生态系统服务价值的测算方法, 以期推动昆虫服务价值评价的研究发展, 同时对生态系统服务价值与政府决策的权衡进行探讨。Abstract: This paper used Costanza and Chee's value substitution method to calculate the value of Catharsius molossus' ecosystem services including dung decomposition, pasture pollution treatment and nitrogen cycle in Henan Province. The economic value of the three services added up to 5.1 billion RMB which indicated that the insects provided substantial value of ecological services. The proposed method to calculate the value of insects' ecosystem services would promote the research on the valuation of insects' services. The paper also discussed the tradeoff between ecosystem services and decision-making of government.
-
南方集体林区以丘陵为主的地形特征决定了林业生产过程中客观上存在着机械替代率低的问题。由于大量农村劳动力进入城市非农就业，林业经营主体需要通过雇佣劳动力以解决劳动力不足的问题[1]；日趋扩大的规模化经营对雇佣劳动力的依赖程度急剧增大，委托代理关系也应运而生。委托代理关系下的营林生产过程往往存在着信息不对称和“道德风险”等问题，导致雇主对劳动过程的监督不到位[2]；而相对粗放的林业生产中，雇主无法准确观测劳动力的投入程度，雇工生产效率较低[3]。因此，在南方集体林区不断推进林业规模化经营的背景下，开展雇工劳动质量对农户营林技术效率的影响研究具有重要的理论和现实意义。国内外学者关于雇佣劳动对农林业生产影响的研究主要集中在农业领域。一种观点认为雇佣劳动对农业生产效率产生了负面影响。雇佣劳动存在“道德风险”，雇佣劳动力的边际产出小于自有劳动力的边际产出，要保证劳动质量必须要付出大量的监督成本。有研究发现：存在雇工行为的家庭农业经营雇主，需要花费10%的劳动力时间来监督雇工，对劳动力在其他用途上的配置产生挤出效应。还有研究表明：团队组织内部监督不完备与激励不足，将会降低生产效率[4]。另一观点则认为雇佣劳动对农业生产效率有正向影响。张五常[5]认为分成制、定租制和混合制等委托-代理机制[3]是有效的。相比之下，目前针对林业雇工研究较少，从雇工视角探讨南方集体林区农户营林生产效率影响因素的研究还未见报道。本研究基于浙江省的实地调研数据，采用计量模型测算样本农户的营林技术效率，分析雇工劳动质量对农户营林技术效率的影响，以期为林业规模化经营中存在的雇工劳动问题提供决策依据。
1. 雇工劳动质量对农户营林技术效率影响的理论机制分析
1.1 研究假说
根据委托-代理理论[3]，存在于雇佣劳动中“道德风险”会导致雇佣劳动力的边际产出小于自有劳动力的边际产出；而相对粗放的林业生产管理会增加劳动力的监督难度，加重林业雇佣劳动中的机会主义行为。即营林生产过程中雇佣劳动会在一定程度上影响技术效率的提高。同时，雇工投入水平的差异对农户营林技术效率的影响也是不同的。雇工投入越多，则家庭自有劳动力越少，劳动监督难度越大，雇工劳动力的道德风险也越大，从而导致劳动生产质量下降，林业生产效率会因此受到影响（落入到生产可能性曲线内部）。就家庭自有劳动力与雇工关系而言，拥有剩余索取权的家庭自有劳动力更多地参与到营林生产中，则边际产出会随之升高，劳动质量会得到更大程度的保证，雇佣劳动对林业生产效率的负面影响程度也相对减弱。因此，可提出以下研究假说：农户营林生产过程中的雇佣劳动会对营林生产造成负面影响，且雇工劳动质量不同，雇佣劳动对农户营林技术效率的影响程度也不同。
1.2 理论模型
技术效率（technical efficiency, TE）是目前测度生产单位效率水平最常见的指标之一，一般用生产单位的实际产出与其理论上所能实现的最大潜在产出的比值来表示，本研究选用技术效率来衡量样本经营农户营林生产效率水平。用Yi表示第i个农户的林业产出，则随机前沿生产函数（stochastic frontier approach，SFA）可表示为：Yi＝f（Xi，β）exp（Vi-Ui）。其中：f（Xi，β）代表生产前沿面，Xi为生产要素的投入，本研究指资本、土地和劳动力的投入；（Vi-Ui）为混合误差，其中Vi表示随机误差，包含测度误差及不可控因素；Ui为技术效率损失，表示农户i的技术非效率项，服从独立的截断正态分布N（mi，δu2），其中mi为数学期望，δu2为方差。
技术无效率函数（mi）可设定为：$ {m_i} = {\delta _0} + \sum\limits_j {{\delta _j}{Z_{ji}} + {\mu _i}} $。其中：i表示第i位样本农户，Zji表示影响技术效率的各外生变量，δj为待估参数，反映各外生变量对技术效率的影响程度，μi为纯随机误差项。
在此基础上求解出技术效率水平：Ti＝E（Yi∣Ui，Xij）/E（Yi∣Ui＝0，Xij）＝exp（-Ui），从而求解出平均技术效率：$ T = \frac{1}{n}\sum\limits_{i = 1}^n {{T_i}} $。其中：n为农户数量，本研究采用最大似然估计法来估计前沿生产函数的参数。
2. 数据来源、实证模型与变量选择
2.1 数据来源
本研究以浙江省南方集体林区为案例点，充分考虑地理位置、自然资源和社会经济发展差异，按农民人均可支配收入分层抽样选取浙江省建德市、开化县2个县（市）作为样本县（市）。2地农民人均可支配收入为23 998和15 736元，存在明显差异，反映出浙江省不同地区经济发展的差异性；森林覆盖率为76.2%和80.9%，均超出浙江省的平均水平，林业代表性较好。采取随机抽样法，在2个样本县（市）随机选取6个乡（镇），其中开化3个乡（镇），建德3个乡（镇）。为保证问卷的质量与信息的真实性，调查采取“一对一”的访谈模式。共得到农户有效样本245户，其中有雇工参与的农户151户，无雇工参与的农户94户。
杉木Cunninghamia lanceolata是样本地区分布最广的树种之一，也是当地农户最主要的林业收入来源之一；因此，本研究以杉木为案例树种进行调研。调查内容包括农户家庭基本情况、林地基本情况、1个营林周期内最大地块上杉木的生产投入和产出情况等。样本分布见表 1。
表 1 农户有效样本分布情况Table 1. Specific distribution of effective sample of farmers县（市） 乡（镇） 村 合计/户 比例/% 开化 华埠 许家源 20 8.2 联丰 20 8.2 池淮 芹源 20 8.2 玉坑 20 8.2 芹阳 泉坑 21 8.6 小桥头 20 8.2 建德 李家 沙墩头 8 3.3 长林 23 9.4 石鼓 9 3.7 龙桥 1 0.4 新桥 2 0.8 李家 2 0.8 建德 大同 上马 1 0.4 小溪源 24 9.8 永平 1 0.4 竹林 1 0.4 竹源 13 5.3 航头 大店口 21 8.6 东村 14 5.7 溪沿 1 0.4 罗源 1 0.4 曹源 2 0.8 总计 245 100.0 2.2 实证模型与变量选择
在测算林业技术效率之前，需要先确定前沿函数的具体形式，常见的前沿函数形式有C—D生产函数和超越对数生产函数。相对于前者，后者不仅形式灵活、易估计、包容性强，而且允许要素间替代弹性可变，没有对技术变化附加任何限制条件，因此近似性更好。考虑到南方集体林区农户拥有林地地块数量普遍较多，而农户林业生产中最重视且存更大雇佣劳动可能性的一般是最大地块，因此，本研究采用农户最大地块1个经营周期内的营林投入产出变量来构建生产函数。超越对数生产函数公式如下：lnYi＝β0+β1lnKi+β2lnLi+β3lnMi+β4（lnKi）2+β5（lnLi）2+β6（lnMi）2+β7lnKilnLi+β8lnMilnMi+β9lnKilnMi+Vi-Ui。其中：Y表示最大地块1个经营周期内的林业产出（主伐量）（m3）；K表示农户在最大地块的1个经营周期内所投入的费用（元）；L表示农户在最大地块的1个经营周期内所投入的劳动投工数量（工）；M表示最大地块面积（hm2）；β0~β9为待估参数。
通过对超越对数生产函数求导可以分别计算出各投入要素的产出弹性。资本投入要素的产出弹性为：η1＝β1+2β4lnK+β7lnL+β9lnM。劳动力投入要素的产出弹性为：η2＝β2+2β5lnL+β7lnK+β8lnM。土地投入要素的产出弹性为：η3＝β3+2β6lnM+β8lnL+β9lnK。农户技术无效率函数表示为：$ {m_i} = {\delta _0} + \sum\limits_{j = 1}^{j = 5} {{\delta _j}{h_j} + \sum\limits_{j = 6}^{j = 16} {{\delta _j}{Z_{ji}} + {\mu _i}} } $。其中hj为关键变量，代表雇工劳动质量，表示农户最大地块1个经营周期内雇佣劳动力的投工质量情况，Zj为控制变量。
为全面测度营林雇工劳动质量问题，本研究构建了雇工年龄比例、雇工性别比例、雇工投工所占比例、受过技术培训的雇工所占比例等具体指标[6]作为关键变量。具体如下。①雇工年龄比例：40岁以下雇工所占比例（h1）、60岁以上雇工所占比例（h2）；②雇工性别比例：男雇工所占比例（h3）；③雇工投工所占比例：总投工中雇工所占比例（h4）；④受过技术培训的雇工所占比例（h5）。
一般研究将影响技术效率的地块特征、农业生产特征、户主特征等外生因素作为控制变量。①地块特征主要包括地块的地理位置（最大地块离家距离）、农户家庭山林总面积、农户家庭林地总块数、最大地块立地质量。地块离家距离越近，越便于林农对林地进行管理，所以，预期最大地块离家距离会对农户营林技术效率产生负向影响。农户家庭山林总面积越大、林地总块数越多意味着农户管理林地的难度也越大，平均到每个地块的管理时间就会越少，因此，预期其会对农户营林技术效率产生负向影响[7]。土地细碎不利于先进机械设备和技术的推广，控制病虫害难度加大，难以实现规模经营[8]，虽然也有研究表明耕地面积与生产率之间的反向关系[9]，但有些研究却表明此种关系并不显著[10]，因此预期其总效应不明确。地块的立地质量对效率的影响是显而易见的，肥沃的土壤相对于贫瘠的土壤更能提高农户的生产效率[11]，因此预期其会对农户营林技术效率会产生正向影响。②农业生产特征包括家庭务农人数、家庭总收入以及农户是否为补贴户。家庭务农人数会直接影响到劳动力要素的投入，因此预期其会对农户营林技术效率产生正向影响。农户家庭总收入的增加会加大农户的林业投资[12]，预期会对农户营林技术效率产生正向影响。农户若为补贴户，林业补贴的增加也会增加农户对林业的资本投入，因此假设其会对农户营林技术效率产生正向影响。③户主特征。包括户主年龄、户主的受教育年限、户主健康状况、户主是否担任过村干部。户主年龄对其技术效率的影响方向取决于该农户是更富有经验还是更守旧[13]；户主受教育年限越长，越能有效利用先进的农业生产技术[14]，因此，假设其会对农户营林技术效率产生正向影响。户主健康状况对农户营林技术效率的影响也是显而易见的：户主健康状况越好，越有利于家庭营林生产劳动和决策，因此假设其会对技术效率产生正向影响。干部身份一方面会带来收入效应，即干部获得先进生产技术和农业生产信息的渠道更多，这会对农户营林技术效率产生正效应；另一方面，干部身份同样存在替代效应，即干部用于家庭经营的时间更少，从事家庭经营的机会成本也较高，这会对农户营林技术效率产生负效应[15]。因此其具体影响尚不可知。具体控制变量如下：Z6表示农户是否为补贴户（0代表否，1代表是）；Z7表示户主是否为村干部（0代表否，1代表是）；Z8为农户家庭务农人数（人）；Z9为户主年龄（岁）；Z10为户主受教育年限（年）；Z11为户主健康状况（1代表好，2代表中，3代表差）；Z12表示农户家庭总收入（元）；Z13为农户家庭山林总面积（hm2）；Z14表示农户经营的山林总块数（块）；Z15表示农户经营山林中最大地块的立地质量（1代表好，2代表中，3代表差）；Z16表示最大地块离家的距离（km）。
3. 结果与分析
3.1 描述性统计分析
由表 2可知：样本农户最大地块整个营林周期内单位面积的平均产出为108.62 m3·hm-2，农户最大地块整个营林周期内单位面积的平均资本投入为8 214.47元·hm-2，单位面积平均劳动力投入为375.23工·hm-2，平均林地投入为2.02 hm2。
表 2 随机前沿生产函数模型变量的描述性统计Table 2. Descriptive statistics of variables in Stochastic Frontier Approach's model统计值变量 最大地块总产出/（m3·hm-2） 最大地块资本投入/（元·hm-2） 最大地块面积/hm2 最大地块劳动力投入/（工·hm-2） 平均值 108.62 8 214.47 2.02 375.23 标准差 94.70 5 187.26 1.14 278.28 技术效率损失模型中（表 3），户主的平均年龄为57.24岁，平均受教育年限为7.20 a，可以看出该地区的劳动力质量较差。88.00%的立地质量为中等及以上，最大地块离家距离平均为1.97 km，表明样本地区的立地质量和交通条件对林业经营相对有利。样本农户家庭户均地块为3.41块，表明样本地区林地细碎化问题并不严重。农户的家庭总收入均值为95 501.74元，说明样本地区当地的经济条件利于林业发展。所有样本农户中仅有27.00%的农户为补贴户，说明国家的林补政策还未真正地惠及该地区。
表 3 技术效率损失模型变量的描述性统计Table 3. Descriptive statistics of variables in the loss of technical efficiency's model变量类型 具体变量 平均值 标准差 最小值 最大值 雇工劳动质量 40岁以下雇工所占比例 0.45 0.24 0.10 1.00 60岁以上雇工所占比例 0.49 0.17 0.01 1.00 男雇工所占比例 0.85 0.16 0.20 1.00 总投工中雇工所占比例 0.72 0.25 0.05 1.00 雇佣的工人中受过技术培训的工人所占比例 0.78 0.42 0.10 1.00 户主特征 户主年龄（岁） 57.24 9.38 27 86 户主教育年限 7.20 3.52 0 16 户主是否为村干部（0代表否，1代表是） 0.33 0.46 0 1 户主健康状况：好 0.84 0.36 0 1 户主健康状况：差 0.12 0.22 0 1 农业生产特征 家庭务农人数（人） 1.18 1.05 0 5 家庭总收入（元） 95 501.74 104 235.30 520 724 652 是否为补贴户（0代表否，1代表是） 0.27 0.44 0 1 地块特征 家庭总地块数 3.41 2.87 0 20 最大地块质量：好 0.54 0.49 0 1 最大地块质量：差 0.12 0.31 0 1 最大地块离家距离（km） 1.97 2.03 0.02 15 山林总面积 3.52 8.26 0.03 96.67 说明：“户主健康状况”参照组为“中”，“最大地块质量”参照组为“中” 为消除因自变量之间多重共线性导致的模型估计结果偏差，在模型估计前对雇工质量各指标进行相关性检验。结果发现（表 4）：雇工质量各指标之间不存在多重共线性问题，各指标可以作为自变量放入模型进行估计。
表 4 雇工质量各指标系数相关矩阵Table 4. Relevance matrix of index coefficients of employee quality40岁以下雇工所占比例 60岁以上雇工所占比例 男雇工所占比例 总投工中雇工所占比例 雇佣的工人中受过技术培训的工人所占比例 40岁以下雇工所占比例 1.00 60岁以上雇工所占比例 0.18 1.00 男雇工所占比例 0.24 0.15 1.00 总投工中雇工所占比例 0.01 0.00 0.33 1.00 雇佣的工人中受过技术培训的工人所占比例 0.24 0.02 0.10 0.23 1.00 3.2 实证结果分析
利用广义似然比（LR）检验可降低对SFA模型的依赖，避免函数形式的误设，从而从设定的待估计模型中筛选出最能拟合样本数据的模型。LR公式可表示为：RL＝-2[lnL（H0）-lnL（H1）] ~χ2（k）。其中：L（H0）和L（H1）分别是零假设H0和备择假设H1下的似然函数值，表示受约束条件的自由度。将LR统计量与临界值进行比较，当LR统计量值大于临界值时拒绝原假设，否则，接受原假设。给出的2个零假设为：（1）规模户和非规模户的前沿面并没有显著的差异，即模型不需要添加是否为规模户的虚拟变量。（2）外生变量对技术效率无任何影响，即模型不需要添加外生变量影响因素。LR验证结果如表 5所示。相对于基准模型，假设1在1%显著性水平上没有被拒绝，而假设2在1%显著性水平上拒绝原假设；说明原假设1对应的模型较好地拟合了样本数据，可作为本研究测度技术效率的主要模型。
表 5 假设检验结果Table 5. Hypothesis test results零假设 LR统计量 自由度 χ2 0.01临界值 结论 H0：不应该设置规模户虚拟变量 0.001 9 20.97 接受 H0：外生变量对技术效率无影响 126.310 14 28.49 拒绝 3.2.1 随机前沿生产函数模型估计结果分析
表 6为随机前沿生产函数模型的估计结果。将表 6的回归系数代入上文生产投入要素产出弹性计算公式中可得到各投入要素的产出弹性。计算得：土地投入要素的产出弹性为2.25，说明样本地区林业生产对土地投入的依赖程度较高，即林业生产中最为稀缺的生产要素是土地，增加土地投入可以大幅度地提高林业产出。资本（-0.09）和劳动力（-0.23）投入要素的产出弹性均为负值，说明目前样本地区林业存在过度投入资本和劳动力的情况，单纯依靠增加林业劳动力和林业资本投入并不会带来林业产出的增加，相反还可能导致林业产出减少。
表 6 随机前沿生产函数模型估计结果Table 6. Estimated results of Stochastic Frontier Approach's model变量 系数 变量 系数 最大地块资本投入 0.589***（0.091） 劳动力投入的平方项 -0.029*（0.016） 最大地块面积（土地投入） 0.746***（0.240） 资本投入×土地投入 0.160**（0.068） 最大地块劳动力投入 0.017（0.125） 土地投入×劳动力投入 0.108*（0.063） 资本投入的平方项 0.051***（0.010） 资本投入×劳动力投入 -0.014（0.035） 土地投入的平方项 -0.410***（0.137） 常数项 -0.368（0.259） 说明：*，**，***分别表示通过10%，5%，1%水平下的显著性检。括号内数值为回归标准误 3.2.2 技术效率损失模型估计结果分析
在245份有效样本农户中，家庭最大地块1个营林周期内有雇佣劳动力的农户有151户，占总体样本的61.63%。总体农户平均技术效率值为0.57，有雇工农户平均技术效率值为0.59，无雇工农户的平均技术效率值为0.76。由表 7可知：如果消除技术效率的损失，总体样本农户的平均技术效率还有43.00%的提升空间。由描述性统计结果可粗略看出无雇工农户的平均技术效率高于有雇工农户的平均技术效率，但雇佣劳动对农户技术效率的具体影响有待进一步计量分析。由样本农户雇工情况对农户营林技术效率损失影响的估计结果（表 8）可知：在控制其他变量不变的情况下，雇工会对农户营林技术效率造成负面影响（P＜0.10）。雇工劳动质量指标中，总投工中雇工所占比例对农户的营林技术效率具有负向影响（P＜0.05）；原因可能是家庭自有劳动力和雇佣劳动力劳动质量存在异质性，雇主对劳动过程的监督很难到位，由此造成总投工中雇工所占比例越高，农户营林技术效率越低。户主年龄对农户的营林技术效率具有正向的影响（P＜0.10）；原因可能是随着户主年龄增大，其营林生产经验越丰富，对家庭营林生产越有利。户主良好的身体状况对农户营林生产技术效率具有正向影响（P＜0.05）；作为家庭最主要的林业劳动力和决策者，户主身体健康程度对林业生产至关重要。家庭总收入对农户营林生产技术效率具有正向影响（P＜0.05）；原因在于农户家庭总收入的增加会减少农户家庭林业生产的资金约束，农户林业投资概率会增大。山林总面积和家庭总地块数都对农户的营林技术效率具有负向影响（P＜0.10）；农户家庭山林总面积越大、家庭总地块数越多，农户管理林地的难度也越大，平均到每个地块的管理时间就会越少，农户无法对林地进行精细化地管理，影响了技术效率的提高。好的地块质量对农户营林生产技术效率具有负向影响（P＜0.10）；这与预期的影响方向相反，可能的原因是：农户会对质量较好的地块相对投入更少的肥料和劳动力等生产要素，因此使立地质量较优的地块产出情况反而不如立地质量较差的地块。
表 7 样本农户营林技术效率总体情况Table 7. Overall situation of technical efficiency of sample farmers in forestry management描述性统计农户类型 平均值 标准差 最小值 最大值 全部农户 0.57 0.20 0.12 0.92 有雇工农户 0.59 0.25 0.06 0.98 无雇工农户 0.76 0.12 0.24 0.92 表 8 雇工情况及雇工劳动质量对农户营林技术效率损失的影响估计结果Table 8. Estimation of the impact of employment and labor quality of employees on technical efficiency of farmers' forestry management变量类型 具体变量 系数 具体变量 系数 雇工情况 是否雇工（0代表否，1代表是） 0.373*（0.212） 40岁以下雇工所占比例 1.139（0.814） 及雇工劳 60岁以上雇工所占比例 -0.474（0.509） 动质量 男雇工所占比例 0.662（0.511） 总投工中雇工所占比例 1.205**（0.538） 雇佣的工人中受过技术培训的工人所占比例 -2.115（1.328） 户主特征 户主年龄（岁） -0.011（0.009） 户主年龄（岁） -0.027*（0.016） 户主教育年限 -0.008（0.022） 户主教育年限 0.052（0.044） 户主是否为村干部（0代表否，1代表是） 0.182（0.160） 户主是否为村干部（0代表否，1代表是） 0.181（0.245） 户主健康状况：好 -0.498**（0.242） 户主健康状况：好 -1.130**（0.489） 户主健康状况：差 -0.180（0.324） 户主健康状况：差 -0.132（0.883） 农业生产 家庭务农人数（人） 0.035（0.078） 家庭务农人数（人） 0.084（0.106） 特征 家庭总收入（元） -0.000**（0.000） 家庭总收入（元） -0.000**（0.000） 是否为补贴户（0代表否，1代表是） 0.333**（0.165） 是否为补贴户（0代表否，1代表是） 0.317（0.259） 地块特征 家庭总地块数 0.035（0.026） 家庭总地块数 0.090*（0.049） 最大地块质量：好 0.221（0.155） 最大地块质量：好 0.640*（0.359） 最大地块质量：差 -0.117（0.241） 最大地块质量：差 0.196（0.460） 最大地块离家距离（km） -0.040（0.037） 最大地块离家距离（km） -0.052（0.061） 山林总面积 -0.005（0.022） 山林总面积 0.042*（0.022） 常数项 1.643***（0.623） 常数项 0.303（1.099） σ2 0.318***（0.070） σ2 0.801***（0.268） γ 0.876***（0.073） γ 0.995***（0.003） 说明：“户主健康状况”参照组为“中”，“最大地块质量”参照组为“中”；*，**，***分别表示通过10%，5%，1%水平下的显著性检验。括号内数值为回归标准误 4. 结论及建议
4.1 结论
本研究发现：样本农户的营林技术效率的平均值为0.57，表明样本农户在当前技术水平下平均56.60%的产出可以通过现有的生产要素组合来获得，样本农户的营林技术效率还有43.40%的提升空间。农户营林生产过程中的雇佣劳动确实会对营林生产造成负面影响，雇工劳动质量不同，对提高农户营林技术效率的影响程度不同；营林生产总投工中雇工所占比例越大，对提高农户营林技术效率的负面影响程度越大。
4.2 建议
建立有效的劳动监督和管理机制。雇佣劳动力在劳动过程中缺乏有效的劳动监督，雇工劳动质量低下，是营林生产技术效率下降的主要原因。因此，农户应根据雇工实际情况建立有效劳动监督机制，在劳动生产可计量的环节尝试使用绩效工资，并根据劳动成果给予一定的激励措施，减少雇工过程中“搭便车”行为的发生，从而提高雇工劳动的质量。在目前林业规模经营日趋普遍的情况下，传统的生产经营和管理方式越来越难以适应，迫切需要建立高效的林业生产管理机制。革新传统林业管理理念，启用具有现代管理才能的人才管理林业经营；在林业生产的各个环节，制定科学合理的管理细则和林业生产流程。
推进适度规模经营，加大林业科技服务投入。解决各林业经营主体营林生产过程中的雇工劳动质量问题的关键是解决农村劳动力不足的问题。因此，推广林业机械化生产、开拓新型经营方式、积极推进林业适度规模经营均能在很大程度上解决上述问题。在地形条件较平缓的地区的林业规模经营户中可以依靠推进农业机械化替代劳动力，解决雇佣劳动所带来的劳动质量问题。同时，在机械替代较为困难的地区，可以尝试开拓林业服务外包、农户之间合作经营和托管经营等新型林业经营方式，促使劳动力要素配置更加专业化，也可以缓解由于雇佣劳动力所带来的劳动质量问题。对规模经营中出现的雇工劳动问题，需要进一步研究以寻找雇工劳动最优的比例，从劳动力层面对农户适度规模经营提出要求。
-
[1] ALCAMO J, BENNETT E M. Ecosystems and Human Well-being:A Frame work for Assessment[M]. Washington D C:Island Press, 2003. [2] BOYD J W, BANZHAF H S. Ecosystem services and government accountability:the need for a new way of judging nature's value[J]. Resources, 2005, 158:16-19. [3] DAILY G C. Nature's Services:Societal Dependence on Natural Ecosystem[M]. Washington D C:Island Press, 1997. [4] HOWARTH R, FARBER S. Accounting for the value of ecosystem services[J]. Ecol Econ, 2002, 41(9):421-429. [5] 赵军, 杨凯.生态系统服务价值评估研究进展[J].生态学报, 2007, 27(1):346-356. ZHAO Jun, YANG Kai. Valuation of ecosystem services:characteristics, issues and prospects[J]. Acta Ecol Sin, 2007, 27(1):346-356. [6] SUTTON P C, COSTANZA R. Global estimates of market and non-market values derived from nighttime satellite imagery, land cover, and ecosystem service valuation[J]. Ecol Econ, 2002, 41(3):509-527. [7] GRêT-REGAMEY A, KYTZIA S. Integrating the valuation of ecosystem services into the input-output economics of an Alpine region[J]. Ecol Econ, 2007, 63(4):786-798. [8] PATTANAYAK S K. Valuing watershed services:concepts and empirics from southeast Asia[J]. Agric Ecos & Environ, 2004, 104(1):171-184. [9] LAL P. Economic valuation of mangroves and decision making in the Pacific[J]. Ocean & Coast Manage, 2003, 46(9/10):823-844. [10] de MENDONCA M J C, SACHSIDA A, LOUREIRO P R A. A study on the valuing of biodiversity:the case of three endangered species in Brazil[J]. Ecol Econ, 2003, 46(1):9-18. [11] 欧阳志云, 王如松.生态系统服务功能、生态价值与机持续发展[J].世界科技研究与发展, 1999, 22(5):45-50. OUYANG Zhiyun, WANG Rusong. Ecosystem services and their economic valuation[J]. World Sci-Tech R & D, 1999, 22(5):45-50. [12] 赵景柱, 徐亚骏, 肖寒, 等.基于可持续发展综合国力的生态系统服务评价研究:13个国家生态系统服务价值的测算[J].系统工程理论与实践, 2003(1):121-127. ZHAO Jingzhu, XU Yajun, XIAO Han, et al. Ecosystem services evaluation based on comprehensive national power for sustainable development:the evaluations on 13 countries[J]. Syst Engineer-theory & Pract, 2003(1):121-127. [13] 赵同谦, 欧阳志云, 王效科, 等.中国陆地地表水生态系统服务功能及其生态经济价值评价[J].自然资源学报, 2003, 18(4):443-452. ZHAO Tongqian, OUYANG Zhiyun, WANG Xiaoke, et al. Ecosystem services and their valuation of terrestrial surface water system in China[J]. J Nat Resour, 2003, 18(4):443-452. [14] 赵同谦, 欧阳志云, 贾良清, 等.中国草地生态系统服务功能间接经济价值评价[J].生态学报, 2004, 24(6):1101-1110. ZHAO Tongqian, OUYANG Zhiyun, JIA Liangqing, et al. Ecosystem services and their valuation of China grassland[J]. Acta Ecol Sin, 2004, 24(6):1101-1110. [15] 赵同谦, 欧阳志云, 郑华, 等.中国森林生态系统服务功能及其价值评价[J].自然资源学报, 2004, 19(4):480-491. ZHAO Tongqian, OUYANG Zhiyun, ZHENG Hua, et al. Forest ecosystem services and their valuation in China[J]. J Nat Resour, 2004, 19(4):480-491. [16] 谢高地, 张钇锂, 鲁春霞, 等.中国自然草地生态系统服务价值[J].自然资源学报, 2001, 16(1):47-53. XIE Gaodi, ZHANG Yili, LU Chunxia, et al. Study on valuation of rangeland ecosystem services of China[J]. J Nat Resour, 2001, 16(1):47-53. [17] 欧阳芳.昆虫的生态服务功能[J].应用昆虫学报, 2013, 50(2):305-310. OUYANG Fang. Insect ecological services[J]. Chin J Appl Entomol, 2013, 50(2):305-310. [18] NEE S. More than meets the eye:earth's real biodiversity is invisible, whether we like it or not[J]. Nature, 2004, 429:804-805. [19] 田华.河南昆虫资源的研究现状及展望[J].贵州农业科学, 2009, 37(9):113-115. TIAN Hua. Advances of studies and prospects on insects resources in Henan Province[J]. Guizhou Agric Sci, 2009, 37(9):113-115. [20] CHEE Y E. An ecological perspective on the valuation of ecosystem services[J]. Biolog Conserv, 2004, 120(4):549-565. [21] 赵景柱, 肖寒, 吴刚.生态系统服务的物质量与价值量评价方法的比较分析[J].应用生态学报, 2000, 11(2):290-292. ZHAO Jingzhu, XIAO Han, WU Gang. Comparison analysis on physical and assessment methods for ecosystem services[J]. Chin J Appl Ecol, 2000, 11(2):290-292. [22] SAS INSTITUTE. SAS/STAT User's Guide, Release 6.12[R]. Cary (NC):SAS Institute, 1996. [23] 白明, 杨星科.蜣螂的生态价值和保护意义[J].应用昆虫学报, 2010, 47(1):39-46. BAI Ming, YANG Xingke. Ecological value and conservation significances of dung beetles[J]. Chin Bull Entomol, 2010, 47(1):39-46. [24] ANDERSON J R, MERRITT R W, LOOMIS E C. The insect-free cattle dropping and its relationship to increased dung fouling of rangeland pastures[J]. J Econ Entomol, 1984, 77(1):133-141. [25] GILLARD P. Coprophagous beetles in pasture ecosystems[J]. J Australian Inst Agric Sci, 1967, 33:30-34. [26] FINCHER G T. The potential value of dung beetles in pasture ecosystems[J]. J Georgia Entomol Soc, 1981, 16:301-316. 期刊类型引用(4)
1. 朱臻，薛家依，宁可. 规模化经营背景下劳动监督对营林质量的影响研究：来自南方集体林区三省规模户的实证数据. 农林经济管理学报. 2021(01): 78-91 . 
百度学术2. 高润，刘洁. 林业经济实现可持续发展的重要性分析. 林产工业. 2021(04): 99-101 . 百度学术3. 肖琦，洪英. 新形势下林业科技对林业发展的重要性分析. 林产工业. 2021(05): 69-71+76 . 百度学术4. 修丕师，张晓梅. 重点国有林区职工林下经济创业潜力及支持路径研究. 林业经济问题. 2020(06): 626-633 . 百度学术其他类型引用(8)
-
-
链接本文:
https://zlxb.zafu.edu.cn/article/doi/10.11833/j.issn.2095-0756.2014.05.018
点击查看大图
计量
- 文章访问数: 3980
- HTML全文浏览量: 674
- PDF下载量: 595
- 被引次数: 12
下载:
