Instructions for Monthly/Seasonal Correlations

Variable, Years and Levels Available TopGo to Page Top

Data is available from Jan 1948 to the present for most variables and is updated the first or 2nd week of the month with the previous month's data. Sigma level variables are available only until Apr 1998, ONLY pressure and air temperature are available for the tropopause level, and ONLY u and v winds, relative humidity and pressure for the surface (.995 sigma). Streamfunction, divergence, vorticity and velocity potential are available ONLY on sigma levels.

Sea Surface Temperature: This is from the NOAA OI SST data (prescribed) over the ocean. Over land it is the skin temperature (prescribed as well). Sea Ice (which varies near land) affects these values and so the plotting range may not be as great over the tropics as desired. Using a lower latitude or specifying the contour range should help.

Index time series have cover a range of years. Some go back as far as 1948 but some don't. Please check the time series for actual data.

How Correlations are Calculated TopGo to Page Top

For seasonal mean correlations, a value is calculated for each year for both the index time series and the variable that is being correlated with. A correlation is then calculated from these seasonal means For monthly correlations, each individual monthly value is used. The seasonal correlation for one month season is identical to the monthly correlation.

The formula for the correlation is...

r=sum(xy)/sqrt[sum(x**2)*(y**2)]

...where x and y have the time means removed. There ARE other ways of calculating correlations.

The time indexes are not defined for some seasons (e.g. PNA in summer). If you choose that season for the correlation, you will get back a blank plot. Likewise, any missing values for the index will not be used the the calculation. Please look at the time series before interpreting the results.

Year Range TopGo to Page Top

Default years are 1948-2011 with the exception of seasons that would begin in 1947 (e.g. DJF) and variables that do not extend over the whole time period. For the former, 1949-2011 is used and for the latter the ending year varies. You can use your own range of years up to the latest possible data. Note that there is no data before January 1948. You should get an error message if you ask for data that is not there. For example, the season DJF 1948 contains December 1947 and therefore will not be calculated. The teleconnection values and the monthly means are usually updated within the next month.

Lead/Lag TopGo to Page Top

The default lead/lag is 0 months. You can specify that the index time series either lead (come before) or lag (come after) the variable being correlated with. You can specify the correlations in months only. So, if you are doing seasonal correlations of FMA, a one month lead would correspond to the index time series of JFM. Please do not specify a lead of more than 12 months. The months used are indicated on the plot.

Contour Options TopGo to Page Top

A desired contour interval and range can be input instead of the default being used. Different plots can be easily compared (and the resulting gifs could be animated). For this option to work, the interval AND the range must be input. There must be at least 2 and less than 33 contours. The contour interval must be positive and the range must go from low to high.

Black and White Option: Black and white contours and shading are available for use on non-color printers.

Atmospheric and Oceanic Indices TopGo to Page Top

Monthly SST and atmospheric indices (e.g. PNA) are obtained from NOAA CPC except where noted below. Please see their site for details on the time series and how they were created. Calculation methods for other time series can be found by following links in the right column. Please reference time-series use in publications!

See Climate Indices Plotting Page for a description and references for the different time series.

Random Calculation: Calculation Details

The time series you get is a randomly calculated "red noise" time series. In this case, the solution is calculated from the integrated stochastic differential equation

dx/dt = (-1/T)x + whitenoise

1/T is set to -1/3 months or -.333 for all cases. x(0) is determined using a changing seed value and the white noise is obtained by sampling from a gaussian distribution. Values are set for all 12 months of 1958-1999 and seasonal values are calculated from that. Code was graciously provided by Cecile Penland of NOAA/PSL. Other definitions of a random time series could have been used.

The intention is to demonstrate how spatial patterns that look "reasonable" will arise during correlation analysis even if the index time series has no physical basis. Please do not try to do Monte Carlo tests by repeatedly running the web page.

Rotating Plots TopGo to Page Top

To rotate map, choose "custom" map projection and then choose either northern or southern polar stereographic projection. For lat range, enter 0 to 90 (for northern hemisphere) or -90 to 0 for southern. For longitude, the center of the longitude range input will be at the bottom of the plot. To center along 0E, choose -180 to 180, for example. To center at 90E, choose -90 to 270. You can plot sectors as well. The longitudes -90 to 90 in the NH will plot the half hemisphere from the US across the Atlantic to Europe.

Instructions for Crossections TopGo to Page Top

  1. Choose a variable with multiple levels (temperature, geopotential heights, omega, humidity, wind, or potential temp).
  2. For 'plot region' choose latitude by height or longitude by height.
  3. For latitude by height plots, choose a latitudinal extent and a range of longitudes (or a single longitude) over which to average.
  4. For longitude by height plots, choose a longitudinal extent and a range of latitudes (or a single latitude) over which to average.
  5. Choose level range (bottom to top). Note that humidity variables only go through 300mb and omega only goes to 100mb.

Plotting Regions TopGo to Page Top

To plot over the dateline, use values from 0 to 720. For example, to plot 180W eastward to 180W, use 180 to 540. Be sure that the western most longitude is less than the eastern., For example, to plot 100W to 70W, use -100 to -70 or 260 to 290 and NOT 100 to 70.

There are 6 custom projections:

Northern Hemisphere: 0-90N, 0-360W in using a polar stereographic projection
Globe: 90S-90N 0E-360W
United States: 20N-65N; 235-285 polar stereographic projection
Tropics: 60S-60N
Tropical Pacific: 35N-35S 100E to 60W
Indio-Pacific (20S-20N,60E-160W)
4 Corner Region US (xxS-xxN,xxE-xxW)

How this Page was Created TopGo to Page Top

The main interface for this page is an HTML form. The data that are input into this form are processed by a Perl script. The script reads the inputs, tests for bad inputs and then executes a FORTRAN code that produces a correlation file (in NetCDF). This file is what is processed by a GrADS script. The GrADS script is run as a batch job with the plot created in a virtual frame buffer (to increase the quality of the GIF output). Plot options are input into the GrADS script at this time. The plot in the frame buffer is converted to GIF and this GIF is displayed as part of a HTML document. The netCDF file and the gif file are kept in a directory where the files are periodically deleted.

Note: Leads and Lags >13 were allowed by the page but not in the instructions. Leads/lags >12 returned correlations a year off. You can now do lags/leads ≤ 24.