MAST-467/667: Introduction to Arctic Oceanography (2021)

NASA's OMG stations 2020

Date: 10/10/2021

Author: Frederike Benz

close all

clear all

%load NASAs OMG data 2020

load output.mat % [NASA (2021)]

% variables

lon = output(:,1); % Longitude

lat = output(:,2); % Latitude

depth = output(:,3);

temp = output(:,4);

sal = output(:,5);

% Disko bay area

ind_disko = find(lat >= 68 & lat <= 72 & lon <= 90 & lon >= 50);

% Variables Disko bay

lon_disko = lon(ind_disko);

lat_disko = lat(ind_disko);

depth_disko = depth(ind_disko);

temp_disko = temp(ind_disko);

sal_disko = sal(ind_disko);

%% MAP Greenland with 2020 OMG stations (red) plus highlighted Disko Bay stations (blue)

clf, axis([-80 -5 55 90])

topoplot(axis,-5:1/8:0),hold on % [Lilly (2021)]

topoplot continents, latratio(55) % [Lilly (2021)]

a = plot(lon.*(-1),lat,'wo','markerfacecolor','r');

b = plot(lon_disko.*(-1),lat_disko,'wo','markerfacecolor','b');

hc=colorbar;hc.Label.String='Bathymetry [km]';

title('NASAs OMG Data 2020 (stations)');

ylabel('Latitude')

xlabel('Longitude')

legend([a b],'Stations','Disko Bay','location','southeast')

set(gca,'fontsize',15,'linewidth',2)

%% Zoom map Disko Bay

clf, axis([-65 -50 66 73])

topoplot(axis,-5:1/8:0) % [Lilly (2021)]

hold on

topoplot continents, latratio(55) % [Lilly (2021)]

b = plot(lon.*(-1),lat,'wo','markerfacecolor','r')

b = 
  Line with properties:

              Color: [1 1 1]
          LineStyle: 'none'
          LineWidth: 0.5000
             Marker: 'o'
         MarkerSize: 6
    MarkerFaceColor: [1 0 0]
              XData: [1×1488442 double]
              YData: [1×1488442 double]
              ZData: [1×0 double]

  Show all properties

c = plot(lon_disko.*(-1),lat_disko,'wo','markerfacecolor','b');

hc=colorbar;hc.Label.String='Bathymetry [km]';

title('Disko Bay - NASAs OMG Data 2020 (stations)');

ylabel('Latitude');

xlabel('Longitude');

set(gca,'fontsize',20,'linewidth',2);

% How many stations airborned XCTD profiles to we have in Disko Bay?

% [lat_disko_station,lon_disko_station,num] = hista(lat_disko,lon_disko,0.1);

% number_stations_disko = length(num)

lat_disko_station = unique(lat_disko);

lon_disko_station = unique(lon_disko);

% The assigned area "Disko Bay" counts 39 stations in 2020.

ind_station1 = find(lon_disko == lon_disko_station(2,1));

% 2. Plot a single profile, that is, temperature T = T(z) and salinity S = S(z)?

lat_profile1 = lat_disko(ind_station1);

lon_profile1 = lon_disko(ind_station1);

depth_profile1 = depth_disko(ind_station1);

temp_profile1 = temp_disko(ind_station1);

sal_profile1 = sal_disko(ind_station1);

clf, d = plot(temp_profile1,depth_profile1,'.','linewidth',2)

d = 
  Line with properties:

              Color: [0 0.4470 0.7410]
          LineStyle: 'none'
          LineWidth: 2
             Marker: '.'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×9520 double]
              YData: [1×9520 double]
              ZData: [1×0 double]

  Show all properties

ylabel('Depth [m]')

xlabel('Temperature [degree C]')

ylim([0 100])

xlim([-1 3])

set(gca,'ydir','reverse')

set(gca,'fontsize',15,'linewidth',2)

clf, d = plot(sal_profile1,depth_profile1,'.','linewidth',2)

d = 
  Line with properties:

              Color: [0 0.4470 0.7410]
          LineStyle: 'none'
          LineWidth: 2
             Marker: '.'
         MarkerSize: 6
    MarkerFaceColor: 'none'
              XData: [1×9520 double]
              YData: [1×9520 double]
              ZData: [1×0 double]

  Show all properties

ylabel('Depth [m]')

xlabel('Salinity')

ylim([0 100])

xlim([20 35])

set(gca,'ydir','reverse')

set(gca,'fontsize',15,'linewidth',2)

References:

NASA (2021), Oceans Melting Greenland, https://omg.jpl.nasa.gov/portal/browse/products/urn:omg:OMG_Ocean_AXCTD_L1

Lilly, J. M. (2021), jLab: A data analysis package for Matlab, v.1.7.1, doi:10.5281/zenodo.4547006, http://www.jmlilly.net/software.