alan robock department of environmental sciences

25
Alan Robock Department of Environmental Sciences Rutgers University, New Brunswick, New Jersey USA [email protected]. edu http://envsci.rutgers.edu/~ robock Climate Dynamics 11:670:461 Lecture 10, 10/6/14

Upload: bert-davis

Post on 02-Jan-2016

30 views

Category:

Documents


1 download

DESCRIPTION

Lecture 10, 10/6/14. Climate Dynamics 11:670:461. Alan Robock Department of Environmental Sciences Rutgers University, New Brunswick, New Jersey USA. [email protected]. http://envsci.rutgers.edu/~ robock. - PowerPoint PPT Presentation

TRANSCRIPT

Page 1: 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

Page 2: Alan Robock Department of Environmental Sciences

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).

Page 3: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

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

Page 4: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

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

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

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.

Page 5: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

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

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

Page 6: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

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

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

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

Page 7: Alan Robock Department of Environmental Sciences

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 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.]

Page 8: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

NOAA Medium Range Forecasts

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

Page 9: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

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

Page 10: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

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

ECMWF forecast skill

Page 11: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

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

Page 12: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

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

Page 13: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

NOAA MJO forecast

from GEFS model

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

Page 14: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

NOAA MJO forecast

from statistical models

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

Page 15: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

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

Page 16: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

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

Page 17: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

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

Page 18: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

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

Page 19: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

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

Page 20: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

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

Page 21: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

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

Page 22: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

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

Page 23: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

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

Page 24: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

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

Page 25: Alan Robock Department of Environmental Sciences

Alan RobockDepartment of Environmental

Sciences

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