Auxiliary Material for Paper 2008GC002010

Depth-shifting cores incompletely recovered from the upper oceanic crust, IODP Hole 1256D

Lisa A. Gilbert
Maritime Studies Program, Williams College and Mystic Seaport, 75 Greenmanville Avenue, Mystic, Connecticut 06355, USA (lgilbert@williams.edu)

Andrea Burke
MIT/WHOI Joint Program, Woods Hole Oceanographic Institution, MS 22, Woods Hole, Massachusetts 02543, USA (aburke@mit.edu)

Gilbert, L. A., and A. Burke (2008), Depth-shifting cores incompletely recovered from the upper oceanic crust, IODP Hole 1256D, Geochem. Geophys. Geosyst., 9, Q08O11
, doi:10.1029/2008GC002010.

Introduction

The auxiliary material includes data used and MATLAB scripts in the main article to depth-shift core data to match log data.  

The data are from ODP Leg 206 and IODP Expeditions 309 and 312, and are publicly available at http://www-odp.tamu.edu/database/.  The only exception is 2008gc002010-txts10.txt, which is the newly shifted data set described in the main article.

For use in MATLAB, the files should be renamed as follows and then loaded into a single directory:
2008gc002010-txts01.txt ->  GRA.dat
2008gc002010-txts02.txt ->  LOG.dat
2008gc002010-txts03.txt ->  PIECES.dat
2008gc002010-txts04.txt ->  MAD.dat
2008gc002010-txts05.txt ->  CORELOGINTEGRATE.m
2008gc002010-txts06.txt ->  DATADJUST.m
2008gc002010-txts07.txt ->  EXAMPLE.m
2008gc002010-txts08.txt ->  FINDNEAREST.m
2008gc002010-txts09.txt ->  MATCHPERCENT.m
2008gc002010-txts10.txt ->  dataset.txt

Data files

Data used by Gilbert and Burke (2008) are given in GRA.txt, log.txt,
pieces.txt, and MAD.txt data files from ODP/IODP core data, in formats similar to Janus: http://www-odp.tamu.edu/database/.  

1. 2008gc002010-txts01.txt -> GRA.dat is a 9-column matrix of GRA density data.
    GRA(:,1) is core number
    GRA(:,2) is section number 
    GRA(:,3) is depth to top of core (mbsf) 
    GRA(:,4) is depth from top of section (cm)
    GRA(:,5) is piece number 
    GRA(:,6) is curated depth (mbsf) 
    GRA(:,7) is density (g/cc) 
    GRA(:,8) is length of core recovered from this section (m) 
    GRA(:,9) is curated bottom of cored interval (mbsf)
  
2. 2008gc002010-txts02.txt -> LOG.dat is a 2-column matrix of downhole logging data.
    LOG(:,1) is depth (mbsf) 
    LOG(:,2) is logging density (g/cc)  

3. 2008gc002010-txts03.txt -> PIECES.dat is a 7-column matrix.
    PIECES(:,1) is core number 
    PIECES(:,2) is section number 
    PIECES(:,3) is piece number 
    PIECES(:,4) is distance below top of section to top of piece (cm) 
    PIECES(:,5) is distance below top of section to bottom of piece (cm) 
    PIECES(:,6) is length of piece (cm) 
    PIECES(:,7) is tide-corrected depth to top of core (mbsf)

4. 2008gc002010-txts04.txt -> DATA.dat is used as DATA in DATADJUST and is a 6-column matrix.
    DATA(:,1) is core number
    DATA(:,2) is section number 
    DATA(:,3) is curated depth (mbsf) 
    DATA(:,4) is physical property data (e.g. density, velocity, porosity)
    DATA(:,5) is piece number 
    DATA(:,6) is the depth of data from top of piece (cm)

MATLAB files

For use in MATLab, following files need to have their extensions 
returned to .m for use in MATLAB: CORELOGINTEGRATE, DATAJUST, FINDNEAREST EXAMPLE, and MATCHPERCENT. 

5. 2008gc002010-txts05.txt -> CORELOGINTEGRATE.m  
% [NEWPIECEDEPTHS, BEFORE, AFTER] = corelogintegrate(GRA, LOG, PIECES, PERCENT)
% 
% This function takes GRA density data, compares it to log density data and
% readjusts the depths of the pieces of the core to minimize the density
% differences. The input variables are: GRA (piece density data), LOG
% (logging density data), PIECES (curated depths of the pieces), and
% PERCENT (decimal percent limit of density to restrict possible logging
% depths).
% -------------------------------------------------------------------------
% INPUTS: GRA, LOG, PIECES
%  GRA is a 9-column matrix of GRA density data from ODP/IODP core data, in
%  a format similar to Janus, http://www-odp.tamu.edu/database/.
%    GRA(:,1) is core number
%    GRA(:,2) is section number 
%    GRA(:,3) is depth to top of core (mbsf) 
%    GRA(:,4) is depth from top of section (cm)
%    GRA(:,5) is piece number 
%    GRA(:,6) is curated depth (mbsf) 
%    GRA(:,7) is density (g/cc) 
%    GRA(:,8) is length of core recovered from this section (m) 
%    GRA(:,9) is curated bottom of cored interval (mbsf)
%  LOG is a 2-column matrix of downhole logging data.
%    LOG(:,1) is depth (mbsf) LOG(:,2) is logging density (g/cc)  
%  PIECES is a 7-column matrix. 
%    PIECES(:,1) is core number 
%    PIECES(:,2) is section number 
%    PIECES(:,3) is piece number 
%    PIECES(:,4) is distance below top of section to top of piece (cm) 
%    PIECES(:,5) is distance below top of section to bottom of piece (cm) 
%    PIECES(:,6) is length of piece (cm) 
%    PIECES(:,7) is tide-corrected depth to top of core (mbsf)
%  PERCENT is user choice of percent match required between core and log
%  data, as a decimal (use 0.2 for 20% match), from the maximum GRA value
%  of a given piece.  Value used in Gilbert and Burke, 2008 is 0.2.
%
% OUTPUT: NEWPIECEDEPTHS,BEFORE,AFTER
%   NEWPIECEDEPTHS is a 4-column matrix. 
%    NEWPIECEDEPTHS(:,1) is core number
%    NEWPIECEDEPTHS(:,2) is section number 
%    NEWPIECEDEPTHS(:,3) is piece number 
%    NEWPIECEDEPTHS(:,4) is new depth to the top of the piece (mbsf)
%   BEFORE is a 4-column matrix representing the data before depth shifting
%    BEFORE(:,1) is curated depth (mbsf) 
%    BEFORE(:,2) is the density value for the core piece
%    BEFORE(:,3) is the nearest logging depth to the curated depth (mbsf)
%    BEFORE(:,4) is the density value at that nearest depth in the log record
%   AFTER is a 4-column matrix representing the data after depth shifting
%    AFTER(:,1) is new depth after shifting of physical property (mbsf) 
%    AFTER(:,2) is the density value for the core piece
%    AFTER(:,3) is the nearest logging depth to the new depth (mbsf)
%    AFTER(:,4) is the density value at that new nearest depth in the log record
%
% ------------------------------------------------------------------------
% *Variations*
%   Other data types (Vp, porosity, etc.) can be substituted as follows:
%   GRA(:,7) can be replaced with data obtained from the recovered core if
%   LOG(:,2) is similarly replaced with the same type of data from the log. 
%   All information in the remaining columns can remain unchanged.
% -------------------------------------------------------------------------

6. 2008gc002010-txts06.txt -> DATADJUST 
% [B,A] = datadjust(DATA, NEWPIECEDEPTHS, LOG)
%
% This script will integrate any physical property data with
% the new piece depths.  
%
% -------------------------------------------------------------------------
% INPUTS: DATA, NEWPIECEDEPTHS, LOG
%  DATA is a 6-column matrix of physical property data from ODP/IODP core data, in
%  a format similar to Janus, http://www-odp.tamu.edu/database/.
%    DATA(:,1) is core number
%    DATA(:,2) is section number 
%    DATA(:,3) is curated depth (mbsf) 
%    DATA(:,4) is physical property data (e.g. density, velocity, porosity)
%    DATA(:,5) is piece number 
%    DATA(:,6) is the depth of data from top of piece (cm)
%
%   NEWPIECEDEPTHS is a 4-column matrix computed by corelogintegrate.m
%    NEWPIECEDEPTHS(:,1) is core number
%    NEWPIECEDEPTHS(:,2) is section number 
%    NEWPIECEDEPTHS(:,3) is piece number 
%    NEWPIECEDEPTHS(:,4) is new depth to the top of the piece (mbsf)
%    
%  LOG is a 2-column matrix of downhole logging data of the physical
%  property found in DATA(:,4)
%    LOG(:,1) is depth (mbsf) 
%    LOG(:,2) is physical property data (see DATA(:,4)) 
%
%  PIECES is a 7-column matrix. 
%    PIECES(:,1) is core number 
%    PIECES(:,2) is section number 
%    PIECES(:,3) is piece number 
%    PIECES(:,4) is distance below top of section to top of piece (cm) 
%    PIECES(:,5) is distance below top of section to bottom of piece (cm) 
%    PIECES(:,6) is length of piece (cm) 
%    PIECES(:,7) is tide-corrected depth to top of core (mbsf)
%
% OUTPUT: B, A
%  B is a 4-column matrix representing the data before depth shifting
%    B(:,1) is curated depth (mbsf) 
%    B(:,2) is physical property data for the core piece
%    B(:,3) is the nearest logging depth to the curated depth (mbsf)
%    B(:,4) is the physical property value at that nearest depth in the log record
%  A is a 4-column matrix representing the data after depth shifting
%    A(:,1) is new depth after shifting of physical property (mbsf) 
%    A(:,2) is physical property data for the core piece
%    A(:,3) is the nearest logging depth to the new depth (mbsf)
%    A(:,4) is the physical property value at that nearest depth in the log record
%
% -------------------------------------------------------------------------

7. 2008gc002010-txts07.txt -> EXAMPLE.m 
% This script runs CORELOGINTEGRATE.m and DATADJUST.m with the data used
% in Gilbert and Burke (2008) as an example.

8. 2008gc002010-txts08.txt ->  FINDNEAREST.m 
% [N] = findnearest(DEPTH, LOGDEPTHS)
%
% This function finds the depth in the logging data (LOGDEPTHS) that is
% closest to the given depth (DEPTH).
%
% -------------------------------------------------------------------------

9. 2008gc002010-txts09.txt -> MATCHPERCENT.m 
% [M] = matchpercent(VALUE, VECTOR, PERCENT)
%
% MATCHPERCENT is called by CORELOGINTEGRATE.
% It takes a single density VALUE (1x6 vector with the following
% columns: 1 depth, 2 density, 3 core, 4 section, 5, piece, 6 depth from
% top of section) from GRA data and compares it to a VECTOR (2 columns: 1
% depth, 2 density) of possible depths and corresponding log density to
% find all of the depths that are within a certain PERCENT (as a decimal)
% of the value.
% [M] is a 8 column matrix: 1 log depth, 2 log density, 3 curated depth,
% 4 GRA density, 5 core, 6 section, 7 piece, 8 distance from top of section
% -------------------------------------------------------------------------

Revised depth dataset

10. 2008gc002010-txts10.txt -> dataset.txt 
Column 1: core 
Column 2: section
Column 3: piece
Column 4: new depth (mbsf) 

After depth-shifting, as described in the manuscript, we this new depth vector results.