Statistics and Visualizations of Profiles in Nares Strait (2016-2021)

Tianyu Zhou, 10/12/2021, at UDel

1. Statistics of maximum sampling depth

Load data located within Nares Strait (pre-selected and stored in a '.mat' file)

clc;
clear;
load('../data/OMG_data_Nares')
profi_num = size(output,1);
year = ones(1,profi_num);
depth_max = ones(1,profi_num);
for n = 1:profi_num
    year(n) = output{n,1};
    data = output{n,5};
    if size(data,2) == 8
        depth_max(n) = max(data(:,2));
    else
        depth_max(n) = max(data(:,3));
    end
end

Make a histogram showing the distribution of maximum sampling depth

f = figure;

set(f,'units','centimeters','position',[1,1,15.5*1.5,9*1.5],'color','w')

FontSize = 20;

histogram(depth_max,12)

ylim([0 17])

text(700,15,['Total profiles: ',num2str(profi_num)],'color','r','FontSize',FontSize)

xlabel('Maximum Sampling Depth (m)','FontWeight','bold')

ylabel('Counts','FontWeight','bold')

grid on

set(gca,'FontSize',FontSize,'tickdir','out')

2. Statistics of sampling year

Make a histogram showing the distributions of years

f = figure;

set(f,'units','centimeters','position',[1,1,15.5*1.5,9*1.5],'color','w')

year_stat = ones(1,6);

for n = 1:6

year_stat(n) = length(find(year==2015+n));

end

b = bar(2016:2021,year_stat);

b.FaceColor = 'flat';

b.CData = jet(6);

ylim([0 22])

text(2015,20,['Total profiles: ',num2str(profi_num)],'color','r','FontSize',FontSize)

set(gca,'FontSize',FontSize,'tickdir','out')

xlabel('Sampling Year','FontWeight','bold')

ylabel('Counts','FontWeight','bold')

grid on

3. Visualizations of all profiles (2016-2021)

Temperature profiles

f = figure;

set(f,'units','centimeters','position',[1,1,9*1.2,12*1.2],'color','w')

FontSize = 15;

hold on

for n = 1:profi_num

data = output{n,5};

if size(data,2) == 8

depth = data(:,2); % m

T = data(:,3); % oC

else

depth = data(:,3); % m

T = data(:,4); % oC

end

% delete the abnormal data and make plots

T(T == -99) = nan;

depth(T == -99) = nan;

plot(T,-depth,'.')

end

xlim([-4 6])

ylim([-1000 0])

box on

grid on

xlabel('Potential Temperature (^{o}C)','FontWeight','bold')

ylabel('Depth (m)','FontWeight','bold')

set(gca,'FontSize',FontSize,'tickdir','out','xaxislocation','top')

Salinity profiles

f = figure;

set(f,'units','centimeters','position',[1,1,9*1.2,12*1.2],'color','w')

FontSize = 15;

hold on

for n = 1:profi_num

data = output{n,5};

% NOTE: a small percentage of probes do not have salinity data

if size(data,2) == 8

depth = data(:,2); % m

S = data(:,5); % ppt

end

% delete the abnormal data and make plots

S(S == -99) = nan;

depth(S == -99) = nan;

plot(S,-depth,'.')

end

xlim([10 40])

ylim([-1000 0])

box on

grid on

xlabel('Salinity (ppt)','FontWeight','bold')

ylabel('Depth (m)','FontWeight','bold')

set(gca,'FontSize',FontSize,'tickdir','out','xaxislocation','top')