Geoscience Reference
In-Depth Information
We dei ne a new variable
newlaminae
that contains the vertical
y
-component
of
laminae
as the i rst column and the years 1 to 22 (counting backwards
in time) as the second colum. h e 22 years are equivalent to the length
of
laminae
. h e function
interp1
is used to interpolate the color-intensity
transects over an evenly-spaced time axis stored in the variable
age
.
newlaminae(:,1) = laminae(:,2);
newlaminae(:,2) = 1 : length(laminae(:,2));
age = interp1(newlaminae(:,1),newlaminae(:,2),CY);
We complete the analysis by plotting the color-intensity curves on both a
length and a time scale for comparison (Fig. 8.8).
subplot(2,1,1), plot(CY,C(:,1),CY,C(:,2),CY,C(:,3))
xlabel('Centimeters'), ylabel('Intensity'), title('Color vs. Length')
subplot(2,1,2), plot(age,C(:,1),age,C(:,2),age,C(:,3))
xlabel('Years'), ylabel('Intensity'), title('Color vs. Age')
h e interpolated color-intensity transects can now be further analyzed
using the time-series analysis tools. h e analysis of a representative red-
color intensity transect across 70-250 varves during the project described
a
b
Fig. 8.8
Color-intensity curves (red, green and blue) plotted against
a
depth and
b
age.