利用MODIS卫星遥感数据反演中国东海叶绿素浓度及SST年际变化的研究报告.docx
《利用MODIS卫星遥感数据反演中国东海叶绿素浓度及SST年际变化的研究报告.docx》由会员分享，可在线阅读，更多相关《利用MODIS卫星遥感数据反演中国东海叶绿素浓度及SST年际变化的研究报告.docx（19页珍藏版）》请在冰豆网上搜索。
利用MODIS卫星遥感数据反演中国东海叶绿素浓度及SST年际变化的研究报告
利用MODIS卫星遥感数据反演中国东海叶绿素浓度及SST年际变化的研究
10海洋科学王雪
摘要：
利用MODIS获取的2003年-2006年卫星数据反演中国东海海域的海表温度〔SST〕和叶绿素a浓度信息，并对其分布规律进展分析。
结果说明：
东海海域的SST和叶绿素a浓度的分布具有明显的分区和季节变化特征。
在年际间，不同海区的SST和叶绿素a浓度均呈现出周期性变化趋势。
在空间上，SST呈现出由近岸向外海递增趋势，南北位差大；叶绿素a浓度那么呈现出由近岸向外海递减的分布趋势。
东海叶绿素a浓度的分布与SST、河口径流、季节等因素有关。
关键词：
东海MODISSST叶绿素a
引言
卫星遥感数据具有周期短、时空分布率高、数据具有可比性，以及卫星遥感能够实现对地球大面积同步观测、卫星遥感平台具有专业化多样化等特点，这使得遥感数据成为全球海洋监测的一个重要数据源，因此遥感技术成为了全球海洋监测的不可替代的技术手段。
在诸多海洋遥感平台所搭载的传感器中，中成像光谱仪〔MODIS〕将时间分布率、空间分布率及光谱分辨率很好地予以结合，每1~2d就能获取分辨率为250~1000m、包括36个波段的全球尺度MODIS数据，其中第8到16共9个波段被广泛应用于海洋遥感。
利用星载或机载遥感器探测海水中的叶绿素光谱辐射，经大气校正后对海水表层叶绿素浓度进展反演，即根据叶绿素的光学特性求的海水中叶绿素浓度的一种方法，该方法使同一时间对大围海域进展叶绿素监测成为可能，因此遥感技术被越来越广泛地应用于海水叶绿素监测。
中国东海是典型的陆架海域,长江等河流入海,使之富集叶绿素、悬移质和黄色物质等影响海洋光学性质的海水组分。
卫星水色遥感是唯一可穿透海水一定深度并获取海洋上层光学参数以及海水组分的技术手段。
海洋叶绿素a〔Chl-a〕是了解海洋中地球生物化学循环的根底和估算海洋生产力的根本指标，也是判断水域的肥瘠程度和评价水域渔业潜在生产力的根本依据。
富集浮游植物的海域是海洋食植动物大密度存在和水产资源丰富存在的根底。
海水叶绿素a浓度的测定不仅与海洋生态系统初级生产力的研究密切相关，而且对于海洋-大气系统中碳循环、环境监测、赤潮灾害监测、海流（上升流、沿岸流）等的研究及渔业管理等都具有重要意义。
海表温度〔SeaSurfaceTemperature，SST〕是影响海洋渔业资源的重要因素之一。
海表温度不仅对海洋鱼类的繁殖、生长等有影响，而且对鱼类的洄游影响也比拟大，从而直接影响渔期的早晚和长短、中心渔场的位置变动以及鱼类的集群程度。
对海表温度的时空分布和变化的研究，是海洋渔业研究的根本容之一。
本研究主要利用MODIS遥感数据，反演东海海区长时间序列的SST和叶绿素a浓度信息，分析研究东海海区SST和表层叶绿素a浓度的分布及时空变化特征，旨在为东海的海洋动力学特征、海洋生产力和海洋生态环境研究提供参考。
1数据来源与处理
1.1数据来源
数据来自ＮＡＳＡ的海洋水色
seadas.gsfc.nasa.gov/
1）在水色首页点击DataArchive→MODIST→L3SMI,选择下载2003年至2006年分辨率为4km的年叶绿素a浓度和年SST数据；
2）在水色首页点击Level3Browser→TerraMODISChlorophyllconcentration〔TerraMODISSeaSurfaceTemperature11μdaytime〕→Seasonalposite，选择下载分辨率为4km的2003年至2006年季叶绿素a浓度〔季SST〕数据；
3）在水色首页点击Level3Browser→TerraMODISChlorophyllconcentration〔TerraMODISSeaSurfaceTemperature11μdaytime〕→Monthlyclimatology,选择下载分辨率为4km的2003年至2006年月叶绿素a浓度〔季SST〕数据。
1.2数据处理
1.2.1叶绿素数据处理
以2003年为例：
clear,clc;
LATLIMS=[2336];
LONLIMS=[117130];
PI=hdfinfo（'T365.L3m_YR_CHL_chlor_a_4km'）;
pin=[];
fork=1:
60,
nm=PI.Attributes（k）.Name;nm（nm==''）='_';
ifisstr（PI.Attributes（k）.Value）,
pin=setfield（pin,nm,PI.Attributes（k）.Value）;
else
pin=setfield（pin,nm,double（PI.Attributes（k）.Value））;
end
end;
lon=[pin.Westernmost_Longitude:
pin.Longitude_Step:
pin.Easternmost_Longitude];
lat=[pin.Northernmost_Latitude:
-pin.Latitude_Step:
pin.Southernmost_Latitude];
[mn,ilt]=min（abs（lat-max（LATLIMS）））;
[mn,ilg]=min（abs（lon-min（LONLIMS）））;
ltlm=fix（diff（LATLIMS）/pin.Latitude_Step）;
lglm=fix（diff（LONLIMS）/pin.Longitude_Step）;
P=hdfread（'T365.L3m_YR_CHL_chlor_a_4km','l3m_data','Index',{[iltilg],[],[ltlmlglm]}）;
P=double（P）;
P（P==-32767）=NaN;
P=（pin.Slope*P+pin.Intercept）;
p=log10（P）;
LT=lat（ilt+[0:
ltlm-1]）;LG=lon（ilg+[0:
lglm-1]）;
[Plg,Plt]=meshgrid（LG,LT）;
m_proj（'lambert','lon',LONLIMS,'lat',LATLIMS）;
m_pcolor（Plg,Plt,p）;shadingflat;
m_gshhs_i（'color','k'）;
m_grid（'linewi',2,'tickdir','out'）;
h=colorbar;
set（get（h,'ylabel'）,'String','Chla（mgm^-3）'）;
set（h,'ytick',log10（[0.10.51351020]）,'yticklabel',[0.10.51351020],'tickdir','out'）;
title（['MODISChla'datestr（datenum（pin.Period_Start_Year,1,0）+pin.Period_Start_Day）'->'...
datestr（datenum（pin.Period_Start_Year,1,0）+pin.Period_End_Day）],...
'fontsize',13,'fontweight','bold'）;
print（gcf,'-dpng','p03'）
1.2.1SST数据处理
clear,clc;
LATLIMS=[2336];
LONLIMS=[117130];
PI=hdfinfo（'T365.L3m_YR_SST_4'）;
pin=[];
fork=1:
60,
nm=PI.Attributes（k）.Name;nm（nm==''）='_';
ifisstr（PI.Attributes（k）.Value）,
pin=setfield（pin,nm,PI.Attributes（k）.Value）;
else
pin=setfield（pin,nm,double（PI.Attributes（k）.Value））;
end
end;
lon=[pin.Westernmost_Longitude:
pin.Longitude_Step:
pin.Easternmost_Longitude];
lat=[pin.Northernmost_Latitude:
-pin.Latitude_Step:
pin.Southernmost_Latitude];
[mn,ilt]=min（abs（lat-max（LATLIMS）））;
[mn,ilg]=min（abs（lon-min（LONLIMS）））;
ltlm=fix（diff（LATLIMS）/pin.Latitude_Step）;
lglm=fix（diff（LONLIMS）/pin.Longitude_Step）;
P=hdfread（'T365.L3m_YR_SST_4','l3m_data','Index',{[iltilg],[],[ltlmlglm]}）;
P=double（P）;
P（P==65535）=NaN;
P=（pin.Slope*P+pin.Intercept）;
LT=lat（ilt+[0:
ltlm-1]）;LG=lon（ilg+[0:
lglm-1]）;
[Plg,Plt]=meshgrid（LG,LT）;
m_proj（'lambert','lon',LONLIMS,'lat',LATLIMS）;
m_pcolor（Plg,Plt,P）;shadingflat;
m_gshhs_i（'color','k'）;
m_grid（'linewi',2,'tickdir','out'）;
h=colorbar;
set（get（h,'ylabel'）,'String','SST（deg-C）'）;
set（h,'ytick',[5101520253035],'yticklabel',[5101520253035],'tickdir','out'）;
title（['MODISSST'datestr（datenum（pin.Period_Start_Year,1,0）+pin.Period_Start_Day）'->'...
datestr（datenum（pin.Period_Start_Year,1,0）+pin.Period_End_Day）],...
'fontsize',13,'fontweight','bold'）;
print（gcf,'-dpng','p4'）
2图形绘制
图1东海叶绿素a浓度的季节分布图〔mg/m3〕
Fig1Seasonalvariationinchloropgyll-aconcentrationinEastChinaSeafrom2003to2006
图2东海SST季节平均分布图〔℃〕
Fig2SeasonalSSTdistributionsinEastChinaSea
图3东海叶绿素a浓度年际变化〔mg/m3〕
Fig3Annualvariationinchloropgyll-aconcentrationinEastChinaSeafrom2003to2006
图4东海SST年际变化（℃）
Fig4AnnualvariationinSSTinEastChinaSea
升级会员 