基于Maxent模型的互花米草潜在生境适宜性分析

第31卷 第11期
V o l .31 N
o .11草 地 学 报
A C T A A G R E S T I A S I N I C A
2023年 11月
N o v . 2023
d o i :10.11733/j
.i s s n .1007-0435.2023.11.030引用格式:方 铧,陈星彤,刘明月,等.基于M a x e n t 模型的互花米草潜在生境适宜性分析[J ].草地学报,2023,31(11):3514-3524
F A N
G
H u a ,C H E NX i n g -t o n g ,L
I U M i n g -y u e ,e t a l .A n a l y s i s o n t h e P o t e n t i a lH a b i t a t S u i t a b i l i t y o f S p
a r t i n a a l t e r -n i f
l o r a B a s e do n M a x e n tM o d e l [J ].A c t aA g r e s t i aS i n i c a ,2023,31(11):3514-3524基于M a x e n t 模型的互花米草潜在生境适宜性分析
方 铧1,
2,陈星彤1,刘明月1,2,3*,张永彬1,苗正红4,满卫东1,2,3,张清文1,2,寇财垚1,2,李 想1,2
(1.华北理工大学矿业工程学院,河北唐山063210;2.唐山市资源与环境遥感重点实验室,河北唐山063210;
3.河北省矿区生态修复产业技术研究院,河北唐山063210;
4.吉林省水利水电勘测设计研究院,吉林长春130033)收稿日期:2023-05-08;修回日期:2023-06-26
基金项目:国家自然科学基金青年基金项目(N o .42101393,N o .41901375);河北省自然科学基金项目(N o .D 2019209322,N o .
D 2022209005);唐山市科技计划重点研发项目(N o .22150221J );吉林省科技发展计划项目(N o .20210203028S F )资助作者简介:方铧(1998-),男,汉族,河北沧州人,硕士研究生,主要从事滨海湿地生态遥感研究,
E -m a i l :f a n g
h @s t u .n c s t .e d u .c n ;*通信作者A u t h o r f o r c o r r e s p o n d e n c e ,E -m a i l :l i u m y
917@n c s t .e d u .c n 摘要:明确影响互花米草(S p a r t i n a a l t e r n i f l o r a )的主导环境因子及其潜在适生区对防控管治互花米草具有重要意义㊂本文以中国大陆滨海地区为研究区域,综合考虑气候㊁地形㊁土壤及海洋环境因子,利用生态位模型(M a x e n t )预测互花米草潜在适生区㊂结果表明:互花米草的分布受海拔影响最为显著,其主要分布于-2~5m 范围内㊂此外,当年平均气温介于14.0ħ~23.5ħ㊁土壤可交换钠盐含量介于26.4%~46.0%㊁气温极差介于14.2ħ~26.
3ħ或海水盐度介于23.0%~28.3%及33.5%~33.8%时,互花米草存在概率均达到60%以上㊂互花米草高度适生区主要分布于江沪浙闽四省市以及莱州湾地区,该区域内适宜互花米草生长的海拔㊁年平均气温㊁气温极差㊁海水盐度分别为-1~0m ,17.9ħ,23.0ħ及27.5%㊂因土壤可交换钠盐含量空间异质性差异较小,故在该区域内对互花米草分布影响不明显㊂此外,大丰麋鹿㊁厦门珍稀海洋物种及漳江口红树林国家级自然保护区面临较高的互花米草入侵风险㊂
关键词:互花米草;环境因子;生态位模型;潜在适生区;入侵风险
中图分类号:Q 948 文献标识码:A 文章编号:1007-0435(2023)11-3514-11
A n a l y s i s o n t h eP o t e n t i a lH a b i t a t S u i t a b i l i t y o
f S p a r t i n a a l t e r n i f l o r a B a s e do n M a x e n tM o d e l
F A N
G
H u a 1,
2,C H E N X i n g -t o n g 1,L I U M i n g -y u e 1,2,
3*,Z HA N G Y o n g -b i n 1,M I A OZ h e n g -h o n g 4
,
MA N W e i -d o n g 1,2,3,Z H A N G Q i n g -w e n 1,2,K O U C a i -y a o 1,2,L IX i a n g
1,2
(1.S c h o o l o fM i n i n g E n g i n e e r i n g ,N o r t hC h i n aU n i v e r s i t y o f S c i e n c e a n dT e c h n o l o g y ,T a n g
s h a n ,H e b e i P r o v i n c e 063210,C h i n a ;2.T a n g s h a nK e y L a b o r a t o r y o fR e s o u r c e s a n dE n v i r o n m e n tR e m o t eS e n s i n g ,T a n g s h a n ,H e b e i P r o v i n c e 063210,C h i n a ;3.H e b e i I n d u s t r i a lT e c h n o l o g y I n s t i t u t e o fM i n eE c o l o g i c a lR e m e d i a t i o n ,T a n g
s h a n ,H e b e i P r o v i n c e 063210,C h i n a ;4.J i l i nH y d r o p o w e rP l a n n i n g &D e s i g n I n s t i t u t e ,C h a n g
c h u n ,J i l i nP r o v i n c e 130033,C h i n a )A b s t r a c t :I
d
e n t i
f y i n
g t
h ed o m
i n a n t e n v i r o n m e n t a l f a c t o r s a n d p o t e n t i a l l y s u i t a b l e a r e a s f o r S p a r t i n a a l t e r -n i f
l o r a p r e v e n t i o na n d c o n t r o l i s c r i t i c a l .I n t h i s s t u d y ,t h e c o a s t a l a r e ao fm a i n l a n dC h i n a i s t a k e na s t h e r e s e a r c h a r e a ,a n d t h e p o t e n t i a l l y s u i t a b l e a r e a s o f S .a l t e r n i f
l o r a a r e p r e d i c t e db y u s i n g t h eM a x e n tm o d -e l ,c o n s i d e r i n g t h e c l i m a t e ,t o p o g r a p h y
,s o i l ,a n dm a r i n e e n v i r o n m e n t a l f a c t o r s .T h e r e s u l t s s h o wt h a t :t h e d i s t r i b u t i o no f S .a l t e r n i f
l o r a i s m o s ts i g n i f i c a n t l y a f f e c t e db y t h ee l e v a t i o n ,m a i n l y d i s t r i b u t e di nt h e r a n g e o f-2~5ma b o v e t h e s e e l e v e l .I n a d d i t i o n ,t h e a n n u a lm e a n t e m p e r a t u r e ,s o i l e x c h a n g
e a b l e s o d i u m s a l t s ,e x t r e m e t e m p e r a t u r e v a r i a t i o n ,a n d s e a w a t e r s a l i n i t y a l s o p l a y e d a n e s s e n t i a l r o l e i n t h e d i s t r i b u t i o n o
f S .a l t e r n i f l o r a .W h e n t h e a n n u a lm e a n t e m p e r a t u r ew a s 14.0ħ~23.5ħ,s o i l e x c h a n
g e a b l e s o d i u ms a l t s c o n t e n tw a s 26.4%~46.0%,t
h e e x t r e m e t e m p e r a t u r e v a r
i a t i o nw a s 14.2ħ~26.3ħ,a n d t h e s e a w a t e r s a -l i n i t y w a s 23.0%~28.3%o r33.5%~33.8%,t h ee x i s t e n c e p r o b a b i l i t y o f S .a l t e r n i f
l o r a w a sa b o v e 60%.T h e h i g h l y s u i t a b l e a r e a s f o r S .a l t e r n i f
l o r a w e r em a i n l y d i s t r i b u t e d i n J i a n g s u ,S h a n g h a i ,Z h e j i a n g ,
第11期方铧等:基于M a x e n t模型的互花米草潜在生境适宜性分析
a n dF u j i a n p r o v i n c e s a n d t h eL a i z h o uB a y a r e a.T h e e l e v a t i o n,a n n u a lm e a n t e m p e r a t u r e,e x t r e m e t e m p e r a-t u r e v a r i a t i o n,a n d s e a w a t e r s a l i n i t y s u i t a
b l e f o r S.a l t e r n i f l o r a g r o w t h i n t h i s r e g i o nw e r e-1~0m,17. 9ħ,23.0ħ,a n d27.5%,r e s p e
c t i v e l y.T h e e f f e c t o f e x c h a n g e a b l e s o
d i u mo n t h
e d i s t r i b u t i o no
f S.a l t e r n i-f l o r a w a s n o t i m m e d i a t e l y a p p a r e n t d u e t o t h e s m a l l s p a t i a l h e t e r o
g e n e i t y o f i t s c o n t e n t i n t
h e s o
i l.I n a d d i-t i o n,D a f e n g M i l uN a t u r eR e s e r v e,X i a m e nR a r eM a r i n e S p e c i e sN a t i o n a lN a t u r eR e s e r v e,a n dZ h a n g
j i a n g k-o u M a n g r o v eF o r e s tN a t i o n a lN a t u r eR e s e r v e f a c e ah i g h r i s ko f S.a l t e r n i f l o r a i n v a s i o n.
K e y w o r d s:S p a r t i n aa l t e r n i f l o r a;E n v i r o n m e n t a l f a c t o r;M a x e n tm o d e l;P o t e n t i a ls u i t a b l ea r e a;I n v a s i o n r i s k
外来生物入侵被认为是导致全球生物多样性下降的第二大因素[1],仅次于栖息地退化[2],对生态环境㊁生物多样性和社会经济造成巨大危害[3]㊂因此研究外来生物的入侵力与入侵风险,预测外来生物的潜在适生区并进行生物防治成为生态学研究热点之一㊂
互花米草(S p a r t i n aa l t e r n i f l o r a)为禾本科米草,属多年生草本植物,原产于大西洋西海岸及墨西哥湾[4],具有广盐性㊁耐淹性和耐低氧能力[5]㊂出于保滩护岸㊁促淤消浪的目的,互花米草被引入中国[6]㊂后经人为引种㊁潮汐和风力等传播方式,互花米草迅速扩散,向南已至广西,向北已达河北南堡等沿海区域[7-9]㊂互花米草引种起到促淤造陆与提高滨海湿地碳汇能力等积极影响的同时,也引发了本土植物数量和面积大规模缩减㊁海鸟生存环境遭到破坏等负面影响[10-12],被国家环境保护总局认定为16种最严重的入侵物种之一[13]㊂按照2022年至2025年互花米草整治需求,自然资源部等部门提出全面遏制互花米草扩散态势,确保滨海湿地生态安全的行动计划㊂已有研究发现,互花米草客观地理分布远未达其最大潜在适生范围,仍存在继续扩散蔓延的可能性㊂科学预测和评价互花米草潜在适生区,明确影响互花米草分布的主导环境因子及其响应特征对修复滨海湿地生态系统具有十分重要的意义㊂
生态位(M a x e n t)模型通过物种已知地理分布点和相关环境因子,探寻在环境因子约束下最大熵时物种分布概率[14]㊂目前,已有研究采用领域模型(D O MA I N)[15]㊁遗传算法(G A R P)和M a x e n t 模型等预测了入侵植物潜在分布[16],并揭示了入侵植物分布对环境因子的响应㊂D O M I N模型预测所需样本多,预测精度受样本偏差影响大[17]㊂G A R P模型假阳性率偏高,预测结果存在物种客观非适生区成为适生区的可能性[18]㊂相比于前两种模型,M a x e n t模型操作相对简单[19],对小样本㊁不规则采样㊁数据偏差小的数据容错性较好,且对物种表现出较高的预测精度,是预测物种潜在适生区的首选模型[20-21]㊂目前,已有学者采用M a x-e n t模型对互花米草潜在适生区进行预测,如张静涵等应用M a x e n t模型预测互花米草在我国的潜在适生区,发现其主要集中在江浙沪及福建㊁广东及广西沿海区域[22]㊂Y u a n等[23]基于I P C C未来气候数据预测互花米草在不同空间尺度下的潜在分布,结果表明在未来气候下,互花米草在中国从海岸中部向海岸两侧蔓延,而在美国从西海岸向东海岸蔓延,且互花米草在入侵过程中对高温干燥环境适应能力越来越强㊂以往利用M a x e n t模型预测互花米草分布相关研究多基于气候因素或结合I P C C未来不同气候模式开展,多因素综合作用下的互花米草潜在分布研究尚且不足㊂
作为模型预测变量,气候㊁海洋环境㊁土壤及地形因子等极大地影响了生境可侵性㊂降水㊁光照与温度差异等气候因子通过影响互花米草建群与扩散来影响其分布[24]㊂地形差异导致互花米草在潮间带的持续淹水时长和淹水深度不同[25],进而影响互花米草入侵能力㊂土壤为植物的生长发育提供物质基础,土壤可交换钠盐含量㊁酸碱度等不同程度地限制或促进植物的生长[26];在一定土壤盐分范围内,土壤盐分越高,潮汐带来的营养元素越丰富[5]㊂此外,海洋环境因素如海水流速㊁海水盐度等也会影响植物的建群概率[27]㊂因此,为理解互花米草入侵的环境影响机制及预测其潜在适生区域,本研究综合考虑气候㊁土壤㊁地形及海洋环境因素,利用刀切法(J a c k k n i f e)结合P e r s o n相关分析进一步筛选预测互花米草适生区域的环境因子,基于M a x e n t模型开展中国大陆沿海地区互花米草适生区预测,探讨互花米草对环境因素的响应特征,为互花米草防控管治和滨海湿地的可持续发展提供科学依据㊂
1材料与方法
1.1研究区概况
本文以中国大陆海岸线为中心,向陆一侧
5153
草地学报第31卷
10k m缓冲区作为陆向边界;向海侧基于E T O P O1数据提取第一条连续的15m等深线为海向边界,海陆边界之间区域作为研究区,地理范围为20ʎ13'~43ʎ26'N,104ʎ26'~123ʎ10'E,海岸线总长度超过18000k m,跨温带季风㊁亚热带季风气候区,夏季高温多雨㊂研究区内秦皇岛㊁汕头及北海的潮汐类型为日潮,其他中国大陆沿海区域潮汐类型皆为半日潮[28];海岸类型多样,除山东半岛㊁辽东半岛部分地区为基岩海岸外,杭州湾以北多为沙质和淤泥质海岸;杭州湾以南地区则以基岩海岸为主,部分地区分布红树林㊁珊瑚等生物海岸[29]㊂研究区内盐生植被类型丰富,以芦苇㊁互花米草㊁红树林㊁碱蓬及短叶茳芏等为主㊂截止2023年,该区域包含滨海湿地㊁珍稀海洋生物㊁候鸟珍禽㊁地质地貌等各类国家级自然保护区及国际重要湿地共30个,对于维持生物多样性与滨海湿地可持续发展具有重要作用㊂
1.2互花米草分布数据
物种分布点的质量和数量直接影响模型的预测精度[30]㊂本文互花米草分布数据来源于2016 2020年野外实地调查,调查时间为每日海水落潮时,采用高精度手持G P S记录互花米草的点位,每个点位相隔50m以上,获取津冀滨海地区44个互花米草分布点;在野外调查的基础上,利用2015年中国科学院互花米草栅格数据集(h t t p://w w w.
i g a d c.c n/)进行栅格转点,同时生成渔网,每个渔网中保留一个互花米草分布点,经哨兵二号遥感影像逐一核对后共获取到473个互花米草有效分布点,与野外实地调查互花米草分布点合并,去除重复记录后获得互花米草分布点共516个㊂将互花米草分布点坐标数据按照物种名㊁经度和纬度的格式保存为.c s v样本数据格式文件㊂
1.3环境因子
基于前人[25,30-31]研究并结合互花米草的生态习性,综合选取具有较强生物学意义的气候㊁土壤㊁地形及海洋环境因子4大类共31个环境变量构建模型预测变量库㊂气候因子源自中国气象数据网站(h t t p://d a t a.c m a.c n)气象站点观测数据集,共14个气候因子;利用联合国粮农组织世界土壤数据库(h t t p s://d a t a.i s r i c.o r g)提取有机碳含量㊁土壤可交换钠盐等11个土壤因子;海洋环境因子源自国家海洋科学数据中心(m d s.n m d i s.o r g.c n)与海洋数据中心(h t t p s://i r i.c o l u m b i a.e d u)的站点观测数据,共获取5个海洋环境因子;海拔数据由来源于中国科学院资源环境科学与数据中心(h t t p s://w w w. r e s d c.c n)的陆地地形数据与来源于N O A A(w w w. n g d c.n o a a.g o v)的海洋地形数据镶嵌生成㊂气候因子与海洋环境因子原始数据为点数据,采用K r i g i n g 插值法将观测数据插值㊂将以上数据掩膜提取为研究区范围边界,统一空间分辨率为500mˑ500m,地理坐标系为C G S_WG S_1984,并转化为A S C I I 格式(表1)
㊂
图1研究区及互花米草分布点
F i g.1 S t u d y a r e a a n dd i s t r i b u t i o no f
S p a r t i n a a l t e r n i f l o r a
1.4环境因子筛选与M a x e n t模型预测
M a x e n t模型通过物种已知地理分布点和相关环境因子预测物种的潜在分布[32]㊂由于参与模型预测的31个环境因子间存在多重共线性,造成M a x e n t模型构建数据过度拟合,影响模拟结果[33]㊂因此,首先利用M a x e n t模型中的刀切法分析各环境因子对模型预测的贡献率,剔除贡献率较小的环境因子;再对剩余环境因子进行P e r s o n相关性分析,当相关性系数绝对值大于0.8时,剔除对模型预测相对贡献率较小的环境因子[34],最终优选9个环境因子作为模型建模参量(表2)㊂
6153
第11期方铧等:基于M a x e n t模型的互花米草潜在生境适宜性分析
表1用于模型分析的环境因子
T a b l e1 E n v i r o n m e n t a l f a c t o r s f o rm o d e l a n a l y s i s
代码C o d e环境因子E n v i r o n m e n t a l f a c t o r s数据集D a t a s e t单位U n i t B i o1日积温ȡ10ħA c c u m u l a t e d t e m p e r a t u r eȡ10ħ气候因子C l i m a t e f a c t o rħB i o2年平均气温A n n u a lm e a n t e m p e r a t u r e气候因子C l i m a t e f a c t o rħB i o3年均月最高温度M e a nm o n t h l y m a x i m u mt e m p e r a t u r e气候因子C l i m a t e f a c t o rħB i o4年均月最低温度A n n u a lm e a nm o n t h l y m i n i m u mt e m p e r a t u r e气候因子C l i m a t e f a c t o rħB i o57月份平均温度A v e r a g e t e m p e r a t u r e i n J u l y气候因子C l i m a t e f a c t o rħB i o61月份平均温度A v e r a g e t e m p e r a t u r e f o r J a n u a r y气候因子C l i m a t e f a c t o rħB i o77月份平均降水量A v e r a g e p r e c i p i t a t i o n i n J u l y气候因子C l i m a t e f a c t o r m m B i o81月份平均降水量A v e r a g e p r e c i p i t a t i o n i nJ a n u a r y气候因子C l i m a t e f a c t o r m m B i o9露点温度T h e d e w p o i n t t e m p e r a t u r e气候因子C l i m a t e f a c t o rħB i o10年平均降水量M e a na n n u a l p r e c i p i t a t i o n气候因子C l i m a t e f a c t o r m m B i o11年平均地面气压A n n u a lm e a n s u r f a c e p r e s s u r e气候因子C l i m a t e f a c t o r P a B i o12年平均风向A n n u a lm e a nw i n dd i r e c t i o n气候因子C l i m a t e f a c t o r-B i o13年平均风速A n n u a lm e a nw i n d s p e e d气候因子C l i m a t e f a c t o r m㊃s-1 B i o14气温极差E x t r e m e t e m p e r a t u r e v a r i a t i o n气候因子C l i m a t e f a c t o rħB i o15海表压力S e a s u r f a c e p r e s s u r e海洋环境因子M a r i n e f a c t o r P a B i o16海水盐度S e a w a t e r s a l i n i t y海洋环境因子M a r i n e f a c t o r% B i o17海水水位S e aw a t e r海洋环境因子M a r i n e f a c t o r m B i o18海表温度S e a s u r f a c e t e m p e r a t u r e s海洋环境因子M a r i n e f a c t o rħB i o19海水流速T h ew a t e r f l o wr a t e海洋环境因子M a r i n e f a c t o r m㊃s-1 B i o20土壤有效含水量A v a i l a b l e s o i lw a t e r c o n t e n t土壤因子S o i l f a c t o r% B i o21土壤电导率S o i l c o n d u c t i v i t y土壤因子S o i l f a c t o r% B i o22土壤可交换钠盐S o i l e x c h a n g e a b l e s o d i u ms a l t s土壤因子S o i l f a c t o r% B i o23硫酸盐含量S u l f a t e c o n t e n t土壤因子S o i l f a c t o r% B i o24碳酸盐或石灰含量C a r b o n a t e o r l i m e c o n t e n t土壤因子S o i l f a c t o r% B i o25交换性盐基C o m m u t a t i v e b a s e土壤因子S o i l f a c t o r% B i o26土壤阳离子交换力S o i l c a t i o ne x c h a n g e c a p a c i t y土壤因子S o i l f a c t o r% B i o27酸碱度P H土壤因子S o i l f a c t o r-B i o28粘土含量C l a y c o n t e n t土壤因子S o i l f a c t o r% B i o29淤泥含量S i l t c o n t e n t土壤因子S o i l f a c t o r% B i o30有机碳含量O r g a n i c c a r b o nc o n t e n t土壤因子S o i l f a c t o r% B i o31海拔E l e v a t i o n地形因子T e r r a i n f a c t o r m
表2M a x e n t建模参量
T a b l e2 M o d e l i n gp a r a m e t e r
代码C o d e环境因子E n v i r o n m e n t a l f a c t o r s数据集D a t a s e t单位U n i t
B i o1日积温ȡ10ħA c c u m u l a t e d t e m p e r a t u r eȡ10ħ气候因子
C l i m a t e f a c t o rħ
B i o2年平均气温A n n u a lm e a n t e m p e r a t u r e气候因子
C l i m a t e f a c t o rħ
B i o57月份平均温度A v e r a g e t e m p e r a t u r e i nJ u l y气候因子
C l i m a t e f a c t o rħ
B i o81月份平均降水量A v e r a g e p r e c i p i t a t i o n i n J a n u a r y气候因子
C l i m a t e f a c t o r m m
B i o14气温极差E x t r e m e t e m p e r a t u r e v a r i a t i o n气候因子
C l i m a t e f a c t o rħ
B i o16海水盐度S e a w a t e r s a l i n i t y海洋环境因子M a r i n e f a c t o r%
B i o22土壤可交换钠盐S o i l e x c h a n g e a b l e s o d i u ms a l t s土壤因子S o i l f a c t o r%
B i o30有机碳含量O r g a n i c c a r b o n c o n t e n t土壤因子S o i l f a c t o r%
B i o31海拔E l e v a t i o n地形因子T e r r a i n f a c t o r m
构建互花米草分布预测M a x e n t模型时,随机选取75%的互花米草分布数据作为训练集,剩余的25%作为测试集,重复运行模型10次,选择最优模拟结果进行可视化栅格数据转化[24]㊂通过刀切法测定各环境因子对模型构建的贡献率,创建环境变量响应曲线获取物种分布概率对各环境变量的响应情况㊂利用受试者工作特征曲线(R O C)评价模型预测精度㊂将R O C曲线下面积值(A U C)作为模型精度评价标准,A U C值越大模型预测精度越高,评价标准为:较差(0.6~0.7)㊁一般(0.7~0.8)㊁良好(0.8~0.9)和优秀(0.9~1.0)[35]㊂
1.5适宜区等级划分
基于自然间断点分级法[36],将互花米草适生等
7153
草 地 学 报
第31卷
级(S u i t a b i l i t y I n d e x ,S I )划分为非适生区(S I <0.05)㊁低度适生区(0.05ɤS I <0.2)㊁中度适生区(0.2ɤS I <
0.6)㊁高度适生区(0.6ɤS I <1)4个等级㊂
2 结果与分析
2.1 模型预测精度评价及制约互花米草分布的主导环境因子
基于物种分布数据和环境变量,通过M a x e n t 模型对互花米草潜在适生区进行模拟,获取最优预测结果R O C 曲线,如图2所示,最优模型训练集和测试集的A U C 值分别达到0.877和0.871,模型预测有很高的准确性,对于互花米草潜在适生区分布具有较好的指示性
㊂
图2 互花米草潜在适生区最优预测结果R O C 曲线F i g
.2 R O Cc u r v e o f t h e p r e d i c t e d r e s u l t s o f p o t e n t i a l h a b i t a t a r e a o f S .a l t e r n i f
l o r a 在使用单一变量构建模型时,海拔㊁年平均气温㊁土壤可交换钠盐及气温极差均具有较高的
A U C 值(
图3)㊂其中,海拔与其他环境因子相比,与互花米草潜在适生区的耦合性更高,对互花米草潜在适生区影响最大㊂年平均气温㊁气温极差与土壤可交换钠盐提供了较大增益,表明其对互花米草分布具有一定的影响㊂7月份平均温度㊁1月份平均降水量㊁有机碳含量㊁日积温(ȡ10ħ)和海水盐度相对其他建模因子对互花米草分布影响作用较小
㊂
图3 各环境因子刀切法(J a c k k n i f e )检验A U C 得分
F i g
.3 J a c k k n i f e t e s tA U Cs c o r e s o f v a r i o u s e n v i r o n m e n t a l f a c t o r s
2.2 互花米草潜在分布对环境因子的响应
为揭示环境因素对互花米草潜在分布的影响,本文选取刀切法评价单一环境变量对互花米草潜在分布重要性,优选5个主导环境因子绘制互花米草存在概率的响应曲线(图4)㊂互花米草存在概率随年平均气温与气温极差的增长总体呈现出先升高后下降的趋势(图4A ,4B )㊂其中,随着年平均气温的增大,互花米草存在概率迅速升高,至17.8ħ处互花米草存在概率达到最大,后缓慢下降又迅速下降(图4A )
;气温极差小于12.8ħ或大于29.9ħ时,互花米草存在概率小于20%(图4B );海水盐度位于23.0%~28.3%及33.5%~33.8%范围内,互花米草存在概率达60%以上㊂其中,互花米草存在概率呈 双峰 状态,在海水盐度23.3%及33.7%处达到峰值(图4C
)㊂随着可交换钠盐含量的升高,互花米草存在概率缓慢升高,当土壤可交换钠盐含量超过43.8%后迅速下降(图4D )㊂互花米草主要生长在海拔-1~2m 的地区,在该区间互花米草存在概率达到85%以上㊂海拔过低导致互花米草过度淹水造成土壤缺氧,限制了互花米草在沿海方向的生长㊂海拔过高,互花米草淹水缺失造成过度蒸腾会限制其沿陆地方向的生长㊂故海拔超过2m 以后,随着海拔的不断升高,互花米草存在概率逐渐减小,直至消失(图4E
)㊂8
153
第11期方 铧等:基于M a x e n t
模型的互花米草潜在生境适宜性分析
图4 存在概率对主要环境因子的响应曲线
F i g .4 R e s p o n s e c u r v e o f e x i s t e n c e p r o b a b i l i t y t
om a i ne n v i r o n m e n t a l v a r i a b l e s 2.3 互花米草在我国潜在适生区及适宜性评价
M a x e n t 预测结果显示:互花米草潜在适生区总体呈现出 中部高,南北低 的分布趋势(图5)㊂其中,互花米草高度适生区主要分布于中国大陆沿海中部偏南的江浙沪沿海和福建北部与南部海湾地区,除浙江台州市至福建福州市局部间断,高度适生区呈现大面积连续分布特征,此外,山东莱州湾是黄河以北唯一的高度适生区㊂中度适生区主要集中在渤海湾㊁胶东半岛㊁福建与广东中部沿海地区及北部湾西部㊂低度适生区主要分布在河北秦皇岛㊁广东深圳以东与江门以西沿海地区㊂非适生区与低度适
生区相接,大面积集中于高纬度与低纬度沿海地区,如辽东湾与雷州半岛沿海地区㊂
本研究所得我国沿海国家自然保护区与国际重要湿地互花米草潜在入侵风险与对应范围内现有互花米
草实际分布趋势较为一致[9,37
]㊂整体上看,高度与中度
适生区范围内的国家级自然保护区与国际重要湿地大量分布于中国中部与东南沿海地区,渤海湾㊁北部湾及雷州半岛范围内的国家级自然保护区与国际重要湿地则主要位于低度与非适生区内㊂进一步研究发现中国沿海地区国家级自然保护区互花米草潜在入侵风险有所差异(图6),其中黄河三角洲国家级自然保护区面积的20.18%为互花米草高度适生区,有关部门仍需加大互花米草防控力度㊂江苏盐城湿地珍禽国家级自然保护区与大丰麋鹿国家级自然保护区内均存在大面积互花米草潜在高度适生区,其互花米草高度入侵风险分别高达
49.99%与77.78%㊂此外,
上海崇明东滩鸟类国家级自然保护区内与九段沙湿地国家级自然保护区内互花米草潜在高度适生区面积分别达到保护区面积的42.00%与
38.54%,中度适生区面积分别达到50.58%与29.85%,
互花米草潜在入侵风险较高㊂福建厦门珍稀海洋物种国家级自然保护区与漳江口红树林国家级自然保护区内面临互花米草入侵压力均相对较高,其潜在高度适生区面积分别达69.97%与78.65%㊂相较于山口红树林国家级自然保护区,北仑河口国家级自然保护区受互花米草高度入侵风险相对较小,仅为35.26%
㊂
图5 互花米草潜在适生分布
F i g .5 P o t e n t i a l S u i t a b i l i t y
G r a d e s o f S .a l t e r n i f
l o r a 9
153
草地学报第31
卷
图6国家级自然保护区及国际重要湿地内互花米草各适生区面积占比
F i g.6 T h e p r o p o r t i o n o f s u i t a b l e a r e a s o f S.a l t e r n i f l o r a i nS t a t e-l e v e lN a t u r eR e s e r v e s a n dW e t l a n d s o f I n t e r n a t i o n a l I m p o r t a n c e
本研究分别在高㊁中㊁低度及非适生区生成86个随机点,保证点与点之间距离大于50m,基于随机点提取影响互花米草分布的五个主导环境变量值,进一步分析不同适生区内互花米草对环境变量响应差异(图7),研究发现:高度与低度适生区内互花米草受年平均气温与气温极差影响较为相似,但高度适生区内互花米草集中分布于年平均温度17.9ħ左右对应区域,而低度适生区内互花米草主要分布于17.5ħ的地区,且在该温度以下互花米草较高度适生区分布更集中;气温极差影响下的高度与低度适生区互花米草主要分布于23.0ħ与24.6ħ气温极差对应区域,高度适生区形成窄波峰,而低度适生区波峰宽度稍大,表明在高度适生区内适宜互花米草密集分布的气温极差范围小于低度适生区㊂中度适生区内互花米草较为均匀的分布在各年平均气温范围,气温极差处于15.0ħ与26.2ħ条件下互花米草分布最集中㊂由于非适生区在中国沿海跨度较大,故非适生区互花米草在年平均气温与气温极差各值域范围内分布较均匀,但在年平均气温处于12.6ħ与23.0ħ及气温极差处于27.7ħ与14.0ħ时互花米草分布更密集㊂土壤可交换钠盐含量在研究区范围内异质性差异较小,含量高的地区主要分布于唐山至潍坊沿海地区,但位于该区域内的互花米草分布与生成随机点数目均较少,因此各分区之间土壤可交换钠盐含量对互花米草分布影响无明显区别㊂不同适生区内,互花米草对海水盐度的敏感程度存在差异㊂高度适生区内互花米草主要分布在海水盐度为27.5%的区域㊂中度与非适生区内互花米草均分布在海水盐度较高的地区㊂低度
适生区内互花米草则广泛分布于各海水盐度范围内㊂互花米草在高度适生区与中度适生区内均在海拔-1~0m范围分布较为集中,但低度适生区与非适生区内海拔跨度大幅增加㊂
3讨论
本文采用M a x e n t模型对外来入侵植物互花米草的潜在适生区进行预测,解析了互花米草分布对气候㊁土壤㊁海洋环境以及地形的响应情况,并预测了互花米草的潜在适生区,为互花米草精准防控提供决策支持㊂本研究表明:互花米草主要分布在莱州湾㊁江浙沪及福建沿海地区,与张丹华等研究结果较为一致,但张丹华预测表明辽宁沿海地区是互花米草中度适生区,而在当前环境背景下本文研究结果中辽宁省为互花米草非适生区,这与互花米草实际在辽宁省的分布较少或缺失有关[8-9,30,38],但不排除气候等环境变化或互花米草驯化向高纬度低海拔地区扩散;陈思明研究结果中互花米草主要适生区域未包含莱州湾地区[36],原因是本研究选取的莱州湾地区互花米草分布点位年份较新,更接近现阶段互花米草分布,且参与建模的环境因子更丰富导致㊂与已有研究对比[22-23,31,36,38],本文构建的研究区符合最新年份互花米草分布范围,且丰富了模型建模因子,综合选取对互花米草影响较大的气候㊁海洋环境㊁土壤以及地形因子进行模型预测,能够准确反映互花米草在我国沿海地区的潜在分布㊂
0253
第11期方 铧等:基于M a x e n t
模型的互花米草潜在生境适宜性分析
图7 不同适生区内互花米草环境变量的比较
F i g .7 C o m p a r i s o no f S .a l t e r n i f
l o r a e n v i r o n m e n t a l v a r i a b l e sb e t w e e nd i f f e r e n t s u i t a b l e a r e a s 本文预测的互花米草高度适生区与客观分布大致相同,主要分布于江苏至福建沿海区域,该区域与
互花米草原产地纬度相同[22]
,具有和互花米草原产
地相似的气候特征与滩涂特征,在该区域内互花米草具有强大的入侵能力;津冀与山东中部沿海地区㊁福建与广东中部沿海地区及广西南部沿海地区是互花米草中度适生地区,然而如天津塘沽和大港等地区大范围填海造田和政府治理等人为干预使互花米
草空间分布不断收缩[35
],导致实际采样点减少,参
与模型预测的分布点减少㊂如无人为干预情况下,互花米草中度适生区有成为高度适生区的可能性㊂研究显示,辽东湾㊁秦皇岛沿海地区为互花米草低度与非适生区,广东除中部沿海地区为中度适生区,其余皆为低度与非适生区㊂导致该现象产生的原因可能是辽宁沿海地区缺少参与模型预测的地理分布点,同时互花米草入侵存在环境阈值效应,极端高温
或极端低温会抑制互花米草的生长[
39
]㊂本研究表明年平均气温与气温极差对互花米草
潜在分布影响较大,互花米草存在概率随年平均气温与气温极差的升高呈现出先升高后降低的趋势,与赵相健等得到的互花米草成功入侵存在气候阈值
效应的研究结果一致[
40]
㊂但相比张丹华研究中适宜互花米草生长的年平均气温跨度较大[
30]
㊂且有研究结果显示互花米草分布面积呈现由低纬度地区向高纬度地区先增加后减少的趋势,然而高纬度地
区种源的种子萌发率普遍高于低纬度地区[41-42]
,这
预示着中纬度沿海地区仍是互花米草入侵的主要区域,但不排除纬度更高的华北与东北沿海地区有成为互花米草高度适生区的可能性㊂考虑到互花米草分布受到海拔的影响,当海拔处于-1~2m 的区域
时,互花米草存在概率最高,而陈思明[36]
预测互花
米草在海拔-50~200m 均有较高分布㊂与以往研究比较,本文结果更符合互花米草生长特性,这是由于海拔基本决定了潮汐的强弱,中潮区潮汐为互花米草繁殖带来充足的营养物质[30
];而处于低海拔的
互花米草频繁受到强潮汐的作用与长时间海水淹
浸,其生长繁殖遭受负面影响,故互花米草向海一侧扩散存在阈值;伴随海拔持续升高,互花米草受到有限潮水作用㊁过度蒸腾及与本地物种竞争等影响,高潮滩互花米草成体生长与种子活性均受到抑制,因
此,互花米草在向陆地入侵时存在高程阈值[
43]
㊂互花米草生长在海陆潮间带,受海水盐度影响明显,适
1
253。