alan robock department of environmental sciences

Alan RobockDepartment of Environmental Sciences

Rutgers University, New Brunswick, New Jersey USA

[email protected]

http://envsci.rutgers.edu/~robock

Climate Dynamics11:670:461

Lecture 10, 10/6/14



Alan RobockDepartment of Environmental

Sciences

Predictability: How can we predict the climate decades into the future when we can’t even predict the weather for next week?

Predictability of the first kind: Predict the future based on initial conditions, with boundary conditions constant. This is limited by the chaotic nature of the atmosphere, which is a physical system with built-in instabilities, in vertical convection (e.g., thunderstorms) and horizontal motion (e.g., baroclinic instability - development of low pressure systems, such as hurricanes and Nor’Easters).


Sciences

Consider a prediction using the above equation of the future state of the variable X, say the surface air temperature. The subscript n indicates the time, say the day, and a is a constant representing the physics of the climate system. X for any day is a times its value on the previous day minus X squared on the previous day.

With such a simple equation, it should be possible to predict X indefinitely into the future. Right?

X0 2.200 2.200 2.210 2.20a 3.930 3.940 3.930 3.93

Precision (decimal places)

3 3 3 2

n (Time Step) Control PhysicsInitial

ConditionsRounding

0 2.200 2.200 2.210 2.201 3.806 3.828 3.801 3.812 0.472 0.429 0.490 0.463 1.632 1.506 1.686 1.604 3.750 3.666 3.783 3.735 0.675 1.004 0.556 0.756 2.197 2.948 1.876 2.397 3.807 2.924 3.853 3.688 0.468 2.971 0.297 0.929 1.620 2.879 1.079 2.7710 3.742 3.055 3.076 3.2111 0.703 2.704 2.627 2.3112 2.269 3.342 3.423 3.7413 3.769 1.999 1.735 0.7114 0.607 3.880 3.808 2.2915 2.017 0.233 0.465 3.7616 3.859 0.864 1.611 0.6417 0.274 2.658 3.736 2.1118 1.002 3.408 0.725 3.8419 2.934 1.813 2.324 0.3520 2.922 3.856 3.732 1.25

Xn+1 = a Xn - Xn2


Sciences

X0 2.200 2.200 2.210 2.20

a 3.930 3.940 3.930 3.93Precision

(decimal places)3 3 3 2


ConditionsRounding

0 2.200 2.200 2.210 2.201 3.806 3.828 3.801 3.812 0.472 0.429 0.490 0.463 1.632 1.506 1.686 1.604 3.750 3.666 3.783 3.735 0.675 1.004 0.556 0.756 2.197 2.948 1.876 2.397 3.807 2.924 3.853 3.688 0.468 2.971 0.297 0.929 1.620 2.879 1.079 2.7710 3.742 3.055 3.076 3.2111 0.703 2.704 2.627 2.3112 2.269 3.342 3.423 3.7413 3.769 1.999 1.735 0.7114 0.607 3.880 3.808 2.2915 2.017 0.233 0.465 3.7616 3.859 0.864 1.611 0.6417 0.274 2.658 3.736 2.1118 1.002 3.408 0.725 3.8419 2.934 1.813 2.324 0.3520 2.922 3.856 3.732 1.25

Xn+1 = a Xn - Xn2

Let’s assume that a is exactly 3.930 and that a prediction with three decimal places is the exact solution.

Then let’s consider three types of errors: imprecise knowledge of the physics of the climate system, imprecise initial conditions, and rounding due to limited computer resources.

This example is from Edward Lorenz.


Sciences

X0 2.200 2.200 2.210 2.20

a 3.930 3.940 3.930 3.93Precision



ConditionsRounding

0 2.200 2.200 2.210 2.201 3.806 3.828 3.801 3.812 0.472 0.429 0.490 0.463 1.632 1.506 1.686 1.604 3.750 3.666 3.783 3.735 0.675 1.004 0.556 0.756 2.197 2.948 1.876 2.397 3.807 2.924 3.853 3.688 0.468 2.971 0.297 0.929 1.620 2.879 1.079 2.7710 3.742 3.055 3.076 3.2111 0.703 2.704 2.627 2.3112 2.269 3.342 3.423 3.7413 3.769 1.999 1.735 0.7114 0.607 3.880 3.808 2.2915 2.017 0.233 0.465 3.7616 3.859 0.864 1.611 0.6417 0.274 2.658 3.736 2.1118 1.002 3.408 0.725 3.8419 2.934 1.813 2.324 0.3520 2.922 3.856 3.732 1.25

Xn+1 = a Xn - Xn2


Sciences

X0 2.200 2.200 2.210 2.20

a 3.930 3.940 3.930 3.93Precision



ConditionsRounding

0 2.200 2.200 2.210 2.201 3.806 3.828 3.801 3.812 0.472 0.429 0.490 0.463 1.632 1.506 1.686 1.604 3.750 3.666 3.783 3.735 0.675 1.004 0.556 0.756 2.197 2.948 1.876 2.397 3.807 2.924 3.853 3.688 0.468 2.971 0.297 0.929 1.620 2.879 1.079 2.7710 3.742 3.055 3.076 3.2111 0.703 2.704 2.627 2.3112 2.269 3.342 3.423 3.7413 3.769 1.999 1.735 0.7114 0.607 3.880 3.808 2.2915 2.017 0.233 0.465 3.7616 3.859 0.864 1.611 0.6417 0.274 2.658 3.736 2.1118 1.002 3.408 0.725 3.8419 2.934 1.813 2.324 0.3520 2.922 3.856 3.732 1.25

Xn+1 = a Xn - Xn2

Xn+1 = a Xn - Xn2

0.00.51.01.52.02.53.03.54.0

0 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20

n

Xn

Control Physics Initial Conditions Rounding


Sciences

Predictability: How can we predict the climate decades into the future when we can’t even predict the weather for next week?Predictability of the second kind: Predict the future based on boundary conditions, independent of initial conditions. If there are slowly-varying (with respect to the atmospheric predictability limit of 2-3 weeks) boundary conditions (e.g., greenhouse gases, stratospheric aerosols, sea surface temperatures, soil moisture, snow cover) that can be predicted, then the envelope of the weather can be predicted. [The first two examples are external to the climate system, and the last three are internal.]


Sciences

NOAA Medium Range Forecasts

http://www.hpc.ncep.noaa.gov/medr/medr.shtml




Sciences

https://climatedataguide.ucar.edu/climate-data/era-interim





Sciences

http://www.ecmwf.int/products/forecasts/d/charts/medium/verification/timeseries/ccafadrian/

ECMWF forecast skill


Sciences

MJO forecast from Dee et al.

(2011).

Correlation of> 0.6 has skill.

Dee, D. P., et al., 2011: The ERA-Interim reanalysis: configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc., 137, 553–597. DOI:10.1002/qj.828


Sciences

NOAA MJO forecasts

(1) Reanalysis 2 (R2) is used since CFS operational forecast utilizes the R2 as initialization data.

(2) The indices for the latest 4 days are calculated using NCEP GDAS (Global Data Analysis System).

(3) The MJO definition used here is identical to the Matt Wheeler's (Wheeler and Hendon 2004), i.e., to represent the MJO, the first two EOFs of combined fields of OLR, u850 and u200 are used. The followings are some details of the forecast models.

(4) CFS operational: this is a 2003 version and two member ensemble mean is used.

(5) GFS offline: this runs exactly the same as CFS operational model (e.g. the same R2 initial data) except that air-sea interactionis not allowed. Four member ensemble mean is used.

(6) GFS operational: the model keeps being updated. Model climatology from the GFS offline model is used. The 11-member ensemble mean is used.

(7) AR: Autoregressive time series model.(8) PCRLAG: Lagged multiple linear regression.

For details please contact to Kyong-Hwan Seo ([email protected]).

http://www.cpc.ncep.noaa.gov/products/people/wd52qz/mjoindex/description_methods_forecasts.html

mailto:[email protected]

http://www.cpc.ncep.noaa.gov/products/people/wd52qz/mjoindex/description_methods_forecasts.html


Sciences

NOAA MJO forecast

from GEFS model

http://www.cpc.ncep.noaa.gov/products/people/wd52qz/mjoindex/index/diagram_40days_forecast_GEFS_membera.gif


Sciences

NOAA MJO forecast

from statistical models

http://www.cpc.ncep.noaa.gov/products/people/wd52qz/mjoindex/index/diagram_40days_forecast.gif


Sciences

Dee, D. P., et al., 2011: The ERA-Interim reanalysis: configuration and performance of the data assimilation system. Q. J. R. Meteorol. Soc., 137, 553–597. DOI:10.1002/qj.828

Persistence

ERA-Interim


Sciences

CONSTRUCTED ANALOG METHOD, Huug van den Dool, http://www.cpc.ncep.noaa.gov/products/people/wd51hd/


Sciences

http://www.ecmwf.int/products/forecasts/d/charts/seasonal/forecast/seasonal_range_forecast/nino_plumes_public_s4!3.4!plumes!201310/


Sciences



Sciences

http://www.ecmwf.int/products/forecasts/d/charts/seasonal/forecast/eurosip/nino_plumes_euro_public!3.4!201309!/


Sciences


alan robock department of environmental sciences

Documents

climate system

climate decades

imprecise initial conditions

boundary conditions

future state

weeks boundary conditions

physical system

previous day