roy troll border art
nefsc banner
Technical Memoranda Reference Documents Classic Publications Contract/Grant Reports
CMER Publications Series Information Links and Contacts Annual/Biennial Lists
Web Manager Email Search Publications Publications Home Site Map
CONTENTS
Introduction
Changes in the Reference Annual Cycles
Changes in the Calculation of Anomalies
Changes in the Calculation of Regional Average Values
Data Tables
References
Northeast Fisheries Science Center Reference Document 04-08

Revised Procedures for Calculating Regional Average Water Properties for Northeast Fisheries Science Center Cruises

by David G. Mountain, Maureen H. Taylor and Cristina Bascuñán

National Marine Fisheries Serv., Woods Hole Lab., 166 Water St., Woods Hole, MA 02543

Print publication date April 2004; web version posted July 28, 2004

Citation: Mountain, D.G.; Taylor, M.H.; Bascuñán, C. 2004. Revised procedures for calculating regional average water properties for Northeast Fisheries Science Center cruises. Northeast Fish. Sci. Cent. Ref. Doc. 04-08; 53 p.

get acrobat reader Download complete PDF/print version


Introduction

Temperature and salinity measurements are made routinely on most Northeast Fisheries Science Center (NEFSC) survey and research cruises.  The area of coverage extends from Cape Hatteras northeastward to the Gulf of Maine, and from near the coast to just seaward of the shelf break.  Since 1991 a report of the water properties has been issued almost every year (e.g., Taylor and Bascuñán, 2003).  This report includes regional average values of surface and bottom water temperature and salinity, and of the anomalies of temperature and salinity for five regions of the shelf (Figure 1).  The anomalies are determined relative to reference annual cycles derived from data collected during the NEFSC Marine Resources Monitoring and Assessment Program (MARMAP) program.  Holzwarth and Mountain (1992) present similar regional values for surface and bottom temperature as measured on the NEFSC spring and fall bottom trawl survey cruises from 1963 to 1990.

The methods used to determine the MARMAP annual cycles are described by Mountain and Holzwarth (1989).  The methods used to calculate the anomalies, the regional average values, and the uncertainty in the regional average values are described by Holzwarth and Mountain (1992).  This report documents changes that have been made to the reference annual cycles and to the methods used to calculate the regional average values.  The revised cycles and methods will be used in the annual summary of water properties beginning with the report for 2003 (Bascuñán and Taylor, 2004).  All previously reported regional average values have been recalculated using the revised annual cycles and methods, and are included in this report.  In most cases differences between the original and the revised regional values are small (< 0.1°C or < 0.1 PSU).


Changes in the Reference Annual Cycles

The NEFSC MARMAP program (1977-1987) made observations of plankton and water properties at approximately 160 standard station locations 4-8 times per year.  Temperature and salinity measurements were made using water bottles with reversing thermometers.  The water samples collected were analyzed for salinity by a laboratory salinometer.  During the last year of the program an electronic conductivity/temperature/depth (CTD) instrument was used.

Because the MARMAP measurements were made at fixed locations over a 10 year period, characteristic annual cycles of the surface and bottom temperature and salinity could be determined for each of the locations.  The annual cycles were derived by sequentially fitting up to three harmonic cycles (periods of one, one-half, and one-third year) to the observations at a location (Mountain and Holzwarth, 1989).  For each harmonic a significance test was performed to determine if the fit was statistically significant.  In some cases no significant harmonic cycle was found and the annual cycle consisted only of the mean of the observations.

In 2002 an ORACLE database was established for NEFSC’s hydrographic data.  As part of that effort the historic hydrographic data was reviewed for quality before it was entered into the database.  This quality control process resulted in some MARMAP data either being changed or deleted from the data set.  Because of these changes the annual cycles based on the MARMAP data were recalculated using the data in the ORACLE database.  In determining the bottom annual cycles, observations had to be no more than 10m or 10% of the water depth (which ever was greater) from the bottom to be included in the calculations as representing a bottom value.  Generally, only minor changes resulted in the annual cycles.  In a few cases no significant annual cycle could be determined where previously one had been, and in a few other cases a cycle was determined where previously one had not been.


Changes in the Calculation of Anomalies

To calculate the anomaly for an observation, the MARMAP stations closest to the observation location are determined.  The expected value of the property (i.e., surface or bottom temperature or salinity) at each of these MARMAP stations is determined from the appropriate annual cycle for the calendar day on which the observation was made.  The MARMAP values are averaged by a weighting inversely proportional to the square of the distance each is from the observation location.  The anomaly is the difference between the observed value and the resulting averaged MARMAP value.

For calculating bottom anomalies, only MARMAP stations with bottom depths similar to the bottom depth of the observation station are included.  The criterion used is that the depth difference has to be smaller than 25m or 25% of the bottom depth of the observation station, whichever is greater.

When calculating salinity anomalies, MARMAP stations that do not have a significant annual cycle are included by using the annual mean value for that station (i.e., the C1 value in Table 1).  For calculating temperature anomalies, MARMAP stations without an annual cycle previously had not been included.  Under the new procedures MARMAP stations without an annual cycle are included in calculating temperature anomalies by using the mean value for that station.  The only MARMAP stations that do not have an annual cycle for bottom temperature are in the deep Gulf of Maine or seaward of the shelf break (>100m depth).  Therefore, this change in the procedure for calculating bottom temperature anomalies only affects observations in those areas.


Changes in the Calculation of Regional Average Values

The regional average values are determined by interpolating the observed values to a fine scale grid, with each grid point representing a known area.  Previously the separation between grid points was 0.2 degrees latitude and 0.25 degrees longitude.  In the revised procedures the separation has been reduced to 0.1 degree in latitude and 0.1 degree in longitude to provide a more accurate representation of an area.  In addition the computer code used to determine if a station was within the boundaries of a region has been modified to correct for an error that in a few, specific situations identified stations close to but south of a regional boundary as being within the region.


Data Tables

Annual Cycles:

The MARMAP annual cycles are calculated as the sum of three harmonics.  For calendar day ‘d’, the value of the parameter ‘V’ is calculated by:

V(d) =  3 Σ Ci + Ai * cos(ri * d) + Bi * cos(ri * d)
             i=1

where  ri = 2*π (i/365)  and the C’s, A’s and B’s are the harmonic coefficients.

Table 1, Table 2, Table 3, and Table 4 contain the coefficients for surface temperature, bottom temperature, surface salinity and bottom salinity, respectively.  The columns in the tables are:

MARMAP standard station number
Latitude
Longitude
Bottom depth
The harmonic coefficients C1, A1, B1, C2, A2, B2, C3, A3, B3
The standard deviation of the original data values from the fitted curve

Regional Property Values:

For each cruise in the NEFSC hydrographic data base Table 5, Table 6, Table 7, Table 8, and Table 9 contain the revised area average surface and bottom temperature and temperature anomaly values for the five regions in Figure 1Table 10, Table 11, Table 12, Table 13, and Table 14 contain the same information for salinity.  In these tables the columns are:

Cruise name
Region  Code ()
Year
Calendar day
Decimal-Year (year + calendar day/365)
For the surface:
  Number of stations in the region
  Value
  Anomaly
  SDV1
  SDV2
  Flagg
For the bottom:
  Number of stations in the region
  Value
  Anomaly
  SDV1
  SDV2
  Flagg
Cruise Purpose Code

The Region Codes are:

1            MAB South
2            MAB North
3            Georges Bank
4            Gulf of Maine West
5            Gulf of Maine East

SDV1 indicates the uncertainty in the regional average anomaly value associated with the uncertainties in the MARMAP annual cycles (See Holzwarth and Mountain, 1990).  SDV2 is the standard deviation of the anomaly values for the stations within the region. A Flagg value of 0 indicates that there were enough observations to calculate an area weighted value for the region.  A Flagg value of 1 indicates that the spatial distribution of the observations was not sufficient to calculate an area average for the region and instead a simple average of the observations was calculated.

The cruise purpose codes are:

  01            BLUEFISH SURVEY
  02            LOBSTER SURVEY
  03            COD SPAWNING SURVEY
  10            NMFS NEFSC BOTTOM TRAWL SURVEY
  16            NMFS ACOUSTICS SURVEY
  20            GLOBEC BROADSCALE SURVEY
  21            GLOBEC PROCESS STUDY
  22            ECOSYSTEM MONITORING
  23            MARMAP
  24            LARVAL HERRING SURVEY
  25            12-MILE DUMPSITE
  27            WARM CORE RING CRUISE
  50            NMFS NEFSC CLAM AND QUAHOG SURVEY
  60            NMFS NEFSC SEA SCALLOP SURVEY
  70            NMFS NEFSC GEAR COMPARISON
  80            FOOD CHAIN DYNAMICS FEEDING ECOLOGY SURVEY
  81            COASTAL OCEAN PROGRAM
  82            APEX PREDATOR SURVEY
  91            BENTHIC HABITAT STUDY
  90            MARINE MAMMAL SURVEY
  93            MISCELLANEOUS NON-RANDOM RESOURCE INVESTIGATION CRUISE

Regional average surface and bottom temperature for the spring and fall bottom trawl surveys from 1963 to 1990 are listed in Table 15 and Table 16.  The columns are the same as those for Tables 5-14 except that cruise name column is omitted.


References

Bascuñán, C. and M.H. Taylor  2004.  Description of  the 2003 oceanographic conditions on the northeast contental shelf.  Northeast Fisheries Science Center Reference Doc.  04-xx;  xx p.  Available from:  National Marine Fisheries Service, 166 Water St., Woods Hole, MA  02543.

Holzwarth, T.J. and D.G. Mountain  1990.   Surface and bottom temperaure distributions from the Northeast Fisheries Center spring and fall bottom trawl survey program, 1963-1987, with addemdum for 1988-1990.  Northeast Fisheries Science Center Reference Doc. 90-03;  77 p.  Available from:  National Marine Fisheries Service, 166 Water St., Woods Hole, MA  02543.

Mountain, D.G. and T.J. Holzwarth.  1989.  Surface and bottom temperature distribution for the northeast continental shelf.  NOAA Tech. Mem. NMFS-F/NEC-73; 32 p.

Taylor, M.H. and C. Bascuñán.  2001.  Description of  the 2000 oceanographic conditions on the northeast contental shelf.  Northeast Fisheries Science Center Reference Doc. 01-01;  93 p.  Available from:  National Marine Fisheries Service, 166 Water St., Woods Hole, MA  02543.