Bred Vectors, Singular Vectors, and Lyapunov Vectors in Simple and Complex Models
|
|
- Mitchell Walton
- 5 years ago
- Views:
Transcription
1 Bred Vectors, Singular Vectors, and Lyapunov Vectors in Simple and Complex Models Adrienne Norwood Advisor: Eugenia Kalnay With special thanks to Drs. Kayo Ide, Brian Hunt, Shu-Chih Yang, and Christopher Wolfe for their guidance and suggestions.
2 Outline Introduction Computation of bred vectors (BVs), singular vectors (SVs), and Lyapunov vectors (LVs) Results Lorenz 1963 model which has 3 degrees of freedom A Fast-Slow coupled Lorenz model developed by Peña and Kalnay (2004) with 9 degrees of freedom A quasi-geostrophic model with 15,015 degrees of freedom, but a single type of instability. SPEEDY, a full atmospheric model with 135,240 degrees of freedom that contains several different types of instabilities. Summary
3 Bred, Singular, and Lyapunov Vectors Instabilities inherent in the atmosphere-ocean system would degrade forecasts even if models and observations were perfect. BVs, SVs, and LVs are three types of vectors frequently used to explore the instabilities of dynamical systems. BVs are by far the easiest (and cheapest) to compute while LVs are the most difficult (and expensive). All perturbations integrated with the tangent linear model (TLM) will, in time, align with the leading LV. Thus linear BVs, those with small amplitudes and short rescaling windows, align with LV1. However, when LV2 grows faster than LV1, BVs align with LV2 (Norwood et al., 2013). SVs are orthonormal and optimized for a particular integration window using a specific norm, but as the window length approaches infinity the final SVs converge to orthonormalized LVs (the basis of the Wolfe and Samelson 2007 method to obtain LVs).
4 Computing Bred Vectors Forecast values Initial random perturbation Bred Vectors ~LLVs Control forcast Integrate nonlinear model, M, to obtain control trajectory, x c. Choose rescaling amplitude, δ 0, and integration window, IW, to target the desired mode of growth (Peña and Kalnay, 2004). Add a small perturbation p x p t i = x c t i + δ 0, p where p is the initial direction of the perturbation and δ p 0 is the amplitude. Integrate x p t i forward using the nonlinear model. The bred vector is bv t i+iw = x p t i+iw x c t i+iw. Rescale to δ 0, and repeat the process. time
5 Computing Singular Vectors For nonlinear model, M, calculate the tangent linear model M ij t i, t i+iw = M i and the adjoint M x ij t i, t i+iw. j Initial singular vectors, ξ j, are the columns of Ξ where M MΞ = ΞS where S is diagonal with entries σ j 2, the squares of the singular values The final singular vectors, η j, are obtained by integrating the TLM forward starting from the initial singular vector. M t i, t i+iw ξ j t i = σ j η j t i+iw ξ 2 M σ 1 η 1 ξ 1 σ 2 η 2
6 Computing Lyapunov Vectors (Wolfe and Samelson, 2007) Fig. 6.7: Schematic of how all perturbations will converge towards the leading Local Lyapunov Vector Extend the integration window for the SVs forward and backward in time until the SVs converge. These asymptotic initial and final SVs shall be denoted as ξ j and η j, respectively. Every perturbation integrated forward in time aligns with the leading LV. Thus φ 1 = < φ 1, η j > η 1. SVs are orthogonal so LV2 has a component in the direction of η 2. Thus φ 2 =< φ 2, η 1 > η 1 + < φ 2, η 2 > η 2, and so on. Every perturbation integrated backward in time aligns with the fastest decaying LV. Thus φ N = < φ N, ξ N > ξ N, where N is the degrees of freedom of the system. Likewise φ N 1 = < φ N 1, ξ N 1 > ξ N 1 + < φ N 1, ξ N > ξ N, and so on random initial perturbations trajectory leading local Lyapunov vector
7 Wolfe and Samelson (WS07) Cont d n φ n = j=1 < η j, φ n > η j (1) N φ n = j=n < ξ j, φ n > ξ j (2) Use (1) to find < φ n, ξ k >, for k n, and (2) to find< φ n, η j >, for k n. Use substitution and patience to find the coefficients of (1) and (2). These Lyapunov vectors (LVs) will be invariant under the tangent linear model, i.e. they can be integrated with the TLM for intervals that are not too long.
8 Results: Lorenz (1963) Model This model has only 3 degrees of freedom. x = σ y x y = ρx y x z = xy βz using the standard parameters of σ = 10, β = 8 3 The integration time step is.01., and ρ = 28.
9 Lorenz (1963) BVs, SVs, and LVs BVs were computed using δ 0 =.1 and IW =.02 units. These vectors are essentially equal to the LV1 and are labeled linear BVs. BVs were also computed using δ 0 = 1 and IW =.08 units. These are less linear than the above and are labeled nonlinear. SVs were computed using integration windows of.02 units.
10 Results with Lorenz (1963) Linear BV Growth LV1 Growth Nonlinear BV Growth FSV1 Growth Red stars indicate periods of fastest growth Fast growth typically signals a regime change (Evans et al., 2004; Norwood et al., 2013) The linear BV is most similar to LV1, with an average correlation (cosine between the vectors) of.996. All of the leading vectors can be used to predict regime changes but FSV1 is the best predictor.
11 Fast-Slow Coupled Model with Extratropics (Peña and Kalnay, 2004) Fast Extratropical Atmosphere Fast Tropical Atmosphere Slow Ocean.08 1 The temporal scaling factor for the ocean, τ o =.1 so the ocean is 10x slower than the other subsystems. Ocean Trajectory The extratropics are like weather noise when studying El Niño, which appears chaotically every 2-7 years.
12 Fast-Slow Coupled Model (FSCM; Peña and Kalnay, 2004) x e = τ e σ y e x e c e Sx t + k 1 y e = τ e ρx e τ e y e τ e x e z e + c e Sy t + k 1 z e = τ e x e y e τ e βz e c t z t x t = σ y t x t c SX + k 2 c e Sx e + k 1 y t = ρx t y t x t z t + c SY + k 2 + c e Sy e + k 1 z t = x t y t βz t + c z Z + c t z e X = τ o σ Y X c x t + k 2 Y = τ o ρx τ o Y τ o SXZ + c y t + k 2 Z = τ o SXY τ o βz c z Z Lower case variables are the fast modes (extratropics and tropics), upper case are the slow modes (ocean). c s are coupling coefficients. k s are uncentering parameters τ and S are temporal and scaling factors, respectively.
13 LV Growth Approximately Corresponds to Particular Subsystems LV1 Growth LV4 Growth LV8 Decay FSV1 Growth Fast Extratropics LV5 Growth LV6 Growth LV9 Decay Fast Tropics LV2 Growth LV3 Growth LV7 Decay Slow Mode BV Growth Slow Ocean
14 Fast-Slow Coupled Model with Convection We accelerate the extratropical atmosphere by setting τ e = 10. This makes the weather noise more like convective noise. The coupling remains the same. The convective and tropical subsystems are weakly coupled and the tropical and ocean subsystems are strongly coupled. Fast Convection Fast Tropical Atmosphere Slow Ocean.08 1
15 LV and SV Growth Rates Approximately Correspond to Particular Subsystems LV1 Growth LV2 Growth FSV1 Growth Fast Convection LV3 Growth LV6 Growth FSV7 Decay Fast Tropics LV7 Growth FSV8 Decay FSV9 Decay Slow Mode BV Growth Slow Ocean
16 A Quasi-Geostrophic Model Based upon Rotunno and Bao (1996), the nondimensional form is q t + ψ x q y ψ y q x = 0 where potential vorticity q = ψ xx + ψ yy + ψ z S z z, ψ is the geopotential, and S z is the stratification parameter. The model has 7 levels on a 65 x 33 grid leading to 15,015 degrees of freedom. This model only exhibits baroclinic instability.
17 BVs for the QG Model Five different BVs with an amplitude of 1 and integration window of 24 hours were computed. Two different random initial BV perturbations on top of one anther.
18 BVs for the QG Model Five different BVs with an amplitude of 1 and integration window of 24 hours all converge to the leading LV. Two different random initial BV perturbations on top of one anther. These two BVS collapse into a single vector, i.e. they converge to the leading LV.
19 BVs for the QG Model Five different BVs with an amplitude of 1 and integration window of 24 hours all converge to the leading LV. Two different random initial BV perturbations on top of one anther. Initial dimension ~ 5 These two BVS collapse into a single vector, i.e. they converge to the leading LV.
20 BVs for the QG Model Five different BVs with an amplitude of 1 and integration window of 24 hours all converge to the leading LV. SVs computed using 4, 5, 6, and 7 day windows remain different. Consequently we cannot compute the LVs using the WS07 algorithm. Two different random initial BV perturbations on top of one anther. Initial dimension ~ 5 These two BVS collapse into a single vector, i.e. they converge to the leading LV. Final dimension = 1
21 What happens with SPEEDY? SPEEDY is a full atmospheric model with weather waves, convective instabilities, and even inertia gravity waves (i.e., Lamb waves) excited by tropical convection. There are 7 levels on a 96x48 grid with 6 variables (2 only at the surface) leading to 135,240 degrees of freedom. Five BVs were computed. With an amplitude of 1 m/s and integration window of 24 hours, this BV targets baroclinic instabilities, stronger in the winter hemisphere than in the summer hemisphere. They do NOT converge to a leading LV Two different BVs on top of one another.
22 SPEEDY Cont d We now reduce the amplitude of the wind to 1 cm/s and use a rescaling window of 6 hours to target convective instabilities. The BVs clearly remain distinct, although there are some unstable regions where they align to a local LV.
23 SPEEDY Cont d If we take a very small amplitude (1 mm/s) and a very short rescaling window (40 minutes), we obtain a leading LV corresponding to a global Lamb Wave, probably triggered by convection in the Warm Pool. There appears to be a leading LV for the SPEEDY model, but it probably would be useless for applications.
24 Summary BVs require the least amount of computational effort, time, and memory BVs can target instabilities through proper choice of perturbation size and integration window. BVs, SVs, and LVs can be used as predictors of regime changes in simple models. LVs can distinguish between various modes of growth, but not as cleanly as BVs, being influenced by more than one mode if there is strong coupling between the two. SVs can only differentiate between more than one mode of growth when the modes are very different in terms of frequency. SVs failed to converge for the QG model. Thus we were unable to use WS07 to compute the LVs for more complex models. BVs identified 3 types of instabilities: baroclinic waves, convection, and global Lamb waves. Contrary to what Toth and Kalnay (1997) hypothesized, there is a global leading LV for a weather model! However, the leading LV of complex systems with several types of instabilities will identify the fastest mode (such as global Lamb waves triggered by tropical convection), but this is not useful for weather forecasting applications.
25 Summary BVs require the least amount of computational effort, time, and memory BVs can target instabilities through proper choice of perturbation size and integration window. BVs, SVs, and LVs can be used as predictors of regime changes in simple models. LVs can distinguish between various modes of growth, but not as cleanly as BVs, being influenced by more than one mode if there is strong coupling between the two. SVs can only differentiate between more than one mode of growth when the modes are very different in terms of frequency. SVs failed to converge for the QG model. Thus we were unable to use WS07 to compute the LVs for more complex models. BVs identified 3 types of instabilities: baroclinic waves, convection, and global Lamb waves. Contrary to what Toth and Kalnay (1997) hypothesized, there is a global leading LV for a weather model! However, the leading LV of complex systems with several types of instabilities will identify the fastest mode (such as global Lamb waves triggered by tropical convection), but this is not useful for weather forecasting applications. Thank you!!
26 Thank you!!
27 Summary Cont d SVs failed to converge for the QG model. Thus we were unable to use WS07 to compute the LVs for more complex models. BVs identified 3 types of instabilities: baroclinic waves, convection, and global Lamb waves. The leading LV of complex systems with several types of instabilities will identify the fastest mode, which may not be useful for applications.
28 Summary Cont d SVs failed to converge for the QG model. Thus we were unable to use WS07 to compute the LVs for more complex models. BVs identified 3 types of instabilities: baroclinic waves, convection, and global Lamb waves. The leading LV of complex systems with several types of instabilities will identify the fastest mode, which may not be useful for applications. Thank you!!
29 Definition of Growth Rates Bred vector (BV) growth rates are defined as where dt is the integration time step. 1 IW dt ln( bv δ 0 ), Singular vector (SV) growth rates are defined as ln σ j dt IW. Lyapunov vector (LV) growth rates are defined as 1 φ n d φ n dt
30 Results with Lorenz (1963) All of the leading vectors can be used to predict regime changes.
31 FSCM with Convection Weak coupling Weaker coupling
32 Fast Mode Vectors The BVs can distinguish between the fastest and slowest modes of growth The FSVs are able to distinguish between the various subsystems and modes of growth, unlike with the previous setup The LVs correspond to particular subsystems.
33 Slow Mode Vectors Each vector is able to capture the slowest mode of growth, but the strong coupling between the tropical and ocean subsystems leads to the tropical subsystem influencing the growth of these vectors as well.
34 Mixed LVs and SVs Fast Tropics Slow Ocean
ABSTRACT. instability properties of dynamical models, Lyapunov vectors (LVs), singular vectors
ABSTRACT Title of Thesis: BRED VECTORS, SINGULAR VECTORS, AND LYAPUNOV VECTORS IN SIMPLE AND COMPLEX MODELS Adrienne Norwood, Doctor of Philosophy, 2015 Directed by: Professor Eugenia Kalnay, Atmospheric
More informationChapter 6: Ensemble Forecasting and Atmospheric Predictability. Introduction
Chapter 6: Ensemble Forecasting and Atmospheric Predictability Introduction Deterministic Chaos (what!?) In 1951 Charney indicated that forecast skill would break down, but he attributed it to model errors
More informationBred vectors: theory and applications in operational forecasting. Eugenia Kalnay Lecture 3 Alghero, May 2008
Bred vectors: theory and applications in operational forecasting. Eugenia Kalnay Lecture 3 Alghero, May 2008 ca. 1974 Central theorem of chaos (Lorenz, 1960s): a) Unstable systems have finite predictability
More informationBred Vectors: A simple tool to understand complex dynamics
Bred Vectors: A simple tool to understand complex dynamics With deep gratitude to Akio Arakawa for all the understanding he has given us Eugenia Kalnay, Shu-Chih Yang, Malaquías Peña, Ming Cai and Matt
More informationWe honor Ed Lorenz ( ) who started the whole new science of predictability
Atmospheric Predictability: From Basic Theory to Forecasting Practice. Eugenia Kalnay Alghero, May 2008, Lecture 1 We honor Ed Lorenz (1917-2008) who started the whole new science of predictability Ed
More informationRelationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF
Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Eugenia Kalnay and Shu-Chih Yang with Alberto Carrasi, Matteo Corazza and Takemasa Miyoshi 4th EnKF Workshop, April 2010 Relationship
More informationEnsemble prediction and strategies for initialization: Tangent Linear and Adjoint Models, Singular Vectors, Lyapunov vectors
Ensemble prediction and strategies for initialization: Tangent Linear and Adjoint Models, Singular Vectors, Lyapunov vectors Eugenia Kalnay Lecture 2 Alghero, May 2008 Elements of Ensemble Forecasting
More informationRelationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF
Relationship between Singular Vectors, Bred Vectors, 4D-Var and EnKF Eugenia Kalnay and Shu-Chih Yang with Alberto Carrasi, Matteo Corazza and Takemasa Miyoshi ECODYC10, Dresden 28 January 2010 Relationship
More informationCoupled ocean-atmosphere ENSO bred vector
Coupled ocean-atmosphere ENSO bred vector Shu-Chih Yang 1,2, Eugenia Kalnay 1, Michele Rienecker 2 and Ming Cai 3 1 ESSIC/AOSC, University of Maryland 2 GMAO, NASA/ Goddard Space Flight Center 3 Dept.
More informationCoupled Ocean-Atmosphere Assimilation
Coupled Ocean-Atmosphere Assimilation Shu-Chih Yang 1, Eugenia Kalnay 2, Joaquim Ballabrera 3, Malaquias Peña 4 1:Department of Atmospheric Sciences, National Central University 2: Department of Atmospheric
More informationData assimilation; comparison of 4D-Var and LETKF smoothers
Data assimilation; comparison of 4D-Var and LETKF smoothers Eugenia Kalnay and many friends University of Maryland CSCAMM DAS13 June 2013 Contents First part: Forecasting the weather - we are really getting
More informationMultiple-scale error growth and data assimilation in convection-resolving models
Numerical Modelling, Predictability and Data Assimilation in Weather, Ocean and Climate A Symposium honouring the legacy of Anna Trevisan Bologna, Italy 17-20 October, 2017 Multiple-scale error growth
More information6.5 Operational ensemble forecasting methods
6.5 Operational ensemble forecasting methods Ensemble forecasting methods differ mostly by the way the initial perturbations are generated, and can be classified into essentially two classes. In the first
More informationComparing Local Ensemble Transform Kalman Filter with 4D-Var in a Quasi-geostrophic model
Comparing Local Ensemble Transform Kalman Filter with 4D-Var in a Quasi-geostrophic model Shu-Chih Yang 1,2, Eugenia Kalnay 1, Matteo Corazza 3, Alberto Carrassi 4 and Takemasa Miyoshi 5 1 University of
More informationPhase space, Tangent-Linear and Adjoint Models, Singular Vectors, Lyapunov Vectors and Normal Modes
Phase space, Tangent-Linear and Adjoint Models, Singular Vectors, Lyapunov Vectors and Normal Modes Assume a phase space of dimension N where Autonomous governing equations with initial state: = is a state
More informationABSTRACT. from a few minutes (e.g. cumulus convection) to years (e.g. El Niño). Fast time
ABSTRACT Title of dissertation: DATA ASSIMILATION EXPERIMENTS WITH A SIMPLE COUPLED OCEAN-ATMOSPHERE MODEL Tamara Singleton, Doctor of Philosophy, 2011 Dissertation co-directed by: Professors Eugenia Kalnay
More informationSeparating fast and slow modes in coupled chaotic systems
Separating fast and slow modes in coupled chaotic systems M. Peña and E. Kalnay Department of Meteorology, University of Maryland, College Park Manuscript submitted to Nonlinear Processes in Geophysics
More informationJ9.4 ERRORS OF THE DAY, BRED VECTORS AND SINGULAR VECTORS: IMPLICATIONS FOR ENSEMBLE FORECASTING AND DATA ASSIMILATION
J9.4 ERRORS OF THE DAY, BRED VECTORS AND SINGULAR VECTORS: IMPLICATIONS FOR ENSEMBLE FORECASTING AND DATA ASSIMILATION Shu-Chih Yang 1, Matteo Corazza and Eugenia Kalnay 1 1 University of Maryland, College
More informationHOW TO IMPROVE FORECASTING WITH BRED VECTORS BY USING THE GEOMETRIC NORM
HOW TO IMPROVE FORECASTING WITH BRED VECTORS BY USING THE GEOMETRIC NORM Diego Pazó Instituto de Física de Cantabria (IFCA), CSIC-UC, Santander (Spain) Santander OUTLINE 1) Lorenz '96 model. 2) Definition
More informationPredictability and data assimilation issues in multiple-scale convection-resolving systems
Predictability and data assimilation issues in multiple-scale convection-resolving systems Francesco UBOLDI Milan, Italy, uboldi@magritte.it MET-NO, Oslo, Norway 23 March 2017 Predictability and data assimilation
More informationLyapunov Vectors and Bred Vectors in a Fast- Slow Model
Lyapuov Vectors ad Bred Vectors i a Fast- Slow Model Adriee Norwood With great thaks to: Eugeia Kalay Kayo Ide Christopher Wolfe ad Yig Zhag April 9 2012 Outlie Wolfe- Samelso (2007) algorithm for fidig
More informationSome ideas for Ensemble Kalman Filter
Some ideas for Ensemble Kalman Filter Former students and Eugenia Kalnay UMCP Acknowledgements: UMD Chaos-Weather Group: Brian Hunt, Istvan Szunyogh, Ed Ott and Jim Yorke, Kayo Ide, and students Former
More informationEnsemble Forecasting at NCEP and the Breeding Method
3297 Ensemble Forecasting at NCEP and the Breeding Method ZOLTAN TOTH* AND EUGENIA KALNAY Environmental Modeling Center, National Centers for Environmental Prediction, Camp Springs, Maryland (Manuscript
More information2. Outline of the MRI-EPS
2. Outline of the MRI-EPS The MRI-EPS includes BGM cycle system running on the MRI supercomputer system, which is developed by using the operational one-month forecasting system by the Climate Prediction
More informationData assimilation for the coupled ocean-atmosphere
GODAE Ocean View/WGNE Workshop 2013 19 March 2013 Data assimilation for the coupled ocean-atmosphere Eugenia Kalnay, Tamara Singleton, Steve Penny, Takemasa Miyoshi, Jim Carton Thanks to the UMD Weather-Chaos
More informationABSTRACT. The theme of my thesis research is to perform breeding experiments with
ABSTRACT Title of Dissertation / Thesis: BRED VECTORS IN THE NASA NSIPP GLOBAL COUPLED MODEL AND THEIR APPLICATION TO COUPLED ENSEMBLE PREDICTIONS AND DATA ASSIMILATION. Shu-Chih Yang, Doctor of Philosophy,
More informationM.Sc. in Meteorology. Numerical Weather Prediction Prof Peter Lynch
M.Sc. in Meteorology UCD Numerical Weather Prediction Prof Peter Lynch Meteorology & Climate Centre School of Mathematical Sciences University College Dublin Second Semester, 2005 2006. In this section
More informationA Truncated Model for Finite Amplitude Baroclinic Waves in a Channel
A Truncated Model for Finite Amplitude Baroclinic Waves in a Channel Zhiming Kuang 1 Introduction To date, studies of finite amplitude baroclinic waves have been mostly numerical. The numerical models,
More informationSensitivity of Forecast Errors to Initial Conditions with a Quasi-Inverse Linear Method
2479 Sensitivity of Forecast Errors to Initial Conditions with a Quasi-Inverse Linear Method ZHAO-XIA PU Lanzhou University, Lanzhon, People s Republic of China, and UCAR Visiting Scientist, National Centers
More informationLocal Predictability of the Performance of an. Ensemble Forecast System
Local Predictability of the Performance of an Ensemble Forecast System Elizabeth Satterfield Istvan Szunyogh University of Maryland, College Park, Maryland To be submitted to JAS Corresponding author address:
More informationAccelerating the spin-up of Ensemble Kalman Filtering
Accelerating the spin-up of Ensemble Kalman Filtering Eugenia Kalnay * and Shu-Chih Yang University of Maryland Abstract A scheme is proposed to improve the performance of the ensemble-based Kalman Filters
More informationTESTING GEOMETRIC BRED VECTORS WITH A MESOSCALE SHORT-RANGE ENSEMBLE PREDICTION SYSTEM OVER THE WESTERN MEDITERRANEAN
TESTING GEOMETRIC BRED VECTORS WITH A MESOSCALE SHORT-RANGE ENSEMBLE PREDICTION SYSTEM OVER THE WESTERN MEDITERRANEAN Martín, A. (1, V. Homar (1, L. Fita (1, C. Primo (2, M. A. Rodríguez (2 and J. M. Gutiérrez
More informationNonlinear fastest growing perturbation and the first kind of predictability
Vol. 44 No. SCIENCE IN CHINA (Series D) December Nonlinear fastest growing perturbation and the first kind of predictability MU Mu ( ) & WANG Jiacheng ( ) LASG, Institute of Atmospheric Physics, Chinese
More informationUse of the breeding technique to estimate the structure of the analysis errors of the day
Nonlinear Processes in Geophysics (2003) 10: 1 11 Nonlinear Processes in Geophysics c European Geosciences Union 2003 Use of the breeding technique to estimate the structure of the analysis errors of the
More informationRISE undergraduates find that regime changes in Lorenz s model are predictable
RISE undergraduates find that regime changes in Lorenz s model are predictable Erin Evans (1), Nadia Bhatti (1), Jacki Kinney (1,4), Lisa Pann (1), Malaquias Peña (2), Shu-Chih Yang (2), Eugenia Kalnay
More informationWork at JMA on Ensemble Forecast Methods
Work at JMA on Ensemble Forecast Methods Hitoshi Sato 1, Shuhei Maeda 1, Akira Ito 1, Masayuki Kyouda 1, Munehiko Yamaguchi 1, Ryota Sakai 1, Yoshimitsu Chikamoto 2, Hitoshi Mukougawa 3 and Takuji Kubota
More informationSTOCHASTICALLY PERTURBED BRED VECTORS IN MULTI-SCALE SYSTEMS BRENT GIGGINS AND GEORG A. GOTTWALD
STOCHASTICALLY PERTURBED BRED VECTORS IN MULTI-SCALE SYSTEMS BRENT GIGGINS AND GEORG A. GOTTWALD School of Mathematics and Statistics, University of Sydney, NSW 2006, Australia Abstract. The breeding method
More informationEvolution of Analysis Error and Adjoint-Based Sensitivities: Implications for Adaptive Observations
1APRIL 2004 KIM ET AL. 795 Evolution of Analysis Error and Adjoint-Based Sensitivities: Implications for Adaptive Observations HYUN MEE KIM Marine Meteorology and Earthquake Research Laboratory, Meteorological
More informationExploring and extending the limits of weather predictability? Antje Weisheimer
Exploring and extending the limits of weather predictability? Antje Weisheimer Arnt Eliassen s legacy for NWP ECMWF is an independent intergovernmental organisation supported by 34 states. ECMWF produces
More informationChapter 2. The continuous equations
Chapter. The continuous equations Fig. 1.: Schematic of a forecast with slowly varying weather-related variations and superimposed high frequency Lamb waves. Note that even though the forecast of the slow
More informationEnsemble Kalman Filter potential
Ensemble Kalman Filter potential Former students (Shu-Chih( Yang, Takemasa Miyoshi, Hong Li, Junjie Liu, Chris Danforth, Ji-Sun Kang, Matt Hoffman), and Eugenia Kalnay University of Maryland Acknowledgements:
More informationClimate Change and Predictability of the Indian Summer Monsoon
Climate Change and Predictability of the Indian Summer Monsoon B. N. Goswami (goswami@tropmet.res.in) Indian Institute of Tropical Meteorology, Pune Annual mean Temp. over India 1875-2004 Kothawale, Roopakum
More informationVortices in accretion discs: formation process and dynamical evolution
Vortices in accretion discs: formation process and dynamical evolution Geoffroy Lesur DAMTP (Cambridge UK) LAOG (Grenoble) John Papaloizou Sijme-Jan Paardekooper Giant vortex in Naruto straight (Japan)
More informationWeight interpolation for efficient data assimilation with the Local Ensemble Transform Kalman Filter
Weight interpolation for efficient data assimilation with the Local Ensemble Transform Kalman Filter Shu-Chih Yang 1,2, Eugenia Kalnay 1,3, Brian Hunt 1,3 and Neill E. Bowler 4 1 Department of Atmospheric
More informationWhat is a Low Order Model?
What is a Low Order Model? t Ψ = NL(Ψ ), where NL is a nonlinear operator (quadratic nonlinearity) N Ψ (x,y,z,...,t)= Ai (t)φ i (x,y,z,...) i=-n da i = N N cijk A j A k + bij A j + f i v i j;k=-n j=-n
More informationContents. Parti Fundamentals. 1. Introduction. 2. The Coriolis Force. Preface Preface of the First Edition
Foreword Preface Preface of the First Edition xiii xv xvii Parti Fundamentals 1. Introduction 1.1 Objective 3 1.2 Importance of Geophysical Fluid Dynamics 4 1.3 Distinguishing Attributes of Geophysical
More informationInitial Uncertainties in the EPS: Singular Vector Perturbations
Initial Uncertainties in the EPS: Singular Vector Perturbations Training Course 2013 Initial Uncertainties in the EPS (I) Training Course 2013 1 / 48 An evolving EPS EPS 1992 2010: initial perturbations
More informationInitial ensemble perturbations - basic concepts
Initial ensemble perturbations - basic concepts Linus Magnusson Acknowledgements: Erland Källén, Martin Leutbecher,, Slide 1 Introduction Perturbed forecast Optimal = analysis + Perturbation Perturbations
More informationINTERCOMPARISON OF THE CANADIAN, ECMWF, AND NCEP ENSEMBLE FORECAST SYSTEMS. Zoltan Toth (3),
INTERCOMPARISON OF THE CANADIAN, ECMWF, AND NCEP ENSEMBLE FORECAST SYSTEMS Zoltan Toth (3), Roberto Buizza (1), Peter Houtekamer (2), Yuejian Zhu (4), Mozheng Wei (5), and Gerard Pellerin (2) (1) : European
More informationEnsemble Assimilation of Global Large-Scale Precipitation
Ensemble Assimilation of Global Large-Scale Precipitation Guo-Yuan Lien 1,2 in collaboration with Eugenia Kalnay 2, Takemasa Miyoshi 1,2 1 RIKEN Advanced Institute for Computational Science 2 University
More informationComparison of 3D-Var and LETKF in an Atmospheric GCM: SPEEDY
Comparison of 3D-Var and LEKF in an Atmospheric GCM: SPEEDY Catherine Sabol Kayo Ide Eugenia Kalnay, akemasa Miyoshi Weather Chaos, UMD 9 April 2012 Outline SPEEDY Formulation Single Observation Eperiments
More information2.5 Shallow water equations, quasigeostrophic filtering, and filtering of inertia-gravity waves
Chapter. The continuous equations φ=gh Φ=gH φ s =gh s Fig..5: Schematic of the shallow water model, a hydrostatic, incompressible fluid with a rigid bottom h s (x,y), a free surface h(x,y,t), and horizontal
More informationWeight interpolation for efficient data assimilation with the Local Ensemble Transform Kalman Filter
QUARTERLY JOURNAL OF THE ROYAL METEOROLOGICAL SOCIETY Published online 17 December 2008 in Wiley InterScience (www.interscience.wiley.com).353 Weight interpolation for efficient data assimilation with
More informationAn implementation of the Local Ensemble Kalman Filter in a quasi geostrophic model and comparison with 3D-Var
Author(s) 2007. This work is licensed under a Creative Commons License. Nonlinear Processes in Geophysics An implementation of the Local Ensemble Kalman Filter in a quasi geostrophic model and comparison
More informationData Assimilation: Finding the Initial Conditions in Large Dynamical Systems. Eric Kostelich Data Mining Seminar, Feb. 6, 2006
Data Assimilation: Finding the Initial Conditions in Large Dynamical Systems Eric Kostelich Data Mining Seminar, Feb. 6, 2006 kostelich@asu.edu Co-Workers Istvan Szunyogh, Gyorgyi Gyarmati, Ed Ott, Brian
More informationThe Local Ensemble Transform Kalman Filter (LETKF) Eric Kostelich. Main topics
The Local Ensemble Transform Kalman Filter (LETKF) Eric Kostelich Arizona State University Co-workers: Istvan Szunyogh, Brian Hunt, Ed Ott, Eugenia Kalnay, Jim Yorke, and many others http://www.weatherchaos.umd.edu
More informationEstimating Lyapunov Exponents from Time Series. Gerrit Ansmann
Estimating Lyapunov Exponents from Time Series Gerrit Ansmann Example: Lorenz system. Two identical Lorenz systems with initial conditions; one is slightly perturbed (10 14 ) at t = 30: 20 x 2 1 0 20 20
More informationEvidence of growing bred vector associated with the tropical intraseasonal oscillation
Click Here for Full Article GEOPHYSICAL RESEARCH LETTERS, VOL. 34, L04806, doi:10.1029/2006gl028450, 2007 Evidence of growing bred vector associated with the tropical intraseasonal oscillation Yoshimitsu
More informationIn two-dimensional barotropic flow, there is an exact relationship between mass
19. Baroclinic Instability In two-dimensional barotropic flow, there is an exact relationship between mass streamfunction ψ and the conserved quantity, vorticity (η) given by η = 2 ψ.the evolution of the
More information+ ω = 0, (1) (b) In geometric height coordinates in the rotating frame of the Earth, momentum balance for an inviscid fluid is given by
Problem Sheet 1: Due Thurs 3rd Feb 1. Primitive equations in different coordinate systems (a) Using Lagrangian considerations and starting from an infinitesimal mass element in cartesian coordinates (x,y,z)
More informationp = ρrt p = ρr d = T( q v ) dp dz = ρg
Chapter 1: Properties of the Atmosphere What are the major chemical components of the atmosphere? Atmospheric Layers and their major characteristics: Troposphere, Stratosphere Mesosphere, Thermosphere
More informationOrganization. I MCMC discussion. I project talks. I Lecture.
Organization I MCMC discussion I project talks. I Lecture. Content I Uncertainty Propagation Overview I Forward-Backward with an Ensemble I Model Reduction (Intro) Uncertainty Propagation in Causal Systems
More informationLateral Boundary Conditions
Lateral Boundary Conditions Introduction For any non-global numerical simulation, the simulation domain is finite. Consequently, some means of handling the outermost extent of the simulation domain its
More informationGFD II: Balance Dynamics ATM S 542
GFD II: Balance Dynamics ATM S 542 DARGAN M. W. FRIERSON UNIVERSITY OF WASHINGTON, DEPARTMENT OF ATMOSPHERIC SCIENCES WEEK 8 SLIDES Eady Model Growth rates (imaginary part of frequency) Stable for large
More informationChapter 3. Stability theory for zonal flows :formulation
Chapter 3. Stability theory for zonal flows :formulation 3.1 Introduction Although flows in the atmosphere and ocean are never strictly zonal major currents are nearly so and the simplifications springing
More informationLocal Ensemble Transform Kalman Filter
Local Ensemble Transform Kalman Filter Brian Hunt 11 June 2013 Review of Notation Forecast model: a known function M on a vector space of model states. Truth: an unknown sequence {x n } of model states
More informationData assimilation : Basics and meteorology
Data assimilation : Basics and meteorology Olivier Talagrand Laboratoire de Météorologie Dynamique, École Normale Supérieure, Paris, France Workshop on Coupled Climate-Economics Modelling and Data Analysis
More informationA review of ensemble forecasting techniques with a focus on tropical cyclone forecasting
Meteorol. Appl. 8, 315 332 (2001) A review of ensemble forecasting techniques with a focus on tropical cyclone forecasting Kevin K W Cheung, Department of Physics and Materials Science, City University
More informationP3.11 A COMPARISON OF AN ENSEMBLE OF POSITIVE/NEGATIVE PAIRS AND A CENTERED SPHERICAL SIMPLEX ENSEMBLE
P3.11 A COMPARISON OF AN ENSEMBLE OF POSITIVE/NEGATIVE PAIRS AND A CENTERED SPHERICAL SIMPLEX ENSEMBLE 1 INTRODUCTION Xuguang Wang* The Pennsylvania State University, University Park, PA Craig H. Bishop
More informationBALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere. Potential temperature θ. Rossby Ertel potential vorticity
BALANCED FLOW: EXAMPLES (PHH lecture 3) Potential Vorticity in the real atmosphere Need to introduce a new measure of the buoyancy Potential temperature θ In a compressible fluid, the relevant measure
More informationObservability, a Problem in Data Assimilation
Observability, Data Assimilation with the Extended Kalman Filter 1 Observability, a Problem in Data Assimilation Chris Danforth Department of Applied Mathematics and Scientific Computation, UMD March 10,
More informationLagrangian Analysis of 2D and 3D Ocean Flows from Eulerian Velocity Data
Flows from Second-year Ph.D. student, Applied Math and Scientific Computing Project Advisor: Kayo Ide Department of Atmospheric and Oceanic Science Center for Scientific Computation and Mathematical Modeling
More informationNonlinear variational inverse problems
6 Nonlinear variational inverse problems This chapter considers highly nonlinear variational inverse problems and their properties. More general inverse formulations for nonlinear dynamical models will
More informationSimple Mathematical, Dynamical Stochastic Models Capturing the Observed Diversity of the El Niño Southern Oscillation (ENSO)
Simple Mathematical, Dynamical Stochastic Models Capturing the Observed Diversity of the El Niño Southern Oscillation (ENSO) Lecture 5: A Simple Stochastic Model for El Niño with Westerly Wind Bursts Andrew
More informationLecture 2 ENSO toy models
Lecture 2 ENSO toy models Eli Tziperman 2.3 A heuristic derivation of a delayed oscillator equation Let us consider first a heuristic derivation of an equation for the sea surface temperature in the East
More informationThe impacts of stochastic noise on climate models
The impacts of stochastic noise on climate models Paul Williams Department of Meteorology, University of Reading, UK The impacts of στοχαστικός noise on climate models Paul Williams Department of Meteorology,
More informationEnKF Localization Techniques and Balance
EnKF Localization Techniques and Balance Steven Greybush Eugenia Kalnay, Kayo Ide, Takemasa Miyoshi, and Brian Hunt Weather Chaos Meeting September 21, 2009 Data Assimilation Equation Scalar form: x a
More informationBalance. in the vertical too
Balance. in the vertical too Gradient wind balance f n Balanced flow (no friction) More complicated (3- way balance), however, better approximation than geostrophic (as allows for centrifugal acceleration
More informationAFRICAN EASTERLY WAVES IN CURRENT AND FUTURE CLIMATES
AFRICAN EASTERLY WAVES IN CURRENT AND FUTURE CLIMATES Victoria Dollar RTG Seminar Research - Spring 2018 April 16, 2018 Victoria Dollar ASU April 16, 2018 1 / 26 Overview Introduction Rossby waves and
More informationthat individual/local amplitudes of Ro can reach O(1).
Supplementary Figure. (a)-(b) As Figures c-d but for Rossby number Ro at the surface, defined as the relative vorticity ζ divided by the Coriolis frequency f. The equatorial band (os-on) is not shown due
More information2D.4 THE STRUCTURE AND SENSITIVITY OF SINGULAR VECTORS ASSOCIATED WITH EXTRATROPICAL TRANSITION OF TROPICAL CYCLONES
2D.4 THE STRUCTURE AND SENSITIVITY OF SINGULAR VECTORS ASSOCIATED WITH EXTRATROPICAL TRANSITION OF TROPICAL CYCLONES Simon T. Lang Karlsruhe Institute of Technology. INTRODUCTION During the extratropical
More informationHigh initial time sensitivity of medium range forecasting observed for a stratospheric sudden warming
GEOPHYSICAL RESEARCH LETTERS, VOL. 37,, doi:10.1029/2010gl044119, 2010 High initial time sensitivity of medium range forecasting observed for a stratospheric sudden warming Yuhji Kuroda 1 Received 27 May
More informationThe general circulation: midlatitude storms
The general circulation: midlatitude storms Motivation for this class Provide understanding basic motions of the atmosphere: Ability to diagnose individual weather systems, and predict how they will change
More informationThe Impact of Background Error on Incomplete Observations for 4D-Var Data Assimilation with the FSU GSM
The Impact of Background Error on Incomplete Observations for 4D-Var Data Assimilation with the FSU GSM I. Michael Navon 1, Dacian N. Daescu 2, and Zhuo Liu 1 1 School of Computational Science and Information
More informationTraveling planetary-scale Rossby waves in the winter stratosphere: The role of tropospheric baroclinic instability
GEOPHYSICAL RESEARCH LETTERS, VOL. 39,, doi:10.1029/2012gl053684, 2012 Traveling planetary-scale Rossby waves in the winter stratosphere: The role of tropospheric baroclinic instability Daniela I. V. Domeisen
More informationGoals of this Chapter
Waves in the Atmosphere and Oceans Restoring Force Conservation of potential temperature in the presence of positive static stability internal gravity waves Conservation of potential vorticity in the presence
More informationIntroduction of Seasonal Forecast Guidance. TCC Training Seminar on Seasonal Prediction Products November 2013
Introduction of Seasonal Forecast Guidance TCC Training Seminar on Seasonal Prediction Products 11-15 November 2013 1 Outline 1. Introduction 2. Regression method Single/Multi regression model Selection
More informationNCEP ENSEMBLE FORECAST SYSTEMS
NCEP ENSEMBLE FORECAST SYSTEMS Zoltan Toth Environmental Modeling Center NOAA/NWS/NCEP Acknowledgements: Y. Zhu, R. Wobus, M. Wei, D. Hou, G. Yuan, L. Holland, J. McQueen, J. Du, B. Zhou, H.-L. Pan, and
More informationInstability of the Chaotic ENSO: The Growth-Phase Predictability Barrier
3613 Instability of the Chaotic ENSO: The Growth-Phase Predictability Barrier ROGER M. SAMELSON College of Oceanic and Atmospheric Sciences, Oregon State University, Corvallis, Oregon ELI TZIPERMAN Weizmann
More informationRevision of TR-09-25: A Hybrid Variational/Ensemble Filter Approach to Data Assimilation
Revision of TR-9-25: A Hybrid Variational/Ensemble ilter Approach to Data Assimilation Adrian Sandu 1 and Haiyan Cheng 1 Computational Science Laboratory Department of Computer Science Virginia Polytechnic
More informationQuantifying Uncertainty through Global and Mesoscale Ensembles
Quantifying Uncertainty through Global and Mesoscale Ensembles Teddy R. Holt Naval Research Laboratory Monterey CA 93943-5502 phone: (831) 656-4740 fax: (831) 656-4769 e-mail: holt@nrlmry.navy.mil Award
More informationHandling nonlinearity in Ensemble Kalman Filter: Experiments with the three-variable Lorenz model
Handling nonlinearity in Ensemble Kalman Filter: Experiments with the three-variable Lorenz model Shu-Chih Yang 1*, Eugenia Kalnay, and Brian Hunt 1. Department of Atmospheric Sciences, National Central
More informationRotating stratified turbulence in the Earth s atmosphere
Rotating stratified turbulence in the Earth s atmosphere Peter Haynes, Centre for Atmospheric Science, DAMTP, University of Cambridge. Outline 1. Introduction 2. Momentum transport in the atmosphere 3.
More informationPredictability of large-scale atmospheric motions: Lyapunov exponents and error dynamics
Accepted in Chaos Predictability of large-scale atmospheric motions: Lyapunov exponents and error dynamics Stéphane Vannitsem Royal Meteorological Institute of Belgium Meteorological and Climatological
More informationCOMPARISON OF THE INFLUENCES OF INITIAL ERRORS AND MODEL PARAMETER ERRORS ON PREDICTABILITY OF NUMERICAL FORECAST
CHINESE JOURNAL OF GEOPHYSICS Vol.51, No.3, 2008, pp: 718 724 COMPARISON OF THE INFLUENCES OF INITIAL ERRORS AND MODEL PARAMETER ERRORS ON PREDICTABILITY OF NUMERICAL FORECAST DING Rui-Qiang, LI Jian-Ping
More informationPart-8c Circulation (Cont)
Part-8c Circulation (Cont) Global Circulation Means of Transfering Heat Easterlies /Westerlies Polar Front Planetary Waves Gravity Waves Mars Circulation Giant Planet Atmospheres Zones and Belts Global
More informationTargeted Observations of Tropical Cyclones Based on the
Targeted Observations of Tropical Cyclones Based on the Adjoint-Derived Sensitivity Steering Vector Chun-Chieh Wu, Po-Hsiung Lin, Jan-Huey Chen, and Kun-Hsuan Chou Department of Atmospheric Sciences, National
More informationMesoscale Atmospheric Systems. Surface fronts and frontogenesis. 06 March 2018 Heini Wernli. 06 March 2018 H. Wernli 1
Mesoscale Atmospheric Systems Surface fronts and frontogenesis 06 March 2018 Heini Wernli 06 March 2018 H. Wernli 1 Temperature (degc) Frontal passage in Mainz on 26 March 2010 06 March 2018 H. Wernli
More informationLocal Ensemble Transform Kalman Filter: An Efficient Scheme for Assimilating Atmospheric Data
Local Ensemble Transform Kalman Filter: An Efficient Scheme for Assimilating Atmospheric Data John Harlim and Brian R. Hunt Department of Mathematics and Institute for Physical Science and Technology University
More informationNecklace solitary waves on bounded domains
Necklace solitary waves on bounded domains with Dima Shpigelman Gadi Fibich Tel Aviv University 1 Physical setup r=(x,y) laser beam z z=0 Kerr medium Intense laser beam that propagates in a transparent
More information