African Easterly Waves, Observations and Theory

Transcription

1 African Easterly Waves, Observations and Theory Nick Hall LEGOS / Univ. Toulouse Stephanie Leroux LTHE / Univ. Grenoble George Kiladis ESRL NOAA, Boulder Chris Thorncroft - SUNY, Albany

2 Three questions How do easterly waves grow? How are they initiated? How do they interact with convection?

3 The birth of hurricane Alberto

4 The West African monsoon TEJ AEJ Hadley Cell E H R Θe

5 Climatology of transient disturbances shading: OLR contours: wind and standard deviation of OLR and vorticity

6 <6-day Perturbation kinetic energy: day filtered (shading) P KE = u 2 + v mb 850 mb (a) mean PKE 700 mb and EOF1 (8.6 %) (b) mean PKE 850 mb and EOF1 (9.1 %) EOF EOF2 (5.1 %) - EOF2 (4.6 %) EOF fig3

7 Zonal wind 600 mb June-Sept NCEP (a) Mean Jun-Sep. non filt. zonal wind 600 mb Climatology (b) Jun-Sep. non filt. zonal wind 600 mb EOF1 (16.2 %) (c) Jun-Sep d filt. zonal wind 600 mb EOF1 (6.6 %) (d) Jun-Sep d filt. zonal wind 600 mb EOF1 (10.2 %) EOF2 (6.6 %) EOF2 (5.7 %) EOF2 (8.5 %) Unfiltered EOFs day filtered EOFs day filtered EOFs fig 1

8 Space-Time Spectrum of JJA Antisymmetric OLR, 15S-15N Wheeler and Kiladis 1999

9 Regression Analysis OLR and 850mb streamfunction regressed against space-time filtered OLR (TD band) at 10N,10W : June-Sept Day -5 OLR anomaly > < 10 W/m2, streamfunction ci=2x105 m2/s

10 Regression Analysis OLR and 850mb streamfunction regressed against space-time filtered OLR (TD band) at 10N,10W : June-Sept Day -5-4 OLR anomaly > < 10 W/m2, streamfunction ci=2x105 m2/s

11 Regression Analysis OLR and 850mb streamfunction regressed against space-time filtered OLR (TD band) at 10N,10W : June-Sept Day OLR anomaly > < 10 W/m2, streamfunction ci=2x105 m2/s

12 Regression Analysis OLR and 850mb streamfunction regressed against space-time filtered OLR (TD band) at 10N,10W : June-Sept Day OLR anomaly > < 10 W/m2, streamfunction ci=2x105 m2/s

13 Regression Analysis OLR and 850mb streamfunction regressed against space-time filtered OLR (TD band) at 10N,10W : June-Sept Day OLR anomaly > < 10 W/m2, streamfunction ci=2x105 m2/s

14 Regression Analysis OLR and 850mb streamfunction regressed against space-time filtered OLR (TD band) at 10N,10W : June-Sept Day OLR anomaly > < 10 W/m2, streamfunction ci=2x105 m2/s

15 Regression Analysis OLR and 850mb streamfunction regressed against space-time filtered OLR (TD band) at 10N,10W : June-Sept Day OLR anomaly > < 10 W/m2, streamfunction ci=2x105 m2/s

16 Regression Analysis OLR and 850mb streamfunction regressed against space-time filtered OLR (TD band) at 10N,10W : June-Sept Day OLR anomaly > < 10 W/m2, streamfunction ci=2x105 m2/s

17 Regression Analysis OLR and 850mb streamfunction regressed against space-time filtered OLR (TD band) at 10N,10W : June-Sept Day OLR anomaly > < 10 W/m2, streamfunction ci=2x105 m2/s

18 Regression Analysis OLR and 850mb streamfunction regressed against space-time filtered OLR (TD band) at 10N,10W : June-Sept Day OLR anomaly > < 10 W/m2, streamfunction ci=2x105 m2/s

19 Regression Analysis OLR and 850mb streamfunction regressed against space-time filtered OLR (TD band) at 10N,10W : June-Sept Day OLR anomaly > < 10 W/m2, streamfunction ci=2x105 m2/s

20

21 NCEP2 2-6d filtered V700mb JJAS day filtered 700 mb meridional wind - NCEP day filtered PKE NCEP2 PKE d d filtered filtered pke700mb JJAS.2006 JJAS P KE = u 2 + v longitude longitude filtered V700mb (m/s) filtered PKE (m^2.s^-2)

22 A brief history of ideas on African transient disturbances In the beginning, there was GATE: Eady/Charney, Burpee,... Thorncroft and Hoskins... etc.

23 A brief history of ideas on African transient disturbances In the beginning, there was GATE: Eady/Charney, Burpee,... Thorncroft and Hoskins... etc.

24 A brief history of ideas on African transient disturbances In the beginning, there was GATE: Eady/Charney, Burpee,... Thorncroft and Hoskins... etc.

25 A brief history of ideas on African transient disturbances In the beginning, there was GATE: Eady/Charney, Burpee,... Thorncroft and Hoskins... etc. then came Stone. Farrel, etc - and their ideas applied to the African easterly jet:

26 A brief history of ideas on African transient disturbances In the beginning, there was GATE: Eady/Charney, Burpee,... Thorncroft and Hoskins... etc. then came Stone. Farrel, etc - and their ideas applied to the African easterly jet:

27 A brief history of ideas on African transient disturbances In the beginning, there was GATE: Eady/Charney, Burpee,... Thorncroft and Hoskins... etc. then came Stone. Farrel, etc - and their ideas applied to the African easterly jet:

28 A brief history of ideas on African transient disturbances In the beginning, there was GATE: Eady/Charney, Burpee,... Thorncroft and Hoskins... etc. then came Stone. Farrel, etc - and their ideas applied to the African easterly jet:

29 A brief history of ideas on African transient disturbances In the beginning, there was GATE: Eady/Charney, Burpee,... Thorncroft and Hoskins... etc. then came Stone. Farrel, etc - and their ideas applied to the African easterly jet: but the AEJ is also short, and variable

30 A brief history of ideas on African transient disturbances In the beginning, there was GATE: Eady/Charney, Burpee,... Thorncroft and Hoskins... etc. then came Stone. Farrel, etc - and their ideas applied to the African easterly jet: but the AEJ is also short, and variable

31 A brief history of ideas on African transient disturbances In the beginning, there was GATE: Eady/Charney, Burpee,... Thorncroft and Hoskins... etc. then came Stone. Farrel, etc - and their ideas applied to the African easterly jet: but the AEJ is also short, and variable

32 A brief history of ideas on African transient disturbances In the beginning, there was GATE: Eady/Charney, Burpee,... Thorncroft and Hoskins... etc. then came Stone. Farrel, etc - and their ideas applied to the African easterly jet: but the AEJ is also short, and variable

33 A brief history of ideas on African transient disturbances In the beginning, there was GATE: Eady/Charney, Burpee,... Thorncroft and Hoskins... etc. then came Stone. Farrel, etc - and their ideas applied to the African easterly jet: but the AEJ is also short, and variable - necessary conditions for African easterly waves? - AEW feedback on AEJ? - it starts to get complicated...

34 Possible Jet - Wave relationships: instability hypothesis simple triggering T T variable jet triggering T T can we see these ideas at work in the observations?

35 Regression of day 600 mb zonal wind (m/s) onto 700 mb PKE index Day -6 Day -4 Day DAY 0 Day+2 Day+4 Day m/s 1 N Σα i X i 1 N Σα2 i αi is the PKE index in the box Xi is the zonal wind fig8

36 E vectors Regression of day E-vectors onto 700 mb PKE index (vectors ~ 1 m 2 /s 2, contours of div E every 10-6 ms -2 ) ) E = (u 2 v 2, u v E vectors converge where easterlies are accelerated by transients. An eastward pointing E vector has north-south elongated eddies. A westward pointing E vector has east-west elongated eddies. A southward pointing E vector has NE-SW oriented eddies

37 E vectors Regression of day E-vectors onto 700 mb PKE index (vectors ~ 1 m 2 /s 2, contours of div E every 10-6 ms -2 ) ) E = (u 2 v 2, u v E vectors converge where easterlies are accelerated by transients. An eastward pointing E vector has north-south elongated eddies. A westward pointing E vector has east-west elongated eddies. A southward pointing E vector has NE-SW oriented eddies

38 E vectors Regression of day E-vectors onto 700 mb PKE index (vectors ~ 1 m 2 /s 2, contours of div E every 10-6 ms -2 ) ) E = (u 2 v 2, u v E vectors converge where easterlies are accelerated by transients. An eastward pointing E vector has north-south elongated eddies. A westward pointing E vector has east-west elongated eddies. A southward pointing E vector has NE-SW oriented eddies

39 E vectors Regression of day E-vectors onto 700 mb PKE index (vectors ~ 1 m 2 /s 2, contours of div E every 10-6 ms -2 ) ) E = (u 2 v 2, u v E vectors converge where easterlies are accelerated by transients. An eastward pointing E vector has north-south elongated eddies. A westward pointing E vector has east-west elongated eddies. A southward pointing E vector has NE-SW oriented eddies

40 PV gradient Regression of day 320K PV gradient onto 700 mb PKE index (contours K m Kg -1 s -1 ) Day -6 Day -4 Day -2 DAY 0 Day +2 Day +4 Day fig10

41 Relationship between convection and transients Regression of day NOAA OLR (W m -2 ) onto 700 mb PKE index Day -6 Day -4 Day -2 DAY 0 Day+2 Day+4 Day fig11

42 Try it the other way round Regressions onto OLR index (a) d filt. OLR Day -6 Day -4 Day -2 DAY 0 Day+2 Day+4 Day PKE (b) d filt. PKE at 700 mb Day -6 Day -4 Day -2 DAY 0 Day Day Day fig12

43 Conclusions so far... - Prior to increased transient activity, the upstream AEJ is reinforced. - Afterwards, the downstream AEJ is shifted to the north. - Eddy momentum forcing from transients is consistent with this AEJ development. - Increased transients over W. Africa are associated with increased convection on intraseasonal timescales. - Increased transient activity tends to follow peaks in convective activity in an upstream location.

44 generalized development Generic equation for the development of a nonlinear system q t = Lq.vq + f introduce a state vector Φ =( v, q,...) this becomes dφ dt = LΦ +Φ QΦ + f where Φ is the diagonal matrix that contains the elements of Φ

45 basic state separation: 1 specify a basic state that is a solution of the equations Φ =Φ 0 + Φ 1 dφ 0 f = f 0 + f 1 dt and LΦ 0 + Φ 0 QΦ 0 + f 0 =0 so =0 dφ dt = LΦ +Φ QΦ + f dφ 1 dt = LΦ 1 + Φ 1 QΦ 0 + Φ 0 QΦ 1 + Φ 1 QΦ 1 + f 1 or to put it another way dφ 1 dt = L 0 Φ 1 + O(Φ 2 1)+f 1 ADVANTAGE: eddy development independent of separation DISADVANTAGE: basic state unrealistic so nonlinear term large - linearization of questionable relevance

46 basic state separation: 2 specify a realistic basic state (for example, the time mean flow) but this time in fact we have and Φ = Φ +Φ f = f + f LΦ + Φ QΦ + f 0 dφ dt =0 mean advection LΦ + Φ QΦ + Φ QΦ + f =0 dφ transient eddy forcing dt = LΦ + Φ QΦ + Φ QΦ + dφ dt = LΦ +Φ QΦ + f [Φ QΦ Φ QΦ ] ADVANTAGE: realistic basic state so meaningful linearization possible DISADVANTAGE: linear development equation not independent of time mean transient forcing + f or to put it another way dφ dt = L meanφ + [ ] O(Φ 2 ) O(Φ 2 ) + f

47 a perturbation model So much for theory. How do we solve these equations? We use a dynamical model. We can appeal to data to deduce the appropriate forcing functions: data dφ dt = LΦ +Φ QΦ + f(t) model dψ dt = LΨ +Ψ QΨ + g We now define g using data, so that if we initialize the model with Φ it will not develop. so g = LΦ Φ QΦ We can easly find g by integrating the unforced model from Φ for just one timestep. From the time-mean budget equation we also see that this definition of g gives g = f + Φ QΦ So this forcing represents the time mean diabatic forcing plus the mean transient eddy forcing. These are the two processes that maintain the time-mean circulation.

48 stability analysis With our data-derived forcing, for small perturbations, integrating dψ dt = LΨ +Ψ QΨ + g is equivalent to integrating dψ dt = L mean Ψ We use the dynamical model to analyse the linear growth problem for normal modes of Lmean L mean e n (x, y, z) =λ n e n (x, y, z) for a single mode λ n = σ + iω Ψ n = e n (x, y, z)e (σ+iω)t in general en is complex so the solution takes the form Ψ =[A(x, y, z) sin ωt + B(x, y, z) cos ωt] e σt

49 response to forcing We can use our dynamical model to find the response to a perturbation forcing f dψ dt = L mean Ψ + f (and if we keep f small the response is linear) Start with another example from African easterly waves. This time we use a (convective) heating anomaly as f, to trigger a response. The response still looks like the normal mode that we found before. But it decays in time.

50 explicit transients: a simple GCM Let s reconsider the definition of our forcing function g. Recall the development equations: data dφ dt = LΦ +Φ QΦ + f(t) model dψ dt = LΨ +Ψ QΨ + g If we set g = f g = LΦ Φ QΦ Φ QΦ then this is the same as setting i.e. we have subtracted out the time-mean forcing due to the transients. Again, we can calculate this forcing by initializing the unforced model from a series of values of and then taking the time-average. Φ Φ If the model is now initialized with it will develop in time. In fact we hope it will develop its own explicit transient activity. And we hope that it will be realistic. But there is no guarrantee that this simple GCM will have a realistic climatology. The only thing that is guarranteed is that: LΨ + Ψ QΨ + Ψ QΨ = LΦ + Φ QΦ + Φ QΦ

51 nudge nudge Another way of forcing a model is to push it towards a desired climatology in a restricted region, and look at the effect on the solution outside that region. This is called nudging. dψ dt = LΨ +Ψ QΨ + g + ( ) Φn Ψ Nudging involves an additional constant forcing term and a damping term. In a linear experiment, the appropriate model is: τ dψ dt = L mean Ψ + ɛ ( ) Φn Φ Ψ τ τ This can be useful technique for diagnosing climate anomalies or simulating other people s GCMs with a simple model.

52 Previous work

53 The African Easterly Jet sector mean 60W-60E slice at Greenwich slice at 12N Thorncroft Which (if any) zonal section is appropriate? - Are these basic states unstable? - Is instability important? - Which structures are the most efficient at extracting energy from this flow? 650 mb zonal wind NCEP JJAS mean 68-98

54 Effect of low level damping on modal growth rate Hall and Sardeshmukh 1998

55 Normal Mode Structures: Zonal Basic State Horizontal streamfunction at 850 mb and 12N Greenwich Slice Growth rate: days-1 Period: 4.88 days

56 Normal Mode Structures: Zonal Basic State Horizontal streamfunction at 850 mb and 12N Sector Greenwich Mean 60W-60E Slice Growth rate: days-1 Period: 4.88 days

57 Normal Mode Structures: Zonal Basic State Horizontal streamfunction at 850 mb and 12N Sector Greenwich Mean 60W-60E Slice Growth rate: days-1 Period: days

58 Normal mode structures: 3D basic state Growth rate: days-1 Period: 5.56 days

59 Normal mode structures: 3D basic state Add damping below 800 mb: 2 days on momentum, 4 days on temperature Growth rate: days-1 Period: days

60 Observations Model

61 Fluxes through the looking glass Low Low High C PV W - + W C - Easterly waves transport heat southwards (away from the Saharan heat source). - But they still have convergent momentum transport, so they deplete the African easterly jet.

62 Eddy Fluxes Maritime decay region Wave maximum region Wave inception region 60-30W 30W E northward temperature flux v't' momentum flux u'v'

63 model 0-30W average regression 10W 17.5W 500mb 500mb 500mb 500mb EQ 30N EQ 30N EQ 30N

64 Some wave properties

65 basic state damped modal structure d-1 Some more basic states.33 Growth rates are given for undamped experiments. Damped modes are all close to neutral Active year.17 Quiet year mb zonal wind ci=2 m/s 850 mb streamfunction and vertical motion

66 Initial value problems - If normal mode solutions are neutral they can tell us about efficient structures. - But we still lack a complete theory for the generation and intermittence of AEW events. heat for one day in an elipse with a cosine squared bell shaped distribution and various vertical profiles

67 Initial value problems - If normal mode solutions are neutral they can tell us about efficient structures. - But we still lack a complete theory for the generation and intermittence of AEW events. heat for one day in an elipse with a cosine squared bell shaped distribution and various vertical profiles

68 An influence function for heating experiments shallow convective profile H 0 =5σ 4

69 Sensitivity to size and profile deep convective deep convective with low level evaporation [ ] H 0 = π(1 σ) sin(π(1 σ)) H 0 = 9π2 3π (1 σ) cos (1 σ) 2(3π + 2) 2 shallow convective radius 5 deg shallow convective radius 7 deg H 0 =5σ 4 H 0 =5σ 4

70 Are triggered waves sensitive to a varying easterly jet? 10-day averaged basic states NCEP II reanalysis ( ) >> 336 jets >> 336 perturbation experiments zonal wind at 650 mb

71 Are triggered waves sensitive to a varying easterly jet? 10-day averaged basic states NCEP II reanalysis ( ) >> 336 jets >> 336 perturbation experiments zonal wind at 650 mb

72 Are triggered waves sensitive to a varying easterly jet? 10-day averaged basic states NCEP II reanalysis ( ) >> 336 jets >> 336 perturbation experiments zonal wind at 650 mb day 2 day 5 day 8 stream-function at 850 mb

73 Are triggered waves sensitive to a varying easterly jet? 10-day averaged basic states NCEP II reanalysis ( ) >> 336 jets >> 336 perturbation experiments zonal wind at 650 mb day 2 day 5 day 8 stream-function at 850 mb

74 Are triggered waves sensitive to a varying easterly jet? 10-day averaged basic states NCEP II reanalysis ( ) >> 336 jets >> 336 perturbation experiments zonal wind at 650 mb day 2 day 5 day 8 stream-function at 850 mb

75 Are triggered waves sensitive to a varying easterly jet? 10-day averaged basic states NCEP II reanalysis ( ) >> 336 jets >> 336 perturbation experiments zonal wind at 650 mb day 2 day 5 day 8 stream-function at 850 mb

76 A measure of the intensity of the wave response: stream-function at 850 mb - day 10 Root mean square stream function anomaly at 850 mb, averaged from day 1 to day 18. How does the basic state influence the intensity of the response?

77 Dependance of wave response on spatial EOFs of basic state 650 mb zonal wind? 30 1st mode of variation (35%) EOF1 (35%) !60! nd mode EOF2 of variation (12%) (12%) !60!

78 Dependance of wave response on spatial EOFs of basic state 650 mb zonal wind? 30 1st mode of variation (35%) EOF1 (35%) !60! st Principal Component !50!100! nd Principal Component bottom 20% weakest wave responses (blue dots) top 20% strongest wave responses (orange dots) The strongest waves are mainly associated with jets that are strong in the south and west. 2nd mode EOF2 of variation (12%) (12%) !60!

79 Composite jets for strong and weak responses zonal wind at 650 mb!4! zonal wind at 650 mb!4!8!8 4!4 0!4!8!10!60!40! !4 0!10!60!40! !4 Composite jet for the top 20% strongest wave responses Composite jet for the bottom 20% weakest wave responses

80 Some more conclusions... - Taking account of the jet entrance/exit structure of a realistic AEJ, even in a linear model analysis can take us a lot further towards explaining the 3-d structure of observed AEW composites. - A modest amount of surface damping with a realistic 3-d basic state produces neutral modes. We therefore need to reexamine the traditional instability hypothesis for the timing of AEW growth. - Instead, we look for finite amplitude perturbation events to initiate AEWs. - The strength of the wave response is sensitive to the position of the finite amplitude trigger - convection over the Darfur region is the most efficient. - The response is also very sensitive to the state of the AEJ, as suspected from the observational analysis, although the details differ. however... - The modelling results shown so far do not take into account two-way wavejet interacitons. - Neither do they take into account tropical - extratropical interactions

81 Use the simple GCM (perpetual JJAS mode) June-September NCEP2 reanalyses Simple GCM Climatological zonal wind (contours 5 m/s)

82 < 120 day transient covariances momentum flux c.i. 10/5 m 2 /s 2 temperature flux c.i. 3/2 K m/s EKE c.i. 50/15 m 2 /s 2 Reanalyses Model

83 zoom over the African easterly jet shading shows standard deviation

84 momentum flux c.i. 0.5/0.3 m 2 /s 2 temperature flux c.i. 0.3/0.2 K m/s EKE c.i. 1/0.25 m 2 /s 2

85 eastward propagating disturbances: reanalysis meridional wind anomaly and day filtered EKE (contours)

86 eastward propagating disturbances: model

87 regression analysis from model <120 day filtered streamfunction regressed onto <10 day filtered meridional wind at base point 9.5N 7.5W

88 standard run transients suppressed north of 30N Contours of zonal wind. EKE shaded in green and yellow transients suppressed south of equator transients suppressed north of 30N and south of equator transients suppressed east of 45E in the NH

89 regression analysis from model day filtered variables regressed onto AEW activity index 650 hpa streamfunction zonal wind

90 regression analysis from reanalysis day filtered variables regressed onto AEW activity index 650 hpa streamfunction zonal wind

91 Last set of conclusions - A convection-free model with explicit transient activity can reproduce observed transient activity, albeit with a reduced amplitude. - It can also produce realistic African easterly waves with an amplitude that is stronger than our simulations triggered by convective heat sources. - The easterly waves it produces again resemble linear modes of the AEJ. - However, they disappear entirely if variability is restricted to the tropical band, shedding further doubt on the idea that small perturbations can grow on the AEJ. - The simulated AEWs also show intermittency, without the need for intermittent upstream heating. - Experiments show that the Atlantic storm track is crucial for the existence of these intermittent AEWs. - In the model, statistical precursors appear to originate in the Atlantic storm track, but it is difficult to corroberate this finding from the observations.