基于水声学方法的天目湖鱼类季节和昼夜空间分布研究

2015-01-18 08:48孙明波谷孝鸿曾庆飞王银平毛志刚谷先坤
生态学报 2015年17期
关键词:天目湖鱼类季节

孙明波, 谷孝鸿, 曾庆飞, 王银平, 毛志刚, 谷先坤

1 中国科学院南京地理与湖泊研究所湖泊与环境国家重点实验室, 南京 210008 2 中国科学院大学, 北京 100049

基于水声学方法的天目湖鱼类季节和昼夜空间分布研究

孙明波1,2, 谷孝鸿1,*, 曾庆飞1, 王银平1,2, 毛志刚1, 谷先坤1,2

1 中国科学院南京地理与湖泊研究所湖泊与环境国家重点实验室, 南京 210008 2 中国科学院大学, 北京 100049

采用水声学方法系统地对天目湖春夏秋冬四季鱼类水平和垂直空间分布进行了调查研究,同时对昼夜探测的差异性进行了比较,并对不同季节间鱼类聚群形态进行了探讨。研究结果表明,天目湖鱼类在不同季节间存在规律性水平迁移,从全湖角度分析,天目湖鱼类资源昼夜间的水平空间分布无明显差异特征,但在春末和夏季的局部区域里,昼夜间鱼类存在一定近岸—远岸的水平迁移;不同季节和昼夜间,鱼类垂直分布差异明显,且在存在温跃层的夏季7月昼间,鱼类密度垂直分布与水温的垂直分布关系密切,温跃层以下的鱼类密度基本为0;天目湖鱼类在四个季节都属于成群分布类型,但季节间鱼类聚群形态不同,在冬季的1月呈现出典型的聚集,相对于昼间,春夏季鱼类在夜间分布更为均匀、分散,且水平探测表明在夜间水体表层存在大量鱼类分布。

水声学; 鱼类; 空间分布; 季节; 昼夜

在自然条件下,鱼类并不是均匀随机的分布在水体中。在水体中不同的时空会形成特定的生境,鱼类也会根据生境因子进行索饵、洄游、栖息、聚群、繁殖等行为,且随鱼类种群不同而各有特点。季节变化和昼夜变化都会引起鱼类的各种垂直和水平迁移行为,如何在湖泊等大尺度水体对鱼类这种季节和昼夜空间分布变化进行研究是传统鱼类调查方法面临的一大难题,随着计算机技术和硬件性能的提高,水声学技术为相关研究提供了一种有效途径。

水声学方法具有快速高效、调查区域广、不损害生物资源、提供连续数据、自然状态下定位鱼类空间分布、准确估算鱼类密度和资源量等优势[1],该方法在国外有广泛应用[2- 6],国内也有了初步的发展[7- 11],而近年来采用水声学方法对鱼类在水体中的季节和昼夜空间分布的研究也得到了一些进展,如在对维多利亚湖的调查中,Albert等[12]发现不同季节中鱼类时空分布有明显的差异,在中心湖区8月份测得的鱼类密度要显著高于2月份。Axenrot和Guillard等[13- 15]在研究中则对鱼类昼夜分布特性进行了探讨,认为夜间鱼类个体较为均匀的分散在水体中,测得的单体目标强度更为准确。Kubecka和Wolter 等[16-17]的研究表明鱼类昼夜间在远岸区和近岸区有明显的水平迁移行为,并且不同种群体现出不同的行为。Frank和Kubecka等[16, 18]的昼夜调查研究中也发现鱼类分布在夜间要比昼间更分散。虽然国外分别进行了部分不同季节和昼夜的鱼类资源水声学调查研究,但在我国湖泊中仍缺乏春夏秋冬四个季节昼夜间的调查,对这方面进行系统研究非常必要,是鱼类行为生态学研究的基础,对水体渔业生产具有现实指导意义。

1 研究区域和研究方法

1.1 研究区域概况

图1 天目湖区域图Fig.1 The map of Lake Tianmu

天目湖(119.419747°E,31.291017°N,图1)位于江苏省溧阳市南部丘陵地区,集水域面积148.5 km2,湖泊面积9.8 km2,最大水深16 m,平均水深6—8 m,南部有平桥河、中田河、下宋河流入。天目湖以放养的鲢(Hypophthalmichthysmolitrix)和鳙(Aristichthysnobilis)为主要优势种[19]。2011—2012年冬季由于鱼类活动能力下降,天目湖采用机动船只拖曳白板树枝的方式将上下游鱼类驱赶至布置在中游的集鱼网箱,然后在赶鱼结束后的冬末春初陆续进行渔业捕捞与放流。2012年春季、夏季和秋季的4月—11月则采用在湖区上中下游布置9个网簖的方式进行捕捞。放流鱼类主要是规格在0.4 kg/尾的鲢和鳙以及规格约0.2 kg/尾的长春鳊(Parabramispekinensis)。

1.2 水声学调查方法

本研究从2011年冬季到2012年秋季共对天目湖进行6次鱼类资源水声学调查(表1),水声学探测仪器为Biosonics DT-X型回声探测仪(半功率波束角6.5°的分裂波束换能器,工作频率为208 kHz)。采用Garmin GPS 17x HVS对位置数据同步采集存储,使用BioSonics Acquisition 6.0软件进行水声学数据采集,脉冲收发周期为0.125 s,脉冲宽度为0.5 ms,目标强度信号数据收集阈值为- 130 dB。使用快艇进行走航调查,并利用Garmin Oregon450导航仪进行航线导航,采取“之”字形航线走航调查,先进行垂直探测,将换能器用铁架固定于船舷,吃水0.4 m,数据收集距离为0.5—20 m,考虑垂直探测受近场和表层鱼类回避行为的影响,本文探索在返航时进行水平探测,并将其与垂直探测拟合计算单位面积鱼类密度[18, 20],水平探测时同样将换能器用铁架固定于船舷,吃水1 m,并给予3°向下的倾角,垂直于船体走航方向探测,数据收集距离为1—30 m。调查航速约为10 km/h。在调查前使用厂家原配的直径36 mm的碳化钨标准球对仪器进行实地校准。

按照Aglen[21]覆盖率公式计算3次水声学调查的覆盖率:

(1)

式中,L(m)为水声学调查走航航程,A(m2)为湖泊水面面积,Dc为水声学调查覆盖率。调查信息如下表1,其中2011年12月,天目湖正在上游进行冬季捕捞赶鱼作业,2012年1月是在捕捞赶鱼作业结束后(鱼类资源调查为不包括集鱼网箱的湖区鱼类资源),2012年3月是在捕捞与放流结束后。

表1 6次鱼类资源水声学调查基本信息

1.3 水声学数据处理和GIS建模

对采集到的水声学数据用BioSonics Visual Analyzer 4.1进行分析,其中单体回声检测识别的设定参数: Echo Threshold(回波阈值)- 60 dB,Correlation Factor(相关系数) 0.9,Min Pulse Width Factor (最小脉冲宽度系数)0.75,Max Pulse Width Factor (最大脉冲宽度系数) 3,End Point Criteria(回波时间计算点)- 12 dB,TVG(Time-Varied Gain, 时变增益)为40lgR。以从下游至上游每300个脉冲发射时间所航行的距离(约110m)作为一个回波数据分析单元。垂直探测波束数据分析范围1—20 m,分析结果包括湖区单位面积鱼类个体数量FPUA(Fish Per Unit Area, 尾/m2)、单位体积鱼类个体数量FPCM(Fish Per Cubic Meter, 尾/m3)、单元起始坐标、单体识别结果、单元平均水深等,另记下每个单元的中心脉冲坐标。水平探测数据分析波束从5 m开始,终止距离根据接近湖岸、船只摆动、水底凸起、水面回波、水底回波等实际情况人工划定,有效分析结果为FPCM,代表上层水体(0—4m)的鱼类密度(FPCM)。

由于鱼类在秋冬季倾向于深水越冬和反应减缓,根据水声学调查结果进行鱼类密度的统计分析,若上层水体垂直探测(1—4m)和水平探测的鱼类密度无显著性差异(P>0.05),只采取垂直探测(1—20m)的分析作为调查结果。若两者差异性显著,采取以下方式进行拟合计算各单元调查结果。

FPUA=FPCMh×h(h≤4)

(2)

FPUA=FPUAv+FPCMh×4 (h>4)

(3)

FPCM=FPUA/h

(4)

式中, FPUAv(尾/m2)为分析单元垂直探测(4—20m)的单位面积鱼类尾数,FPCMh(尾/m3)为分析单元水平探测单位体积鱼类尾数,h(m)为分析单元的平均水深。

采用ArcGIS10.0软件进行鱼类资源分布的建模[9],将分析计算出的FPUA、单元中心坐标GPS数据导入ArcGIS平台,采用IDW方法进行栅格插值运算[22-23]。设定栅格大小为0.0001°,大小约为11.5 m×9.5 m,设定幂值p=2,做出鱼类分布模型图。

另在天目湖下游大坝处的深水区设置监测区域P(图1),该区域为1000个脉冲的航行距离(约360m),采用垂直探测的数据研究该区域鱼类垂直密度变化,由于上下层水体升温和降温不一致,在深水水体中水温和溶解氧都会在垂直空间呈现一定的变化,本研究在1月、5月、7月对1号点的垂直温度和溶解氧采用水质分析仪YSI进行实时监测,12月、3月、10月采用天目湖同期监测资料,以期探讨鱼类密度垂直分布与温度、溶解氧的关系。

1.4 使用的经验公式和分析软件

所进行的数据描述性统计、正态性检验、非参数检验、相关性分析均采用SPSS 17.0。鱼类水平分布建模采用ArcGIS10.0。

2 结果与分析

2.1 天目湖鱼类水平分布

图2 天目湖鱼类资源水平分布Fig.2 The horizontal distribution of fish resource in Lake Tianmu

本研究对天目湖进行了春夏秋冬四个季节的鱼类资源水声学调查,鱼类资源在水平空间的分布随季节变化明显,采用GIS对其建立分布模型如图2。鱼类水平分布在一定程度上受到不同季节的渔业活动影响,2011年的冬季12月份天目湖进行赶鱼作业,受其影响鱼类呈现逃逸到中上游趋势,导致下游接近大坝处的鱼类资源量较低;而2012年的冬季1月份赶鱼结束后鱼类资源量在下游处的分布明显较高,应是鱼类返回下游的深水区越冬;2012年的春季初3月份鱼类仍较集中分布在下游的深水区;2012年春末5月份和夏季7月份鱼类则明显向中上游扩散,较为均匀的分散在湖区中,上游鱼类密度偏高;2012年秋季10月份鱼类在湖区仍分散分布,但相对7月份鱼类资源量在下游深水区分布偏高。

由此可得出,天目湖鱼类在春夏秋冬四个季节存在一定规律的水平迁移,春季鱼类由下游深水区向上游迁移,夏季鱼类较为均匀的分布在湖区中,秋季鱼类向下游深水区迁移,冬季鱼类在下游深水区越冬。

在3、5、7、10月的鱼类资源水平分布图(图2)中,鱼类资源昼夜分布在全湖水平空间上无明显差异特征,而在天目湖发电站的取水口处发现昼间鱼类资源量高于夜间。对近岸区(距岸5 m以内的区域)的鱼类密度(FPCM)进行分析,5月和7月夜间近岸区单位体积计的鱼类密度显著(P<0.05)低于昼间,而在其他月份无显著性(P>0.05)差异,这表明在春末和夏季天目湖鱼类在昼夜间存在一定的近岸—远岸水平迁移。

2.2 天目湖鱼类垂直分布

比较水体上层(0—4 m)鱼类密度与中下层(4 m以下)鱼类密度(图3),在冬季的12月和1月昼间,中下层鱼类密度显著(P<0.05)高于上层;春季3月的昼夜探测结果都为中下层鱼类密度与上层鱼类密度无显著差异(P>0.05);春季5月和夏季7月的昼间中下层鱼类密度显著(P<0.05)高于上层,但夜间则相反;秋季的10月昼夜探测都为中下层鱼类密度显著(P<0.05)高于上层。

对天目湖最深处大坝前P号点监测(图4),夏季(2012年7月)和春末(2012年5月)时,在天目湖的深水区存在明显的温跃层,温跃层下温度和溶解氧都较低,而其他月份这种现象不明显,冬季各水层水温和溶解氧变化不大。对温跃层明显的夏季7月份进行鱼类垂直分布与水温的关系研究,如图5所示,鱼类密度垂直分布与水温、溶解氧密切相关,水深7m以下鱼类密度基本为0,而这部分水层也完全处于温跃层之下。

2.3 天目湖鱼类聚群形态

6次调查探测到的鱼类密度变异系数都大于1,根据孙儒泳[24]种群分布型的标准划分,各调查阶段鱼类都属于成群分布类型。2011年12月昼、2012年1月昼、2012年3月昼(夜)、2012年5月昼(夜)、2012年7月昼(夜)、2012年10月昼(夜)调查实测水温为11.5、4.1、14.8(13.2)、21.4(20.9)、31.6(31.4)、18.7(18.3℃),对相应月份的的回波图像分析,鱼类随着温度变化在空间上展现出不同的聚群行为,在温度较低的冬季12月鱼类有少量聚群现象(图6A),而在温度最低的冬季1月出现了鱼群的典型团聚(图6B),主要出现在4—7m的中下层,并且在1月鱼类密度的变异系数也最大,而后的春季3月和5月、夏季7月和秋季10月随着温度上升鱼类无较多成群团聚出现,仅有少量如图6C的聚群存在,在5月和7月鱼类相对均匀的分散在水体中(如图6D)。在2012年的春夏季的调查中,对昼夜探测得到回波图像分析,鱼类在夜间时更为均匀地分散在水体中,在敞水区分布较多(图6E),且水平探测表明在夜间存在大量鱼类分散在表层水体中(图6F),鱼类在夜间更为活跃。

图3 天目湖鱼类在不同水层的垂直分布比较 Fig.3 The comparison of fish vertical distribution between different water layers in Lake Tianmu

图4 水温、溶解氧与水深的关系Fig.4 The relationships between water depth and water temperature, DO

图5 7月鱼类密度与水温、溶解氧的关系Fig.5 The relationships between fish density and water temperature, DO in July

3 讨论

3.1 天目湖鱼类水平分布变化及影响因素

图6 鱼类成群分布形态Fig.6 Fish aggregations of different seasons

3.2 天目湖鱼类垂直分布变化及影响因素

在秋季10月和冬季12月、1月,天目湖中下层鱼类密度都要显著(P<0.05)高于上层,这种垂直分布同样体现了鲢鳙深水越冬的生态习性;但在春夏季天目湖鱼类的垂直分布较为复杂,初春3月昼夜调查的结果表明鱼类密度在上层和中下层无显著差异(P>0.05),可能是由于3月份刚放养了大量幼鱼,其在垂直水层间的游动较为活跃所致;鲢鳙在春末、盛夏、秋初索饵行为非常活跃,喜集群游至水域的中上层,尤其是水质较肥的明水区[25-26],但本研究中春末5月和盛夏7月的昼间中下层鱼类密度却显著(P<0.05)高于上层,这种不一致性可能是由于天目湖属于旅游景区,昼间来往游船较多和水上娱乐项目造成的,另外调查中船只对鱼类垂直分布的干扰因素也不容忽视,这些影响因子都降低了天目湖鱼类在上层的分布,夜间的调查表明上层鱼类密度要高于中下层则呈现出了天目湖以鲢鳙为主鱼类的生态习性[25-26]。已有的研究[29]表明天目湖存在典型的温跃层现象且随季节变化明显,夏季尤为强盛,冬季基本消失。本研究中也发现了这种规律,在水层垂直水温变化明显的7月,通过对实时同地监测的鱼类密度垂直分布分析,鱼类密度垂直分布与水体垂直水温变化密切相关,温跃层以下鱼类密度基本为0,可能是鱼类对这一区域的低温和低溶解氧的回避行为,这与Drastik[28]的研究结果具有一致性。

3.3 天目湖鱼类聚群形态及对鱼类水声学调查方法的影响

根据变异系数划分,天目湖鱼类在空间上为成群分布,但这种区域成群可能是受多种环境因子影响的结果,并不都是鱼类主动聚群。通过春夏秋冬四个季节不同温度下的调查探测,本研究监测到了随着温度降低而出现的鱼类空间聚群的不同现象(图6),在冬季鱼类出现了明显的团聚,且在1月份团聚程度最高,这种聚群与谭细畅和Churnside等[30-31]发现的鱼类主动聚群回波图像具有相似性,应该是鱼类随着温度降低逐渐出现的鱼类聚群行为,调查的回波图像也显示随着春夏季温度升高,这种团聚逐渐消失。春夏季的昼夜调查回波图像表明鱼类在夜间更为分散,在这种状态下更适合鱼类个体水声学探测[13- 15],但春夏季在水体表层分布着大量鱼类,建议鱼类资源水声学调查同时进行垂直探测和水平探测,若只进行垂直探测则可能会因鱼类对船只的回避行为造成评估的资源量偏低[18]。

[1] Simmonds J, MacLennan D. Fisheries Acoustics: Theory and Practice. 2nd ed. Oxford: Blackwell Publishing, 2005: 1- 8.

[2] Auvinen H, Jurvelius J. Comparison of pelagic vendace (Coregonusalbula) stock density estimation methods in a lake. Fisheries Research, 1994, 19(1/2): 31- 50.

[3] Shirakihara K, Yoshida M, Nishino M, Takao Y, Sawada K. Acoustic evaluation of the vertical distribution of dwarf ayuPlecoglossusaltivelisaltivelisin Lake Biwa. Fisheries Science, 2001, 67(3): 430- 435.

[4] Wanzenbock J, Mehner T, Schulz M, Gassner H, Winfield I J. Quality assurance of hydroacoustic surveys: the repeatability of fish-abundance and biomass estimates in lakes within and between hydroacoustic systems. ICES Journal of Marine Science, 2003, 60(3): 486- 492.

[5] Fabi G, Sala A. An assessment of biomass and diel activity of fish at an artificial reef (Adriatic sea) using a stationary hydroacoustic technique. ICES Journal of Marine Science, 2002, 59(2): 411- 420.

[6] Hughes S. A mobile horizontal hydroacoustic fisheries survey of the River Thames, United Kingdom. Fisheries Research, 1998, 35(1/2): 91- 97.

[7] 王珂, 段辛斌, 刘绍平, 李志华, 陈大庆. 三峡库区大宁河鱼类的时空分布特征. 水生生物学报, 2009, 33(3): 516- 521.

[8] 张慧杰, 杨德国, 危起伟, 杜浩, 张辉, 陈细华. 葛洲坝至古老背江段鱼类的水声学调查. 长江流域资源与环境, 2007, 16(1): 86- 91.

[9] 陶江平, 艾为明, 龚昱田, 陈志俭, 李尚鲁, 谢起浪, 陈少波, 周志明. 采用渔业声学方法和GIS模型对楠溪江鱼类资源量及空间分布的评估. 生态学报, 2010, 30(11): 2992- 3000.

[10] 谭细畅, 史建全, 张宏, 陶江平, 杨建新, 祁洪芳, 李新辉. EY60回声探测仪在青海湖鱼类资源量评估中的应用. 湖泊科学, 2009, 21(6): 865- 872.

[11] 陈国宝, 李永振, 赵宪勇, 陈毓桢, 金显仕. 南海北部海域重要经济鱼类资源声学评估. 中国水产科学, 2005, 12(4): 445- 451.

[12] Getabu, A, Tumwebaze R, MacLennan D N. Spatial distribution and temporal changes in the fish populations of Lake Victoria. Aquatic Living Resources, 2003, 16(3): 159- 165.

[13] Axenrot T, Didrikas T, Danielsson C, Hansson S. Diel patterns in pelagic fish behaviour and distribution observed from a stationary, bottom-mounted, and upward-facing transducer. ICES Journal of Marine Science, 2004, 61(7): 1100- 1104.

[14] Guillard J, Lebourges-Dhaussy A, Brehmer P. Simultaneous Sv and TS measurements on Young-of-the-Year (YOY) freshwater fish using three frequencies. ICES Journal of Marine Science, 2004, 61(2): 267- 273.

[15] Freon P, Gerlotto F, Soria M. Diel variability of school structure with special reference to transition periods. ICES Journal of Marine Science, 1996, 53(2): 459- 464.

[17] Wolter C, Freyhof J. Diel distribution patterns of fishes in a temperate large lowland river. Journal of Fish Biology, 2004, 64(3): 632- 642.

[18] Knudsen F R, Sgrov H. Benefits from horizontal beaming during acoustic survey: application to three Norwegian lakes. Fisheries Research, 2002, 56(2): 205- 211.

[19] 孙明波, 谷孝鸿, 曾庆飞, 王银平, 毛志刚, 谷先坤. 基于水声学方法的天目湖鱼类资源捕捞与放流的生态监测. 生态学报, 2013, 33(23): 7553- 7562.

[21] Aglen A. Random Errors of Acoustic Fish Abundance Estimates in Relation to the Survey Grid Density Applied. Norway: Symposium on Fisheries Acoustics, 1983: 293- 298.

[22] Petitgas P. Geostatistics for fish stock assessments: a review and an acoustic application. ICES Journal of Marine Science, 1993, 50(3): 285- 298.

[23] Li X, Cheng G D, Lu L. Comparison of spatial interpolation methods. Advance in Earth Sciences, 2000, 15(3): 260- 265.

[24] 孙儒泳. 动物生态学原理 (第三版). 北京: 北京师范大学出版社, 1987: 283- 295.

[25] 李林春. 鱼类养殖生物学. 北京: 中国农业科学技术出版社, 2007: 251- 252.

[26] 倪勇, 朱成德. 太湖鱼类志. 上海: 上海科学技术出版社, 2005: 134- 136.

[27] 高永霞, 朱广伟, 贺冉冉, 王芳. 天目湖水质演变及富营养化状况研究. 环境科学, 2009, 30(3): 673- 679.

[29] 张运林, 陈伟民, 杨顶田, 黄文钰, 江晶. 天目湖热力学状况的监测与分析. 水科学进展, 2004, 15(1): 61- 67.

[30] 谭细畅, 李新辉, 林建志, 周东华, 高翔, 李嘉敏. 基于水声学探测的两个广东鲂产卵群体繁殖生态的差异性. 生态学报, 2009, 29(4): 1756- 1762.

[31] Churnside J H, Demer D A, Mahmoudi B. A comparison of lidar and echosounder measurements of fish schools in the Gulf of Mexico. ICES Journal of Marine Science, 2003, 60(1): 147- 154.

Seasonal and daily spatial distribution of fish in Lake Tianmu based on the hydroacoustic method

SUN Mingbo1,2, GU Xiaohong1,*, ZENG Qingfei1, WANG Yinping1,2, MAO Zhigang1, GU Xiankun1,2

1StateKeyLaboratoryofLakeScienceandEnvironment,NanjingInstituteofGeography&Limnology,ChineseAcademyofSciences,Nanjing210008,China2UniversityofChineseAcademyofSciences,Beijing100049,China

The seasonal and diurnal changes of horizontal and vertical spatial distribution of fish were studied by applying the hydroacoustic method in Lake Tianmu. The results showed that there was regular horizontal migration of fish between different seasons and no significant difference in horizontal distribution between day and night in Lake Tianmu. However, there existed a slight diurnal change in the inshore—offshore horizontal migration in some parts of the lake area during the late spring and in the summer. Significant differences in vertical distribution were found between seasons and between day and night. The vertical distribution of fish density was significantly affected by the water temperature stratification in July when there was a thermocline, and the fish density in hypoliminion was almost zero. The fish in Lake Tianmu was distributed in groups in the four seasons, but the typical cluster form occurred only in January. Relative to the daytime, the fish distribution at night was more uniform and scattered. Moreover, horizontal survey showed that a large number of fish was spread in the upper water layer at night.

hydroacoustic; fish; spatial distribution; seasons; day and night

国家自然科学基金(31270506); 国家“十二五”科技计划(2012BAD25B06/07); 江苏省自然科学基金(BK20131059)

2014- 01- 01;

日期:2014- 11- 03

10.5846/stxb201401010001

*通讯作者Corresponding author.E-mail: xhgu@niglas.ac.cn

孙明波, 谷孝鸿, 曾庆飞, 王银平, 毛志刚, 谷先坤.基于水声学方法的天目湖鱼类季节和昼夜空间分布研究.生态学报,2015,35(17):5597- 5605.

Sun M B, Gu X H, Zeng Q F, Wang Y P, Mao Z G, Gu X K.Seasonal and daily spatial distribution of fish in Lake Tianmu based on the hydroacoustic method.Acta Ecologica Sinica,2015,35(17):5597- 5605.

猜你喜欢
天目湖鱼类季节
基于MFCC和ResNet的鱼类行为识别
天目湖日月潭演绎的佳话
鱼类运动会
我喜欢的季节7
季节蠕变
季节的变换
TЫCяЧEЛETний КУЛьTУрнЫй гороДоК
花的季节
天目湖砂锅鱼头