%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% clc clear all load Figure1_data.mat %%%% tide variables:ptide ntide lf-tide; AMOC variable:amoc; global reanalysis:gb_ssha; downscaled reanalysis:ds_ssha %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% figure(1) %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%10 10 h1=subplot('position',[0.1 0.6 0.8 0.3]) set(h1,'position',[0.1 0.6 0.8 0.3]) a1=area(1912:1:2023,ptide); set(a1,'facecolor','r','edgecolor','r') hold on a2=area(1912:1:2023,ntide); set(a2,'facecolor','b','edgecolor','b') xlim([1910 2024]) set(gca,'linewidth',1,'fontsize',8,'Ytick',[-80:40:80],'ylim',[-88 88],'xtick',[1910:10:2021],'xlim',[1910 2021],'tickdir','out','YMinorTick','off','XMinorTick','off') box on hold on plot(1912:2023,lf_tide,'k','linewidth',2) hold on [AX,H1,H2]=plotyy(1961:2022,gb_ssha*1000,1961:2020,amoc(4:end)-mean(amoc(4:end))) set(H1,'Color','m','linewidth',2,'Marker','.','Markersize',18) set(H2,'color',[255 150 0]/255,'linewidth',2) set(AX(1),'YColor','k','linewidth',1,'fontsize',11,'Ytick',[-100:50:100],'ylim',[-100 100],'xtick',[1910:10:2020],'xlim',[1910 2024],'tickdir','out') set(AX(2),'YColor',[255 150 0]/255,'linewidth',1,'fontsize',11,'Ytick',[-4:2:4],'ylim',[-4 4],'xtick',[1910:10:2020],'xlim',[1910 2024],'tickdir','out') hold on plot(1961:2022,ds_ssha*1000,'color',[0.15 0.66 0.47],'linewidth',2,'Marker','.','Markersize',18) hold on lin(1:62)=-300; plot(1961:2022,lin,'color',[255 150 0]/255,'linewidth',2) clear lin hold on text(1930,114,'U.S. Northeast Coast sea level timeseries Versus AMOC index','fontsize',12) text(1902,114,'A','fontsize',12,'fontweight','bold') l1=legend('TG','TG','low pass filtered TG','global reanalysis','downscaling reanalysis','AMOC(reanalysis)') set(l1,'edgecolor','w','orientation','horizontal') hold on text(1912,-61,'Cor(TG, global reanalysis): 0.72','fontsize',11,'color',[0.15 0.66 0.47]) hold on text(1912,-70,'Cor(TG, downscaling reanalysis): 0.65','fontsize',11,'color',[0.15 0.66 0.47]) hold on text(1912,-79,'Cor(low pass filtered TG obs, AMOC): -0.73','fontsize',11,'color',[0.15 0.66 0.47]) xlabel('Year','fontsize',11) ylabel('sea level (mm)','fontsize',11) text(2031,-20,'AMOC (Sv)','fontsize',11,'rotation',90,'color',[255 150 0]/255) % [R,P]=corrcoef(tide(50:111),gb_ssha(1:62)); % [R,P]=corrcoef(tide(50:111),ds_ssha(1:62)); % [R,P]=corrcoef(lf_tide(50:109),amoc(4:63)); %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% h2=subplot('position',[0.1 0.3 0.22 0.22]) set(h2,'position',[0.1 0.3 0.22 0.22]) bb(1,:)=[102 102 255]/255; bb(2,:)=[0 0 255]/255; bb(3,:)=[0 64 255]/255; bb(4,:)=[0 128 255]/255; bb(5,:)=[0 204 255]/255; bb(6,:)=[0 255 255]/255; bb(7,:)=[200 255 255]/255; bb(8,:)=[240 255 255]/255; bb(9,:)=[255 240 240]/255; bb(10,:)=[255 210 210]/255; bb(11,:)=[255 166 166]/255; bb(12,:)=[255 111 111]/255; bb(13,:)=[255 55 55]/255; bb(14,:)=[255 0 0]/255; bb(15,:)=[190 0 0]/255; bb(16,:)=[255 51 153]/255; colormap(bb) m_proj('miller','lon',[260 360],'lat',[20 68]); hold on m_coast('color',[0 0 0]); hold on m_contourf(sea_level_lon,sea_level_lat,sea_level_reg,[-40:5:40],'linestyle','none') m_grid('tickdir','out','linestyle','none','fontsize',10,'linewidth',1,'xtick',[260:30:360],'ytick',[20:20:90],'Xminortick','on','Yminortick','on'); % set(findobj('tag','m_grid_color'),'facecolor','none') caxis([-40 40]) hold on m_coast('patch',[0 0 0]); hold on m_text(264,70.9,'Sea level during -1*AMOC','fontsize',11,'color','k') m_text(236,70.9,'B','fontsize',12,'fontweight','bold') %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% h3=subplot('position',[0.38 0.3 0.22 0.22]) set(h3,'position',[0.38 0.3 0.22 0.22]) m_proj('miller','lon',[260 360],'lat',[20 68]); hold on m_coast('color',[0 0 0]); hold on colormap(bb) m_contourf(sea_level_lon,sea_level_lat,sea_level_steric,[-40:5:40],'linestyle','none') m_grid('tickdir','out','linestyle','none','fontsize',10,'linewidth',1,'xtick',[260:30:360],'ytick',[20:20:90],'Xminortick','on','Yminortick','on'); % set(findobj('tag','m_grid_color'),'facecolor','none') caxis([-40 40]) hold on m_coast('patch',[0 0 0]); hold on m_text(275,70.9,'Steric component','fontsize',11,'color','k') m_text(242,70.9,'C','fontsize',12,'fontweight','bold') %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% h4=subplot('position',[0.66 0.3 0.22 0.22]) set(h4,'position',[0.66 0.3 0.22 0.22]) m_proj('miller','lon',[260 360],'lat',[20 68]); hold on m_coast('color',[0 0 0]); hold on colormap(bb) m_contourf(sea_level_lon,sea_level_lat,sea_level_mass,[-40:5:40],'linestyle','none') m_grid('tickdir','out','linestyle','none','fontsize',10,'linewidth',1,'xtick',[260:30:360],'ytick',[20:20:90],'Xminortick','on','Yminortick','on'); % set(findobj('tag','m_grid_color'),'facecolor','none') caxis([-40 40]) hold on m_coast('patch',[0 0 0]); hold on m_text(275,70.9,'Mass redistribution','fontsize',11,'color','k') m_text(242,70.9,'D','fontsize',12,'fontweight','bold') ch1=colorbar('horizontal') set(ch1,'xtick',[-40:5:40],'xticklabel',{'-40','','-30','','-20','','-10','','0','','10','','20','','30','','40'},'position',[0.23 0.29 0.5 0.013],'fontsize',10,'TickLength',[0.026]) clear ans bb print('-depsc2', 'Figure1_in_article.eps')