福建闽江水口水库凤眼莲空间分布特征及其动态*

2012-12-11 06:42潘文斌
湖泊科学 2012年3期
关键词:凤眼莲支流水口

陈 潇,潘文斌,王 牧

(福州大学环境与资源学院,福州350108)

福建闽江水口水库凤眼莲空间分布特征及其动态*

陈 潇,潘文斌**,王 牧

(福州大学环境与资源学院,福州350108)

为研究闽江水口水库凤眼莲空间分布特征及动态变化,根据2001-2009年水口库区无云覆盖的24幅卫星遥感影像,构建研究区域凤眼莲干流-支流分布的空间结构,并通过建立端元的混合象元线性分解模型,对福建省水口水库不同时间、不同区域的凤眼莲分布特征进行系统分析,结果表明,随着时间的推移,凤眼莲的暴发区逐渐向下游迁移,同时,支流对研究区域凤眼莲的面积有重要影响,特别是河流较长、周边有人口居住的支流是凤眼莲源头治理的重点;各种人为或自然原因造成的水体富营养化、河流流速变缓等是凤眼莲泛滥的重要原因.建议针对研究区域上游及其周边进行养殖污染、工业污染整治,强化对污染源的监督管理,推进生态环境建设与保护,建设生态功能保护区等.本研究可为了解水口水库的凤眼莲时空分布规律,进而制定相应的预防与控制管理措施提供必不可少的数据基础和理论依据.

凤眼莲;水口水库;空间分布特征及其动态

凤眼莲,拉丁名Eichhornia crassipes(Water Hyacinth),原产于南美的巴西东北部,属雨久花科,它还有许多俗名,如水葫芦、凤眼蓝、水凤仙、水风信子等.凤眼莲现广泛分布于北美、亚洲、大洋洲和非洲等众多国家,是目前世界上危害最严重的水生漂浮植物,被列为世界十大害草之一[1].1901年凤眼莲作为观赏植物引入中国,1950s开始在我国南方地区作为畜禽饲料栽培,由于繁殖速度极快,又几乎没有竞争对手和天敌,易在生长区内形成优势物种,目前已在我国南方多处水域发生凤眼莲种群大暴发事件.福建位于全国东南部,其温暖潮湿的气候利于外来物种(凤眼莲)生长繁殖.近年来凤眼莲在江河、水库泛滥成灾,致使水质恶化,水生动植物大量死亡,并严重影响行洪排涝、航运、灌溉、发电、水产养殖及旅游,给社会、经济带来巨大损失[2-3].本文研究区水口库区则是福建的重灾区之一.

1 材料和方法

1.1 研究区域概况

水口水库位于闽江中下游的闽清县雄江镇下蹼村(26°18'7″~26°38'25″N,118°10'44″~118°48'48″E),库区跨越闽清县、古田县、尤溪县和南平市延平区.水口水库1993年3月开始蓄水,1996年11月竣工,为河道型水库,总库容26×108m3,库区主河道长94 km,面积约99.6 km2,水库最大坝高101 m,坝顶高程74 m,是目前中国东南地区最大的水力发电站[13].2002年起闽江上游出现少雨干旱、气温偏高的气候现象,上游来水量减少,汛期水口大坝也基本没有开闸泄洪,因此绝大部分库区上游漂流而下的凤眼莲被拦截在库内兹生蔓延,且水口库区的水流缓慢,给凤眼莲提供了快速生长的环境[3].

1.2 技术路线

研究基本思路和技术路线如图1.

图1 技术路线图Fig.1 Technology roadmapping

1.3 数据收集与样点布设

共收集了2001-2009年水口水库天气状况良好的卫星遥感影像24幅,其中LANDSAT TM/ETM+为15幅(2001年1月2日、2001年5月23日、2002年3月7日、2002年5月26日、2003年2月6日、2003年3月26日、2006年6月22日、2006年9月26日、2006年11月5日、2007年1月8日、2007年5月8日、2007年7月27日、2007年9月13日、2009年6月6日、2009年7月16日),CEBERS CCD为9幅(2001年2月16日、2005年10月21日、2003年12月18日、2004年10月22日、2004年12月16日、2006年1月7日、2006年6月12日、2006年8月29日、2007年8月13日).同时收集水口水库相关的新闻资料,对比当地拍摄的照片,发现水口库区内漂浮的都是凤眼莲,因此,本研究假设影像中河道内发现的水生植被均为凤眼莲.

运用ERDAS对卫星影像进行几何纠正及研究区域的水体提取.基于研究区域卫星影像的不同地物(水体和凤眼莲)的光谱特征,选择432波段组合作为凤眼莲信息提取的模式.将水口水库出现过凤眼莲的地点设为观察区,共将水口库区分为8个支流汇水区,并将干流分成7段(图2),其中G开头的点代表干流的观察区,Z开头的点代表支流的观察区.

图2 研究区子区布设Fig.2 Map of observation zones in the Shuikou Reservoir

1.4 混合像元分解

通过归一化光谱,运用约束线性光谱混合分析方法计算研究区域每种端元比例[14],公式如下:

2 结果与分析

2.1 总体趋势分析

将2001-2009年研究区域划分的7段干流分别统计发现,随着时间的推移,凤眼莲的暴发区逐渐向下游迁移.2003、2004年暴发区多在流域上游,而2007、2009年暴发区在流域下游(图3).从总体上看,上游凤眼莲暴发多在2006年前,而下游多在2006年后.

图3 研究区域凤眼莲暴发区分布(a、b:CEBERS卫星遥感影像;c、d:Landsat TM卫星遥感影像;影像中红色部分为凤眼莲分布区域;影像a缺少该时间水口库区坝前的影像)Fig.3 The distribution of Eichhornia crassipes in survey region(a,b:CEBERS satellite imagery;c,d:Landsat TMsatellite imagery)

从8个支流情况看,支流1为凤眼莲最初出现的地点,且每年都有出现,2003年前面积较小,但2003年12月之后面积增加,从面积计算结果看出,大部分时间该支流上的凤眼莲总面积都在105m2以上;支流3(尤溪)在凤眼莲大暴发的2003、2004、2009年出现数量较多;支流5(武步溪)分别在2003年和2004年的12月份、2007年的1月和5月份及2009年分布面积较大,2006年虽有出现,但分布面积很小;在研究范围内,支流7(古田溪)凤眼莲于2003、2004、2006、2007、2009年均有出现,且从研究数据看出2003-2004年的暴发地区从上游逐渐向下游迁移,2006年后该支流的凤眼莲仅在Z7-3出现,但除了2007年的5-9月,其分布面积均在9×104m2以下;支流8仅在2007年5月份和2009年凤眼莲占据河道面积较大;支流2、支流4、支流6出现的次数不多且占据面积小.从以上支流情况看出,该研究区域除了支流1为源头外,河流较长,河岸有居民居住的支流(如支流3、5、7)也对凤眼莲的分布面积有重要影响.此外,干流中一些与支流交汇处较接近的点也是受灾频率较多的地区,根据金樑等的研究表明,流域中成灾的凤眼莲并不是当地自然产生,而是从黄浦江的上游与其支流随水流而下并在某处聚集形成[4].由此看出,干流中的凤眼莲大部分从支流漂流而至.因此,在河道中,支流是凤眼莲源头治理的重点.

有学者[22]将同伴的修改意见按意义单元切分,意义单元是指关于单个的同伴写作的自包含信息; 然后制定编码方案,对评论片段进行分类,为了测量评判者间信度,两个编码者独立编码评论,每个编码者先分别编码部分意义单元(大约是总意义单元的15%),确定编码者的信度(包括评价和范围两个维度),通过讨论解决编码者之间的不一致。

2004年10月和12月对比,支流1凤眼莲面积变化不大;支流3凤眼莲数量急剧增加;支流5凤眼莲数量增加,在12月占据1/2河道;支流7凤眼莲面积急剧减少,仅Z7-3点在12月有极少量凤眼莲;干流2、干流3和干流5段凤眼莲面积增加.

2006年支流1的凤眼莲所占面积较多,特别是点Z1-2在6-9月占据整个河道,1月和11月份其它污染物占据整个河道,凤眼莲所占面积较少;点Z1-1在1月份有少量凤眼莲,但6月份消失,8月份再次出现,且面积逐月增加;Z1-3点位于干流交界口,河水具有流动性,凤眼莲随水流漂向下游,因此数量相对较少.Z4点在6月份有大量凤眼莲聚集在岸边.G3-2点在6月12日的影像上有一定量的凤眼莲附着在岸边.除此之外,2006年干流总体趋势是1-9月凤眼莲面积较少,11月份部分观察点的凤眼莲面积大幅度增加(图4、图5).

图4 2006年干流凤眼莲面积Fig.4 Area of Eichhornia crassipes in mainstream in 2006

图5 2006年支流凤眼莲面积Fig.5 Area of Eichhornia crassipes in tributaries in 2006

2007年凤眼莲的面积从总体趋势上看,5月份凤眼莲较多,5-7月下降,8-9月略微增加.特别是点G7-1、Z7-3和Z5-2,在5月份占据整个河道,至7月份骤减.根据新华社报道①http://fj.news.163.com/07/0418/11/3CC17DD1004400BA.html.,2007年3月,水口水电站的坝头上凤眼莲泛滥,5月份针对凤眼莲进行了人工清理,因此出现5-7月份河道中凤眼莲面积下降,8-9月份面积略微上升的现象.以支流的总体情况看,除点Z5-2、Z7-3和Z8-2,其余观察点1-5月份凤眼莲面积减少,7-9月份面积逐渐增加(图6);以干流的总体情况看,流域上游(干流1~干流5)在1月份凤眼莲面积较5月份多,下游(干流6~干流7)1月份面积较5月份少(图7),说明凤眼莲具有漂流特性,从上游漂流至下游.

2009年河道中凤眼莲所占面积相对往年更多,特别是支流1、支流3、支流8和干流6、干流7.在6月份点G6-1和点G7-2凤眼莲占据整个河道,受灾水域长度近2 km,船舶通行严重受阻,是凤眼莲泛滥最为严重的一次,但到7月份,除点Z3-3、点Z5-2和支流8,凤眼莲面积总体上大量减少,因此认为在6月份凤眼莲暴发后,受灾区经过人工清理.由于2009年入春来,气温偏高,使得原本生长快的凤眼莲加速繁殖.导致2009年是2004年以来情况最严重的一年.

研究区域中凤眼莲面积在不同年份相同时间段有一定差异:2004年10月约是2005年10月的1.3倍,2007年8月约是2006年8月的1.65倍以上,2007年9月约是2006年9月的1.67倍,2009年7月约是2007年7月的2倍.

总体来看,2001-2009年、2003年12月和2004年的凤眼莲所占面积较多,2005年和2006年除个别地区,凤眼莲总体面积不多,2007年凤眼莲面积再次上涨,至2009年凤眼莲面积加倍上涨,是2004年来情况最严重的一年.

图6 2007年干流凤眼莲面积Fig.6 Area of Eichhornia crassipes in mainstream in 2007

图7 2007年支流凤眼莲面积Fig.7 Area of Eichhornia crassipes in tributaries in 2007

2.2 研究区域水质状况

据调查,研究区域上游的沙溪河段是闽江流域工业最为集中的区域,沿途污染源密集,接纳大量畜禽养殖、生活和工业废水,成为流域内污染最严重的河段.根据2005年福建省政府转发的关于闽江流域水环境保护规划文件显示,2004年沙溪流域优于Ⅲ类水质比例和水环境功能达标率均为61.5%,主要超标项目有氨氮、五日生化需氧量(BOD5)、溶解氧.研究区域水质环境功能均为Ⅲ类标准,为工业、渔业、景观娱乐用水(表1)[15].从南平市环境保护局发布的南平市水质状况看出,2006年前研究区域及其上游水质基本未达标,2006年后水质大部分达标(表2).

表1 研究区域水环境功能区划Tab.1 Regionization of Water Environmental Function in survey region

表2 南平市水域水质状况Tab.2 Water quality conditions of rivers in Nanping City

2.3 研究区域凤眼莲频发成因分析

考虑到2006年前上游污染较严重,导致上游水质恶化,这可能是上游凤眼莲暴发多在2006年前的一个原因.研究发现,水环境中的营养水平对凤眼莲的生长繁殖有明显影响,当水体中有足够的氮、磷养分,并具备合适的光照,温度和氧气等外因时,就会引发凤眼莲的快速生长[16-17].并且,凤眼莲的耐受性强,在水质很差的环境下(如工业废水)也能正常生长,其在污染水体中的生长速度为天然水体中的2倍,并能有效地竞争并排斥其他植物物种,形成大面积的纯种群[18-19].2006年后,凤眼莲多在下游暴发,可能是上游的凤眼莲随着水流往下漂流的原因.

资料显示①《福建省人民政府办公厅转发省环保局关于闽江流域水环境保护规划的通知》.,闽江流域有多个大中型水电站,还有大量无序开发的小水电站,根据水利部门提供的信息,至2004年初,南平市已建成各类水电站七百多座.过度和不合理的水电开发使自然急流变为人工平湖,河流流速变缓,自净能力降低,污染物淤积,水质变差,库区富营养化程度加剧.邬红娟等[20]认为湖泊、河流和水库的水滞留时间对浮游植物生物量和生产力有很大影响,当流速超过浮游植物繁殖速率,生物量的积累和生产力就会流失.再加上研究区域多为渔业养殖用水,河道间有较多的养殖网箱,这不仅促进水流变缓,还有利于凤眼莲在网箱周围生长.金樑等以上海黄浦江为研究区,发现凤眼莲在水体中的固着机制与河流的整齐度以及河道内其它水生植物的分布存在很大关系,河岸不整齐,近岸边生有许多水生植物的河道和流水较慢的小河流易于成为凤眼莲生长和繁殖的发源地[4,21].

对照实际情况,观察凤眼莲出现频率较多的地点(表3),发现这些观察点的共性:①沿岸多为村庄或乡镇,常有居民生活污水和农业污水排入河流,且多处河道进行渔业养殖,使得水体氮、磷元素过多,水体富营养化.根据蔡雷鸣[13]2001年与2004年对水口库区水质指标的监测结果显示,2004年水中的总氮、总磷含量均超过Ⅲ类标准.同时结合福州监测站2006-2009年对水口库区雄江镇的水质监测结果,2007-2009年个别月份的总磷含量超标.由此可见,库区部分地区已呈现富营养化趋势.②河道结构较为复杂,河流沿岸角落处容易成为死角,河水在此处不流动.

表3 2001-2009年各观察点出现凤眼莲的频率Tab.3 The frequency of occurrence of Eichhornia crassipes from 2001 to 2009

2001年1月起,研究区域内就有个别地区存在凤眼莲,但数量极少.2001年5月水口坝前已覆盖有少量凤眼莲.2003年12月研究区域凤眼莲大面积暴发,多段河道被凤眼莲占据.此后,2003、2004、2007、2009年均为凤眼莲扩散较严重的年份,特别是2004年和2009年.中国气象科学数据库中获得的南平市气象资料表明,2003、2004、2007、2009年正好是年降雨量相对较少的年份(表4),由于降雨量小、水流量小、流速较缓,有利于凤眼莲固着生长,并且水体富营养化浓度升高,导致凤眼莲快速生长.同时2003、2004、2007、2009年也是日照时数或年平均气温相对高的年份,由于凤眼莲性喜光,适应于热带、亚热带气候,在光照强、温度为25 ~35℃气候条件下时,植株成长迅速[15,22].

表4 2001-2009年研究区气象数据Tab.4 Meteorological data in survey region from 2001 to 2009

3 总结与讨论

本文运用2001-2009年的卫星遥感影像,对比分析得出水口水库中凤眼莲的时空分布情况:

1)2001年水口水库已有少量凤眼莲存在,2003年12月凤眼莲已在水库库区大量暴发,其中2004年、2007年、2009年凤眼莲占据河道面积较多,特别是2004年和2009年.

2)每年凤眼莲数量较多的时间段出现在5、6月份和9-12月份.原因为5、6月份凤眼莲大量繁殖,影响水上交通,当地政府和环保部门针对其进行人工清理,使之数量减少.但由于清理未彻底,剩余的凤眼莲继续生长,使9-12月份凤眼莲再次泛滥.

3)随着时间的推移,凤眼莲的暴发区逐渐向下游迁移.从总体上看,上游凤眼莲暴发多在2006年前,下游凤眼莲暴发多在2006年后.

4)支流应该是凤眼莲源头治理的重点.支流1是凤眼莲最先出现的地点,也是出现次数最多的地点,可以将其认为是凤眼莲最初产生点之一,首先要对支流1的凤眼莲进行重点清除.除此之外,河流较长、河岸边有居民居住的支流,如支流3(尤溪)、支流5(武步溪)和支流7(古田溪)凤眼莲出现的频率高、面积大,对研究区域凤眼莲的面积有重要影响,因此,应将支流3、5、7视为凤眼莲重要来源,采取相应措施.

5)凤眼莲泛滥有多方面原因,除了客观因素(闽江上游出现罕见的少雨干早、气温偏高的气候现象),沿江一带养殖业的发展、工业废水、人类生活污水的排放造成水体氮、磷元素过多以及过度和不合理的水电开发使自然急流变为人工平湖,河流流速变缓等都是重要原因.

若要预防凤眼莲年年泛滥,应每年都对凤眼莲积极清理,特别是支流和易形成死角的岸边,建议在前一年冬天开始作业,可达到事半功倍的效果.同时对沿岸排放的污水进行达标处理,加强研究源头综合治理措施,推广库区上游低污染的生产方式.根据福建省人民政府的要求,2010年闽江南平段水域功能区达标率应在95%以上,2020年全流域水质全面稳定达到或优于功能分区的环境质量标准.为达到研究区域水环境的保护,应针对研究区域上游及其周边,建设生态功能保护区,力争使流域的源头区、重要水源涵养区、水土保持重点预防保护区,重点监督区和重要渔业水域等重要生态功能区的生态系统和生态功能得到保护与恢复;在切实抓好现有自然保护区建设与管理的同时,建设一批新的自然保护区,使各类良好自然生态系统及重要物种得到有效保护;建立、健全生态环境保护监管体系,使生态环境保护措施得到有效执行,重点资源开发区的各类开发活动严格按规划进行,生态环境破坏恢复率有较大幅度提高;加强生态示范区和生态农业县建设,并继续开展城镇环境综合整治,进一步加快能源结构调整和工业污染源治理[22],重点完成以下工程:

1)养殖污染整治:严格落实畜禽养殖“两禁区”规定,全面完成禁养区内养殖场搬迁拆除工作,加快禁养区外规模养殖场的污染治理,特别是研究区域上游的畜禽污染综合整治,并加强先进技术推广,防止出现污染反弹.此外,研究区域多处进行渔业养殖,应加强重点渔场及其他渔业水域等重要水生资源繁育区的保护,严格渔业资源开发的生态环境保护监管.

2)工业污染整治:结合污染减排,加强工业污染整治,加快推进工业开发区、重点行业企业污染治理,针对研究区域上游沙溪流域为闽江流域工业最为集中的区域,应加强工业污染排放监控,强化对污染源的监督管理,并停止一切产生严重环境污染的工程项目建设.

3)生态环境建设与保护:以保障流域生态功能为目标,推进矿山整顿、水土流失治理、生态林保护和自然保护区建设,整顿规范河砂开采,落实水电站生态下泄流量要求,对于布局不合理,影响生态平衡的水库进行整改或拆除.加强库区水葫芦及漂浮垃圾的清理工作,开展小流域污染整治,大力推进农村环境整治,引导和推动县市及乡村开展生态创建工作.

[1]汪凤娣.外来入侵物种凤眼莲的危害及防治对策.黑龙江环境通报,2003,27(3):21-23.

[2]任明迅.水生植物凤眼莲的种内基因型差异与入侵模式[学位论文].北京:中国科学院植物研究所,2006.

[3]庄先林,郑 树,吴世铭等.福建省水口库区水葫芦专项整治研究.海峡科学,2007,(7):76-77.

[4]金 樑,王晓娟,高 雷等.上海市凤眼莲种群的时空分布及控制对策.生态学杂志,2005,24(12):1454-1459.

[5]王 云,周国庆,周国峰等.洞庭湖区水葫芦分布及危害的评价.草业科学,2008,25(3):97-101.

[6]王 瑞,王印政,万方浩.外来入侵植物一年蓬在中国的时空扩散动态及其潜在分布区预测.生态学杂志,2010,(6):1068-1074.

[7]Cavalli RM,Laneve G,Fusilli L et al.Remote sensing water observation for supporting Lake Victoria weed management.Journal of Environmental Management,2009,90:2199-2211.

[8]章 莹,卢剑波.外来入侵物种互花米草(Spartina alterniflora)及凤眼莲(Eichhornia crassipes)的遥感监测研究进展.科技通报,2010,26(1):130-137.

[9]沈永明,刘咏梅,陈全站.江苏沿海互花米草盐沼扩展过程的遥感分析.植物资源与环境学报,2002,11(2):33-38.

[10]Everitt JH,Yang C,Escobar DE et al.Using remote sensing and spatial information technologies to detect and map two aquatic macrophytes.Journal of Aquatic Plant Management,1999,37:71-80.

[11]Venugopal G.Monitoring the effects of biological control of water hyacinths using remotely sensed data:A case study of Bangalore,India.Singapore Journal of Tropical Geography,1998,19(1):91-105.

[12]袁信芳,施华宏,王晓蓉.太湖着生藻类的时空分布特征.农业环境科学学报,2006,25(4):1035-1040.

[13]蔡雷鸣.福建闽江水口库区飘浮植物覆盖对水体环境的影响.湖泊科学,2006,18(3):250-254.

[14]Adams JB,Smith MO,Johnson PE.Spectral mixture modeling:A new analysis of rock and soil types at the Viking Landerl site.Journal of Geophysical Research,1985,91(B8):8098-8112.

[15]程永隆等.福建省水(环境)功能区划.福建:福建省水利厅,福建省环境保护局,2004.

[16]杨红玉,张江山,仝 川等.水环境中氮、磷营养水平对水葫芦生长的影响.安徽农学通报,2008,14(23):67-71.

[17]张秋芳,刘 波,官雪芳等.水葫芦对氮磷营养的需求及分布规律初探.中国农学通报,2007,23(5):387-389.

[18]齐玉梅,高伟生.凤眼莲净化水质及其后处理工艺探讨.环境科学进展,1999,7(2):136-140.

[19]杨凤辉.上海黄浦江凤眼莲灾害的发生机理及控制对策初探.复旦学报,2002,41(6):599-603.

[20]邬红娟,郭生练.水库水文情势与浮游植物群落结构.水科学进展,2001,12(1):51-55.

[21]孙 超,陈振楼,毕春娟等.黄浦江上游水源保护区不同河段水葫芦最佳治理对策.杂草科学,2008,(3):14-17.

[22]新华社.全国生态环境保护纲要[EB/OL].http://www.people.com.cn/GB/channel1/907/20001222/358911.html.

Spatial distribution characteristics and dynamics of Eichhornia crassipes in the Shuikou Reservoir,Fujian Province

CHEN Xiao,PAN Wenbin&WANG Mu
(College of Environment and Resources,Fuzhou University,Fuzhou 350108,P.R.China)

In order to study the spatial distribution characteristics and dynamics of Eichhornia crassipes in the Shuikou Reservoir,this study,using 24 unclouded remote sensing satellite images from 2001 to 2009 which covered the Shuikou Reservoir region,reconstructed the spatial structure of Eichhornia crassipes distributed in the mainstream and tributaries,and analyzed the situation of Eichhornia crassipes in the Shuikou Reservoir at different times and different observation zones according to the Linear Spectral Mixture Analysis Method.The results indicate that the outbreak area of Eichhornia crassipes moves from the upper reaches to the lower reaches gradually as time goes on.And some long tributaries with high riparian populations played an important role in the contribution of Eichhornia crassipes in this survey region.The significant factors about the fast growth of Eichhornia crassipes included water eutrophication,low velocity,etc..It is suggested that local government should take actions to fight against the pollution from the industry,aquaculture,and poultry farming in this area,enhance the supervision and protection of ecological environment,and establish nature reserves.This study can enhance the understanding of the temporal and spatial distributions of Eichhornia crassipes,and then provide fundamental data and theory for managing and controlling the nuisance plants.

Eichhornia crassipes;Shuikou Reservoir;spatial distribution characteristics and dynamics

http://www.jlakes.org.E-mail:jlakes@niglas.ac.cn

ⓒ2012 by Journal of Lake Sciences*2011-06-09收稿;2011-09-25 收修改稿.陈潇,女,1986 年生,硕士研究生;E-mail:amethyst_cx@126.com.

** 通信作者;E-mail:WenbinPan@fzu.edu.cn.

猜你喜欢
凤眼莲支流水口
水口水电站7号机组完成增容改造和并网发电
水生植物组合培植对水体中重金属富集效果研究
老政法把水葫芦变“宝葫芦”
老政法把水葫芦变“宝葫芦”
水口、水口环境与风水(下篇)
千年沧桑话水口
光影视界
阿什河某支流河岸带草本植物重要性初步研究
金沙江支流东川玉碑地遗址
基于富营养化水体修复的凤眼莲放养及采收条件研究