CN116307872B - 生态优先目标下草地生态承载力指数和载畜量获取方法 - Google Patents
生态优先目标下草地生态承载力指数和载畜量获取方法 Download PDFInfo
- Publication number
- CN116307872B CN116307872B CN202310210619.4A CN202310210619A CN116307872B CN 116307872 B CN116307872 B CN 116307872B CN 202310210619 A CN202310210619 A CN 202310210619A CN 116307872 B CN116307872 B CN 116307872B
- Authority
- CN
- China
- Prior art keywords
- grazing
- grassland
- ecological
- bearing capacity
- index
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 244000144972 livestock Species 0.000 title claims abstract description 61
- 238000000034 method Methods 0.000 title claims abstract description 33
- 238000009304 pastoral farming Methods 0.000 claims abstract description 161
- 241001465754 Metazoa Species 0.000 claims abstract description 51
- 241000894007 species Species 0.000 claims abstract description 50
- 238000002474 experimental method Methods 0.000 claims abstract description 32
- 241001494479 Pecora Species 0.000 claims description 38
- 244000025254 Cannabis sativa Species 0.000 claims description 17
- 238000004364 calculation method Methods 0.000 claims description 16
- 238000010606 normalization Methods 0.000 claims description 7
- 230000008569 process Effects 0.000 claims description 6
- 238000012545 processing Methods 0.000 claims description 4
- 230000008859 change Effects 0.000 claims description 3
- 230000001419 dependent effect Effects 0.000 claims description 3
- 238000010586 diagram Methods 0.000 claims description 3
- 238000004519 manufacturing process Methods 0.000 description 5
- 238000011156 evaluation Methods 0.000 description 3
- 239000004459 forage Substances 0.000 description 3
- 230000015556 catabolic process Effects 0.000 description 2
- 238000006731 degradation reaction Methods 0.000 description 2
- 238000009826 distribution Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000009313 farming Methods 0.000 description 2
- 241000270955 Achnatherum calamagrostis Species 0.000 description 1
- 241000234642 Festuca Species 0.000 description 1
- 230000006978 adaptation Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000013459 approach Methods 0.000 description 1
- 230000033228 biological regulation Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 230000008094 contradictory effect Effects 0.000 description 1
- 238000013480 data collection Methods 0.000 description 1
- 238000013075 data extraction Methods 0.000 description 1
- 230000003247 decreasing effect Effects 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 238000013401 experimental design Methods 0.000 description 1
- 230000005764 inhibitory process Effects 0.000 description 1
- 238000011835 investigation Methods 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000012423 maintenance Methods 0.000 description 1
- 239000000463 material Substances 0.000 description 1
- 238000010197 meta-analysis Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 238000012544 monitoring process Methods 0.000 description 1
- 238000005457 optimization Methods 0.000 description 1
- 238000012946 outsourcing Methods 0.000 description 1
- 238000009331 sowing Methods 0.000 description 1
- 238000007619 statistical method Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000009897 systematic effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q10/00—Administration; Management
- G06Q10/06—Resources, workflows, human or project management; Enterprise or organisation planning; Enterprise or organisation modelling
- G06Q10/063—Operations research, analysis or management
- G06Q10/0639—Performance analysis of employees; Performance analysis of enterprise or organisation operations
- G06Q10/06393—Score-carding, benchmarking or key performance indicator [KPI] analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F16/00—Information retrieval; Database structures therefor; File system structures therefor
- G06F16/20—Information retrieval; Database structures therefor; File system structures therefor of structured data, e.g. relational data
- G06F16/29—Geographical information databases
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06F—ELECTRIC DIGITAL DATA PROCESSING
- G06F17/00—Digital computing or data processing equipment or methods, specially adapted for specific functions
- G06F17/10—Complex mathematical operations
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06Q—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES; SYSTEMS OR METHODS SPECIALLY ADAPTED FOR ADMINISTRATIVE, COMMERCIAL, FINANCIAL, MANAGERIAL OR SUPERVISORY PURPOSES, NOT OTHERWISE PROVIDED FOR
- G06Q50/00—Information and communication technology [ICT] specially adapted for implementation of business processes of specific business sectors, e.g. utilities or tourism
- G06Q50/02—Agriculture; Fishing; Forestry; Mining
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A40/00—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production
- Y02A40/70—Adaptation technologies in agriculture, forestry, livestock or agroalimentary production in livestock or poultry
Landscapes
- Engineering & Computer Science (AREA)
- Business, Economics & Management (AREA)
- Physics & Mathematics (AREA)
- Human Resources & Organizations (AREA)
- Theoretical Computer Science (AREA)
- General Physics & Mathematics (AREA)
- Data Mining & Analysis (AREA)
- Databases & Information Systems (AREA)
- Economics (AREA)
- Strategic Management (AREA)
- Marketing (AREA)
- Educational Administration (AREA)
- Tourism & Hospitality (AREA)
- Mathematical Physics (AREA)
- General Business, Economics & Management (AREA)
- General Engineering & Computer Science (AREA)
- Entrepreneurship & Innovation (AREA)
- Development Economics (AREA)
- General Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Marine Sciences & Fisheries (AREA)
- Mining & Mineral Resources (AREA)
- Health & Medical Sciences (AREA)
- Quality & Reliability (AREA)
- Primary Health Care (AREA)
- Agronomy & Crop Science (AREA)
- Remote Sensing (AREA)
- Animal Husbandry (AREA)
- Game Theory and Decision Science (AREA)
- Algebra (AREA)
- Computational Mathematics (AREA)
- Mathematical Analysis (AREA)
- Mathematical Optimization (AREA)
- Operations Research (AREA)
- Pure & Applied Mathematics (AREA)
- Software Systems (AREA)
- Management, Administration, Business Operations System, And Electronic Commerce (AREA)
Abstract
本发明公开了一种生态优先目标下草地生态承载力指数和载畜量获取方法。所述草地生态承载力指数获取方法包括：通过放牧梯度实验获得区域不同草地类型梯度的放牧强度下不同放牧小区的放牧数据及物种多样性数据；计算不同同一草地类型内不同放牧小区的基础载畜压力指数，并将基础载畜压力指数和归一化后的物种多样性数据进行曲线拟合，根据拟合曲线的极大值确定不同类型草地的的生态承载力指数。依据区域不同类型草地的的生态承载力指数，结合各草地类型面积、单位面积产草量，汇算获取生态优先目标下区域的生态载畜量。本发明获取的草地生态承载能力指数以生态保护优先，可更准确地评价保护地区域草地的适宜承载能力。
Description
技术领域
本发明涉及草地载畜能力评估方法的技术领域。
背景技术
现有技术中通常以草地不发生退化的约束下，其产草量所能承载的最大羊单位(SU)数确定草地的适宜载畜量。即在进行草地畜牧管理时，以畜牧生产能力作为优先的管理目标，其采用的管理指标一般为草地载畜压力指数Ip,即以Ip为1作为草地畜牧生产的管理目标，来经营调控生产畜牧产品的草地上放养的牲畜数量。作为一种衡量草地放牧对生态系统植被生产力的影响及草畜矛盾特征的指数，草地载畜压力指数的计算方式为：Ip＝Cs/Cl，其中，Cs为草地现实载畜压力,即单位面积草地实际承载的羊单位数量(SU/hm2)，其计算方式为：Cs＝(Cn×(1+Ch)×Gt)/(Ar×365)，其中，Cn为年末家畜存栏数，按标准羊单位计算，其他牲畜换算为标准羊单位，Ch为家畜出栏率,根据目标区域的实际情况而定；Ar为草地面积(hm2)，Gt为草地放牧时间；Cl为草地理论承载力，即单位面积草地在不发生退化前提下，适宜承载的羊单位(SU/hm2)，其计算方式为：Cl＝(Ym×Ut×Co×Ha)/(Sf×Gt)，其中，Ym为单位面积草地的产草量(kg/hm2)；Ut为牧草利用率，Co为草地可利用率，Ha为草地可食牧草比率；Gt为草地放牧时间(d)；Sf为一个羊单位家畜的日食量，按照每羊单位每天采食1.8kg干草计算。
在以草地载畜压力指数Ip为1作为草地畜牧的管理标准时，其实际对应了以下的畜牧管理目标：当天然草地和其他途径如人工草地、农作物、外购等提供的饲草料总量等于牲畜所需的饲草料总量时，草地生长的牧草恰好足够供给家畜生长，即达到生产力管理下的草畜平衡。可以看出，该标准并未考虑到随着草地长期处于Ip为1的放牧强度条件下时，草地的生态系统可能出现功能减弱、物种多样性下降等问题。因而原有管理方式不能长期维持草地最佳的物种多样性水平、生态服务等最优生态状况。同时完全禁牧的草地管理方式又不符合生态系统中度干扰理论，即只有在适度干扰和人类生产活动影响下生态系统才能保持最大活力，保有最优生态功能和物种多样性状态。如何准确确定草地维持最优的生物多样性和生态状态的放牧牲畜数量，进而管理以生态优先作为目标的草地区域尚未有科学系统的技术方法。
已有的草地载畜能力评估管理技术方法尤其不适用于需要以保护生物多样性、维持或恢复优良生态状态为主要目标的区域，如生态红线区、国家公园或自然保护区。因此，亟需获得一种新的以生态保护、生物多样性最优作为优先管理目标的草地载畜能力和区域载畜量的评价指标和方法，以对国家公园或自然保护区等具有生态保护要求的区域进行畜牧管理。
发明内容
针对现有技术的缺陷，本发明的目的在于提出一种以生态保护作为优先管理目标的、可更科学准确地评价草地的生态承载能力和区域生态载畜量的草地生态承载力指数和载畜量的获取方法。
本发明的技术方案如下：
一种生态优先目标下草地生态承载力指数的获取方法，其包括：
S1核算需要进行草地生态承载力评估的区域内的草地类型，在不同草地类型上各设置1个放牧梯度实验区，在各放牧梯度实验区内布设若干放牧小区，在放牧小区上按梯度的放牧强度进行放牧实验，即放牧梯度实验，和/或，通过技术资料获得需要通过所述放牧梯度实验获得的数据；
S2通过所述放牧梯度实验和/或所述技术资料获得放牧强度变化对各草地类型的生态影响的数据，包括记录放牧过程的放牧数据及放牧后各放牧小区的物种多样性数据；
S3基于所述放牧数据，计算各草地类型的放牧梯度实验区下不同放牧小区的载畜压力指数，即基础载畜压力指数Ip(i)；
S4对不同草地类型下各放牧小区的物种多样性数据进行归一化处理，得到归一化后的物种多样性数据；
S5对同一草地类型下不同放牧小区的基础载畜压力指数Ip(i)和所述归一化后的物种多样性数据进行曲线拟合，根据拟合曲线确定不同类型草地的草地生态承载力指数Ip-s(i)。
以上技术方案以生态优先作为约束，可获得更准确且更能实现可持续畜牧的畜牧方案。
根据本发明的一些优选实施方式，所述放牧梯度实验数据包括：放牧区域的草地类型、经纬度、海拔(m)、多年平均产草量(kg/hm2)，放牧区面积(hm2)、放牧持续时间(天/年)、放牧强度(只/hm2)、放牧区域牲畜种类及其数量。
根据本发明的一些优选实施方式，所述物种多样性数据选自丰富度指数、均匀度指数、优势度指数中的一种或多种。
根据本发明的一些优选实施方式，所述放牧数据中，不同种类的牲畜数据均换算为以羊单位为标准单位的标准化数据。
根据本发明的一些优选实施方式，所述基础载畜压力指数Ip(i)通过以下计算模型得到：
Ip(i)＝Cs(i)/Cl(i)
Cs(i)＝(Cn(i)×Gt(i))/(Ar(i)×365)
Cl(i)＝(Ym(i)×Ut(i)×Co(i)×Ha(i))/(Sf×Gt(i))
其中，Cs(i)为第i个放牧小区的草地现实载畜压力、即该放牧小区的单位面积实际承载的羊单位数量(SU/hm2)，Cl(i)为第i个放牧小区的草地理论承载力、即该放牧小区的单位面积草地适宜承载的羊单位数量(SU/hm2)；Cn(i)为第i个放牧小区实际承载的羊单位数量(SU)；Ar(i)为第i个放牧小区的草地面积(hm2)；Gt(i)为第i个放牧小区的实际放牧时间；Ym(i)为第i个放牧小区的单位面积干草产草量(kg/hm2)；Ut(i)为第i个放牧小区的牧草利用率；Co(i)为第i个放牧小区的草地可利用率；Ha(i)为第i个放牧小区的草地可食牧草比率；Sf为一个羊单位家畜的日食量。
以上优选实施方式可通过草地载畜压力指数Ip量化不同放牧梯度，将放牧牲畜的采食压力量化为标准的可跨区域比较的草地载畜压力指数。
在更具体的一些实施方式中，草地放牧时间Gt(i)可按全年放牧365d计算；一个羊单位家畜的日食量Sf可按照每个羊单位每天采食1.8kg干草计算。
根据本发明的一些优选实施方式，所述归一化处理通过以下计算模型实现：
其中，表示归一化后的物种多样性数据，xj表示第j个放牧小区的物种多样性指标；MaxXj和MinXj分别表示第j个放牧小区对应的草地类型下全部放牧小区中物种多样性指标的最大值和最小值，xjmean表示其物种多样性指标的平均值，xjstd表示其物种多样性指标的标准差。
本发明通过如以上的归一化处理，可使在同一草地类型下不同多样性本底水平的区域获得的物种多样性数据实现跨区域可比较。
根据本发明的一些优选实施方式，所述曲线拟合包括：以同一草地类型下不同放牧小区的基础载畜压力指数Ip(i)作为自变量X，该放牧小区的归一化后的物种多样性数据作为因变量Y，作出各归一化后的物种多样性数据随不同草地类型的基础载畜压力指数Ip(i)变化的趋散点图，并利用lm()和/或poly()函数进行曲线拟合，获得曲线对应的一元多项式回归方程。
根据本发明的一些优选实施方式，所述不同类型草地的草地生态承载力指数Ip-s(i)为通过optimize()函数求取得到的所述拟合曲线的极大值。
本发明进一步提供了一种生态优先目标下草地生态载畜量的获取方法，其包括以上步骤S1-S5，及：
S6：根据以下计算模型获得含有不同草地类型的待评估区域的草地生态载畜量：
Cp-s＝ΣCp-s(i)
Cp-s(i)＝Ci·Art(i)
Ci＝Ip-s(i)·Clt(i)
其中，Cp-s为生态优先目标下的草地生态载畜量(SU)；Cp-s(i)为生态优先目标下该区域内第i种草地类型的草地生态载畜量；Ci为生态优先目标下该区域内第i种草地类型的草地生态承载力(SU/hm2)；Clt(i)为第i种草地类型的草地理论承载力、即该草地类型的单位面积草地适宜承载的羊单位数量(SU/hm2)；Art(i)为该区域第i种草地类型的草地面积(hm2)。
以上获取方法可分草地类型汇算草地生态承载力Ci，通过汇算各草地类型的草地生态载畜量Cp-s(i)获取目标区域总草地生态载畜量Cp-s，获得准确的生态载畜量。
基于本发明的草地生态承载力指数获取方法可获得准确衡量以生物多样性保护为先的草地适宜的生态载畜量。在草地生态承载力指数获取中，可根据中度干扰理论(在适度的干扰下，物种多样性最高)，设计放牧梯度实验，综合考虑了各草地类型的特点，得到各草地类型的草地生态承载力指数，进而基于各草地类型的生态承载力指数汇算得到区域草地生态载畜量。
本发明以各草地类型下的物种多样性最优作为约束条件，可获取以生态保护为优先目标的管理区域的最适宜的载畜规模或载畜数量，特别是可以不同草地类型下多样性拐点出现时的草地载畜压力指数作为汇算依据，进一步在管理中获得草地放牧承载力的关键参数，并基于获取的关键参数，分草地类型汇算目标区域生态优先管理目标下草地的生态承载力和生态载畜量。
附图说明
图1为本发明的草地生态承载力指数和草地生态载畜量获取方法的流程示意图。
图2为本发明所述放牧梯度实验的一种具体实施例。
图3为实施例1中待评估区域的草地类型空间分布图。
图4为实施例1中求取拐点的示例图。
具体实施方式
以下结合实施例和附图对本发明进行详细描述，但需要理解的是，所述实施例和附图仅用于对本发明进行示例性的描述，而并不能对本发明的保护范围构成任何限制。所有包含在本发明的发明宗旨范围内的合理的变换和组合均落入本发明的保护范围。
参照附图1，根据本发明的技术方案，在一些具体实施方式中，草地生态承载力指数获取方法可包括：
S1识别区域内重要草地/植被类型，在各典型草地类型设置放牧梯度实验和/或通过文献资料获取对应数据；
S2通过放牧梯度实验和/或文献资料进行数据收集和提取，获得畜牧对放牧区域的生态影响的数据，即基础实验数据，包括放牧数据和放牧后各草地类型的物质多样性数据；
其中，放牧梯度实验为在梯度的放牧强度下进行的放牧实验，其一种具体的实施方式如：在具有不同草地类型的放牧区域内设计梯度的放牧强度，在毗邻的不同放牧小区之间以铁丝网或其它围栏方式隔绝牲畜的自行跨越迁移。对特定载畜数量的牲畜进行长期限定采食区域的放牧实验，记录放牧数据并最终通过样方调查等统计方法获取放牧梯度实验各放牧小区的生物多样性变化情况。
更具体的，所述放牧数据可包括：放牧区域的草地类型、经纬度、海拔(m)、多年平均产草量(kg/hm2)，放牧小区面积(hm2)、放牧持续时间(天/年)、放牧强度(只/hm2)、放牧牲畜种类及其数量。
所述物种多样性数据优选为多样性指数，包括如丰富度指数、均匀度指数、优势度指数等指数中的一种或多种，且不同种类的牲畜可按照NY/T 635-2015《天然草地合理载畜量计算》中提出的不同种类牲畜折算为羊单位进行换算的方式，转换为标准羊单位/hm2。
更具体的，参照附图2，在一种具体实施例中，所述放牧梯度实验可设计如：
按评估区域内的N个典型草地类型划定放牧梯度实验区N个，每块1km×1km，各实验区内划定边长为150m的正方形实验小区即放牧小区36个，每个实验小区建立牲畜不能翻越的围栏边界，内部可设置一个10×10m的禁牧空白对照方格，以0羊单位放牧压力为起始，设立12个放牧压力梯度，随机分配施行至36个实验小区中(每个放牧梯度同时施行于3个小区作为重复)。在各放牧小区内设置的1个禁牧空白对照方格边界建立牲畜无法翻越的围栏，确保内部植被不被采食，以该对照方格内测定获得草地年产草量作为无采食情况下全年总产草量。
根据以上实验设计，结合各小区放牧牲畜产生的梯度载畜压力，可以测算各实验小区内的基础草地载畜压力指数Ip(i)。同时，监测中获取各实验小区内的物种多样性数据，用于后续分析。
另一方面，若通过文献资料提取上述所需数据用于后续计算，则在缺乏放牧对照样方的情况下，可用放牧强度为0的放牧小区的产草量数据或临近禁牧区域产草量数据替代。
S3基于所述基础实验数据，计算并提取出不同放牧强度下的不同类型的草地的载畜压力指数，获得基础载畜压力指数Ip(i)；
进一步的，所述基础载畜压力指数的获取如下：
Ip(i)＝Cs(i)/Cl(i)
Cs(i)＝(Cn(i)×Gt(i))/(Ar(i)×365)
Cl(i)＝(Ym(i)×Ut(i)×Co(i)×Ha(i))/(Sf×Gt(i))
其中，Cs(i)为第i个放牧小区的草地现实载畜压力、即该放牧小区的单位面积实际承载的羊单位数量，Cl(i)为第i个放牧小区的草地理论承载力、即该放牧小区的单位面积草地适宜承载的羊单位数量；Cn(i)为第i个放牧小区实际承载的羊单位数量；Ar(i)为第i个放牧小区的草地面积；Gt(i)为第i个放牧小区的实际放牧时间，其可按全年放牧356d计算；Ym(i)为第i个放牧小区的单位面积干草产草量；Ut(i)为第i个放牧小区的牧草利用率；Co(i)为第i个放牧小区的草地可利用率；Ha(i)为第i个放牧小区的草地可食牧草比率；Sf为一个羊单位家畜的日食量，其可按照每个羊单位每天采食1.8kg干草计算。
S4对同一草地类型下各放牧小区的物种多样性数据进行归一化处理，得到归一化后的如物种多样性数据；
进一步的，所述归一化处理可通过如下的处理模型实现：
其中，表示归一化后的物种多样性数据，xj表示第j个放牧小区的物种多样性指标；MaxXj和MinXj分别表示第j个放牧小区对应的草地类型下全部放牧小区中物种多样性指标的最大值和最小值，xjmean表示其物种多样性指标的平均值，xjstd表示其物种多样性指标的标准差。
S5根据所述基础载畜压力指数Ip和所述归一化后的物种多样性数据，确定不同类型草地的草地生态承载力指数Ip-s(i)；
进一步的，在一些具体实施方式中，所述草地生态承载力指数Ip-s(i)通过拐点数据确定，如下：
以所得基础载畜压力指数Ip(i)作为自变量X，所得归一化后的物种多样性数据作为因变量Y，作出各归一化后的物种多样性数据随不同放牧强度下的基础载畜压力指数Ip(i)的变化的趋散点图，并利用如lm()和/或poly()函数进行曲线拟合，获得曲线对应的一元多项式回归方程；其后利用如optimize()函数求取该方程的极大值，即拐点值，以该拐点值作为类型i的草地的所述草地生态承载力指数Ip-s(i)。
S6根据所述不同类型草地的草地生态承载力指数Ip-s(i)及其对应的草地基础指数确定不同类型草地的草地生态承载力Ci及待评算区域的总草地生态承载量指数Cp-s；其中，所述草地基础指数包括不同草地类型的产草量均值及年均放牧时间；
进一步的，所述不同草地类型的草地生态承载力Ci，不同草地类型的草地生态载畜量Cp-s(i)及待评估区域的总草地生态载畜量Cp-s确定方式如下：
Cp-s＝ΣCp-s(i)
Cp-s(i)＝Ci·Art(i)
Ci＝Ip-s(i)·Clt(i)
其中，Cp-s为生态优先目标下的草地生态载畜量(SU)；Cp-s(i)为生态优先目标下该区域内第i种草地类型的草地生态载畜量；Ci为生态优先目标下该区域内第i种草地类型的草地生态承载力(SU/hm2)；Clt(i)为第i种草地类型的草地理论承载力、即该草地类型的单位面积草地适宜承载的羊单位数量(SU/hm2)；Art(i)为该区域第i种草地类型的草地面积(hm2)。
本发明以该草地生态承载力Ci及目标区域的总的草地生态载畜量Cp-s，而非直接的以区域草地载畜压力指数Ip＝1作为区域畜牧管理的管理指标，可得到对区域生态状况维持提升角度来说更佳的管理方案，从而通过放牧压力调控管理获取并维持区域最优的生态状态。
实施例1
待进行草地生态承载力评估的区域具有如附图3所示的草地类型空间分布，其主要草地类型包括高寒草原和高寒草甸，其中高寒草原占该区域总面积的36％，主要分布在中部、北部地区；高寒草甸占该区域总面积的52％，主要集中连片分布在南部，并覆盖中部部分区域。该区域内草地型多样，其中覆盖面积最大的草地型为紫花针茅草地，占总面积的26％，主要分布在东北部；其次为紫花针茅+杂类草草地型、高山嵩草草地型、藏嵩草草地型和藏嵩草+苔草草地型等，分别占总面积的13％、12％、11％和9％。
在该区域内的典型草地类型高寒灌丛、高寒草原和高寒草甸上设置放牧梯度实验或收集已有放牧梯度实验数据资料，根据本发明具体实施方式提供的方法获取各草地类型在生态优先目标下的生态承载力指数Ip-s(i)，即拟合基础草地载畜压力指数-归一化多样性指数的曲线并求取拐点(其中，在高寒草甸区域求取拐点的过程如附图4所示)，其后分别计算得到各个草地类型的生态优先目标下的草地生态承载力Ci(SU/hm2)，结合各草地类型面积和多年平均产草量，得到各个草地类型在生态优先目标下的草地生态载畜量Cp-s(i)(SU)，进而在区域上通过汇算各草地类型草地生态载畜量求得该区域的总生态载畜量Cp-s，作为区域牲畜放牧管理的依据，如在偏离生态载畜量的地方，通过增加或减少区域的载畜压力管理区域生态载畜量，从而维持区域生态状态和物种多样性处于最优状态。
依据本方案提出的生态优先目标下的草地生态承载力指数(即Ips)，计算得到该区域生态优先目标下平均草地生态承载力为0.05SU/hm2，全区域在生态优先目标下可承载的牲畜数量为9.72万羊单位。其中，高寒灌丛可养畜1.29万羊单位，高寒草原可养畜1.3万羊单位，高寒草甸可养畜7.13万羊单位，具体计算情况如下表1所示。
表1.生态优先目标下草地适宜载畜量汇算
而依据以往的生产优先目标的草地理论承载力指数(即Ip＝1)计算可得，维持草地不退化的生产优先目标下，平均草地承载力为0.12SU/hm2，全县可养畜26.24万羊单位。其中，高寒灌丛可养畜0.81万羊单位，高寒草原可养畜2.27万羊单位，高寒草甸可养畜23.15万羊单位。
可以看出，在生态优先管理目标约束下，为维持区域物种多样性最优状态，相比原有生产可持续(Ip＝1)状态，区域放牧牲畜数量需要减少58％。
实施例2
根据具体实施方式提供的计算方法，通过文献资料数据荟萃分析的方式，获得分布于我国的主要草地类型的草地生态承载力指数Ip-s，如下表2所示。
表2.不同草地类型的指数对比
通过以上结果可知，温性草甸草原类草地在载畜载畜压力指数为0.50时，可使草地物种多样性维持在最优状态；高寒草甸类草地在载畜载畜压力指数为0.31时，可使草地物种多样性维持在最优状态；高寒草原类草地在载畜载畜压力指数为0.57时，可使草地物种多样性维持在最优状态；高寒灌丛类草地在载畜载畜压力指数为1.59时，可使草地物种多样性维持在最优状态；高寒混播人工草地在载畜载畜压力指数为0.40时，可使草地物种多样性维持在最优状态；温性荒漠草原类草地在载畜载畜压力指数为0.00时，可使草地物种多样性维持在最优状态；典型温性草原类草地在载畜载畜压力指数为0.13时，可使草地物种多样性维持在最优状态。
以上实施例仅是本发明的优选实施方式，本发明的保护范围并不仅局限于上述实施例。凡属于本发明思路下的技术方案均属于本发明的保护范围。应该指出，对于本技术领域的普通技术人员来说，在不脱离本发明原理的前提下的改进和润饰，这些改进和润饰也应视为本发明的保护范围。
Claims (8)
1.一种生态优先目标下草地生态承载力指数的获取方法，其特征在于，其包括：
S1核算需要进行草地生态承载力评估的区域内的草地类型，在不同草地类型上各设置1个放牧梯度实验区，在各放牧梯度实验区内布设若干放牧小区，在放牧小区上按梯度的放牧强度进行放牧实验，即放牧梯度实验，和/或，通过技术资料获得需要通过所述放牧梯度实验获得的数据；
S2通过所述放牧梯度实验和/或所述技术资料获得放牧强度变化对各草地类型的生态影响的数据，包括记录放牧过程的放牧数据及放牧后各放牧小区的物种多样性数据；
S3基于所述放牧数据，计算各草地类型下不同放牧小区的载畜压力指数，即基础载畜压力指数Ip(i)；
S4对不同草地类型下各放牧小区的物种多样性数据进行归一化处理，得到归一化后的物种多样性数据；
S5对同一草地类型下不同放牧小区的基础载畜压力指数Ip(i)和所述归一化后的物种多样性数据进行曲线拟合，根据拟合曲线确定不同类型草地的草地生态承载力指数Ip-s(i)；
其中，所述归一化处理通过以下计算模型实现：
其中，表示归一化后的物种多样性数据，xj表示第j个放牧小区的物种多样性指标；MaxXj和MinXj分别表示第j个放牧小区对应的草地类型下全部放牧小区中物种多样性指标的最大值和最小值，xjmean表示其物种多样性指标的平均值，xjstd表示其物种多样性指标的标准差；
所述曲线拟合包括：以同一草地类型下不同放牧小区的基础载畜压力指数Ip(i)作为自变量X，该放牧小区的归一化后的物种多样性数据作为因变量Y，同一草地类型下作出各归一化后的物种多样性数据随基础载畜压力指数Ip(i)变化的趋散点图，并利用lm()和/或poly()函数进行曲线拟合，获得曲线对应的一元多项式回归方程。
2.根据权利要求1所述的生态优先目标下草地生态承载力指数的获取方法，其特征在于，所述放牧数据包括：草地类型、经纬度、海拔、多年平均产草量，草地放牧区域面积、放牧持续时间、放牧强度、牲畜种类及其数量。
3.根据权利要求1所述的生态优先目标下草地生态承载力指数的获取方法，其特征在于，所述物种多样性数据选自丰富度指数、均匀度指数、优势度指数中的一种或多种。
4.根据权利要求1所述的生态优先目标下草地生态承载力指数的获取方法，其特征在于，所述放牧数据中，不同种类的牲畜数据均换算为以羊单位为标准单位的标准化数据。
5.根据权利要求1所述的生态优先目标下草地生态承载力指数的获取方法，其特征在于，所述基础载畜压力指数Ip(i)通过以下计算模型得到：
Ip(i)＝Cs(i)/Cl(i)
Cs(i)＝(Cn(i)×Gt(i))/(Ar(i)×365)
Cl(i)＝(Ym(i)×Ut(i)×Co(i)×Ha(i))/(Sf×Gt(i))
其中，Cs(i)为第i个放牧小区的草地现实载畜压力、即该放牧小区的单位面积实际承载的羊单位数量，Cl(i)为第i个放牧小区的草地理论承载力、即该放牧小区的单位面积草地适宜承载的羊单位数量；Cn(i)为第i个放牧小区实际承载的羊单位数量；Ar(i)为第i个放牧小区的草地面积；Gt(i)为第i个放牧小区的实际放牧时间；Ym(i)为第i个放牧小区的单位面积干草产草量；Ut(i)为第i个放牧小区的牧草利用率；Co(i)为第i个放牧小区的草地可利用率；Ha(i)为第i个放牧小区的草地可食牧草比率；Sf为一个羊单位家畜的日食量。
6.根据权利要求1所述的生态优先目标下草地生态承载力指数的获取方法，其特征在于，所述不同类型草地的草地生态承载力指数Ip-s(i)为通过optimize()函数求取得到的所述拟合曲线的极大值对应的横坐标值。
7.根据权利要求1-6中任一项所述的生态优先目标下草地生态承载力指数的获取方法，其特征在于，所述放牧梯度实验还包括在不同草地类型上设置的若干放牧区域之外，另外设置至少一个禁牧的空白对比区域，根据该空白对比区域内测得的草地年产草量计算该草地类型的全年总产草量。
8.一种生态优先目标下草地生态载畜量的获取方法，其特征在于，其包括权利要求1-7中任一项所述的S1-S5，及以下的步骤S6：
S6根据以下计算模型获得含有不同草地类型的待评估区域的草地生态载畜量：
Cp-s＝ΣCp-s(i)
Cp-s(i)＝Ci·Art(i)
Ci＝Ip-s(i)·Clt(i)
其中，Cp-s为生态优先目标下的草地生态载畜量；Cp-s(i)为生态优先目标下该区域内第i种草地类型的草地生态载畜量；Ci为生态优先目标下该区域内第i种草地类型的草地生态承载力；Clt(i)为第i种草地类型的草地理论承载力、即该草地类型的单位面积草地适宜承载的羊单位数量；Art(i)为该区域第i种草地类型的草地面积。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310210619.4A CN116307872B (zh) | 2023-03-07 | 2023-03-07 | 生态优先目标下草地生态承载力指数和载畜量获取方法 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202310210619.4A CN116307872B (zh) | 2023-03-07 | 2023-03-07 | 生态优先目标下草地生态承载力指数和载畜量获取方法 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN116307872A CN116307872A (zh) | 2023-06-23 |
| CN116307872B true CN116307872B (zh) | 2023-10-03 |
Family
ID=86800816
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202310210619.4A Active CN116307872B (zh) | 2023-03-07 | 2023-03-07 | 生态优先目标下草地生态承载力指数和载畜量获取方法 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN116307872B (zh) |
Families Citing this family (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN116882694A (zh) * | 2023-07-20 | 2023-10-13 | 中国科学院植物研究所 | 一种高频轮牧方法 |
| CN116957207B (zh) * | 2023-09-19 | 2023-11-28 | 成都信息工程大学 | 一种实时数据采集的草地生态健康监测方法 |
Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105117968A (zh) * | 2015-08-11 | 2015-12-02 | 南京大学连云港高新技术研究院 | 海岸带生态承载力分析方法 |
| CN105893774A (zh) * | 2016-04-25 | 2016-08-24 | 环境保护部南京环境科学研究所 | 生物多样性保护优先区域绿色发展生态安全评价指数的计算机建模方法 |
| CN106599576A (zh) * | 2016-12-13 | 2017-04-26 | 中国科学院、水利部成都山地灾害与环境研究所 | 基于结构动力学的草地生态承载力定量方法 |
| CN113011235A (zh) * | 2019-12-21 | 2021-06-22 | 中国林业科学研究院资源信息研究所 | 基于cpugpu异构平台的多源遥感草地草畜平衡评估方法 |
| CN113887974A (zh) * | 2021-10-11 | 2022-01-04 | 华中师范大学 | 一种基于多源时空数据集成的城市资源环境承载力综合评价方法 |
-
2023
- 2023-03-07 CN CN202310210619.4A patent/CN116307872B/zh active Active
Patent Citations (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN105117968A (zh) * | 2015-08-11 | 2015-12-02 | 南京大学连云港高新技术研究院 | 海岸带生态承载力分析方法 |
| CN105893774A (zh) * | 2016-04-25 | 2016-08-24 | 环境保护部南京环境科学研究所 | 生物多样性保护优先区域绿色发展生态安全评价指数的计算机建模方法 |
| CN106599576A (zh) * | 2016-12-13 | 2017-04-26 | 中国科学院、水利部成都山地灾害与环境研究所 | 基于结构动力学的草地生态承载力定量方法 |
| CN113011235A (zh) * | 2019-12-21 | 2021-06-22 | 中国林业科学研究院资源信息研究所 | 基于cpugpu异构平台的多源遥感草地草畜平衡评估方法 |
| CN113887974A (zh) * | 2021-10-11 | 2022-01-04 | 华中师范大学 | 一种基于多源时空数据集成的城市资源环境承载力综合评价方法 |
Also Published As
| Publication number | Publication date |
|---|---|
| CN116307872A (zh) | 2023-06-23 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN116307872B (zh) | 生态优先目标下草地生态承载力指数和载畜量获取方法 | |
| Euclides et al. | Maintaining post‐grazing sward height of Panicum maximum (cv. Mombaça) at 50 cm led to higher animal performance compared with post‐grazing height of 30 cm | |
| Euclides et al. | Grazing intensity affects forage accumulation and persistence of Marandu palisadegrass in the Brazilian savannah | |
| Rurinda et al. | Managing soil fertility to adapt to rainfall variability in smallholder cropping systems in Zimbabwe | |
| Tamele et al. | Optimal defoliation management of brachiaria grass–forage peanut for balanced pasture establishment | |
| Leeuw et al. | Sheep-grazing with different foraging efficiences in a Dutch mixed grassland | |
| McRoberts et al. | Urea and composted cattle manure affect forage yield and nutritive value in sandy soils of south‐central Vietnam | |
| Scheneiter et al. | Tiller demography in tall fescue (Festuca arundinacea) swards as influenced by nitrogen fertilization, sowing method and grazing management | |
| Ren et al. | Grazing effects on herbage nutritive values depend on precipitation and growing season in Inner Mongolian grassland | |
| Pereira et al. | Grazing management and tussock distribution in elephant grass | |
| Paraiso et al. | Herbage accumulation, nutritive value, and organic reserves of continuously stocked ‘Ipyporã’and ‘Mulato II’brachiariagrasses | |
| Sternberg et al. | Effects of grazing on Bituminaria bituminosa (L) Stirton: a potential forage crop in Mediterranean grasslands | |
| Li et al. | The long-term effects of thinning and mixing on species and structural diversity of Chinese fir plantations | |
| Zhang et al. | Evaluation of maize varieties for mechanical grain harvesting in mid‐latitude region, China | |
| Rebollo et al. | Refuge effects of a cactus in grazed short‐grass steppe | |
| Miqueloto et al. | Herbage utilisation efficiency of continuously stocked pastures during periods of restricted pasture growth | |
| Antunes et al. | Intermittent stocking strategies for the management of Marandu palisade grass in the Brazilian Cerrado biome | |
| Santos et al. | Forage accumulation in brachiaria grass under continuous grazing with single or variable height during the seasons of the year | |
| Interrante et al. | Bahiagrass tiller dynamics in response to defoliation management | |
| Kuinchtner et al. | Thermal sum (degrees-days) rest period as a grazing management tool in natural grasslands: Effects on animal performance and forage production | |
| Biligetu et al. | Quantifying regrowth characteristics of three bromegrass (Bromus) species in response to defoliation at different developmental stages | |
| Alberto Depablos Alviarez et al. | Managing “Marandu” palisadegrass and calopo pastures based on light interception | |
| Pluess | Meta-analysis reveals microevolution in grassland plant species under contrasting management | |
| Gatti et al. | Dynamics of structural traits in two competing C3 grass species: influence of neighbours and nitrogen | |
| Song et al. | Change in seed dormancy status controls seasonal timing of seed germination in Korean pine (Pinus koraiensis) |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| PB01 | Publication | ||
| PB01 | Publication | ||
| SE01 | Entry into force of request for substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |