UNIVERSITY OF IBADAN LIBRARY ? I ) UNIVERSITY OF IBADAN THIS THESIS SUBMITTED BY ...... .r:.r;.~.. :-!~~~.~~~~¥.f?:.?~~~~ . WAS ACCEPTED IN PARTIAL FULFILMENT OF THE REQUIREMENTS FOR THE DEGREE OF DOCTOR OF PHILOSOPHY IN THE FACULTY OF SCIENCE OF THIS UNIVERSITY THE EFFECTIVE DATE OF THE AWARD IS ................7..t.h...·.a..r..c..h.. , 1983 . ....~~O·J-!~~ . ~~~~~ DATE. ·~ECRETARY ~ POSTGRADUATE SCHOOL UNIVERSITY OF IBADAN LIBRARY INITIAL STAGES AND SOME CHARACTERISTICS OF WEST AFRICAN LINE SQUALLS by JOHN AKINTAYO ADEDOYIN B.Sc.(HONS) PHYSICS (IBADAN) \ I A THESIS IN THE DEPARTMENT OF PHYSICS .> SUBMITTED TO THE FACULTY OF SCIENCE IN PARTIAL FULFIU1ENT OF THE REQUIREMENTS FOR TIlE DEGREE OF DOCTOR OF PHILOSOPHY OF THE UNIVERSITY OF IBADAN MARCH, 1983 UNIVERSITY OF IBADAN LIBRARY , . . , fBAOAN UNIVERSITV LIBRARY NACoC. MLAORCK.. ,U ) j ,. ; UNIVERSITY OF IBADAN LIBRARY -2 - DEDICATION To God for being, as ever, faithful to His promises, my father and late mother without whose prayers and support, this work would have been impossible. UNIVERSITY OF IBADAN LIBRARY - 3 - ABSTRACT A review of West African synoptic weather pattern reveals that the sub-region experiences a special kind of atmospheric disturbance - the line squall - whenever the south-westerlies cover, approximately, the whole of Nigeria. Various methods that have been used to study the squalls (i.e. observational investigations, satellite investi- gations and modelling) have not been very successful in isolating the trigger mechanism of the phenomenon. It is been proposed that line squalls are initiated through the amplification (with time) of any wave-like perturbation along the surface of discontinuity between the south-westerlies and the north-easterlies. The amplifying perturbations could block the 650 mb. mid-tropospheric jet which further distorts the 'bump' formed by the undulating perturbation. This distortion forces the south- westerlies further up and they could condense. The precipitates fall into the underlying, dry jet and some of them evaporate; the latent heat of evaporation being supplied by the jet. The iet, now cooler, sinks. On sinking, the jet could hit the surface of the earth on which it forms the squall front and crawls; thereby lifting the south-westerlies ahead of it. The cycle of condensation, evapo- ration and sinking then continunes. A gravity-wave model of this mechanism is presented through the solution of a frequency equation with the aid of a two-layer atmospheric model. The solution is an eigenvalue problem from which UNIVERSITY OF IBADAN LIBRARY - 4 - many modes of different growth rates and phase velocities could be obtained. Some of these phase velocities ,,;111 be complex - the real part representing the pt::1sevelocity Hhile the imaginary part represents the ampLf.f LcatLon , Waves 'td.ththe largest amplifications (Le. the largest imaginary part) are those th'lt could possibly block the 650 mb. mid-tropospheric jet and trigger off line squalls. Among others, this proposal on the trigger mechanism of line squalls is able to explain: (i) the preference of highlands as places of origin of line squalls, (ii) the close association between the speeds of propagation of line squalls and the mid-tropospheric jet and (iii) the observed overturning of the atmosphere after the passage of line squalls. UNIVERSITY OF IBADAN LIBRARY - 5 - ACKNO\VLEDGEMENT I am grateful to my supervisor, Professor F.B.A. Giwa, who not only laid the foundation of Meteorology in the Department of Physics, University of Ibadan, Nigeria, but also sugp,ested this project and painstakingly guided me from the beginning to the end. I am indebted to the following people for their assistance, in various ways, during and after the trying days of 1979: Professor O. Awe, Dr. C.E.F. Oni and Dr. I.A. Babalola. The University of Ibadan contributed greatly to this project by granting me graduate assistantship for some part of this course and later on, a bursary award. Finally, I most sincerely thank my wife, Bola, for her patience, loyalty and co-operation throughout the period of this work. UNIVERSITY OF IBADAN LIBRARY - 6 - CERTIFICATION This is to certify that the work described in tbis thesis was carried out under my supervision by John Akintayo Adedoyin in the Department of Physics, University of Ibadan, Ibadan, Nigeria. F.B.A. Giwa, B.Sc., Ph.D., D.I.C. Professor in the Department of Physics and Director of the Computing Centre University of Ibadan Ibadan, Nigeria. UNIVERSITY OF IBADAN LIBRARY -7 - CONTENTS Page Title ... 1 Abstract 3 Acknowledgement 5 Certification 6 List of plates 10 List of figures 11 Notations and symbols 12 CHAPTER 1 I N T ROD U C T ION 14 ... CHAPTER 2 A REVIEW OF LINE SQUALL OBSERVATIONS 2-1 West African synoptic weather pattern 20 2-2 Observational investigations of tropical line squall 25 2-3 Satellite investigations of line squalls 41 CHAPTER 3 A REVIEW OF LINE SQUALL MODELS 3-1 Mathematical models 50 3-1-1 Gravity wave models 50 3-1-2 Non-linear steady state model 51 3-2 Numerical models 52 3-2-1 Two-dimensional models 53 3-2-2 Three-dimensional models 58 UNIVERSITY OF IBADAN LIBRARY - 8 - Page 3-3 Other theoretical models 60 3-3-1 Plume model 60 3-3-2 Bubble-theory model 61 3-4 Syntheses of line squalls 61 CHAPTER 4 EVOLUTION OF LINE SQUALL 4-1 Necessary atmospheric features in the development of line squall 69 4-2 Initial stages and trigger mechanism of line squall 73 4-3 Fully-developed line squall 74 CHAPTER 5 FUNDAMENTAL EQUATIONS 5-1 Basic and perturbation equations 84 5-2 Frequency equation 87 5-3 Effect of atmospheric stability on frequency equation 90 CHAPTER 6 SOLUTION TO THE FREQUENCY EQUATION 6-1 Boundary conditions 94 6-2 Mathematical derivations 98 6-3 Special case 102 UNIVERSITY OF IBADAN LIBRARY - 9 - CHAPTER 7 Page SUMMARY AND CONCLUSION 104 LIST OF REFERENCES 115 ._. _ ._.-. UNIVERSITY OF IBADAN LIBRARY - 10 - LIST OF PLATES PLATES TITLE PAGE l(a) Line squall at ~innaon 12 September 1977. View to the E at 1555 GMT showing advancing line of cumulonimbus 29 l(b) Line squall at Minna on 12 September 1977 View to the NE at 1611 GMT showing roll cloud ••• 30 l(c) Line squall at }unna on 12 September 1977. View to SSE at 1613 GMT showing front of cumulonimbus belt and edge of rain area 31 1(d) Line squall at Hinna on 12 September 1977. View to W at 1648 GMT showing individual cumulonimbus 32 2(a) The 'Fortune' squall of 5 September 1974. Visible SMS - I image of West Africa at 0300 43 2(b) The 'Fortune' squall of 5 September 1974. Visible SMS - I image of West Africa at 0930 44 2(c) The 'Fortune' squall of 5 September 1974. Visible SMS - I image of West Africa at 1100 45 2(d) The 'Fortune' squall ot 5 September 1974. Visible SMS - I image of West Africa at 1230 46 2(e) The 'Fortune' squall of 5 September 1974. Visible SMS - I image of West Africa at 1400 47 UNIVERSITY OF IBADAN LIBRARY - 11 - LIST OF FIGURES FIG. TITLE PAGE 1. Three disturbance lines in West Africa 16 2. North-South section of the prevat l.mg winds over ~-1estAfrica 21 3. Monthly frequency and diurnal variation of squalls over 15 mls at various Nigerian stations 24 4. Number of disturbance lines \olhichgenerated and decayed in each 50 square over the contient and ocean during GATE 35 5. Time~height section of wind component over a GATE ship (2 - 17 September 1974) 37 6. Tephigram showing T9 Td for varous post-squall soundings 39 7. Two types of arc lines in satellite imagery 48 8. Idealized history of a line squall 63 9. A sketch of vertical profile of wind over West Africa between April and September 70 10. The process of formation of line squalls 75 11. Schematic croaa-eect.Ion through a class of squall system 80 12. A section of the interface between the air masses in a two-layer atmospheric model 95 UNIVERSITY OF IBADAN LIBRARY - 12 - NOTATIONS AND SYMBOLS u x - component of velocity v y - component of velocity w vertical velocity p pressure .!2E.. w Dt p density of air g acceleration due to gravity v coefficient of viscosity f Corioli's parameter s twice the angular velocity of the earth L latent heat of vaporisation H scale height r condensed water vapour per unit mass of air parcel q specific humidity C specific heat at constant pressurep Cv specific heat at constant volume e potential temperature UNIVERSITY OF IBADAN LIBRARY - 13 - c - C P v K Cp T temperature z vertical height a frequency of waves wave number along x-direction B wave number along y-direction R gas constant UNIVERSITY OF IBADAN LIBRARY - 14 - CHAPTER 1 INTRODUCTION The weather pattern in West Africa could be classified into two main seasons: the wet and the dry. One of the distinguishing factors between the seasons is the amount of rainfall. Almost all the rain- fall within the sub-region are recorded during the wet season which is approximately from April to September. About half of the total annual rainfall recorded in any particular location within West Africa is due to the precipitation which usually accompanies isolated thermal convection in the atmosphere (Obasi, 1976). Such pockets of convective processes are referred to as LOCAL CONVECTIVE STORMS. As expected, the frequency and intensity of local convective storms vary from place to place; hence, total annual rainfall is not uniform in the sub-region. Apart from local convective storms, there exists a special kind of storm that travels and thus covers a large area of land whenever it occurs. These 'travelling' storms are also accompanied by heavy rain- fall. More than half of the total annual rainfall in West Africa is attributable to these storms which are kno~m as LINE SQUALLS. The extent of 'travel' of line squalls (about two to three nights of 'travelling' at an average rate of 15 m/s (Fortune, 1977» puts them in the class of mesoscale systems. UNIVERSITY OF IBADAN LIBRARY - 15 - r Line squalls can be defined in various ways depending on the mode of identification and place of occurrence. However, the definition given by Zipser (1977) seems all-encompassing and therefore adequate. According to him, line squalls are: cumulonimbus clouds, organised in linear fashion, associated with a pseudo-cold front (squall front) at the surface. propagating with considerable speed with respect to the ambient low-level air~ in the general direction of the squall wind in the cold air behind the squall front. Line squalls have been referred to by various names at different times and places. Hamilton and Archbold (1945) called them DISTURBANCE LINES (D.L. 's) because the squalls occur along a fictitious bow-shaped line that has N-W to S-E orientation (Fig. 1). Other names are: chubascos (in Central America), haboobs (in Sudan) and sumatras (in Malaysia). These various names indicate the geographical spread of the occurrence of line squalls. Thus, they are not entirely a local West African (or tropical) sub-synoptic phenomenon. Marriot (1892), Prohaska (1907) and Browning and Ludlam (1962) reported the existence of this kind of storms in the mid-latitudes. Although mid-latitude and tropical line squalls propagate at speeds close to the speed of the mid-tropospheric winds in their respective regions (ioe about 15 m/s in West Africa), they differ in structure because the anvil in the latter extends behind rather than in front of the convective elements UNIVERSITY OF IBADAN LIBRARY 16 ~D-Ll --- - - -D-L 2--- -._ .D._L._3 .-. o 480 960kms I I Fig. 7.' Three dis turbance lines (DL1,DL 2, DL3) in West Africa (the fict itious bow-shaped lines propagate westwards across West Africa) (trom Garnier, 1967-) UNIVERSITY OF IBADAN LIBRARY - 17 - (Fortune, 1977; Mansfie1d~ 1977; Zipser, 1977). This difference in structure is a direct consequence of the variations in wind profile between the mid-latitudes and the tropics. In the mid-latitudes, wind speeds increase almost steadily with height and as such, the high wind speed aloft makes the updraught air, which originates ahead of the storms, flow out at upper levels whereas in the tropics, there is a low-level jet rather than strong upper level winds (Bolton, 1981). The importance of line squalls lies principally in their contri- bution to the economic growth of West Africa (vis-a-vis agriculture). This economic importance is most appreciated if it is realised that the total annual rainfall in the Sahel region of West Africa is due, mainly, to these storms. This implies that years of little squall activities are usually dry in the Sahel. The consequent drought leads to loss of life (both human and animal) and poor agricultural output. Apart from their ,economic and social implications, hazards in the field of aviation is another reason why studies on the initiation, maturity and dissipation of West African line squalls have engaged the attention of researchers for almost four decades now (Hamilton and Archbold, 1945; Dhonneur, 1970; LeRoux, 1976; Fortune, 1977; Okulaja, 1978; Bolton, 1981). On the global scale, the contribution of the phenomenon to the overall transfer of mass, energy and momentum is of interest to meteorologists. In chap~ers 2 and 3, reviews of some of UNIVERSITY OF IBADAN LIBRARY - 18 - the methods that have been employed in the study of line squalls are presented. Broadly, these methods could be classified into three: obse~1ational investigations, satellite investigations and modelling. Whichever method is used to study line squalls however, complete theories on them must account for their growth, the fully-developed stage and dissipation. Furthermore, such theories must account for their orientation, direction of movement and the almost-uniform speed of propagation that is common to West African line squalls. It is also desirable to explain the limitation of line squalls to a perticu1ar season of the year and the similarities between a local convective storm and line squall. Almost all features of the fully-developed (or mature) line squall have been documented (Hamilton and Archbold, 1945; Mansfield, 1977; Bolton, 1981). Outside West Africa, but within the tropics, the following works on line squalls are notable: Zipser (1969, 1977), Belts, Grover and Moncrieff (1976) and Miller and Betts (1977) while the description of the Wokingham storm by Browning and Ludlam (1962) gave a good picture of mid-latitude 'travelling' storms. Data from satellite photographs have also added to our knowledge of the matured line squall. In spite of these works, however, there is still some doubt as to the meteorological factors on which the origin of line squalls are UNIVERSITY OF IBADAN LIBRARY - 19 - dependent. More precisely, the trigger mechanism which at verious times have been linked with insolation, orography and synoptic conver- gence, is still not fully understood. This study has, therefore, been directed towards the ~-::'1.:!s.ttaigaels in the development of line squalls; paying special attention to the trigger mechanism. In chapter 4, the mechanism of initiation of line squalls, along with a description of the atmospheric features during the matured stages, are described. The process of initiation of line squalls is mathematically modelled in chapters 5 and 6. A frequency equation is derived in chapter 5 from the basic set of hydrodynamical equations of motion. This frequency equation is solved in chapter 6. Since line squall is not a single mode but rather an amplifying patch of convection that consists of many modes of different growth rates and phase velocities (Bolton, 1981), solutions to the frequency equation are expected to reveal modes with the greatest amplifications. It is proposed that such modes play vital roles in the initiation of line squalls. UNIVERSITY OF IBADAN LIBRARY - 20 - CHAPTER 2 A REVIEW OF LINE SQUALL OBSERVATIONS 2-1 West African synoptic weather pattern As a background to understanding the various stages in the development and study of West African line squalls, a brief review of West African synoptic weather pattern is essential. The air masses and the prevailing winds in West Africa are shown in figure 2. One of these air masses is labelled the MONSOON. According to Hamilton and Archbold (1945), the monsoon winds originate off the coast of South Africa and blow south--easterly for some time before they veer to their direction of south-west in which they arrive at the coast of West Africa. The long track of these winds over the Atlantic ocean accounts for their considerable moisture content by the time they arrive at the southern part of the sub-region. On the other hand, HARMATTAN WINDS blow across the northern part of West Africa. In contrast to the monsoon south-westerlies, this air mass is very dry because of its long track over the Sahara desert. The north-easterlies and the south-westerlies meet at a zone of strong convergence called the Inter-Tropical Convergence Zone (ITCZ) (Fig. 2). Other names for this zone are: Inter-Tropical Front (ITF) and Inter-Tropical Discontinuity (ITD). The ITCZ moves North or South in sympathy (but with a phase lag of about one and a half months (Bolton, 1981» with the motion of the sun between the tropics of Cancer and UNIVERSITY OF IBADAN LIBRARY 21 metres 6000 Westerlies Easterlies 3000 Harmattan (north ecsterly) Monsoon (sout h westerly) Zone A Zone B Zone C 1Z'N-- - --- - -- - -8°N- - - --- -- - 6N (January) (July) 2d'N- - - - - - - - - - - - - - - -6°N Ag·2: North - South secti on of the prevailing over West Africa (from Hamilton and Archbold, 7945) UNIVERSITY OF IBADAN LIBRARY - 22 - Capricorn. This zone is farthest North (about 20oN) around July and at its southermost extent (about 8oN) around January. When the zone is far North, West Africa is covered by the monsoon south-westerlies; this is the wet season. Towards the end of the year when the zone is far South, West Africa experiences the dry season because the sub-region is almost entirely covered by the dry north-easterly aarmattan winds. In discussing the characteristic weather pattern of West Africa, Hamilton and Archbold (1945) divided the sub-region into four weather zones that run East-West (Fig. 2). Balogun (1974) adopted the suggestions of Hamilton and Archbold (1945) and Walker (1959) and added a fifth zone to the four previously mentioned. Zone A of Hamilton and Archbold is the region north of the rTCZ while Zone B extends from the ITCZ to about 150 km southwards. Zone C, with a width of about 500 km lies south of Zone B. Zone D~ which is over the land between July and September, lies south of Zone C. In other months of the year, Zone D is over the Atlantic ocean. Within Zone D, there is an isothermal (or inversion) layer between 800 mb. and 850 mb. pressure levels. As we shall discuss later, the existence of this inversion might have something to do with the observed dissipation of line squalls around the coasts of West Africa. Zone C is distinguished from other zones because it is UNIVERSITY OF IBADAN LIBRARY - 23 - periodically traversed from East to West by D.L.'s (Hamilton and Archbold, 1945). According to Garnier (1967), the movement of the ITez controls both the number and duration of the weather types experienced in different parts of West Africa. Broadly, the resulting pattern is a latitudinal one. For instance, Ibadan (70 26' N 30 54¥ E) experiences weather conditions in the following way: Zone A late December to part of January Zone B February and part, or all of March Zone e April (or part of March) to about mid or late July Zone D late July and part, or all of August Zone e late August to end of October or early November Zone B briefly in November and early December Zone A late December to part of January Since Zone e is, as said earlier, traversed by D.L.'s, a station like Ibadan experiences intense D.L. activities twice in a year as the break-down of weather conditions given above shows. This explains the two peaks observed on the chart of monthly frequency of line squalls at some Nigerian stations (Fig. 3). The peaks coincide with the periods when the stations experience Zone e weather conditions. In section 4-1, it is shown that all the essential atmospheric features that favour the development of persistent storms are present at such stations during UNIVERSITY OF IBADAN LIBRARY 24 >. u sc 0" .•e.... >. x..:- co ~ A 0 0 c o +'" d 20 \... d > d 10 \c... ::J a 0 '=O"'-:-4""'-r-;-o....,.T'<"><""T'" 1-:::,.....-norr'="T:7-r= 300 air ~ ~ 400 1.00 500 600 600 700 700 800 800 Ambient ~d la 900 900 Squall front 100km Fig.11: Schematic cross section through a class of squall system. All flow is relative to the squall line which is moving from right to left. Circled numbers are typical values of Etw In ll( (from Zipser, 1977) UNIVERSITY OF IBADAN LIBRARY - 81 - The extent and persistence of line squall systems suggest the possibility of organised ascents in the upper half of the tropo- sphere. These ascents would be very vital to the sustenance of the storm. But ascents would normally give rise to downdraughts. t-fansfield(1977), Miller and Betts (1977), Zipser (1977) and Bolton (1981) have all observed two scales of down draughts in line squall systems. According to Zipser (1977) one of these downd raught s could be located within the lowest 200-400 m behind the squall front. Its origin may be the ambient cloud-layer air (900-800 mb).lifted into the cumulonimbus and sinking after becoming negatively bouyant or from higher up in the cumulonimbus where downdraughts ncrmally originate from water loading and entrainment of air with low wet- bulb potential temperature, 6w , (Zipser, 1977). The other scale of downdraughts has lower ew . Zipser (1977) postulated that their origin must therefore be from areas of low ew in the ambient winds. This source must be typically above 750 mb. These downdraughts approach the line squall from the front and the rear (Fig. 11). While overtaking from the rear is more easily understood because of desceuding 650 mb. jet, the entry from the front is more difficult to visualise. However, the picture becomes clearer if it is realised that along the squall front, cumulobimbus UNIVERSITY OF IBADAN LIBRARY - 82 - clouds are not continuous in distance ar in time. Fujita (1963) noticed mesoscale high pressure regions immediately following line squall passage as well as the occasional presence of a mesoscale low-pressure area some distance behind line squalls. Although this observations was over the temperate zone, Zipser (1977) has also observed mesohighs and mesolows within the tr6pics. A large depth of cold saturated air combine with large liquid water content (Sanders and Emanuel, 1977) to cause a rise in pressure in regions immediately following line squall passage. Mesolows exist at a distance of about 100 km behind the squall front. Data compiled by Aspliden et al. (1976)(Fig. 4) show that line squalls scarcely propagate beyond the coasts of West Africa. Hence, it would be right to conclude that line squalls decay around the coastal areas. This decay might be as a result of extensive convective activities around the coasts such that squalls propagating into these areas die out. l\nother reason for the decay might be the fact that the 650 mb. jet now impinges on the ocean surface rather than the solid earth on which it could spread out and form squall fronts which force the south-westerlies to rise. The layer of inversion (Fig. 2) observed around the coasts could also contribute to the decay of line squalls by inhibiting convective processes. This layer acts like a 'lid' that prevents upward growth of clouds(Balogun, 1974). UNIVERSITY OF IBADAN LIBRARY - 83 - This 'lid' is never taken off unless the force of convection is enought to destabilise the inversion layer. In the absence of precipitation from convective elements, line squalls would dissipate completely. Attention can now be focused on a gravity wave model for the initial stages in the development of line squalls. UNIVERSITY OF IBADAN LIBRARY - 84 - CHAPTER 5 FUNDAMENTAL EQUATIONS 5-1 Basic and perturbation equations The basic equations in the theory of atmospheric circulations are given below in notational forms. All symbols have their usual meaning. However, for the sake of clarity, the meanings of these notations and symbols have been given. au + 2 at \7(~)2 + fAu + \p7P+ g = v\72u ·.... 5.1 ...Q.. (Lr + Dt C e) = 0 ·.... 5.2 P Dp + Dt p\7u = 0 ·.... 5.3 DDt (r+q) = 0 5.4 Te = (p rPo ·.... 5.5 5.1 is the momentum equation with viscosity term,. 5.2 is the energy equation and 5.3 is the equation for continuity of mass. 5.4 represents the conservation of water substance and 5.5 is the definition of potential temperature. The Corio lis parameter, f, in the momentum equation could be neglected because the square of the time scale of line squalls is much less than -2f • All molecular effects (e.g. viscosity) are also neglected because the scale of motion is much more than molecular dimenSions. UNIVERSITY OF IBADAN LIBRARY - 85 - If DDtr = Q, 5.1 - 5.5, after sinp1ification, become: · . • •. 5.6 • • • •• 5. 7 aw + aw + dW + waa\z01+ 1. ap + = 0at uax Vay P dZ g • ••.•. 5.8 • . • •• 5.9 · . . .. 5.10 19..+Js+~+~+Q= 0 at dX dy dZ · . . .. 5.11 ••••• 5.12 The set of hydrodynamica1 equations 5.6 - 5.12 consists of non- linear equations and, hence, they are difficult to solve. The perturbation method reduces them to a linear form. The basic assumption of the perturbation method is that atmospheric motions are made up of small perturbations superimposed on a basic atmospheric state. Fundamentally, the basic motion as well as the perturbations must separately and jointly satisfy the hydrodynamical equations. Thus, the atmospheric variables in 5.6 - S.12'could be expressed as linear combinations of their basic UNIVERSITY OF IBADAN LIBRARY - 86 - ;>- and perturbation values as: u = u+s u' v = v+s v' w = w e = es+ e' • . . •. 5.13 q = qs+ q' p = ps+ p' where subscript's' denotes the basic values and the primed variables are the perturbation quantities. The set of equations 5.6 to 5.12 could be converted into their equivalents in the pressure co-ordinates. As stated earlier, an advantage of the pressure co-ordinates is the availability of more data (on atmospheric parameters) along isobaric surfaces. Further, the frequency equation being sought appears neater in the pressure co- ordinates. (The equivalence of the frequency equation in z-coordinates have been documented by many authors e.g. Giwa (1965)). The set of values in 5.13 could be incorporated, along with the co-ordinate conversion, into the hydrodynamical equations to give, after lineari- sation: UNIVERSITY OF IBADAN LIBRARY - 87 - ;:.- au' au' dUS 1 p' at + us -a+x wap+- P = 0s ax av' + u -d-V+' 1 ap'at s dX psay = 0 __ l_(dw_+ Us aw) +1:. ap' + R Ps3 a ax Ps ap p T' = 0 s . • .. 5.14 au' + av'+ 3y aw = 0dX dP -a+e' u -ae' ae -A s ax's a+x OJdP - +LQ c = 0p T' = Pe,(~)K Po 5.14 is a set of perturbation equations which could be used to derive a frequency equations. 5-2 Frequency equation The perturbations superimposed on the hydrodynamical equations constitute the disturbances in the atmosphere. Since atmospheric disturbances are wave-like in nature, these perturbations could be expressed in the form: = i(o t + (Xx + By)u' u(p) e 0 = i(o t + (Xx + By)v' v(p)e 0 etc. UNIVERSITY OF IBADAN LIBRARY - 88 - ,,t- ::.~ S'Uch that -a-u='at ia u'0 au' ax = iau' and au'ay = iBu' Le. a at = ia0 a = ia •.•.•• a· 5.15ax a ay = i8 If it is assumed that density does not vary along the horizontal (i.e. x direction) in the mean atmosphere, and so, if G' = =-aG'ax and 1 ap' = aaGy'p ay s SiRce, as mentioned in chapter 4, the horizontal extent of the line squall phenomenon is much more than its vertical extent, the hydrostatic relation could be used in this mathematical formulation. UNIVERSITY OF IBADAN LIBRARY - 89 - In effect, the vertical acceleration terms will be neglected in the set of equations in 5.14. If 5.15 is incorporated into 5.14 and tracer element, ~, placed on the latent heat term, the following set of equations is obtained: ia* au u' + wap+s iaG' = 0 ·... 5.16 iav' + i8G' ::: 0 ·... 5.17 aG' +-RT'ap Ps = 0 •••••5.18 iau' + i8v' + aawp = 0 ·... 5.19 iaS' + ass LQwap+cll ::: 0 ·... 5.20 p P T' = e' (-p!!·l ..... 5.21 0 Eliminating u', Vi, 6', T' and G' from the set of equations 5.16 - 5.21, the result is: 2'- 2 ~l+ _R_Ts[.~LgP s ~a_L_~]\1w ar 2a a 2 p2 = 0 •••• 5.22p g R6sCs p ap Cp dZ R2T ~LPs ar) aTIf r = -~2L(lC --R6- ap ~, ••.•. 5.23g p s 5.22 becomes 2 at (a2~+ + 82 ap2 p2 )w =2 0 .... 5.24s a * a = a + u a0 s UNIVERSITY OF IBADAN LIBRARY - 90 - 5.24 is the frequency equation being sought and its solution will be discussed in the next chapter. 5-3 Effects of atmospheric stability on the frequency equation The measure of atmospheric stability in 5.24 is r and as defined in 5.23 2 r R T J..lLP dr dT= __ s [-.-8.(1 s 2 C )- ~ dP +~P s dZg If Y represents the environmental lapse rate and rd denotes the dry adiabatic lapse rate t J..lLP '\ r = ___ s ~) _ yJ 5.25R9s dP In a dry atmosphere, the latent heat term in 5.25 will be non- existent because there would not be any condensation i.e. u = 0 Thus, 2 R rd = --TT f rd-y] g where the subscript '2' is used to distinguish the atmospheric stability within such a dry environment. In the context of West Africa, 5.26 would represent the stability condition within the north-easterlies. In a moist environment, J..l= 1 because of the possibility of condensation during a convective process. This possibility would then justify the retention of the latent heat term. The term UNIVERSITY OF IBADAN LIBRARY -- - 91 -LP [ 1 __ 8 ~fd Re op s in 5.25 would then represent a modification of the dry adiabatic lapse rate due to precipitation. If a layer of saturated air is assumed, this term will be the wet adiabatic lapse rate, fw • LP i.e. rw = fd[ 1 - R6s aor?l s Since the south-westerlies are moist, the stability of this air mass could be written as: • . • . • • .. 5.27 The south-westerlies have been shown in section 4-1 to be convectively unstable and as such, the environmental lapse rate, y, within this air mass lies between r w and fd• This inequality implies that fl (5.27) will be negative. If such a negative value of r is used in the solution of 5.24, the condi- tion for obtaining a sinusoidal function of time for any perturba- tion will be o2 < O. UNIVERSITY OF IBADAN LIBRARY - 92 - This implies that the frequency of the waves generated is purely imaginary and such waves are known to amplify, without any propaga- tion, with time. In practice, an environmental lapse rate that is greater than the dry adiabatic lapse rate is a sign of instability. Such an environment cannot be maintained for along period because mechanical mixing would occur and thereby stabilise the atmosphere. This means that the environmental lapse rate, even in dry air, is scarcely larger than rd and if it does, it is for a short period. Therefore, r2 (5.26), which is a measure of the stability of the north-easterlies, will be positive. If such a positive value of r is used in the solution of 5.24, the condition for obtaining a sinusoidal function of time would be The implication of this is that the frequency ~f the waves generated will be real. Such waves propagate~ without any amplification, with time. Given a situation, therefore, where there are two air masses, whose separate stability conditions give complex and real values for wave frequencies respectively, one on top of the other (e.g. the north-easterlies on the south-westerlies) the waves generated UNIVERSITY OF IBADAN LIBRARY - 93 - along the surface of discontinuity would both propagate as well as a amplify with time. The velocities of such waves would be complex; the real part represents the phase velocity while the imaginary part is a measure of the amplification of the wave. UNIVERSITY OF IBADAN LIBRARY - 94 - CHAPTER 6 SOLUTION TO THE FREQL~NCY EQUATION Equation 5.24 could be solved by assuming a single-layer model for the atmosphere. This kind of solution had been discussed by Ogura and Charney (1962), Hiller (1974) and Bolton (1981). However, since line squalls require two layers of air with different static stabilities for their initiation ann sustenance, it would be more realistic to solve equation 5.24 with the aid of a two-layer model of the -atmosphere. 6-1 Boundary conditions Figure 12 shows a section of the interface between the air masses in a two-layer atmospheric model. If wand u are the vertical and horizontal components of velocity respectively, the velocity component perpendicular to the surface is wcose - usine • This value must be.continuous if no air parcels accumulate at the interface. For small € (fig. 12), aax£ = tane ~ eine. Also, coae ~ I Therefore, wcose - usine a£~Tjl-U- ax UNIVERSITY OF IBADAN LIBRARY 95 w "k1----~u I I I I I I I I I I € I Fig. 12 : .A section of the interface between the arr masses In a two-layer atmospheric model UNIVERSITY OF IBADAN LIBRARY - 96 - w = DDet: (fig. 12) i.e. w'" .... 6.1 = iooe: + a.use: (using 5.15) i.e. w '"ioe: • • • • • • • • • •• 6.2 From 6.1, 3e: ae:w - uax '"at = ioo£ but w - ae:uax is continuous at the interface, hence e: is continuous. If e: is continuous. then (from 6.2) iW(1 is continuous • w = H W (discussions Hitlr:Gi~~ra1.,981) P where H is scale height. Since w o is continuous, then aapW is continuous. Infinite pressure gradient are never observed at the interface between air masses. If ever they occur, it is only for a very short time because air would quickly rush from a region of high pressure to that of low pressure. So, pressure is necessarily continuous at the interface. UNIVERSITY OF IBADAN LIBRARY - 97 - i.e. p(z+£) is continuous. p(z+£) = p(z) + £~Clz ClsP(z)= ps(z) + p'(z) + £ Clz Ps (z) is a constant and the underlined term is a product of small quantities which makes it negligible; therefore, ClP(z ) p'(z) + £ ~z is continoous • By definition, w = .QE.Dt Using 5.15 and the hydrostatic relation, w = Lop ' - Pgw i.e. wia = p ' -~ ia . • . • .. 6.3 It has been proved that ClPs(z) p'(z) + c dZ is continuous. So. using 6.2, dP (z ) p' (z) + 1~a- --Cl~:--=- p'(z) - i~ is continuous. From 6.3, this means w a is continuous. For'the scale of motion consLdered in this study, the surface of the earth is usually assumed to be flat. Because of the flatness, UNIVERSITY OF IBADAN LIBRARY - 98 - the vertical component of velocity is taken to be zero at the surface i.e. w = o. To summarise, the following terms are continuous at the lllter- face between the air masses: W a and R W a dWa p 2 2 ap-g(a. +S ) At the surface, w = 0 i.e. R W aa p 2 2 -da=Wp 0 and at theg(a. +S ) top, W = 0 (this is a free surface condition). 6-2 Mathematical derivations To simplify mathematical derivations in our two-layer atmospheric model, parameters like density, temperature and horizontal velocity that should normally vary with height are taken to be constant withic each layer. Subscripts '0', 'B' and 'T' indicate the values of these parameters at the surface, boundary between the air masses and the top of the atmosphere respectively while subscripts 'I' and '2' are used to identify atmospheric parameters within the lower and top layers respectively. We recall equation 5.24 i.e. 2 2 2 ~+ =gf- (a +S ) ap2 p2 " W = 0 . s a "- UNIVERSITY OF IBADAN LIBRARY - 99 - The solution to this differential equation is of the form m m w AP 1 + BP 2 . . .. 6.4 where A and B are constants and are the roots of o .•6.5 From 6.5t let J !- gr(a 2+82)ml = ~ + 2 (J and If A = w = APl/2+A t Bpl/2-A (from 6.4). At the boundary between the air massest it was shown that H w (J P and w (J are continuous. At the surface, w = 0 and at the top w = O. Applying these conditions to 5.24 with the following substitutions: UNIVERSITY OF IBADAN LIBRARY - 100 - Ql'= 2 Q2) g(a + ~ and 2 2 two equations are obtained: g(a +8 ) 1.2 I 1 P 21.2 P -{[ HB-Q2 (A 2+ -2)] -l HB-Q2(- - A )] (-) T }- B2 2 P 02 ·Az = 0 ••. 6.6B . •. 6.7 6.6 and 6.7 constitute a set of simultaneous equations in Al and A2• For non-trivial solutions, the determinant of the coefficients must vanish. If PT 21.2and R2 = (p-) ,the value of this determinant is B 2 °2 - -2- = 0 . • •• 6.8 UNIVERSITY OF IBADAN LIBRARY - 101 - Recalling from chapter 5, al = a + aU0 l °2 = a +0 au2• Substituting for al and a2 in 6.8, "0 1 + Rl 1 + R2{2aul[2:>"(11 _ R ) + 1]- 2au2[ 2:>"(21 ~.R ) + 1) }+l 2 2 2 1+R 2 2 1+R {a "i [2>'1(1-R 11) + 1] - a "z 2 [2k2 (l-R} + 1]} = o • • •• 6.9 At a glance, 6.9 looks like an ordinary quadratic equation in ao but a closer observation reveals that the coefficients are not constants. This means that no analytic solution exists for 6.9. Thus, the equation reduces to an eigenvalue problem to determine ao such that there are non-zero solutions. Some of the eigenvalues obtained in the solution of 6.9 will be complex and as said earlier, the imaginary part of such complex values would represent the amplification of any wave-like pertur- bation. The modes with the largest amplifications (i.e. the largest imaginary parts) are those that are proposed as being responsible for the onset of line squalls. These modes are those UNIVERSITY OF IBADAN LIBRARY - 102 - that will likely awplify up to the extent of blocking the 650 mh. mid-tropospheric jet as explained in the proposal for the initiation of line squalls (section 4-2). A single-layer model could be considered as a special case of equation 6.9. 6-3 Special .case In a one-layer model, ul = u2 = u (say) and ~l:::":2 = " (say). If these are incorporated into 6.9, •.. 6.10 The conditions for 6.10 to be true are: (i) = o a Le. oa- = +u • This represents the velocity of the ambient wind without any perturbation whatsoever. (ii) .\= O. This imlies a2 = 4gra2 For any single layer of air, this condition represents a phase velocity of a :::: a u ± 2/gf UNIVERSITY OF IBADAN LIBRARY - 103 - where 21gf 'represents' the perturbation on the ambient wind u. As considered in section 5-3, if r is positive, the perturbation propagates, without any amplification, through the ambient winds. On the other hand, if r is negative, the perturbation amplifies with time without any propagation. These conditions depict wave-like perturbations in dry and moist air respectively. (Hi) = o. Substituting for R1 and 2HRO2 ' (P.-!')2A 1 + 2A -0::: ...... 6.11 Ps 1 - 2A - 2HO"T Equation 6.11 had been obtained by Mi11ier (1974) while solving a frequency equation for the speed of propagation of external gravity waves. These are 'free-surface' waves which occur in one- layer models through the use of particular boundary conditions. Simplified cases of these waves in hydrostatic, neutral and non- hydrostatic, neutral systems have been discussed by }til1er (1974). UNIVERSITY OF IBADAN LIBRARY - 104 - CHAPTER 7 SUffi