Line.Sta Algorithms

Click for larger ImageCorrected Conversion Algorithm for CalCOFI Station Grid And Their Implementation in Several Computer Languages by Edward Weber & Thomas Moorem, NOAA Fisheries, 2013 CalCOFI Report v 54.

Link to downloadable code (Java, Perl, Python & R) file presented in the CalCOFI Report vol 54, pages 102-106.

Matlab scripts by Robert Thombley, SIO-CalCOFI with Sep2014 error correction by Augusto Valencia, UABC; based on Weber & Moor 2013.

Link to MatLab downloadable code authored by Robert Thombley, SIO-CalCOFI

(Abstract from pdf) “Converting between geographic coordinates in latitude and longitude and the line and station sampling pattern of the California Cooperative Fisheries Investigations (CalCOFI) program is a commonly required task for conducting research on the California Current ecosystem. This note presents several corrections and clarifications to the previously published algorithms for performing these conversions. We include computer code to implement the algorithms in Java™1, Perl, Python, and R. We note that freely available code to conduct the conversions in Fortran, Matlab®2, JavaScript™, and Visual Basic®6 has previously been published, and an online conversion tool is also available. A future version of the PROJ.4 cartographic projections library will also include support for CalCOFI conversions, thereby allowing for convenient conversions using the GRASS GIS, PostGIS, Python, Perl, R, and many other programs and programming languages.”

For individual language – Java, Perl, Python, or R – implementations & downloads, see sections below

Continue reading “Line.Sta Algorithms”

Bottle Sampling Depths

Bottle Sampling Depths, based on downcast profiles

here

The CalCOFI CTD-rosette is equipped with a Sea-Bird Electronic carousel water sampler (SBE 32), a computer-driven, electro-magnetically-released latch system. The 24 ten-liter plastic (PVC) bottles, equipped with epoxy-coated springs & Viton (non-toxic) O-rings, connect to 24 individual triggers by lanyards which keep the bottle ends open. During the downcast, profiles of different sensor measurements vs depth are displayed real-time on a computer screen. Based on the chlorophyll maximum & mixed layer depths, bottles are closed at specific depths to isolate the seawater. The 10 meter bottle spacing shifts up or down (see table below) to resolve steep gradient features such as chlorophyll, oxygen, nitrite maxima and shallow salinity minimum. Salinity, oxygen and nutrients samples are analyzed at-sea for all depths sampled. Chlorophyll-a and phaeopigments samples from the top 200 meters, bottom depth permitting, are also extracted for 24hrs and analyzed at-sea. Most CTD-rosette casts sample 20 depths to a maximum of 515 meters, bottom depth permitting. Occasionally, additional bottle depths or multiple bottles are tripped at the same depth to provide extra water for ancillary projects or primary productivity incubations. Two basin stations, off Santa Monica & Santa Barbara, are sampled beyond 515m to within 10m of bottom. Wire-length permitting, a 3500m deep cast is performed at sta 90.90.

Continue reading “Bottle Sampling Depths”

CTD Setup & Diagnostics

Seabird CTD-Carousel Setup & Deck Unit Diagnostics

Deck Unit 00000000
Deck Unit LED Panel Numbers

SIO-CalCOFI’s method of terminating the sea cable uses a custom 4-pin pigtail to terminate a multi-strand conductive wire that allows each individual conductor to be used in any combination. Terminating the sea cable using this technique may be viewed in the CalCOFI Handbook: CTD Termination.

Conventional two-pin Seabird pigtail terminations are standard on NOAA vessels and are usually performed by the ET. So interfacing the CTD deck unit only requires the attachment of the sea cable to the deck unit. This is done by connecting a two-conductor wire, such as Seabird part no. 31371 or a coaxial cable equipped with MS3106A connector, from the winch junction box in the CTD lab to the Sea Cable port on the back of the deck unit (see the deck unit section below).
Once the sea cable is connected to the Seabird 9plus and Seabird 11plus deck unit, the deck unit can be powered on. If cabled properly and all electronics are functioning, the numeric LED panel of the deck unit should display non-zero numbers. These values are different for each channel which can be selected by rotating the WORD SELECT dial to the right of the LED panel.
If the LED panel displays “0.0.0.0.0.0.0.0” then the deck unit and CTD “fish” are not handshaking correctly. Turn the deck unit off and wait 60 seconds before touching any exposed wire. Typically, with only two pins, the problem is reversed polarity so switching the wires in the junction box will correct this issue. Please note that powering on the deck unit with reversed polarity may blow the fuse(s). So after switching the wires and before re-powering on the deck unit, check both the 1/2A and 2A fuses on the back of the deck unit & replace if necessary. It is prudent to keep a good supply of these fuses on hand. Power on the deck unit and hopefully, the LED panel has numbers.

Please note in the photos, yellow tape is used to keep the deck unit Signal Source set to Fish. If the switch is toggled to Tape, the deck unit will not communicate with the CTD. The yellow tape prevents the switch from changing from Fish to Tape unintentionally.

Continue reading “CTD Setup & Diagnostics”

Distance From Shore

Distance = Distance from shore in nautical miles. The CalCOFI database includes a calculated distance-from-shore value used primarily for vertical section plots published in the hydrographic data reports. CalCOFI line coastline intercepts, digitized from Google maps, are tabulated below; the latitudes and longitudes are the origins (zero starting point) for vertical section plots.

Line Sta (if calculated) Latitude N Longitude W
 93.3, 93  26.2  32.978982  -117.272386
 90.0  27.6  33.497430  -117.741800
 86.7, 87  32.1  33.917500  -118.430800
 83.3, 83  39.1  34.273350  -119.308200
 80.0  50.3  34.474500  -120.476500
 76.7, 77  47.23  35.148000  -120.653000
 73.3, 73  49.6  35.653333  -121.221667
 70.0  50.9  36.196667  -121.723333
 66.7, 67  46.9  36.903333  -121.845000
 63.3, 63  49.3  37.398333  -122.425000
 60.0  49.0  37.988333  -122.815000

Basic Source Code (translated from F77 Fortran)

Public Static Sub OLDDS(LAT1!, Lon1!, D!, Liine$)
‘C
‘C — ROUTINE TO RETURN DISTANCE IN NAUTICAL MILES GIVEN TWO POSITIONS
‘C    IN DECIMAL DEGREES.  USES A METHOD THAT WILL MATCH THE RESULTS
‘C    PRODUCED IN THE OLD VERTICAL SECTION PROGRAM.
‘C
‘C  -> LAT0  — LATITUDE OF FIRST POINT. (DECIMAL DEGREES. SOUTH NEGATIVE.)
‘C  -> LON0  — LONGITUDE OF FIRST POINT. (WEST NEGATIVE)
‘C  -> LAT1  — LAT. OF SECOND POINT.
‘C  -> LON1  — LONG. OF SECOND POINT.
‘C  <- D     — DISTANCE (NAUTICAL MILES) BETWEEN POINTS.
‘C
  Dim COSD#
‘      Real LAT1, LON1, LAT2, LON2, D
‘      Real L, DLO, COSD, PI, F, A
      PI# = 3.141592654
      F# = PI# / 180
‘C
‘      INTRINSIC COS,ABS,SQRT
‘C
      ‘COSD#(A#) = Cos(A# * F#)
      If frmSelData.txtDataIn(6) = “” Then
        LineIn! = (Val(Liine$))
      Else
        LineIn! = Val(frmSelData.txtDataIn(6))
      End If
      D! = 0
      Dist2Orig$ = “”
190
    Select Case CInt(LineIn!)
    Case 93, 94 ‘, 93.3
            Lat0! = 32.969    ‘93.3 26.2 Google projection 2010; old value 32.95298 from Ralf & Nav Program
            Lon0! = -117.2696 ‘Google projection 2010; old value -117.2664
    Case 90, 91
            Lat0! = 33.49743   ‘90.0 27.6  Google checked 2010    ‘4969634422
            Lon0! = -117.7418      ‘745017037
    Case 87, 86 ‘, 86.7
            Lat0! = 33.9175    ‘86.7 32.1 Google projection 2010; old .89568      ‘918195424
            Lon0! = -118.4308  ‘Google projection 2010; old .419      ‘431669307
    Case 83, 84  ‘, 83.3
            Lat0! = 34.27335    ‘83.3 39.1  Google checked 2010  ‘2739356345
            Lon0! = -119.3082     ‘310048755
    Case 80, 79
            Lat0! = 34.4745     ‘80.0 50.3 Google projection 2010; old .4631       ‘4823190269
            Lon0! = -120.4765   ‘Google projection 2010; old .4739     ‘485679527
    Case 77, 76  ‘, 76.7
            Lat0! = 35.148     ‘76.7 47.23  Google checked 2010  ‘1446809515
            Lon0! = -120.653     ‘650653673
    Case 73  ‘, 73.3
            Lat0! = 35.6533333333
            Lon0! = -121.221666667
    Case 70
            Lat0! = 36.1966666667
            Lon0! = -121.723333333
    Case 67  ‘, 66.7
            Lat0! = 36.9033333333
            Lon0! = -121.845
    Case 63
            Lat0! = 37.3983333333
            Lon0! = -122.425
    Case 60  ‘, 57.8, 57.9
            Lat0! = 37.9883333333
            Lon0! = -122.815
    Case Else
     Select Case LineIn!
      Case 91.7, 91.6, 91.8, 92
            Lat0! = 33.25773      ‘from Ralf & Nav Program
            Lon0! = -117.4385
      Case 88.5, 88.6, 88.4, 88
            Lat0! = 33.69985      ‘from Ralf & Nav Program
            Lon0! = -118.0522
      Case 85.4, 85.5, 85.6
            Lat0! = 34.01588      ‘from Ralf & Nav Program
            Lon0! = -118.8243
      Case 82, 81.9, 82.1
            Lat0! = 34.435        ‘from Ralf & Nav Program
            Lon0! = -119.955
      Case 81.7, 81.8, 81
            Lat0! = 34.4174       ‘from Ralf & Nav Program
            Lon0! = -119.8033
      Case 78.5, 78.4, 78.6
            Lat0! = 34.77699      ‘from Ralf & Nav Program
            Lon0! = -120.6302
      Case Else
        If LineIn! <> Val(Liine$) And LineIn! <> Val(StndL$) Then
          LineIn! = Val(StndL$)
          GoTo 190
        Else
          Exit Sub
        End If
      End Select
  End Select
‘C
      Lon0! = Lon0! * -1
      L! = Abs(LAT1! – Lat0!)
      DLO! = Abs(Lon1! – Lon0!)
      If (DLO! > 180) Then DLO! = 360 – DLO!
‘C
      L! = L! * 60
      DLO! = DLO! * 60
      DLO! = DLO! * Cos(F# * ((Lat0! + LAT1!) / 2))
      D! = (Sqr(L! * L! + DLO! * DLO!)) * -1
      Dist2Orig$ = Format(D!, “#0.00”)
      If D! = 0 Then Dist2Orig$ = “”

End Sub

Underway Data Processing

Underway Data Processing Methods

1.       Navigate to Underway Data folder and join into a single file, do so using Advanced File Joiner or similar:  http://download.cnet.com/Advanced-File-Joiner/3000-2094_4-169639.html



a.     CalCOFI collects its own data while aboard SIO ships (New Horizon).  These are saved in the Events folder of the appropriate cruise.  A single file name is similar to: 0901NHMETEV_06012009.csv.  Although each filename is labeled by date, date data my not correlate.  This is of no particular consequence, just be aware of it.



b.     When a cruise is aboard a NOAA ship, it is normal to receive a CD-ROM at the cruise with a variety of SCS data.  These data formats are highly variable by ship and by year.  It will usually take a combination of files to find all that you will require (i.e., one file type for TSG, one for fluorescence, one for lat long).  Just merge each file type, then cut and paste columns into the correct order.  Make all SCS data available on the web.



c.     The resultant  file becomes the RAW UNDERWAY DATA for that cruise.  This file contains all unedited data from the Shipboard Computer System (SCS).  This file may include navigational, meteorological, oceanographic, and other mechanical data.  In your local working directory  (e.g., ../UnderwayData/2009/0911NH) save raw data in tab delimited ASCII format and name as YYMMSS_UnderwayRaw.txt  (e.g., 0901NH_RawUnderway.txt)



d.      *Links to SCS information

                               i.    http://oceanexplorer.noaa.gov/technology/tools/scs/scs.html

                               ii.    http://www.unols.org/meetings/2006/200610inm/SessionIIIa/SessionIIIa_Shields.pdf

                               iii.    http://sdsioa.ucsd.edu/



2.       Open merged file and crop data columns to include the following, in order and with the following column headers:



a.     DATE_UTC, TIME_UTC, LA, LO, ST, SA, FL, FI  (Date(UTC), Time(UTC), Lat, Lon, Temp, Salinity, Fluor, Pump)



                                          i.    *Note – There is the possibility that power gets toggled on/off to underway sensors during the cruise.  This results in changes in the total number of data columns (i.e., the first days may have 100 columns of data and the next few only 98).  Attention  is needed to make sure that when cropping data that the correct data is sub sampled.



                                         ii.     **Note – Different ships have different underway column headers/variable designations/configurations.  Consult Underway Data Formats (link) for more information.



3.        Delete rows with:



a.     pump speed not within than manufacturer’s specifications



                                          i.     Doubling the flow rate from what we recommend will have at least one main problem. This flow rate will probably be too high for the conductivity cell and overtime you can see a degradation of the conductivity readings due to this flow rate. What can happen is that the high flow rate may slowly remove the platinum on the electrode within the conductivity cell, especially in areas with a lots of debris and particulate in the water. This will change the area of the conductivity cell which can cause the readings to be high of correct. I am unsure how long this process would take but it is something we have seen before. I really suggest using the recommended flow rate for this instrument. ~SBE Technical Support



b.     no spatial coordinates



c.     erroneous values



                                          i.    *Note – There is the possibility that one or more columns may be missing some or all of the data.   There are several workarounds depending on the missing variable. Consult WTF Moments (link) for more information.



                                         ii.    **Note – With upwards of 100,000 rows of data it is inevitable that most bad values will go unnoticed at this point.  Subsequent data processing steps (i.e., plotting) will make it easier to find errors and to remove them.  Just try to find the large and obvious discrepancies. 



4.    Save cropped file as 0901NH_UnderwayCrop.txt  in same folder as the raw file.



5.    It is necessary to compare and calculate correction regressions based on the comparison of CTD to Underway Data.  Navigate to the final CTD Downcast .csv file folder and join all csvs into a single file (as in Step 1).  These csv files are the ones with IEH bottle data added to them.  In your local working directory  (e.g., ../UnderwayData/2009/0901NH) save raw data in tab delimited ASCII format and name as YYMMSS_CTDd.txt (e.g., 0901NH_CTDd.txt).  If the YYMM_CTDd.csv file is available open in Excel and crop columns accordingly (#7).  



6.    Open file with MS Excel and select the first column ‘A’ and Text-to-Column the Data with comma delimiter.



a.   *Note – DO NOT save this file as YYMMSS_CTDd.txt while in Excel.  A host of things can go wrong.  Data will be Cut and Pasted into a new file shortly.  If you must you can save it as some tempfile.xlsx.



7.        Crop data columns to include the following, in order and with the following column headers:



a.      Date_Time_UTC, Lat_Dec, Lon_Dec, Temp0, Temp1, BTL_Depth, BTL_Salt, BTL_Chl



                                          i.    *Note – If Underway Data was only available with Date and Time in PST select Date_Time_PST instead.



8.    Sort Ascending (Smallest to Largest) based on BTL_Depth and Delete all depths greater than 2m, all rows with no BTL_Depth values, and all Header rows except for one which needs to be cut and pasted back at the top after we are finished sorting.



9.    Sort Ascending based on Date_Time.



10.  Delete any rows that are missing any information.



11.  Delete any duplicate stations.  You can tell if two subsequent rows are from the same station by the time, delete  the deeper.



12.  Make sure Header row is at the top.



13.  After the Temp 1 column, insert 2 additional blank columns.



a.   Average the two Temp columns into the first blank column, then copy and paste ‘values only’ into the second.  Delete the Temp0, Temp1, and first column.  Header the remaining column as TempAvg.



14.      Add three extra headers/columns  at the end of the existing



a.      …, ST, SA, FL, FI (…, T, S, Fluro, Pump)



15.  Open a previously completed copy of YYMMSS_CTDvsUnderway.xlsx.  Select, Copy, and Paste all rows and columns from the sorted and corrected …CTDd.txt file into the recently opened xlsx file.  Save and RENAME the xlsx file accordingly (i.e., NEW YYMMSS ).  Make sure no old data lingers in the spreadsheet as this would cause errors in the fit of the regression.



a.   *Note – We reuse the xlsx file because it already has multiple graphs set up.  You will of course need to resize the data sets for each graph as necessary, and rename the graph titles.  But if you followed the above steps the columns will be in the right places and everything should fall into place.



16.  Open the …UnderwayCrop.txt in MS Excel, Tab delimited.  Arrange both files in Tile View.



17.  Using this new …CTDvsUnderway.xlsx file’s Date_Time as a guide, scroll down the …UnderwayCrop.txt and find a corresponding date and time.  There will be more than one row of data for the time you are trying to match.  Try not to bias the selection of the Underway data and just pick the first one you see. Select the four columns, Temp, Salinity, Fluor, Pump, from UnderwayCrop and Copy them, and Paste them to the xlsx file in the last four columns ST, SA, etc.  Do this for all the rows in the CTDvsUnderway file.



a.     Because the number of rows will change from cruise to cruise, the graph data source will need to be adjusted accordingly.  This is important as the linear regression equations shown on the graphs are used to correct the Underway Data in subsequent steps.  Save file.



18.  For export to MATLAB it is necessary to Cut the xlsx data (numbers no graph info) and paste it into a new UltraEdit .txt file. Name it accordingly (YYMMSS_CTDvsUnderway.txt).  This is done to remove any .xlsx formatting that may cause MATLAB to choke on import.



19.  Start Matlab and load a previously completed .m file as a template (e.g., 0901NH.m).  The .m file will have incorrect variable names and notes pertaining to that specific cruise.  Change and update accordingly.  Also make sure you are working in your current underway folder.



a.     Import UnderwayCrop.txt and follow along the notes in the script for error checking and saving data as a .mat file.



b.      Import …CTDvsUnderway.txt.  There will be an error with this as the Date_Time will not ingest properly.  Ignore it, the date is not important for plotting  unless GPS has failed and you need the data and time to correlate with the Event Log GPS to get some coordinates to make these maps.



c.      An important step to be aware of is toward the end of the script it asks you to use dlmwrite function to export the final underway data product.



20.     The template m-file will guide you through the entire process.  When you get to the regression corrections take the values from the xlsx file graphs that you generated earlier.  When you make the Matlab plots they should print into your working folder.



21.     Open a previous yymmplot.ai file in Adobe Illustrator and the first plot from Matlab (i.e., t0901NHpath.eps).  Save the .ai file with the current cruise YYMM.



a.     In dealing with the .eps file from Matlab you will find that all we really need is the color path and the gradient bar frame and numbers.  Ungroup and delete all you don’t need (coastline, grid, tick marks, gradient color, etc.)



b.     Since the path is thousands of small circles, select them all, group them, then copy them to the clipboard. Delete the preexisting path in the .ai file and make sure you paste the new one in the correct .ai layer and resize to fit the station layout / cruise path.  This isn’t an exacting process, just make it fit best you can.  Make sure to lock/unlock layers as you go to prevent the pasting / editing in the wrong place.



c.     Go back to the eps and copy the gradient bar box, tickmarks, and numbers.  Since each cruise’s graph has a different range of values these gradient bars change.  The colors stay the same (if you followed the Matlab m-file correctly) but we need to make the new numbers fit the template.   Delete the old box and numbers and paste in the new and resize accordingly.  Resizing will skew the font so once everything is in the right place check the font setting and set to 18pt Myriad Pro with 100% font stretching in both directions.  You will also have to change the decimals and numbers a bit as Matlab tend to separate the numbers after the decimal in eps format.  Just delete the end and retype the number using the left half.



d.     Unlock and Edit  the Title layer to the correct cruise designation and dates and title.



e.     In Illustrator to save as pdf “save a copy”, to save as png Use Export (white background, medium resolution, anti aliased) .  Open each png in Photoshop and shrink to 75 pixels keep proportionality and rename with a         ” _t” at the end for the web thumb.

31.   Upload all data and images to calcofi.com FTP Underway_Data directory.

32.   Edit calcofi.com Joomla site with new underway data accordingly.

 

Primary Productivity Methods

PRIMARY PRODUCTIVITY

 


OVERVIEW: Primary production is estimated from 14C uptake using a simulated in situ technique in which the assimilation of dissolved inorganic carbon by phytoplankton yields a measure of the rate of photosynthetic primary production in the euphotic zone.


 

1. Principle

Seawater samples are incubated with a radioactive substrate to determine the incorporation of inorganic carbon into particulate organic carbon due to photosynthesis at selected light levels.  The data have units of mg-carbon per m3 per half day.

2. Productivity Cleaning Procedures

2.1.

Micro-90 Cleaning solution is diluted to 2% solution using de-ionized water (DW).  Hydrochloric acid (HCl) Trace Metal Grade, Fisher Scientific, solution (1.2M) diluted with DW. Acid-washing of Teflon should be done with great care as Teflon is porous to HCl which can compromise dilute basic stock solutions of 14C -bicarbonate.

2.2.

250 ml polycarbonate incubation bottles are filled to capacity with 2% MICRO for 3 days with the cap on in an inverted position. Next, rinse all Micro away and then rinse down the walls with 20 -30mls 10%HCl and recap and shake  to acid rinse inside bottle.  This should be left overnight 12-16 hours. The acid is removed by rinsing the bottles three times with milliQ clean water before air drying.

2.3.

10 liter rosette sample bottles are cleaned with a 2% MICRO soak for 3 days, rinsed with de-ionized water and then dipped in 10% metals free HCl.  Caps, special coated springs and valve assemblies are also cleaned with a 2% MICRO soak for 3 days and then rinsed with de-ionized water and dried. 

2.4.

All lab ware to be used is cleaned in this manner.

 

3. Preparation  of Isotope Stock

3.1.

To prevent contamination of self or solutions, work with the isotope stock is performed wearing vinyl gloves.

3.2.

A solution of 0.3 g of Na2CO3 anhydrous (ALDRICH 20,442-0, 99.995%) per liter Milli-Q filtered DW in a Micro cleaned 1 liter Teflon bottle to yield a concentration of 2.8 mM Na2CO3.   This solution is filtered through 0.2µM Nucleapore filter to remove particulate carbonate.

3.3.

Concentrated stock, 50ml of NaH-14CO3 (~50-57 mCi mmole; MP Biomedicals LLC.) was diluted with 350 ml the 2.8mM Na2CO3 solution in productivity-cleaned 1 liter polycarbonate graduated cylinder.  It has become necessary to pH this up with an ultra clean 1N NaOH solution to raise the pH to ~10.

3.4.

Specific activity can be checked by diluting the above made solution to working concentrations, ie 50-200µl added to 250ml polycarbonate centrifuge bottle and measuring out triplicate 1ml portions into beta ethanolamine spiked (1.5%v/v) Ecolume scintillation cocktail.

3.5.

 To check for  14C-organic carbon contamination another working aliquot of 200µl can be placed into a scintillation vial and acidified with 0.5ml 10% HCl and placed on a shaker overnight. This is done in the hood as it liberates 14C-CO2. The acidified dpm should be <0.0001% of the total dpm of the 14C preparation.

 

4. Incubation Systems: situ incubation techniques

4.1.

 Incubation apparatus consists of seawater-cooled, temperature monitored incubator tubes wrapped with neutral-density screens which simulate in situ light levels.    

 

 

4.2.

Six incubation depths are selected, they represent 56, 30, 10, 3, 1  and ~0.3 % light level.  These values are estimated using a wand type PAR meter after cleaning tubes and screens covering them.  The near surface light level is reduced to 56% using common plastic screening to prevent a lense effect and subsequent cooking of the surface samples. 

 

5. Sampling

5.1.

Primary productivity samples are taken each day shortly before local apparent noon (LAN).  Light penetration was estimated from the Secchi depth (Using the definition that the 1% light level is three times the Secchi depth).  The depths with ambient light intensities corresponding to light levels simulated by near surface and the on-deck incubators were identified and sampled on the rosette up-cast.  Extra bottles were tripped in addition to the usual 20 levels sampled in the combined rosette-productivity cast in order to maintain the normal sampling depth resolution.

5.2.

 


5.3.

Using a dark sleeve to subdue the light, water samples are transferred to the incubation bottles (250 ml polycarbonate bottles) and stored in a dark box until inoculation. 

 

Triplicate samples (two light and one dark control) were drawn from each productivity sample depth.

 

6. Isotope Addition and Sample Incubation

6.1.

Samples are inoculated with 50-200 µl of 14C as NaHCO3 stock solution of sodium carbonate (Fitzwater et al., 1982).

6.2.

Samples are incubated from LAN to civil twilight in a surface seawater-cooled incubators with neutral-density screens which simulate in situ light levels, corresponding to those from which samples were taken (see 4.2).

6.3.

At civil twilight the incubation is terminated and the time noted.  Sea state and safety is the only exception accepted to delay the end time.

 

7. Filtration

7.1.

At the end of the incubation, all bottles have subsamples of 10mls removed for DO14C analysis. The LTER DOC filtrate apparatus consists of a plexi-glass filtration manifold to hold up to 18 scintillation vials over which syringe needles with 0.45um equivalent micro-syringe filters can be passed through stoppers with 25 ml syringe bodies serving as filter funnels.  The exception to this is dark bottles are only sampled for DO14C on two each high and low chlorophyll stations.

7.2.

Additionally, from dark bottles a 1ml sample is placed into beta mercapto-ethanol spiked (1.5%v/v) Ecolume scintillation cocktail  to determine the specific radioactivity in the samples.   These values are used to calculate an average cruise value after removing outliers.

7.3.

Finally the samples are filtered onto Millipore HA filters and placed in scintillation vials.  One half ml of 10% HCl was added to each sample.  The samples are then allowed to sit, without a cap, at room temperature for at least 3 hours (after Lean and Burnison, 1979).

 

8. 14C Sample Processing

8.1.

After addition of 10mls of Ecolume cocktail, vials are tightly capped and mixed before vials are counted for up to 10 minutes each for 14C  on a Beckman 6100LC liquid scintillation counter set to 1.0% counting precision.

8.2.

Data is captured to a flat file using Beckman data capture software for Windows in ASCII format.  This format is then used to integrate productivity depths into the CalCOFI data processing flow.

 

9. Calculations

Data is presented as mean mg Carbon assimilated per meter cubed of seawater for one half light day. 

mgC/m3 per one half light day =  ((Sampledpm Blankdpm) x W)/R, where

 

   W = 25200 = 12,000 x A  x FT  x 1.05

   12,000= molecular weight of carbon in milligrams

   A = carbonate alkalinity (milliequivalents/liter)

   FT = Total carbon dioxide content/ carbonate alkalinity

   1.05 is the 14C isotope fractionation factor, reflecting preferential use of C12 over C14 by a factor of 5%

   R = dpm added to sample (µCi/200µl x 2.2 x 106)

 

To better understand this equation and variables see Strickland and Parsons (1968).

10. Equipment/Supplies

·         10 liter pri. prod. cleaned sampling bottles

·         Secchi disk

·         Re-pipet dispensers for delivering 20µl, 200µl, 0.5ml

·         Pipets able to measure 1ml and 10ml

·         250 ml polycarbonate centrifuge bottles

·         liquid scintillation counting (LSC) vials

·         Seawater plumbed incubation rack with neutral density screening.

·         Par meter, wand type (Biospherical Instruments)

·          14C sodium bicarbonate stock solution (MP Biomedicals, LLC)

·         Millipore Type HA filters (Fisher Scientific)

·         vacuum filtration system including separate device for DOC filtrate capture

·         Polycarbonate centrifuge bottles

·         Teflon laboratory wares

·         vortex mixer

·         liquid scintillation counter (LS 6000LC Beckman Instruments, Inc.)

11. Reagents

·         Milli-Q filtration/anion exchange water purifier

·         Micro-90 Cleaning solution, Cole Palmer Instrument Co.

·         HCl for trace metal analysis (Fisher Chemical)

·         Na2CO3 (99.995%) Aldrich Chemical

·         NaH-14CO3 solution (cat #17441H MP Biomedicals, LLC.)

·         2-amino ethanol (ethanolamine)  ACS grade

·         Aquasol-II (Dupont)

·         Ecolume (MP Biomedicals, LLC.)

 

12. Re-count check

14C scintillation counts were checked for accuracy by re-counting an entire cruise (n>200) of vials 9 months after original counting.  Depletion due to half life was ignored due to the long half life of 14C.  Results for samples greater than 1000dpms were averaged resulting in a return of counts equal to 101.3%.  Efficiencies had a similar recount statistic of 100.9%.  The exercise lead to evaluating cruise counts where the source of some replicate inconsistency was the result of chemiluminescence problems in which the counter displays a “lumex %”. It is important to monitor for higher lumex numbers which result in elevated counts due to a chemiluminescent reaction. Samples were dark adapted and recounted resulting in much better replicates.

13. References

 

·         Fitzwater, S. E., G. A. Knauer and J. H. Martin, 1982.  Metal contamination and its effect on primary production measurements.  Limnol. Oceanogr., 27: 544-551.

·         Lean, D. R. S. and B. K. Burnison, 1979.  An evaluation of errors in the 14C method of primary production measurement.  Limnol. Oceanogr., 24: 917-928.

·         Steeman-Nielsen, E. (1951). “Measurement of production of organic matter in sea by means of carbon-14”. Nature 267: 684–685.

·         Strickland, J. D. H. and Parsons, T. R. 1968.  A Practical Handbook of Seawater Analysis pp. 267-278.

 

Methods Timeline

CalCOFI Time Series data are over 67 years old. New adopted practices, methods, software and hardware are thoroughly tested to maintain dataset continuity as the program & science evolves. Core measurements are maintained and many new measurements added. CTD temperature sensors, for example, provide data at a much higher resolution than a 20 bottle hydrocast equipped with reversing thermometers pre-9308 (Aug 1993).

Although the CTDs casts on CalCOFI started in 1990, CTD-rosette casts did not replace bottles-on-wire hydro casts completely until Aug 1993 (9308). In 2004, LTER joined the CalCOFI cruises, expanding the seawater analyses, adding new measurements.

Changes of standard practices, methods, software, & equipment will be tabulated here, particularly those affecting the hydrographic data or other data products.

 Cruise/Date

 Measurement, 
 Equipment,  Method

 Changes    

Comments

CalCOFI Methods Timeline
  Methods Timeline Started  

This page is ‘dynamic’ so new items will be added whenever new methods, measurements, equipment are implemented. Older changes will be added as time-permits and they are remembered.

1708SR CCE-LTER ISUS upgraded to fw v3 ISUS upgraded to firmware v3 Seabird acquired Satlantic and took several month to get ISUS service online. Once done, they upgraded the CCE-LTER MBARI-ISUS v2 to firmware v3. Allows easy in-house or at-sea calibration.
1708SR Oxygen Stable Isotope Seawater samples collected for OSI 5ml of seawater collected from two depths (10m & 50m) on 18 stations
1701RL
1704SH
LARS rosette 24-btl LARS rosette used Even though RV Reuben Lasker & FSV Bell M Shimada do not have a LARS CTD deployment system, the height clearance allows the use of the LARS rosette. SIO-CalCOFI will use this rosette as primary unless the height clearance becomes a problem, such as on RV Ocean Starr.
1611SR First RV Sally Ride SIO RV Sally Ride used for CalCOFI RV Sally Ride used for the first time to do CalCOFI station work; acoustics operational but not calibrated
1611SR RINKO III Optode RINKO III Oxygen sensor deployed Optode used on non-basin stations deeper than 500m, where altimeter was not needed
1611SR Seawater samples Domoic Acid & Th-234 Collected ancillary seawater & deploy insitu pumping system
pre-1611SR ISUS (Frank’s) Repaired Peter Frank’s ISUS repaired & upgraded to firmware v3 ISUS was not working so it was sent to Satlantic who repaired it and upgraded the firmware to v3
1609SR New research vessel SIO RV Sally Ride test for CalCOFI Shakedown cruise for RV Sally Ride to do CalCOFI – new LARS, epoxy-coated rosette deployed
1609SR LARS epoxy-coated rosette First use of new 24-bottle rosette 24-btl white epoxy coated aluminum rosette deployed with stainless steel LARS support
1607OS Seawater Samples Domoic Acid samples collected Ancillary seawater sample collection
1604SH Nutrient Analysis New (returning) in-house nutrient chemist Nutrients run by SIO-CalCOFI analyst DGS
1507OC Underway Measurements RV Oceanus MET system logs data into individual data files; combined data are unavailable unless merged manually Unlike other UNOLS vessels, particularly SIO RV New Horizon, the underway TSG & meterological data are logged at 1Hz into individual files. A combinded data file will have to be generated post-cruise by merging individual sensor data files using the common UTC timestamp. Software pending…
1504NH Underway Measurements PCO2 and pH measurements added to CCE-LTER suite

 

1504NH Nutrient Analysis Standard matrix; sample vials

Standards now prepared in low nutrient seawater (collected from the end of lines 93.3 or 90.0 and processed with UV light, filtration and aged before use).  New 30 mL polypropylene tubes in use.

06 Aug 2014 calcofi.com calcofi.com 2014 online

Developmental site web.calcofi.com replaces Joomla 1.5 calcofi.com; Joomla 3+ fully implemented. calcofi.com old site moved to old.calcofi.com

1402SH Nutrient Analysis New in-house Nutrient Technician

Nutrients run by SIO-CalCOFI analyst LJE.

1402SH NCOG DNA/RNA NCOG sample collection started NCOG DNA/RNA samples collected at ~4 depths (10m, chl max, 170m, 515m)
1311NH DIC measurments Sampling expands

Dissovled inorganic carbon samples are now drawn at more locations along the cruise track with 14 profile and 8 additional surface water stations per cruise.

1211NH Rinko Oxygen Optode RINKO III Optode deployed RINKO III Oxygen Optode sensor replaced 2nd SBE43 oxygen sensors on stations where the altimeter was not needed – non-basin stations deeper that 500m.
01Mar2013 web.calcofi.com website update Upgrading calcofi.com to latest Joomla

web.calcofi.com started to migrate calcofi.com 1.5 to 3.0; to improve responsiveness with new dynamic templates – desktop, tablet, and smartphone auto-formatting. Security improvements, auto-updating, jQuery & other new features of version 3.+

1207OS Nutrient Analysis New in-house Nutrient Technician

Nutrients run by former ODF chemist, now a SIO-CalCOFI analyst MTM.

01May2012 dev.calcofi.com blog started SIO-CalCOFI Technical Group developmental blog started.

For metadata: documenting changes in software, data-processing methods, formats, products, & practices.

1203SH Nutrient Analysis Instrument and analyst change, new Seal QuAAtro Analyzer purchased.

Transitioned to in-house nutrient analysis on new Seal QuAAtro Nutrient Analyzer, nutrients run by SIO-CalCOFI analyst DNF.  Replacing sample analysis by ODF chemist on an AA3 analyzer.

Standards now prepared in artificial seawater.

NO3 data calculated using regression coefficients from ISUS voltage vs NO3 for 1202NH & 1203SH look nearly identical.

1203SH atsea.calcofi.com blog online Setup to keep notes during quarterly cruises, particularly those related to data, equipment, & generally noteworthy.

Connected-linked to CalCOFI’s Twitter feed

1202NH Nutrient Analysis  Last cruise where nutrients were run by ODF chemist on an AA3 analyzer.

 

1202NH+ Data Processing IEH retired as primary data processing file format. Sta.csvs and casts.csv adopted. After parallel-processing (IEH & csv file format) 1104, 1108, 1110 cruise data. Sta.csvs & casts.csv data processing is being used to process, merge, quality-control, publish all hydrographic data. IEH-format is a data product along with the hydrographic database.
1108 All measurements & data products Hydrographic data are processed using both old & new methods and compared.

The cruise is once again processed twice, in parallel, using the IEH old-school method and the new csv-format method. GTool development & refinement continues as sta.csvs & casts.csv formats are improved. Final data products are compared for agreement. This is the first cruise to generate data products using new csv-database data processing methods not based on IEHs.

1104 GTool Matlab program developed SIO-CalCOFI Technical Team implement a new graphical matlab tool to point-check CalCOFI hydrographic data

MGS, in cooperation with the CCTG, develops GTool to replace Andyplots, our legacy method to visually inspect all bottle data on one figure. CTD continuous data are plotted with bottle parameters to cross-check and eliminate fliers.

1104 All measurements; data products Hydrographic data are processed using both old & new methods and compared.

The cruise is processed twice, in parallel, using the IEH old-school method and the new csv-format method. GTool, a matlab plotting & data-quality control program script, is developed to ingest CTD & sta.csvs to graphically assess CTD & bottle data-quality. Final data products are compared for agreement. This should be the last cruise processed using 00/20/22 & IEHs.

1101 Reported standard level O2 CTD primary or secondary oxygen sensor bottle–corrected measurements are reported instead of interpolated standard level (ISL) oxygen values Pretains to final, standard-level bottle data reported in all data products – with the reliability of the SBE 43 oxygen sensors and consistently high data quality when calibrated bi-yearly (twice/year) and bottle-corrected. Bottle-corrected CTD primary or secondary oxygen values (whichever sensor is performing better that cast) are reported instead of interpolated (calculated) standard level bottle oxygen values.
1101  Data storage format, ieh phase out begins; all measurements affected Sta.csv & casts.csv implementation begins, replacing 00/20/22 & IEH data process methodology. IEH-method is used as CSV-method is developed.

Bottle data & CTD data processing up til now have been relatively separate although data-quality cross-check are common. Other than CTD temperatures, all hydrographic data reported in IEHs or hydrographic database are bottle sample measurements. In 2011, our new data processing strategy merges the bottle & CTD starting at seawater sample collection. CESL generates individual sta.csvs combining CTD & bottle data immediately after seawater sample collection. Casts.csv is also generated and contains station specific information, replacing the Station Cast Description & Weather files (stacst & weather).

1101  DECODR DECODR migration to the new sta.csv and casts.csv data-processing scheme in full development

DECODR (Data Entry & Compile Oceanographic Data Reports) program modules are adapted to process either old or new hydrographic data formats. Sta.csvs & casts.csv adaptations are implemented and different data products are generated directly or by generating an IEH first then data products.

0901  DIC measurements resume DIC (dissolve inorganic carbon) sampling resume after a hiatus.

DIC, aka Keeling, samples were commonly taken on two CalCOFI stations but the new sampling scheme covers several additional stations and depths.

2008  Nutrient Analysis Began ammonium analysis

Nutrient analysis expanded to include ammonium.

0701  Salinity measurement SalReCap – Salinity Record & Capture program introduced, replacing the PSal (DOS shell executable).

SalReCap, Windows program developed by SIO-CalCOFI, replaces ODF-developed PSal, a DOS executable. Immediate comparisons to CTD salinities during analysis becomes available.

Apr-05  Primary Productivity New specific activity procedure The procedure for calculating cruise 14C specific activity for the productivity assay was changed to reflect daily 14C additions averaged over the course of a cruise.  The six dark bottles have one milliliter removed, added to ethanolamine spiked scintillation cocktail for counting.  Previously specific activity was calculated for a batch of stock and used for the entire cruise.  The new method served to check pipeting, inoculation amounts and any changes in volatile 14C stocks.
0411  Ancillary measurements such as HPLC, DOM, Size fractionations, Epifluorescence LTER affiliation on CalCOFI cruises begins

CCE-LTER (California Current Ecosystem Long Term Ecological Research Site) is established and LTER participation on all CalCOFIs begins. A large suite of additional measurements, particularly seawater samples from the rosette, and PRPOOS vertical plankton tow are added to CalCOFI hydrography.

0310  Oxygen measurement Autotitrator replaces manual modified Winkler titrations

After substantial assay-comparative testing, an ODF-developed oxygen auto-titrator is used to titrate discrete, rosette-bottle, seawater oxygen samples at-sea. New oxygen data format & data processing module is introduced in DECODR.

0204  CalCOFI Data Report Last bound hardcopy of the bi-annual CalCOFI Hydrographic Data Report is published, CC Reference 03-01 31 July 2003

CalCOFI Hydrographic Data Reports will continue being published electronically as pdfs for global distribution. Bound hardcopies, sent my mail, to different libraries & institutions will stop.

0104  Chlorophyll measurement FLog, chlorophyll data logging fluorometer program introduced

Prior to FLog, our 24+hr acetone cold-extracted chlorophyll samples were manually logged by hand to a sample form then key-entered into CODES. FLog records fluorometer values automatically in a data-processing friendly format so no transcription is required.

Jun 2001 Seasave for Windows Seabird releases Windows versions Seasoft/Seasave & SBE Data Processing software for Windows 98/NT  released by Seabird
2000  CTD Fluorometer Seapoint Fluorometer

Seapoint Chlorophyll Fluorometer in use.

9807 (Jul 1998) CTD Seasave Con file readable by Seasave v5 con files prior to 9807 will not open it Windows SBE Data Processing software v5 or v7
9308  Temperature measurement 4sec ave (prior to bottle closure) CTD temperatures replace reversing thermometer measurements

CTD temperatures are the only measurement entering the hydrographic timeseries (until 2011 when CTD standard level temperature, salinity, & oxygen replace interpolated values).

9308  Seabird CTD-Rosette w/ 24-10L Niskin-type PVC Bottles Replace 20 bottle hydro cast & 6 bottle prodo cast

10L Niskin-like bottle constructed by Research support have plastic-coated springs, nylon lanyards, Viton or nitrile o-rings. No metal, rubber, or latex come incontact with seawater samples. Bottles are disassembeld &  ‘productivity-cleaned’ between cruises.

1993  CTD Fluorometer SeaTech Fluorometer

1993-2000

1973  Chlorophyll measurement Began

Discrete chlorophyll analysis was added to the hydrographic dataset.

1961  Nutrient Analysis Expands beyond phosphate

Nutrient analysis expanded to include silicate, nitrate and nitrite.

     

 

 

Chlorophyll Methods

Chlorophyll Determination

 


SUMMARY: Chlorophyll a is extracted in an acetone solution. Chlorophyll and phaeopigments are then measured fluorometrically using an acidification technique.


 

1. Principle

Seawater samples of a known volume are filtered (< 10 psi) onto GF/F filters. These filters are then placed into 10ml screw-top culture tubes containing 8.0ml of 90% acetone. After a period of 24 to 48 hours, the fluorescence of the samples is read on a fluorometer. Then samples are acidified to degrade the chlorophyll to phaeopigments (i.e. phaeophytin) and a second reading is taken. The readings prior to and after acidification are then used to calculate concentrations of both chlorophyll a and ‘phaeopigment’. The method used today is based on those developed by Yentsch and Menzel (1963), Holm-Hansen et al. (1965) and Lorenzen (1967).  Note that concentrations of ‘phaeopigments’ are not a good measure of Chl a degradation products present in the sample since Chl b present in the sample will be measured as ‘phaeopigments’.

2. Sample Drawing

2.1.

Chlorophyll bottles should be rinsed three times with sample prior to filling. The bottles are calibrated for volume, so the sample drawer must insure that air bubbles are not clinging to the sides of the bottle and it is filled completely.  The sensitivity of the fluorometric method allows for sample bottles of ~50 to 250 ml.

3. Sample Filtration

3.1.

Check that the filtration funnels are well seated on the base, and be sure that the filters (Whatman GF/F) are in place. Improperly placed filters or loose funnels will result in loss of sample. The chlorophyll samples are volumetric and should sample loss occur, replace the filter with a new one and redraw the sample.

3.2.

Turn on the vacuum pump, pour the sample into the filter funnel, and open the valve. Check the vacuum pressure to see that it does not exceed 10 psi or ~500mm Hg. Generally samples are filtered in such a way as to insure that the deepest samples (i.e. those typically containing less chlorophyll) are filtered at the same manifold positions in each time. When a shallow cast is performed and a reduced number of samples is taken, it is advisable to filter them on positions typically used for those approximate depths. This reduces the potential for contaminating filter funnels used for filtering deep samples that in general contain low levels of chlorophyll.

3.3.

When a sample has finished filtering, turn off the valve; once all the samples have filtered, turn off the pump; use designated sample forceps to pick off the filter and place it in the appropriate numbered tube containing 8 mL 90% acetone. Make sure that the filter is completely submerged in the acetone.

3.4.

Cap TIGHTLY but be aware that tube tops can break off, then place the sample tubes in a rack. The sample rack is then placed in a refrigerator and the filtration time is recorded.

4. Standardization of Fluorometer

4.1.

A commercially available chlorophyll standard (e.g. Anacystis nidulans, Sigma Aldrich) should be used to calibrate the fluorometer, preferably before and after each cruise. The Chl a standard is dissolved in 100% acetone to yield approximately 0.1mg-Chl per ml solution. 1ml of this solution can then be diluted in 100ml 100% acetone and read in a spectrophotometer at 664nm. A second reading at 750nm is also recorded as a blank value to correct for sample turbidity. The remainder can be aliquotted into cryo tubes and stored in liquid N2 for future use. Chl a standards such stored are stable for years. The initial dilution is made with 100% acetone because it stores better in liquid N2 than those made with 90%. However, since 90% acetone is used for the extraction, it is also used for dilutions when generating a standard curve.

4.2

The concentration (mg l-1) of the standard is determined by the following equation:

Chl a =           Equation1

A664 = absorption at 664nm

A750 = absorption at 750nm

E = Extinction coefficient (100% acetone = 88.15, 90% acetone = 87.67) from Jefferies and Humphrey (1975)

l = cuvette path length (cm)

4.3.

A series of dilutions using 90% acetone (N> 5) are then made and read, recording both Rb and Ra values. Blank values should be subtracted from the Rb and Ra prior to performing calculations. If using a fluorometer with multiple sensitivity and range settings such as a Turner model 10, then the proper blank value must be subtracted for readings taken at a given setting.

4.4.

A calibration factor (F) must be calculated for each fluorometer. It is the slope of the line resulting from plotting the fluorometer reading (x-axis) vs. chlorophyll concentration (y-axis). This line is forced through zero. An acidification coefficient (τ) is the average acid ratio (Rb/Ra) for the pure chlorophyll standards used in the calibration.

4.5.

Calculating chlorophyll and phaeopigment concentration in a sample is accomplished by using the following equations (Knap et al., 1996):

Chl (µg/l) =  Equation2

Phaeo (µg/l) =   Equation3

F = Linear calibration factor (see 4.4)

τ = Average acid ratio (Rb/Ra) – Note that these are actually corrected values, with the blank readings already subtracted.

Ve = Volume of extract (ml)

Vf = Volume of sample filtered (l)

S = Sensitivity setting of fluorometer (Applicable to Turner model 10. If using a model 10AU or another fluorometer, use a value of “1”)

Rng = Range setting of fluorometer (Applicable to Turner model 10. If using a model 10AU or another fluorometer, use a value of “1”)

There are variations of this equation that can be used and other factors that can affect chlorophyll measurements. More detailed descriptions can be obtained in Strickland and Parsons (1968) and Holm-Hansen and Riemann (1978).

Note: After a cruise, the fluorometer is calibrated again and the calibration factors and average acid ratios obtained from pre and post-cruise calibrations are averaged for final data processing.

5. Reading Samples on the Fluorometer

5.1.

The fluorometer should be allowed to warm up for approximately 1/2 hour before using it. Samples must extract in acetone for at least 24 hours prior to reading on the fluorometer and should be read before 48 hours.

5.2.

Samples must be at room temperature prior to reading. One hour before samples are to be read, they should be removed from the refrigerator and allowed to warm up in a dark place.

5.3.

A blank tube containing the same acetone batch used for the extractions should be prepared and read prior to reading samples. This blank should be read before and after every sample run and after door setting have been changed (Turner model 10 fluorometer)

5.4.

A coproporphyrin standard should be read prior to reading samples (D’Sa et al., 1997). While not used in any calculation, it is useful to monitor the performance of the fluorometer over time between calibrations. Significant changes in coproporphyrin readings may indicate a problem with the fluorometer.

5.5.

Remove the filter, shake the sample to insure that it is well mixed, and use a Kimwipe to remove fingerprints from the exterior of the tube prior to running samples.

5.6.

Read the sample and record the number (Rb). Add 100µl of 10% HCl and wait approximately 30 seconds for the number to stabilize and record the value (Ra).

6. Equipment/Supplies

· Whatman 25mm GF/F filters (Fisher Scientific)

· Volumetric sample bottles (~130-150ml)

· Vacuum filtration apparatus with vacuum pump capable of maintaining 10 p.s.i

· Fluorometer and proper filter kit for measuring chlorophyll a/phaeophytin with acidification method (Turner model 10AU uses a 10-037R optical kit).

· Pipet (or re-pipet) capable of delivering 100µl.

· Personal protection equipment (PPE) consisting of gloves and safety glasses.

· Kimwipes or equivalent laboratory wipes.

· 10ml screw-top sample tubes (Fisher Scientific)

· Two sets of forceps (one for sample manipulation and one for replacing clean filters)

· Assorted laboratory glassware, including volumetric flasks for diluting calibration standards

7. Reagents

· Milli-Q or equivalent polished water source.

· HPLC-grade or equivalent low-fluorescing acetone. Note that volume is not conserved when preparing solution of water and acetone. The addition of 413ml Milli-Q water to 3800ml of acetone results in 4130ml of 90% acetone.

· 10% HCl solution

· Chlorophyll a (Sigma Aldrich catalog number C6144)

· Coproporphyrin III tetramethyl ester (Sigma Aldrich catalog number C7157)

8. References

· D’Sa, E.J., Lohrenz, S.E, Asper, V.L., and Walters, R.A. (1997). Time Series Measurements of Chlorophyll Fluorescence in the Oceanic Bottom Boundary Layer with a Multisensor Fiber-Optic Fluorometer. , 167: 889–896. DOI: 10.1175/1520-0426(1997)0142.0.CO;2

· Holm_Hansen, O., Lorenzen, C.J., Holms, R.W., Strickland, J.D.H. (1965). Fluorometric Determination of Chlorophyll. J. Cons.perm.int Explor. Mer. 30: 3-15.

· Holm-Hansen, O., and B. Riemann. (1978). Chlorophyll a determination: improvements in methodology. Oikos, 30: 438-447.

· Jeffery, S.W. and Humphrey, G.F. (1975). New spectrophotometric equations for determining chlorophylls a, b, c1, and c2 in higher plants, algae and natural phytoplankton. Biochem. Physiol. Pflanz. 167: 191-194.

· Knap, A., A. Michaels, A. Close, H. Ducklow and A. Dickson (eds.). (1996). Protocols for the Joint Global Ocean Flux Study (JGOFS) Core Measurements.
JGOFS Report Nr. 19, vi+170 pp. Reprint of the IOC Manuals and Guides No. 29, UNESCO 1994.

· Lorenzen, C. J. (1967) Determination of chlorophylls and phaeopigments: spectrophotometric equations. Limnol. Oceanogr. 12: 343–346.

· Strickland J. D. H., Parsons T. R., (1968). A practical handbook of seawater analysis. Pigment analysis, Bull. Fish. Res. Bd. Canada, 167.

· Yentsch, C.S., Menzel, D.W. (1963). A method for the determination of phytoplankton chlorophyll and phaeophytin by fluorescence. Deep-Sea Res. 10: 221-231

Dissolved Oxygen

DISSOLVED OXYGEN

 


SUMMARY: The amount of dissolved oxygen in seawater is measured using the Carpenter modification of the Winkler method. Carpenters modification (1965) was designed to increase the accuracy of the original method devised by Winkler in 1889. Using Carpenters modification, the significant loss of iodine is reduced and air oxidation of iodide is minimized. Rather than using the visible color of the iodine-starch complex as an indicator of the titration end-point, we use an automated titrator that measures the absorption of ultraviolet light by the tri-iodide ion, which is centered at a wavelength of 350 nm.


 

1. Principle

Manganous chloride solution is added to a known quantity of seawater and is immediately followed by the addition of sodium hydroxide iodide solution.  Manganous hydroxide is oxidized by the dissolved oxygen in the seawater sample and precipitates forming hydrated tetravalent oxides of manganese.

Mn+2 + 2OH ———————————————-> Mn(OH)2 (solid)
2Mn(OH)2 + O2 ——————————————–> 2MnO(OH)2 (solid)

Upon acidification of the sample, the manganese hydroxides dissolve and the tetravalent manganese in MnO(OH)2 acts as an oxidizing agent, setting free iodine from the iodide ions.

Mn(OH)2 + 2H+ ——————————————-> Mn+2 + 2H2O
MnO(OH)2 +4H+ +2I ————————————> Mn+2 +I2 +3H2O

The liberated iodine, equivalent to the dissolved oxygen present in the sample, is then titrated with a standardized sodium thiosulfate solution and the dissolved oxygen present in the sample is calculated.  The reaction is as follows:

I2 + 2S2O3 ————————————————–> 2I + S4O6

2. Reagent Preparation

2.1.  The manganous chloride solution (3M) is prepared by dissolving 600g of reagent grade manganous chloride tetrahydrate, MnCl2•4H20, in Milli-Q water to a final volume of 1 liter.  This solution is then filtered using 47mm glass fiber filters (Whatman GF/F).

2.2  The sodium hydroxide (8M)-sodium iodide (4M) solution is prepared by first dissolving 600g sodium iodide (NaI) in approximately 600ml Milli-Q water.    After the NaI is dissolved, 320g of NaOH is added slowly (caution-the solution will get hot) and the volume is adjusted to 1 liter with Milli-Q.  The solution is then filtered through a GF/F.

2.3.  Sulfuric acid solution (10N) is prepared by slowly adding 280ml of reagent grade concentrated sulfuric acid, H2SO4, to 770ml of Milli-Q water.   This should be prepared with caution as is gets very hot.

2.4.     The sodium thiosulfate solution (0.2N) is prepared by dissolving 50g sodium thiosulfate pentahydrate (Na2S2O3.5H20) and 0.1g anhydrous sodium carbonate (Na2CO3) in Milli-Q water to a final volume of 1 liter.  This solution is prepared approximately 2 weeks before use and stored in an amber glass bottles.

2.5.  The potassium iodate standard (0.0100N) is prepared by first drying potassium iodate(KIO3) in a drying oven for approximately one hour.  Once the KIO3 is dried, carefully measure out 0.3567g, using a 5-place balance, and dissolve in Milli-Q water to a final volume of 1 liter.

3. Sample Drawing

3.1.     Oxygen samples are always drawn first from the Niskin bottles and should be drawn as soon as possible to avoid contamination from atmospheric oxygen.  Approximately 6 inches of Tygon tubing connected to a temperature probe via a y-connector is slipped onto the discharge valve of the Niskin.

3.2.  A calibrated volumetric flask is rinsed three times, then with the seawater still flowing, the end of the tygon tubing is placed into the flask nearing the bottom.  The sample is then overflowed with twice the sample volume while making sure that there are no bubbles in the tubing during the overflow process.  The tubing is then carefully removed from the sample flask to prevent the influx of bubbles.

3.3. Immediately after drawing the sample, 1 ml of manganous chloride solution is added into the flask.  This is followed by the addition of 1ml of sodium hydroxide-sodium iodide solution to the sample.  Both dispensers should be purged to remove air bubbles prior to the addition of these reagents.

3.4.  The stopper is then carefully placed in the bottle to avoid the trapping of air and the temperature of the seawater at the time of sample draw is recorded.

3.5.  After all samples are drawn, they are shaken vigorously to disperse the precipitate uniformly through the flask.  This process is repeated again after the precipitate has settled to the bottom of the flask, or after at least ten minutes.
 

4. Standardization of thiosulfate

4.1.  Proper care in the setup of the auto-titrator is required before running blanks, standards and samples.  The UV lamp is turned on at least 30 minutes prior to the run and should have a stable voltage of 2.4-2.5 volts for the run.  The dosimat tubing is carefully purged so that the lines are completely free of air bubbles, and the water bath inside the auto-titrator is clean and filled with Milli-Q water prior to analysis.  Make sure, prior to your first run, any thiosulfide solution is rinsed off the tip of the line with Milli-Q after purging.

4.2.  Using a Metrohm 655 Dosimat, dispense 10 ml of the standard potassium iodate solution into a clean oxygen flask.  Add a stir bar and rinse down the sides with a small amount of Milli-Q water.  Add 1 ml of the 10N sulfuric acid solution and swirl to ensure the solution is well mixed before adding the pickling reagents.

4.3.  Add 1 ml sodium hydroxide-sodium iodide solution to the acidified flask and swirl gently.  Then add 1ml of Manganous chloride solution, swirl gently, and fill the solution to the neck of the flask with Milli-Q water. 

4.4.  The UV detector on the auto-titrator measures the transmission of ultra-violet light through the standard (as well as seawater sample and blank) as a Metrohm 665 Dosimat dispenses thiosulfate at increasingly slower rates.  The endpoint is reached when no further change in absorption is detected by the detector.  At this point all of the iodine has been consumed.
 

5. Blank determination

5.1.  Using the Metrohm 655 Dosimat, dispense 1 ml of the standard Potassium Iodate solution into a clean oxygen flask.  Add a stir bar and rinse down the sides with a small amount of Milli-Q water.  Add 1 ml of the 10N sulfuric acid solution and swirl to ensure the solution is well mixed before adding the pickling reagents.

5.2.  Add 1 ml sodium hydroxide-sodium iodide solution to the acicified flask and swirl gently.  Then add 1ml of Manganous chloride solution, swirl gently, and fill the solution to the neck of the flask with Milli-Q water.  The solution is then titrated to the end-point as described in section 4.3 above.

5.3.  A second 1 ml aliquot is added to the same solution which is then titrated to a second end-point.  The difference between the first and second titration is used as the reagent blank.
 

6. Sample analysis

6.1.  Samples are analyzed after all of the precipitate settles to the bottom of the flask, after the second shake.  The top of the flask is wiped with a kimwipe to remove moisture containing excess reagent around the stopper and then the stopper is carefully removed.

6.2.  1ml of 10N sulfuric acid is added to the sample and a stir bar is placed inside the flask.  The flask is then secured inside the clean water bath.  The tip of the thiosulfate dispenser is placed inside the sample flask and the automated titration can begin with the use of an auto-titrator program.

 

7. Certified Standard Comparison

7.1 Presented here are the results of a comparison of several certified standards with a CalCOFI prepared standard (weighed and diluted up to set normality in 1 liter volumetric)standard.   Three certified iodate solutions, (Table 1) were tested using the same 0.2N thiosulphate solution.   Iodate concentrations are back calculated from defined thiosulphate concentrations for comparison purposes.   It is noteworthy that the Acculute and Fisher solutions had to be diluted before use.  The same volumetric was used for all dilutions, the same Dosimat and piston was used for all titrations, extensive rinsing between uses.  Titration n=9 as dictated by the maximum number of samples that could be acquired from the 100ml portion for the OSIL standard.  CalCOFI values are a result of cruise and shore based titrations to demonstrate an integrated real use sampling.  Differences between all standards represent less than   +/-  0.5 percent of the signal. 


 


Average[KIO3]


mM


 


STDEV


mM


 


% Diff from CC


       mM


CalCOFI iodate 12/16/2012 used on 1301


 

 

 


10.015722


 


0.009609


 


        n.a.


 

 

 

 

 

 


OSIL iodate oxygen standard


 

 

 

 


9.972671


 


0.013310


 


99.57


 

 

 

 

 

 


Acculute iodate Lot 00702


 

 

 

 


10.049259


 


0.007684


 


100.33


 

 

 

 

 

 


Fisher iodate Lot 125009


 

 

 

 


10.044804


 


0.007650


 


100.29

 7.2 These results are consistent with previously published methods comparison for precision of oxygen measurements1 (WOCE Report 73/91) and previous replicate analysis to verify precision with auto-titration methods2.  Although the Ocean Data Facility (ODF) performed the 1991 comparisons for Scripps and was different than CalCOFI, the Winkler techniques used then and now are the same.  Presently, CalCOFI methods employ a UV end point auto-titrator made by ODF that produces a precision of 0.005-.01 ml/L.  The difference between high and low samples in replicate analysis equals ~0.5% with a standard deviation typically under 0.010 ml/L.  See references:

  1. oceaninformatics.ucsd.edu/calcofi/docs/CCConf06Wolgast.ppt
  2. odf.ucsd.edu/index.php?id=56

8. Calculation and Expression of Results

8.1.  The auto-titrator uses a UV detector that detects changes in voltage as thiosulfate is added to the sample.  The volume of thiosulfate added is recorded at an endpoint once there is no change in voltage.  The end point is determined by a least squares fit using a group of data points just prior to the end point, where the slope of the titration curve is steep, and a group of data points just after the endpoint, where the slope of the curve is close to zero.  The intersection of the two lines is taken as the endpoint.

8.2.  The calculation of dissolved oxygen follows the same principle oulined by Carpenter (1965).  Our results are expressed in mL/L.
 

                       (R-Rblk)VIO3*NIO3*E         DOreg

O2(ml/L)=      —————————  –  ———

                       (Rstd-Rblk)*(Vb-Vreg)         Vb

Where:
R= Sample titration (mL)
Rblk = Blank value (mL)
VIO3 = Volume of KIO3 standard (mL)
NIO3 – Normality of KIO3 standard
E = 5,598 mL O2/equivalent
Rstd = Volume used to titrate standard
Vb – Volume of sample bottle
Vreg = Volume of reagents
DOreg = Oxygen added in reagents
 

9. Equipment/Supplies

  • Volumetrically calibrated 100ml Glass Erlenmeyer flasks with paired ground glass stoppers
  • 3 – 1ml Brinkman reagent dispensers
  • Tygon tubing
  • Fisher Scientific UV Longwave Pencil Lamp, 365 nm and power supply, 115 VAC
  • 2 Metrohm Dosimat 665 Automatic Burets
  • 10ml Metrohm Dosimat Exchange unit
  • 1ml Metrohm Dosimat Exchange unit
  • Metrohm Dosimat Keypad
  • Spare 1ml and 10ml dispensor pistons for the Metrohm Dosimat Exchange units
  • Scripps/STS Auto-titrator Unit and Software
  • PC Computer
  • Waterproof sampling thermometer
  • Concentrated Sulfuric Acid, H2SO4, ACS Grade
  • Manganous Chloride Tetrahydrate, MnCl2•4H20, ACS Grade
  • Sodium Hydroxide, NaOH, ACS Grade
  • Sodium Thiosulfate, Na2S2O3•5H20, ACS Grade
  • Potassium Iodate, Dry high purity KIO3, Alfa Aesar
  • Granular Sodium Iodide, EMD Chemicals via VWR
  • Magnetic Stir Bars

 

10. References

  • Anderson, G. C., compiler, 1971. “Oxygen Analysis,” Marine Technician’s Handbook, SIO Ref. No. 71-8, Sea Grant Pub. No. 9.
  • Carpenter, J. H., 1965. The Chesapeake Bay Institute technique for the Winkler dissolved oxygen method.  Limnol. Oceanogr., 10: 141-143.
  • Culberson, C. H. 1991. Dissolved oxygen. WHP Operations and Methods — July 1991.
  • Parsons, T. R., Y. Maita, C. M. Lalli, 1984. A Manual of Chemical and Biological Methods for Seawater Analysis. Pergamon Press Ltd., 3-28.
  • oceaninformatics.ucsd.edu/calcofi/docs/CCConf06Wolgast.ppt

Salinity

Salinity Determination


OVERVIEW: This procedure describes the method for the determination of seawater salinity using a Guildline Portasal™ Salinometer (Model 8410).
 

1. Principle

During a cruise a 24 bottle rosette is used to collect seawater salinity samples of ~225 mL from all bottles closed from pre-designated depths.  A Guildline Instruments Portasal™ Salinometer (8410A) makes the precise conductivity comparisons between the water samples and a reference water standard.  From these comparisons salinities are then calculated and logged using PC based software that averages data that meet replicate criteria. Concurrent with the water sampling, a Sea-Bird Electronics CTD (SBE 911 plus) profiles in situ data.  Data processing software is used to compare the bottle salts to in situ CTD measurements.  This is useful to confirm bottle closure, monitor CTD performance, and to select the best salinity data for future comparison.  

2. Sample Bottles

Salinity bottles are collected in ~250 ml borosilicate KIMAX-35 bottles.  These square cross sectioned bottles have screw tops and lids with separate plastic thimbles to prevent leakage and evaporation. Numbered stackable divider boxes hold 24 numbered bottles to match the 24 Niskins on the rosette water sampler. Bottles are filled with previously unused sample water to limit salt crystallization and to pre-leech silicate from the glass.  

3. Sample Collection

Numbered salt bottles are drawn from corresponding Niskin bottles and filled to the shoulder after three ~40 ml rinses. The last fill is done without interruption until overflowing; then ~10 ml is poured out over the thimble and the bottle is stoppered and capped. Samples are left to equilibrate to room temperate for >8 hours prior to measurement. Before being placed in the salinometer samples are: gently inverted 3 times to remove any possible stratification, wiped around the bottom of the cap to draw out as much water as possible that may be trapped under and around the cap, wiped around the thimble and threaded neck once the cap is removed to eliminate any excess water or salt precipitate, and flushed by dumping out ~10 ml into a collection bucket.  

4. Guildline Portasal™ Salinometer (8410A)

Guildline Portasal Salinometer model 8410A is used to make the precise conductivity comparisons between the water samples and reference water standards.  Two of these instruments are taken on each cruise.  Guildline specifications state an accuracy of  ±0.003 PSU (same set point temperature as standization and within -2°C and +4°C of ambient), and a precision of 0.0003 PSU.  

3. configuring the PORTASAL™

Bath temperature must be within 2 degrees below and 4 degrees above ambient. To check bath temperature set point press the TSET key. Compare set point to actual by pressing ENTER for actual bath temperature. Set point and actual temperature must be within 0.02°C. Press UP ARROW key to view temperatures of the dual bath thermistors (TH1 and TH2). TH1 and TH2 must agree to within 0.04°C.  

The Portasal must be powered up for >3 hours to ensure bath temperature regulation has begun before Reference values can be calibrated. While the the FUNCTION switch is in STDBY pressing REF key will display alternating Reference readings plus (+), minus (-), and Reference readings. After several cycles, when + and – values are within 1 unit of each other, pressing the COND key will initiate salinometer self-calibration.  

With the FUNCTION switch set to ZERO, pressing the COND key will initiate the ZERO calibration process. When satisfied that the displayed value is stable, press the ZERO key. When the subsequent displayed number is stable (but not necessarily zero) press the COND key. The display should then read 0.00000. Setting the FUNCTION switch back to STDBY will ready the Portasal for Standardization.  

4. SubStandard Preparation

Deep sea water is collected on cruises from bottles tripped at depths greater than 300m in 10 liter polyethtylene jerry jugs. The seawater is then transfered to 50 liter carboys at the lab for filtration. During the transfer 10 ml hypochlorate (bleach) is added to inhibit biological activity.   Concentrated high salinity seawater is prepared by evaporating deep sea water (>300m) by heating the seawater at about 90°C in an oven overnight. A volume of 1800 ml will evaporate down to about 1000 ml in 24 hours, final volume can be adjusted with DI water, this will give a salinity of around 60psu. Filter the concentrate through GF/F filter before using it for adjusting the sub-standard. Filter the deep seawater through GF/F to 46 liter mark. Add 10 ml laundry bleach to inhibit biological activity. Measure the salinity of the filtered water. Adjustment of the substandard salinity is usually done in two steps. For the first adjustment, add concentrate to the filtered sea water to adjust to the desired salinity using the relationship:

(vol. sw)(sal. sw) + (vol. conc.)(sal. conc.) = (vol. sw + vol. conc.)( final sal. sub)

The salinity of the concentrate will be too high to measure on the salinometer. If you have prepared the concentrate as above you can do your initial adjustment using an estimated salinity of 60 for the concentrate. The final salinity of the sub should be close to your standard value ± 0.015 PSU. For the initial adjustment it works well to add about 90% of the concentrate volume calculated above. Be sure you have thoroughly mixed the solution by stirring gently. A motorized paddle wheel is used to thoroughly stir the mixture. Measure the salinity of the sub after the first adjustment. Using the equation above you can now calculate the apparent salinity of the concentrate, by using the salinity of the concentrate as an unknown, along with the known values for the initial salinity of the seawater, and the adjusted value of the substandard as just measured. Applying the calculated salinity value for the concentrate, use the equation again to calculate a new volume of concentrate to be added to get to the final desired salinity of the substandard. If your initial concentrate addition was to large you may have to dilute with DI water using zero for salinity concentrate in the equation. When you have arrived at the desired salinity for the sub-standard add a layer dimethylpolysiloxane to prevent evaporation resulting in salinity changes. 
 

6. IAPSO Standard Seawater

Cut and Paste from WHP Operations and Methods July 1991 

Standard Seawater prepared by Ocean Scienti c International Ltd. (OSI) is the recognized standard for the calibration of instruments measuring conductivity (salinity). This water is natural surface water which has been collected in the North Atlantic and carefully ltered and diluted with distilled water to yield seawater with a conductivity ratio near unity and a salinity near 35. The seawater is sealed in glass vials and labeled with the date, batch number, K15 value and chlorinity. Mantyla (1987) has compared batches of this seawater (P-29 to P-102) and found inter-batch di erences as great as about 0.003 in some of the older batches. Since batch P- 93 however, inter-batch di erences have not exceeded about 0.001.

It is recommended that a single batch of SSW be used during each cruise and that it be identi ed in the cruise report. It is the responsibility of the salinity analyst and chief scientist to ensure that the quality of the salinity observations are the highest attainable. It follows, therefore, that the inter-batch di fferences described by Mantyla (1987) should be used to correct the nal salinities before they are reported. It has also been noted that the conductivity of some Standard Seawater changes with time. Vials more than 4-5 years old should be compared with fresher Standards to determine possible changes in conductivity due to aging.

7. Wolga Water

In addition to using IAPSO SSW, a batch of aforementioned substandard water has been bottled and sealed in glass vials.  This is done to decrease dependence on IAPSO standard while maintaining acceptable standardization accuracies.

8. Standardizing the Portasal with substandard

Before standardization can begin the conductivity cell must be flushed rigorously as it had been filled with detergent or 50% ethanol between uses. To flush, attach the peristaltic pump to the Portasal sampling tube. Turn on both the Portasal and peristaltic pumps, and open the valve on the on the substandard feed line. Flush repeatedly (~10 flushes) until the cell is free of contamination. It is convenient to do this during Portasal configuration. Begin the salt program on the salt computer, and complete requisite data fields.

With the conductivity cell full of substandard and a stable value being displayed turn the FUNCTION switch to READ. When displayed value is stable press STD key. When the Portasal display reads STD STANDARDIZE press ENTER key. The display will show the conductivity of the standard being used. Press ENTER key to proceed (if values of substandard need to be adjusted use ARROW keys.) The BATCH # will be displayed and should match what is being used. If values need to be changed modify with the ARROW keys, pressing the ENTER key will ready Portasal for standardization. The display will read ENTER WHEN READY, press ENTER key. The Portasal will begin measurement and display a substandard value. When stable press the COND key. This will terminate the standardization and the displayed conductivity ratio should match that of your sub-standard within ± 0.00001. Pressing the ENTER key on the salt computer will save it to file.

At the beginning of each sample run the Portasal is standardized with substandard.  However, every other day IAPSO and Wolga standards are also sampled.  All three standards are run before and after the sample run in order to ascertain instrument drift.  

9. Salt Data ACQUISITION Program

The SIO-CalCOFI-authored conductivity recording software, PSal.exe, is a Windows-based data acquisition program that records conductivity values from discreet salinity samples. Averaging five or more conductivity readings from the Guildline Portasal, the operator saves a stable reading. The flow cell is flushed and refilled and additional readings are performed. A pair of conductivity values which agree within 0.0010 standard deviation are recorded and salinity calculated. The software compares this calculated bottle measurement to the matching CTD salinity when available. This comparison is a good indicator of sensor performance and bottle sampling accuracy.
The data are saved as salt run files which may contain more than one station. Keeping the salt runs combined is necessary to calculate the drift-over-time calculated from any change in the end standard or substandard reading. The data are also save in a single, combined, database-friendly csv for database processing (in development).

10. CTD Equipment and Data Processing

A SBE 9ll+ CTD is equipped with dual conductivity and temperature sensors.  These are routinely calibrated (~6 months), but because of superior stability the pressure sensor in calibrated much less frequently (~24 months).

The SBE 11plus Deck Unit applies a real time alignment correction to conductivity during data collection, and the accompanying software generates a marker file for each profile containing CTD data at the instant of each Niskin closure.
Seasoft-processed CTD salinities are processed and imported into station csvs for comparison to bottle salinities.

11. DECODR: Data Entry Compiler & Output Data Reports

CalCOFI-authored Windows program will process the bottle sample conductivity data files generated by the SalReCap.exe data acquisition program. Using beginning and end standards, a linear drift is applied to each conductivity pair before salinity (PSU) is calculated. DECODR combines all stations selected into a single text file for data-quality assessment. When CTD salinities for matching depths are available, DECODR compares the bottle & CTD salinities, calculating & reporting the differences. DECODR, as an option, will save the bottle salinity values into each station’s csv. DECODR can generate data products using bottle or CTD data such as IEH or data reports.