CVDP logo

Go to:

Diagnostics Plots

Methodology and Definitions

General Notes

  • ts (CESM equivalent = ts) is surface ("skin") temperature and is used in lieu of sea surface temperatures (SSTs), tas (TREFHT) is 2m air temperature.
  • pr (PRECC+PRECL or PRECT) is total precipitation , psl (PSL) is sea level pressure, sic (aice) is sea ice concentration, msftmz (MOC) is the ocean meridional overturning mass streamfunction, and zos (SSH) is the sea surface height above geoid.
  • The annual cycle is removed prior to every calculation by subtracting the long-term monthly means. Exception: The annual cycle is not removed for mean spatial maps.
  • The long-term trend removal is performed after the annual cycle is removed (if necessary) and set by the user for 1) data specified in the input namelist_obs file (usually observational data) and 2) data specified in the namelist file (usually model data). Description of method used in this comparison for specified namelist_obs data: A linear trend is removed by month. Description of method used in this comparison for specified namelist data: A linear trend is removed by month. The trend removal method is performed for every metric other than Climatological Averages and Linear Trends.
  • Area-averages are always based on cosine of latitude weighting.
  • For visual clarity, the Y-axis may differ amongst individual panels on a particular plot.
  • Power Spectra: The best-fit first-order Markov red noise spectrum (red curve) and its 95% (blue curve) and 99% (green curve) confidence bounds are shown on each panel. Top X-axis shows the period (in years), and the bottom X-axis shows the frequency (cycles/mo). If calculated, the observational spectrum is overlaid in gray on each model spectrum. The spectra are displayed in variance-conserving form.
  • Wavelet Analysis: A wavelet transform is computed using a Morlet wavelet with a wavenumber of 6. Areas significant at the 95% based on a chi-square test are stippled and the “cone of influence” is hatched. See Torrence, C. and G. P. Compo, 1998: A Practical Guide to Wavelet Analysis. Bull. Amer. Meteor. Soc., 79, 61-78. doi: 10.1175/1520-0477(1998)079<0061:APGTWA>2.0.CO;2.
  • Metrics Tables: Area-weighted pattern correlations and rms differences are calculated between observations and each model simulation (regridded to match the observational grid) for 11 climate metrics. The Total Score column shows the average of the 11 pattern correlations (Z-transformed) and rms differences. The following domains are used to compute the pattern correlations and rms differences: Means, standard deviations, ENSO and PDV: Global; AMV: 63°S:65°N; El Nino and La Nina Hovmöllers: entire longitude/temporal range shown, NAM/NAO (20:90°N) and SAM (20:90°S).
  • EOF significance: If an eigenvalue is significantly seperated from neighboring values a star is appended to the percent variance explained on the plot. Significance is calculated following North et al. (MWR, 1982).
  • Ice Extent: Any grid cell defined as having a value at or above 15% is assumed to be 100% ice covered. The area of these grid cells is summed to create ice extent. For data sets with a northern hemisphere pole hole the CVDP looks for an attribute named pole_hole_area that is attached to the AREA variable in the netCDF file. The format of pole_hole_area should be start month (YYYYMM), end month (YYYYMM), area value, (repeat as necessary). If pole_hole_area is detected the CVDP will add the area value to the calculated northern hemisphere ice extent timeseries from the specified start month to the specified end month, and a * is placed after the dataset name in the ice extent plots.
  • The ensemble mean of a metric is the average of the individual members from that ensemble for that metric.
  • For more information on observational datasets and climate indices, see the Climate Data Guide.

Modes of Variability

ENSO Spatial Composites The normalized December nino3.4 timeseries is used to composite all years greater than 1 standard deviation (El Niño) and all years less that -1 standard deviation (La Niña). The number of El Niño/La Niña events composited is shown in the right subtitle. The December nino3.4 timeseries is based on the December values of the monthly nino3.4 time series smoothed with a 3-point binomial filter. TS/TAS/PSL composites: Temperatures are color shaded and in units of Celsius. Sea level pressure is contoured from -16 to 16hPa by 2hPa; negative contours are dashed. PR composites: Precipitation is color shaded and is in units of mm/day. See Deser, C., A. S. Phillips, R. A. Tomas, Y. Okumura, M. A. Alexander, A. Capotondi, J. D. Scott, Y. -O. Kwon, and M. Ohba, 2012: ENSO and Pacific Decadal Variability in Community Climate System Model Version 4. Journal of Climate, 25, 2622-2651, doi: 10.1175/JCLI-D-11-00301.1.
ENSO Hovmöllers A 1-2-1 running average is applied to the monthly nino3.4 timeseries, and then December values are selected and normalized. Meridional averages are calculated by averaging from 3°S:3°N, and spatial patterns are formed by compositing from Jan year 0 to May year 2 where the nino3.4 (1-2-1 weighted) December (year 0) index is greater than 1 (El Niño) and those years where the index is less than -1 (La Niña). See Deser, C., A. S. Phillips, R. A. Tomas, Y. Okumura, M. A. Alexander, A. Capotondi, J. D. Scott, Y. -O. Kwon, and M. Ohba, 2012: ENSO and Pacific Decadal Variability in Community Climate System Model Version 4. Journal of Climate, 25, 2622-2651, doi: 10.1175/JCLI-D-11-00301.1.
PDV (Pacific Decadal Variability) Monthly index timeseries defined as the leading principal component (PC) of North Pacific (20:70°N, 110°E:100°W) area-weighted SST anomalies. Pattern created by regressing SST anomalies (in Celsius) at each grid box onto the normalized PC timeseries. Low pass-filtered timeseries (black curve) is based on a a 61-month running mean. See Deser, C., M. A. Alexander, S. -P. Xie, and A. S. Phillips, 2010: Sea surface temperature variability: patterns and mechanisms. Ann. Rev. Mar. Sci., 2010.2, 115-143, doi:10.1146/annurev-marine-120408-151453. Also see Mantua, N. J., S. R. Hare, Y. Zhang, J. M. Wallace, and R. Francis, 1997: A Pacific interdecadal climate oscillation with impacts on salmon production. Bull. Amer. Met. Soc., 1069-1079. For more information on the PDO see the Climate Data Guide.
IPV (Standard) (Interdecadal Pacific Variability) Standard monthly index timeseries defined as the leading principal component (PC) of 13yr low pass filtered Pacific (40°S:60°N, 110°E:70°W) area-weighted SST anomalies. Pattern created by regressing SST anomalies (in Celsius) at each grid box onto the normalized PC timeseries. At least 40 years of data are required for the standard IPV index to be calculated. See Meehl, G.A. and A. Hu, 2007: Megadroughts in the Indian Monsoon Region and Southwest North America and a Mechanism for Associated Multidecadal Pacific Sea Surface Temperature Anomalies, Journal of Climate, 19, 1605-1623, doi: 10.1175/JCLI3675.1.
IPV (Henley) (Interdecadal Pacific Variability) Standard monthly index timeseries defined as the 10°S:10°N, 170°E:90°W area-weighted SST anomalies minus the average of the area-weighted anomalies from 25°:45°N, 140°E:145°W and 15°:50°S, 150°E:160°W. Pattern created by regressing resulting SST anomaly timeseries onto the SST anomalies at each grid box. See Henley, B.J., J. Gergis, D.J. Karoly, S. Power, J. Kennedy and C.K. Folland, 2015: A Tripole Index for the Interdecadal Pacific Oscillation, Clim. Dyn., 45, 3077-3090, doi: 10.1007/s00382-015-2525-1.
AMV (Atlantic Multidecadal Variability) Monthly index timeseries defined as area-weighted SST* anomalies averaged over the North Atlantic (0:60°N, 80°W:0°E), where SST* denotes that the global (60°S:60°N) mean SST anomaly has been removed at each timestep. Pattern created by regressing SST* anomalies onto the index timeseries and smoothing with a 9-point spatial filter. Low pass-filtered timeseries (black curve superimposed on the monthly timeseries) is based on a a 61-month running mean. Based on Trenberth, K. E., and D. J. Shea, 2006: Atlantic hurricanes and natural variability in 2005, Geophys. Res. Lett., 33, L12704, doi:10.1029/2006GL026894. Low-pass filtered regression maps ("Regr LP") use a 10-year running mean on both the index timeseries and the field being regressed. For more information on the AMO see the Climate Data Guide.
AMOC The Atlantic Meridional Overturning Circulation (AMOC) is defined as the oceanic meridional mass transport (Sv) in the Atlantic sector. To compute AMOC, we follow the methods of Danabasoglu et al. (2012). Here we use annual averages of the AMOC, weighted by the cosine of the latitude and vertical extent of each model layer. Areas in which AMOC variance is low (standard deviation < 1e-6 Sv) are set to missing values for clarity. The leading EOF and associated principal component (PC) timeseries are computed over the Atlantic basin from 33°S to 90°N. The AMOC patterns are created by regressing the AMOC anomalies (in Sv) onto the normalized PC timeseries. The SST/TAS patterns associated with AMOC variations are created by regressing TAS/SST anomalies (in Celsius) at each grid box over the globe onto the normalized AMOC PC timeseries. A 15-point low-pass Lanczos filter is applied to the AMOC PC (and AMO) timeseries prior to computing lead/lag correlations, with a minimum of 90 years of data required. See Danabasoglu, G., S. G. Yeager, Y. -O. Kwon, J. J. Tribbia, A. S. Phillips, and J. W. Hurrell, 2012. Variability of the Atlantic Meridional Overturning Circulation in CCSM4. J. Climate, 25, 5153-5172, doi: 10.1175/JCLI-D-11-00463.1.

For CCSM4 and CESM1, the MOC variable is read in, the Eulerian Mean, Eddy-Induced and Submeso components are summed, and the Atlantic Ocean + Mediterranean Sea + Labrador Sea + GIN Sea + Arctic Ocean + Hudson Bay transport region is selected. For CCSM2 and CCSM3 the same transport region is selected but only the Eulerian Mean component is used as that is all that is available. For CMIP5 (CMIP6) data the msftmyz (msftmz) variable is read in and the atlantic_arctic_ocean basin is used. For CMIP3 data, the stfmmc variable is read in and the atlantic_ocean geo_region is used.
NAM (Northern Annular Mode) Seasonal/annual PSL averages are formed, square root of the cosine of the latitude weighting is applied, and then the leading EOF and associated principal component (PC) timeseries are computed over 20:90°N, 0:360°E. Pattern created by regressing global PSL anomalies (in hPa) onto normalized PC timeseries. Based on Hurrell, J. W., and C. Deser, 2009: North Atlantic climate variability: The role of the North Atlantic Oscillation. J. Mar. Syst., 78, 28-41, doi:10.1016/j.jmarsys.2008.11.026. Also see Thompson, D. W. J., and J. M. Wallace, 2000: Annular modes in the extratropical circulation. Part I: Month-to-month variability. Journal of Climate, 13, 1000-1016.
NAO (North Atlantic Oscillation) Seasonal/annual PSL averages are formed, square root of the cosine of latitude weighting is applied, and then the leading EOF and associated principal component (PC) timeseries are computed over 20:80°N, 90°W:40°E. Pattern created by regressing global PSL anomalies (in hPa) onto normalized PC timeseries. Based on Hurrell, J. W. and C. Deser, 2009: North Atlantic climate variability: The role of the North Atlantic Oscillation. J. Mar. Syst., 78, 28-41, doi:10.1016/j.jmarsys.2008.11.026. For more information on the NAO see the Climate Data Guide.
SAM/PSA1/PSA2 (Southern Annular Mode, Pacific South American Patterns 1/2) Seasonal/annual PSL averages are formed, square root of the cosine of latitude weighting is applied, and then the 1st (SAM), 2nd (PSA1) and 3rd (PSA2) EOFs and associated principal component (PC) timeseries are computed over 20:90°S, 0:360°E. Patterns created by regressing global PSL anomalies (in hPa) onto normalized PC timeseries. SAM calculation based on Thompson, D. W. J. and J.M. Wallace, 2000: Annular modes in the extratropical circulation. Part I: Month-to-month variability. Journal of Climate, 13, 1000-1016.
PNA/NPO (Pacific North American Pattern, North Pacific Oscillation) Seasonal/annual PSL averages are formed, the square root of the cosine of the latitude weighting is applied, and then the 1st (PNA) and 2nd (NPO) EOFs and associated principal component (PC) timeseries are computed over 20:85°N, 120°E:120°W. Patterns created by regressing global PSL anomalies (in hPa) onto normalized PC timeseries.
SST Regressions SST anomalies (in Celsius) at each grid box are regressed upon the normalized atmospheric mode timeseries.
TAS Regressions TAS anomalies (in Celsius) at each grid box are regressed upon the normalized atmospheric mode timeseries.
PR Regressions PR anomalies (in mm/day) at each grid box are regressed upon the normalized atmospheric mode timeseries.

Climate Indices

Atlantic SST Meridional Mode Defined as the difference between area-averaged SST anomalies (in Celsius) computed over 5:15°N, 20:50°W and area-averaged SST anomalies computed over 5:15°S, 20°W:10°E. See Doi, T., T. Tozuka and T. Yamagata (2009), Interannual variability of the Guinea Dome and its possible link with the Atlantic Meridional Mode. Climate Dynamics, 33, 985-998, doi:10.1007/s00382-009-0574-z.
Atlantic Niño SST Area-averaged SST anomalies (in Celsius) computed over 3°S:3°N, 20°W:0°E. See Zebiak, S. E., (1993): Air–sea interaction in the equatorial Atlantic region. Journal of Climate, 6, 1567–1586.
North Atlantic SST Area-averaged SST anomalies (in Celsius) computed over 0°:60°N, 90°W:0°E.
Tropical North Atlantic SST Area-averaged SST anomalies (in Celsius) computed over 5.5:23.5°N, 15:57.5°W. See Enfield, D.B., A.M. Mestas, D.A. Mayer, and L. Cid-Serrano (1999), How ubiquitous is the dipole relationship in tropical Atlantic sea surface temperatures?, JGR-O, 104, 7841-7848.
Tropical South Atlantic SST Area-averaged SST anomalies (in Celsius) computed over 0:20°S, 30°W:10°E. See Enfield, D.B., A.M. Mestas, D.A. Mayer, and L. Cid-Serrano (1999), How ubiquitous is the dipole relationship in tropical Atlantic sea surface temperatures?, JGR-O, 104, 7841-7848.
niño1+2 Area-averaged SST anomalies (in Celsius) computed over 0:10°S, 80:90°W.
niño3 Area-averaged SST anomalies (in Celsius) computed over 5°S:5°N, 90:150°W. See Trenberth, K. E. (1997) The Definition of El Niño. Bulletin of the American Meteorological Society, 78, 2771-2777.
niño3.4 Area-averaged SST anomalies (in Celsius) computed over 5°S:5°N, 120:170°W. See Trenberth, K. E. (1997) The Definition of El Niño. Bulletin of the American Meteorological Society, 78, 2771-2777.
niño4 Area-averaged SST anomalies (in Celsius) computed over 5°S:5°N, 160°E:150°W.
NPI (North Pacific PSL Index) Winter (December-March) average PSL anomalies (in hPa) area-averaged over 30°:65°N, 160°E:140°W. Based on Trenberth, K. E. and J. W. Hurrell, 1994: Decadal atmosphere-ocean variations in the Pacific, Climate Dynamics, 9, 303-319.
North Pacific Meridional Mode See Chiang, J. C. H. and D. J. Vimont, (2004): Analogous Pacific and Atlantic Meridional Modes of Tropical Atmosphere–Ocean Variability. Journal of Climate, 17, 4143–4158, doi: 10.1175/jcli4953.1.
South Pacific Meridional Mode See Zhang, H., A. Clement and P. Di Nezio, (2014): The South Pacific Meridional Mode: A Mechanism for ENSO-like Variability. Journal of Climate, 27, 769–783, doi: 10.1175/jcli-d-13-00082.1.
Indian Ocean SST Dipole Defined as the difference between area-averaged SST anomalies (in Celsius) computed over 10°S:10°N, 50:70°E and area-averaged SST anomalies computed over 0:10°S, 90:110°E. See: Saji N.H., Goswami B.N., Vinayachandran P.N., Yamagata T., 1999: A dipole mode in the tropical Indian Ocean, Nature, 401, 360-363.
Tropical Indian Ocean SST Area-averaged SST anomalies (in Celsius) computed over 15°S:15°N, 40:110°E.
Southern Ocean SST Area-averaged SST anomalies (in Celsius) computed over 50°:70°S, 0:360°E.

Created Fri Apr 19 22:55:12 MDT 2024

CVDP Version 5.9.0