JPRS ID: 8888 USSR REPORT METEORLOGY AND HYDROLOGY NO.11, OCTOBER 1979

Document Type: 
Collection: 
Document Number (FOIA) /ESDN (CREST): 
CIA-RDP82-00850R000200040057-3
Release Decision: 
RIF
Original Classification: 
U
Document Page Count: 
186
Document Creation Date: 
November 1, 2016
Sequence Number: 
57
Case Number: 
Content Type: 
REPORTS
File: 
AttachmentSize
PDF icon CIA-RDP82-00850R000200040057-3.pdf10.08 MB
Body: 
APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 ~ ~ ~ T ~aG~r N0. 11, N.0"~EM~ER 3979 2$ JRNURI7Y ~980 t FQUO ) 1 aF ~ APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 I~OR OFI~IC'IAI. US1~: ONI.Y JPRS L/8~888 28 January 1980 ~ IJSSR Re ort p METEOROLOGY AND HYDROLOGY No. 11, November 1978 _ ~g f$ FOREIGN BROADCAST INFORMATIOI~ SERViCE FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 NOTE JPRS ~ublications contaiti information primarily from foreign newspapers, periodicals and books, but also from news agency transmissions and broadcasts. Materials from foreign-lan~uage sources are translated; those from English-language sources are transcribed or reprinted, with the original phrasing and other characteristics retained. - Headlines, editorial reports, and material enclosed in brackets [J are supplied by .TPRS. Processing indicators such as [Text] or [ExcerptJ in tha first line of each item, or following the last line of a brief, indicate how the original information was processed. Where no processing indicator is given, the infor- _ mation was summarized or extracted. ~ Unfamiliar names rendered phonetically or transliterated are enclosed in parentheses. Words or names preceded by a ques- _ tion mark and enclosed in parentheses were not clear in the original but have been supplied as appc�opriate in context. Other unattributed parenthetical notes with in the body of an _ item originate with the source. Times within items are as given by source. - The contents of this publication in no way represent the poli- cies, views or attitudes of the U.S. Government. For further information on report content call (703) 351-2938 (ec~nomic); 3468 (political, sociological, military); 2726 (life sciences); 2725 (physical sciences). COPYRIGHT LAWS AND REGULATIONS G~VERNING OWNERSHIP OF Mr1TERIALS REPRODUCED HEREIN REQUIRE THAT DISSEMINATION OF THIS PUBLICATION BE RESTRICTED FOR OFFICIAL USE ONLY. - APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL US~ ONLY JPRS L/8888 28 January 1980 USSR REPORT ~ METEOROLOGY AND HYDROLOGY No. 11, November 1979 Selected articles from the Russian-language journal METEOROLOGIYA I GIDROLOGIYA, Moscow. CONTENTS PAGE - Moistening of the Continents and Intensity of Summer Monsoonal Circulation (G. P. Kurbatkin, et al.) 1 Construction of a Model of the AtmosphPric Boundary Layer tor the Equatorial Zone (Ye. M. Dobryshman) 10 Small Oscillations of the Polytropic Atmosphere and the Filtering Rale of the Hydrostatic Approximation (V. M. Kadyshnikov) 2!? Correlation Between Minimum Pressure ~nd Maximum Wind Velocity in Trbpical Cyclones (V. M. Radikevich, G. G. Tarakanov) 37 Possible Mechanism of Transfer of Disturbances from the Lower Thermosphere into the Meso-Stratosphere (V. I. Bekaryukov, et al.) 48 - Characteristic Diurnal Varialtions of Winds ~n the Upper Mesopause Region - Over Central Europe and Eastern Siberia (R. Schminder, et al.) 58 Evaluation of Errors in Computing Effective :Zadiation ~ (A. I. Budagovskiy, L. Ya. Dzhogan) G4 ~ -a- CIII -USSR- 335&TFOUO] FOR OFFICIAL U5E ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 I FOR OFFICIAL USE ONLY CONTENTS (Continued) page � Optical Characteristics of the Atmosphere in the Tropical Zone of _ the Atlantic Ocean (V. N. Adnashkin, et al.) 72 Present Status of Research on Sea Surface Temperature (F. S. Terziyev, et al.) 82 Maximum Possible Heights of Wind Waves in the Oceans and Seas (G. V. Matushevskiy) 93 Evaluation of Accuracy in Determining Water Discharge in Streams (G. V. Zheleznyakov, B. B. Danilevich) 99 Status and Prospects for the Development of Agroclimatic Investiga- - tions - (I. G. Gringof, Yu. Z~ Chirkov) 104 Degree of Activity of Winter Wheat During Winter Thaws (I. v. svisyuk) 114 - Scientific Center of Soviet Hydrology (Sixtieth Anniversary of the State Hydrological Institute) (V. I. Korzun) 120 Patterns in the Redistribution of Ice in the Waters of the Foreign Arctic (V. I. Smirnov) 135 The Magnus Effect for a Spherical Particle During Deta~ament from a ~ Solid Surface (N. N. Grishin) 141 natermination of Mean Annual Runoff from Slopes When Taking Antierosion Measures (A. I. Gonchar, I. K. Sribnyy) 145 Ensuring the Uniformity of Measurements in the Spstem Operated by the State Committee on Hydrometeorology and Environmental Monitoring (G. N. Kondrashov, L. V. Selivanov) 149 Review of Monograph by V. R. Alekseyev: NALEDI I NALEDNYYE PROTSE�SY (VOPROSY KLASSIFIKATSII I TERMINOLOGII) (Ice Encrustations and I;.e - Encrustation Processes (Problems in Classification and Terminology)), Novosibirsk, Nauka, 1978, 192 pages (B. M. Krivonosov) 154 b FOR OFFICIAL USE ONLY > APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY . CONTENYS (Continued) Page Sixtieth $irthday of Sergey Konstantinovich Cherkavskiy 158 Seventieth Birthday of Andrey Anisimovich Glomozda .............o.... 161 Sixtieth ~irthday of Konstantin Petrovich Vasil~yev 165 = Sixtieth Birthday of Yuriy Ivanovich Chirkov 168 Conferences, Meetings and Seminars fM. A. Butuzova, et al.) 171 Notes from Abroad (V. I. Silkin) 177 - c FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY . ~ PUBLICATION DATA _ English title . METEOROLOGY AND HYDROLOGY Russian title . MET~EOROLOGIYA I GIDROLOGIYA Author (s) ; Editor (s) . Ye. I. Tolstikov Publishing House ; Gidrometeoizdat Place of Publication , Moscow Date of Publication ; November 1979 Signed to press ' ; 23 Oct 79 Copies ~ 3870 COPYRIGHT : "Meteorologiya i gidrologiya", 1979 ci ~ FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY J ~ UDC 551.(513:571) MOIST'ENING OF TI~ CONTINENTS AND INTENSITY OF SUMMER MONSOONAL CIRCULATION Moscow METEOROLOGIYA I GIDROLOGIYA in Russian No I1, Nov 79 pp 5-11 [Article by Correspending Member USSR Academy of Sciences G. P. Kurbatkin, Professor S. Manabe and Doctor D. G. Hahn, Computation Center Siberian De- partment USSR Academy of Sciences and Geophysical Fluid Dynamics Laboratory, Princeton, New Jersey, suhmi.tted for publication 2~ June 1979] Abstract: Using the ~pecr_ral model of the atm~sphere developed at the Geophysical Fluid Dynamics Laboratory (Princeton, New Jersey), a study was made of the influ- - ence of changes in moistening of the con- tinents on the intensity of swmmer mon- soonal clrculation in the middle latitudes. The mc~del includes the annual cycle of cli- mate, the hydrology of the atmosphere and continents. An analysis of the numerical experiments indicated that the drying out of the continents can lead to a decrease of precipitation not only over the contin- ents, but also over the ocean; drying-out � of the continents simultaneously can inten- sify planetary summer monsoonal circulation in th~ middle latitudes, which can be an ; i~nportant condition in the annual cycle of _ ~ climate for summer radiation heating o~ - the ocean. jText] Without allo~vance for the annual variation of solar radiation it is evidently impossible to detect and understand the rela~ive importance _ of different interacting physical processes determining stable and un- stable "weather systems" and in the long run the forming climates. The annual cycle of climate can be quantitatively explained by solution of the problem of intEraction between the atmosphere, continents and ocean in accordance with the characteristic times of the processes participatin~ in this interaction. But how and why the physical processes forming the annual cycle of climate interact with one another to produce climatic, 1 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR UFFICIAL USE ONLY seasonal and year-to-year fluctuations is almost unluiown. For example, one of the unstudied aspects of this prohlem is the influence of changes in t:lie mo:istening of the continents on the intensity of summer atmoapheric ~�Ir~~iilr~t lun. As is well known, planetary monsoonal circulation is secondary circulacion of the convective type. It is not described by the stream function, but is manifested only in the wind velocity potential field a small-scale com- ponent of large-scale horizontal motion of the atmosphere, not subject to direct instrumental measurements. At the present time it is unknown to what extent there is a change in summer planetary monsoonal circ.ulation in the middle latitudes from y~ar to year and what the reasons for these changes are. It is postulated that the reason for its changes may be a change in the moistening of the continents. However, at the present time not even mean seasonal maps of moistening of the continents for each year are being compiled. We have only mean long-term (climatic) seasonal maps of moistening of the continents [1]. Moreover, at the present time mean seasonal maps of the wind velocity potential field are also not being com- piled. This does not make it possible to judge the nature of the changes in summer monsoonal circulation of each year. - Thus, at the present time it is virtually impossible to investigate this ~~rol~Lem by means of a diagnostic analysis of observational data. For studying summer planetary monsoonal circulation in the middle lati- tudes we used a complex spectral model of general circulation of the at- mosphere developed at the Fluid Dynamics Laboratory located at Prince- - ton University in the United States [2]. , Spectral models differ from grid (finite-difference) models in that the dynamic variables in them are represented by a synthesis of a finite sum - of spherical harmonics, and not by the values in a grid of discrete points. The equations of the model predict the spectral components, and not the variables at the grid points. The predicte~ variables in this model are the following: d~ of the stream function, ~ of wind velocity potential (~7 2 is the horizontal Laplace operator), temperature, mixture ratio and logarithm of pressure at the level of the earth's surface. These variables are scalars (the spectral representations of vector values introduce singularities at the poles). In tlie model use was made of a hydrostatic approximation and in the vertical direction use is made of the d coordinate, equal to the ratio of pressure to pressure at the earth's surface in order to introduce topographic ef- fects. SimPle liorizontal viscosity was introduced by attenuation of the model variables by a constant multiplied by 'd 4 of the hydrothermodynamic elements. A simple diffusion scheme wi~ch the mixing length is used vertic- ally. 2 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY The model has nine vertical levels (d = 0.025, 0.095, 0.205, 0.350, 0.515, 0.680, 0.830, 0.940, 0.99~). These were selected in such a way as to de- scribe the lower stratosphere and the Ekman boundary layer. The liorizontal rliomboidal resolution of the model is M= 21. Tiie prognostic equations Qf the model are integrated in time using a quasi-implicit scheme: the linear and nonlinear components of the trends are split and integrated in time implicitly and explicitly respectively. Use is made of temporal smoothing � with oC = 0.01 in each time interval. Ttie model includes the annual cycle, the hydrology of the atnosphere and continents. In order to compute tlie flux of solar radiation tliere is stip- ulation of the seasonal ctiange in insolation at the upper boundary of the model atmosphere. The attenuation of sular radiation and the transfer of long-wave radiation emitted by the earth and atmosphere are computed tak- ing ir~to account the effects of clouds, water vapor, carbon dioxide and - ozone. The carbon dioxide mixing ratio is everywhere assumed to be con- _ stant. The zonally homogeneous distribution of ozone is stipulated as a function of latitud e, altitude and season. The tirie-dependent spatial dis- tribution of water vapor is fosnd as a result of integration in.time for the prognostic equa tion for water vapor, includ ing: three-dimensional advec- tion of water vapor, vertical mixing of water vapor in tlie planetary boun- dary layer, evaporation, nonc~nvective condensation and moist convective ada~tation. In computing radiation fluxes an allowance is made for the time-variable distribution of cloud cover at three levels in dependence on the change in water vapor and air temperature. Ttie temperature of the earth's surface over the continents is determined by tlie boundary condition expressing the accumulation of heat in the soil (tliat is, the fluxes of solar and long-wave radiation and the turbulent fluxes of apparent and latent heat locally together are equal to zero). Over the oceanic part the seasonal change in temperature of the ocean sur- face is stipulated. It is determined by interpolation in time between the four observed distr ibutions of the mean monthly temperature fields of the - ocean surface. In o rder to compute the descendin g flux of sol.ar radiation the albedo of the earth's surface i.s stipulated as a function of latitude over the ocean and as a function of latitude and longitude over the con- tinents; in places where as a result of the computations there is a repro- duction of snow cover or sea ice, the albedo is replaced by higher values. The rates of change in moistening of the continents and thickness of the snow cover are determined by the budget of water, snow and heat at the land surface. Numerical integration in time was carried out for two years and eight model montlis. In this control variant the distribution of moistening of the con- tinents was determined at each moment in time from the budget of evapora- tion, precipitation, melting of snow and river runoff. The model quite precisely reproduces climate and its seasonal changes. Fig- ure 1 shocas the field of atmospheric pressure at sea level (a) and the field of moistening of the continents (b) reproduce3 by this model in the 3 - FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 r�uK ur�r'iC1AL USE UNLY control variant, averaged for July and August (here the moistening values beyond the continents must not be taken into account). Ttiey can be compar- ed with the corresponding fields constructed on the basis of long-term ob- servations [1, 3], which, unfortunately, are not cited here in order to save space. Q) :1: ::i,~:..;.': ' ~ - ~ I~. I , I ~ I I ~ ~ i ~ \ . I i: I . - . \ r: \ \ :i: i : ` \ \ ~ ~ 1 . ':I~�'.~~ I� ~ :i . �`I'. � ~.i r'� ' . ' � ~~I : - ; : \ : , ; ~ . . + \ \ ~T; :rl;:::.~.:::::,::~::�`:r:;�>::. ~ Z, S ~3?6 - Fig. 1. Atmospheric nressure at sea level (a) and moistening of the contin- ents (b), reproduced by Geophysical Hy3rodynamics Laborator}r model (Prince- ton, United States). Isolines are drawn each 10 mb (a) and each 4 cm (b). 1) less than 1000 mb, 2) 1000-1020 mb, 3) more than 1020 mb, 4) less than - 2 cm, 5) 2-10 cm, 6) more than ld cm From an analysis of these maps it is possible to n.ote a coincidence of the principal regions of high and low pressure. However, they are more strong- ly expressed in the model climate than in the real climate. In an analysis of the maps for moistening of the cpntinents it is possible to note a 4 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY coincidence in the control variant with the real climate of extensive re- gions of quite high moistening ~more than 10 cm, shaded regions) in the tropical zone of Africa, in India, on the islands of Indonesia and in the middl.e latituaes of the F.urasian continent. Q~ f / i ~ i /`'i ~ . / , ; ~ / / ~tY . ~ ~ '~j i~ ' ~ ~ j' i~` ' / / / ~ `~i' / I ' ~ y/ ~ ' i/ 0 / ~ _ i ~ / ' / / ' ! , , , . . . ; i , _ , ~ ~ 1' ; ~ i. ~ / ' i / ' ~ ~ : ~i i i 'i!.',~i/, /Y ~ _ . ^ - ~ / ~1, '/'j/,~y', i 4~,~~'~ ~ ; l.~i/~' f/ 'il I i `.~~r y ~ ~ ~ ~ % . , t ~ ~ , Y ~ ~ y, ; ~ 'i / . / 1 ~'>/~y ~ ; / / /j: �"V~ .i . ,'?j' % ~ / ~ ~ ~S, _ _ ~ I ' ~ i . ~A77`-~-:�' : ~ j ~ ' ~ ~ i l ~ / i , / 'Y ~ ~ i. - " . , , ~ ' ~ - ~ / : ~ ~ / ~ . ~ j ~ / . i _ ~ ; Fig. 2. Difference in rate of evaporation (a) and in rate of precipitation (b) between "Arid" and control variants. Of the two simglest nwnerical experiments with a change in moistening of the continents total moistening of the continents or their total des- sication it was the second which could be favorable for a quantitative , study of the lLydro*_hermodynamic interactions between the atmosphere and the continents in the middle latitudes. Therefore, after carrying out and anal- ~zing the control experiment the last three model months (June, July and August) were recalcuiated witfi the regtriction that the continents remain- ed completely waterless during the entire three-month summer period ("Arid" vartant). - 5 _ FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-00850R040240040057-3 FOR OFFICIAL USE ONLY ' tri/tym cm/day , q~0 6) Q) - a4 ' 1~ ~ '/n ~ ~ ~ ~ i r ~ ~ 4r~, i ; , ~ ~ ~`'ii~.i v ~ ~ ' 4S'GUI. ~ 0 45 Kt m. S' ~ 4S G u!. N D 45~0. ur $ _ Fig. 3. Total quantity of precipitation per day, averaged for July-August over the continents (a) and over the ocean (b). ~ Q) i ~ � ; ~ i, ; . - ' / o $ ` _ ~ ~ ~ 'r' ~ ! ` ~ ~ I ,S~ f ; ~ ~ , ~ ~ / , , i - ' ~ ~ ~ ' . ' . O ~ ~ Q-~ : ~7 ~ I ~ ii, - � O I I ~y ' ~ ~ / ~ % Q /'i ~ ~ o f ~1" t!/ Q , /r ~ ~i / / (wr_ ~ : i / ~ _ ~ ~ J~ ~ i ~ / ~ ~ 11~ / - ~ 'i. ~ r - i , ~ ~ ~ - , ' - --r-- 0 . ~ ~ , j" ~i. ~ t~ i � ' /',I / / . ~ ~ I,' ~ / / ~ � ~ 1/ . ' ~ ~ . I. . i / / / i,1; ,%i / � - / / / ~ ' i . ,./i I ~ / � ~ ' i i ~ Fig. 4. Difference between "Arid" and control variants of horizontally smoothed vertical velocities in isobaric coordinate system at S00-mb level (a) and atmospheric pressure at sea level (b), aver~~ed t'or July-August. � 6 ; FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000240040057-3 - FOR OFFICIAL USL ONLY '['fie dessication of the continents r,aturally led to a decrease in evapora- ` r~Ion Erom the surface of the cont~nents and 3ccordingly to a decrease in ~~rrc~f~~itriC[~~n, rigure 2.i show~ the differences in the rnte of evap~ration b~t;wc~~n ~I~~~ "Ar(c1" .u~d c~unlrc~l vurlr~nLH, wi~~rrut~ rig. 2b ~~I~uw4 HimLl~tr cllf- ferences in the rate of precipitatioti. (In Fig. 2 shaded regions corres- ' pond to ne~ative differences, whereas regions without shading correspond _ to positive differences; the isulines are drawn at a logarithmic sca2e: t0.2 _ and fl cm/day) . _ However, precipitation decreased not only over the continents, but also ~ver the oceans. (Fi.gure 3 shows the total quantity of precipitation per day, averaged during July-August over the continents (a) and over the ocean (b); _ 1) control variant, 2) "Arid"). This is associated with an intensification of descending vextical movements over the ocean in the "Arid" variant. Fig- ure 4a shows the differences in vertical velocities ir. an isobaric coordin- ate system at the 500-mb level between the "Arid" and control variants; the unshaded areas show regions of inter.sification of descending (or weakening of ascending) vertical movements. Figure 4a shows an intensification of as- cending (or weakening of descending) movements over the continents in the - "Arid" variant (shading; the isolines are drawn each 20 mb/day). _ Tl~e intensification of ascending vertical currents over the continents in ' the "Arid" variant caused a pressure decrease over the continents as a re- sult of an increase in the temperature of the continental surfaces (and the lower troposphere over them) due to a decrease in the heat loss on evapora- ti.oti and a simultaneous decrease in the cloud cover, increasing the radia- tion heat influx. Figure 4b shows the difference in atmospheric pressure at _ sea level between the "Arid" and control va-riants (shading region~ of pressure decrease; isolines drawn each 4 mb). The intensification of ascending vertical currents over the continents in t11is "Arid" variant must not only intensify the descending vertical cur- rents over the ocean and the inflow of moister air onro the continents from the ocean (intensification of planetary monsoonal circulation in the middle latitudes), but also somewhat increase precipitation over the con- ' tinents. However, a decrease in evaporation ovzr the continents led to both an increase in air temperature and to a decrease in specific humidity and both factors decreased relative humidity to such an extent that this ef�- fect exceeded the preceding effect of a possib"le increase in precipitation over the continents as a result of an intensirication of ascending ver- tical currents over the continents and led to a general decrease in pre- _ cipitation. Figure 5 shows a map of the influence of drying-out of the continents on the intensification and weakening of different interacting processes. This numerical experiment demonstrated that the dry~ing-out of the continents - can lead to a decrease in precipitation not only over the continents, but ~ 11so over the ocean; the drying-out of the continents can simultaneou~ly inr.ensify the planetary summer monsoonal circulation in the middle latitudes, 7 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-00850R040240040057-3 FOR OFFICIAL USE ONLY which can be an important condition in the annual cycle of climate for sum- _ mer radiation heating of the ocean. � ~ . : - / ~ i-T'/~ ~ - CXOpOClil b UCAOp~NUA C - /AIIIMOC~EpHOE (JIlO/IfHfle~ TeMnepamypa Bo3ayxa nobepxHOCmu KoHmuHeH-. ~ ~ Ha ypOdHB MOJ~R ' Ha ypoBHe Mopa Z %~;moB u oKear+vi - //.iii~..'i'i ii i ri}~ /iiii;~~"~ ,i . Bme?raHUe Bo~dyxa Ho noy- A6conwmHOR OnanrHOCme murrzHmbi OcneBcmGue do3ayxa g MyccoNxou uupxynxuuu ~ ~ BOCXOaAUSUP dBUM(eNlIR ObAQ9Hb1U 0lIiHOCU/OPllbl.'OA B/IQMfNOCT - bo3ayxQ 5 / ~ j o0 : 3a xa j ir CKOpocmn oca8rro0 6 crropocme ocaaKOB - ; ~ ~ r 06waA cKOpncmb oca8xo8 i ii~i Fig. 5. Diagram of influence of drying-out of the continents on intensif- ication (unshaded rectangles) and weakening (shaded rectangles) of differ- ent interacting processes. KEY: 1. Atmospheric pressure at sea level 7_. Air temperature at sea level 3. Rate of evaporation from surface of continents and ocean 4. Inflow of air onto continents as a result of monsoonal circulation 5. Ascending air movements 6. Rate of precipitation 7. Cloud cover 8. Absolute humidity 9. Relative humidity 10. Rate of precipitation ' 11. Total rate of precipitation BIBLIOGRAPHY 1. Budyko, M. I., ATLAS TEPLOVOGO BALANSA (Heat Balance Atlas), Moscow. Gidrometeoizdat, 1963. 2. Manabe, S., Ha~in, D. G., Holloway, S. L., "Climate Simulations With ~GFDL Spectral Models of the Atmosphere: Effect of Spectral Truncation," PROC. OF GARP CONFERENCE ON CLIMATE MODELS, Washington, April 1978. 8 FOR OFFICIAL USE ONLY - APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY - 3. Taljaard, J. J., Loon van H., Crutcher, H. L., Jenne, R. L., "Climat~e - of the Upper Air: Part I. Southern Hemi~phere Temperatures, Dew Point~s and Heights at Selected Pressures," NAVAIR 50-IC-55, Superintendent of _ Documents, Washington, D. C., Vol 1, 1969. 9 . FOR OFFTCIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY ~ UDC 551.510.522(-062.4) CONSTRUCTION OF A MODEL OF THE ATMOSPHERIC BOUNDARY LAYER FOR THE EQUATORIAL ZONE Moscow METEOROLOGIYA I GIDROLOGIYA in Russian No 11, Nov 79 pp 12-22 [Article by Professor Ye. M. Dobryshman, Institute of Atmospheric Physics, submitted for publication 22 May 1979] Abstract: The author analyzes the difficulties involved in constructing a model of the plan- etary boundary l~yer for a narrow equatorial zone. The well-known Kuo model of a plane boundary layer is refined and a model of the three-dimensional boundary layer is construct- ed within the framework of the linear theory. In the latter case the vertical velocity max- imwn is in.the layer 1.5-3.0 km, which agrees ~~ith observational data: the angle between the outer zonal flow and the wind velocity vector near the surface is less (31�) than accord- ing to the classical Ekman formula (45�). ~ [Text] The classical formula for the thickness of the planetary boundary layer (PBL), determined from the ratio of the viscosity coefficient v to the main Coriolis parameter .~1, loses sense at the equator. This occurs because the values ~,1 = 2 ulsin ~ with 0 become equal to zero, which reflects, in particular, the fact of an impossibility of using the quasi- geostrophic approximation in a narrow equatorial zone with a width of ap- - proximately 500 km on each side of the equator [3, 4] (~1 = 7.29�10-5 sec'1 is the angular velocity~of the earth's rotatio~, ~ is geographic latitude). In the equatorial zone the relationships between the wind field and the . pressure gradient are extremely complex, not fitting within the framework - of such a simple linear operator as the operator of geostrophic corres- pondence. The fundamental difficulties in formulating a model of the boun- dary layer (BL) for the equatorial zone are as follows: first, it is neces- sary to be able to "splice" the model with the Ekman model of the PBL; sec- ond, the model must take into account all three wind velocity components in the equatorial zone. The latter circumstance indicates the necessity for drawing upon models of the three-dimensional BL, which, as is well known, 10 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY is more complex, since these have a number of specif.tc peculiarities [11]. In this article we will refine one of the models of a plane BL and will examine a variant of a model of a spatial BL. In attempts to formulate a PBL model for the low latitudes using both anal- - ytical methods [13] and numerical methods [10], emphasis is on allowance for the vertical wind velocity component w-- one of the principal factors in virtually any circulation mechanism in the equatorial zone. "Qualitative" considerations on the role of w can lead to contradictory results. For ex- ample, if we use as a point of departure the general idea of a PBL in which the wind changes with altitude in such a way that there is "pumping" of air, it must be assumed that on the BL boundary near the equator there is a pos- itive vertical velocity component (it gives rise to or favors the convection process). On the other hand, a very simple model of a one-dimensional BL [13], not taking into account the principal peculiarities in the dynamics of the equatorial atmosphere, leads to a solution with general descending . ' movement. The idea behin~ the formulation of such models is as follows. We will examine a constant zonal flow in the free atmosphc:re (u = Up = const). The system of equations for the BL is taken in the form a?, a~ � av v a~ v av am 7�' az az, ~ w ds + 2~� ro u=~ az-; dy -F- ~s = 1 - Here z is the vertical coordinate; y is the horizontal cornponPnt, reckoned from the equator to the north; r~ = 6.37�106 m is the earth's mean radius. System (1) must describe the process in the three-dimensional boundary layer because all three wind velocity components are present. In order to reduce - _ the problem to a plane BL system (1) is integrated for y from the equator - (y = 0) to some value y= L(N 500-1,000 km), where the BL structure is close to the structure of the Ekman BL. With y= L the zonal component of wind velocity outside the BL can be determined approximately from the geo- strophic relationship. By postulating the nondependence of each of the aver- aged components on y, we find that the nonlinear terms of the equations of system (1), integrated for y, remain unchanged in form. Thus, for the values c t ~ t - u- ~ udy; v- L f vdy; w= wdy 0 0 0 we obtain the system of equations (the lines over u, v and w have been dropp- ed) da d"u dv 2 w L d~ v dur ~ ~ 2 ) . w dz dz1 ' ce dz r~ a tt='~ dz~ ' L ds ' The ~X-parameter has been introduced for the sake of universality. Its value is dependent on the hypothesis determining the "averaged" value of the zon- al velocity component. In [13] it is assumed t?~at OC = 1. A more accurate value is oC= 1/2; in this case the procedure of removing the mean values from t}ie integral sign is identical for all the terma in yystem (1). 11 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY If we exclude u and v, we obtain one nonlinear equation for w in the form d~ ri~tr~ rf~ w d d~w d~ w , ~ dz~ ~.w d:~ ~is~ ~ w ~is I.~ ds~ - v ds~ 1~ ( 3) It must be solved with the following boundary conditions: ~ ~ dw d~w this is e uivalent to u= v= z-~' dz - dz~~ Q Lt = 0~� z oo v 0 = = ~s = ~ ~4~ - u= U~ u' - W~ 2 w jo UW = ~ - where we use Woo to denote the w value when z--~ oo ; this value is deter- mined frora the expression ~ a Va, v d~ w l ~ 5) ~ ro ds3 s~m d~~ - dz= I Z (The operator d~ w d= w v--~~-- dz3 dz= is typical for houndary layer prohlems [11]). Equation (3) cannot be pre- :~isely integrated; it is reduced to a nonlinear integrodifferential equa- tion ~ f se~d: � dw-~~i~:-C,~e`. dz-}-C. (6) When C2 =~0 and C1 = 80~? equation (6) has the partial solution w= -Sv/z.] It can be seen from this form of the formula that when w< 0 equation (6) ~ is conveniently integrated approximately, for example, by iterations. In the case of westerly flow (u~ 0=~' du/dz> 0 in the lower layer of the atmosphere) the turbulence is compensated by the vertical flow u, as fol- lows from the first equaCion of system (2) or (1). This means that with u) 0 we should have w< 0. As a first approximation it is proposed in [13] that it be asstnned that - w=-wp = const, which after substitution into system (Z) reduces it to an easily integrable linear equation. The general integral for the deriv- ed equation has the form _w~= _ _ - w_(C, + C., z) e � -I- Cq (;a z-F~ C:, ZZ� . If there is rigorous satisfaction of all the conditions (4), then C4 = C5 = 0; C3 =-W and for the remaining two ~rbitrary constants a contradic- tory system is obtained. Accordingly, Kuo [13] satisfies only one of the conditions and a"correction" is introduced into the second. As a result, for w a solution is obtained which in the adopted notations can be written as follows: 1z FQR OF~'ICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFF'ICIAL USE ONLY _ ~n s ~c+ _ W�� ~ ~ Zl e ~ - 1 . ~ / ~~1 In thia case w~ must be expressed through the parameters of the problem, nut lepending on the zero approximation wp), and specifically U~, 2 w/r0. Using the methods of dimensionality theory [1], it can be shown that W oo is not uniquely determined through these parameters. In ac- _ tuality, the combination i U1_3pV^p(2l' lp $ u' o ~ ~ for any p has the dimensionality of velocity, which in this case can be interpreted as w~. In [13] a value was adopted for w~ corresponding to _ - p= 1/4 in formula (8). This gives ' (9) W~, _ [ ro' U~ ~z]' � Assuming 10 m2 �sec-1; U ov = 4 m�sec'1 we obtain [,1~ N 10-z m� sec 1. It is easy to see that s p= ~ W~ =~~U~~)z ~o ~ 4~ 10-s m�sec'2; . when . ~ 6 when p= b W~ = YU~ 1' v~ 2ro ~ 5� 1(1-~ m� sec 1. With other "reasonable" assumptions the W~ value is unnaturally small (when p= 1/2 W~ ^~2�10-5 m�sec l, when p= 1 Wa,^~10-1U m�sec 1, etc.). - The continuous curves in Fig. 1 reproduce the results from [13] for u, v, w; the last two functions are easily found from system (2) after substit- ution there of the zero approximation for w: u,,, ~ i Q�o 1 - ~ (10) n- U~ 1- e v= L r`~ U~, ze . ~ i ~ The structure of these formulas is similar to the Ekman formulas for the PBL. However, it must be remembered that they were derived at the expense of "coarsening" of both the physical model and the mathematical methods for analysis of the initial system of equations. In addition, the model does not take into account the dynamics of the processes transpiring in the equatorial zone, it is unsuitable with Ua, ~ 0-- an easterly flow, which is typical for the equatorial zone [6, 12]. The Kuo model is not 13 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY ciCed for critical purposes, but as an illustration of the difficulties arising in formulating a BL model in the low latitudes in general and in the equatorial zone in particular. - 8 _ y Iv y I ~ 1 . , , o ~ ` ~ 0 1 - Fig. 1. Vertical profiles of wind velocity components in dimensionless form using the Kuo model [13~. The dashed curve shows the proper value w for small z. . The model can be made more precise, for example, by introducing a parametric dependence on y and stipulating the zero approximation~for w not in the farm of a constant, but selecting a function of altitude which is more real for the lower layer, for example, a power-law function. We will assume 2Q1 = - ?Qlo ( h i~~ ~ 11 ~ \ / where h is the characteristic height whicui must be determined proceeding ' on the basis of dimensionality considerations. As in the case nf velocity, the representation of h is unambiguous through the parameters of the prob- lem 2~�J /rp; v; U ~,p . Specifically for any p ~ - p ,,,+2 p - -'-3 p h-r f~ 1 U� . ~ , ~ ~ The h values "reasonable" from the point of view of interpret~~tion of BI, theory are obtained, for example, for . _ - - - - . p-- 4 h-Y 2~W " U~ 10~~t; for p=- ~ 1i=Y 2~ "a ~ W U:,o ^'0.3�103 m. Nonallowance for the Rossby parameter (p = 0) leads to the un- expected result: h= y/U~ 2 m. It is probably possible to interpret this result as follows: this is the characteristic thickness of the Prandtl boundary layer, which is unrelated to the macrocirculation mechanism. Here tne scales are completely different: W^~0(u). ~ Then, the parameter r) 1. In actuality, from the continuity equation _ 14 FOR OFFICIAL USE ONLY ' i APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 - FOR OFFICIAL USE ONLY _ : ~ _ - ` dy dz 0 and by virtue of the condition v~ z~0 = 0 for the function w the point z~ c) must be zero of a higher order than for v; > r~1. Substituting (11) inro the first equation of system (1), we obtain z wo sl~r - - h; z~ aZ = Y as u= C, + CZ ~ e Y h~ ~~+r) dx. _ o Satisfying the boundary conditions2 we find u_in the form ar~� z i t, 1 Y C vh~ (I -r r) � r�4- l, (j2) - tt = Voo (y) ~ , ~~r+l) where y(p,`q) and r(p) are the incomplete and complete Euler gamma func- _ tions respectively [2]. Substituting (11) and (12) into the second equation of system (1), we ob- - tain for v a linear inhomogeneous equation whose solution can be repre- sented in the form - m~ v=[ V�~)') - A~Y)] u(Y~ z) t~ e ~ I~~ (rt 1) dzi ~13) [10, (Y ) z~ wo Z.+t z e�~~Hj hr Z dz.,, X Q ~ z, Y ) 0 where, for the sake of brevity, we have introduced the notations v~ (Y) = v (y, z) ~ W t�o si~'1 i~ wn s~+l (14 ) A(Y) = bf e Y(r+ 1) ?i= dZ~ f Q(zz, Y) e v h~ r+ 1 dzz~ Q(z, Y) - 2~' a u~Y~ z). (It is easy to confirm that all the integrals entering into (13) and (14) converge: wp, h, r- 1> 0) . Differentiating (13) for y and then integrating for z, we obtain a more pre- cise formula for w(z) than (11): v� (y)-A' (Y) r wozi+i 1 1 da _ w=- r 1 ,~T`vh'(r+l)' r+l~ . (r+l) (15) . . mo ~~1 Z , mo t2+~ 's r' e ~ hr ~r t 1~~' ~Q (~y ~ Y) e v hr (r+1) dia dzZ dz,. ~ 15 - - FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 rux ~rrt~ttu, u~r, uNLx An anal.ysis of formulas (12), (13) and (14) shows that with small z uNz; vN z; w N z2. Tficrefore, evidently, as a zero approximation it i~ naCural to take w=-wp(z/h)2, that is, assume r= 2. In this case the gamma func- tLons ;uil.t be oE the parameter 1/3. Strictly speaking, the dependence of the functions on y is not "purely" parametric the derivatives of y enter into the expressior., for (15). As Ucao and Voo it is possible to take a suitable solution corresponding to a stationary model of circulation; in the simplest case it is possible to obtain a solution, for example, from [3]: (1a, = uo - cu ~ , _ -V~ _ j 01 uo ~ 1 - r _ 0 0 o~y (up < 0 is easterly flow). As can be seen from the solution found (13)-(15), allowance for the nonlin- ear terms, even in very approximate form, by the stipulation of w in the _ form (11), which actually brought about linearization, to some degree made ir possible to reflect the interrelationship between the boundary lay- er and the external circulation mechanism. Now we will proceed to an examination of another model in which the spatial - structure of the BL is manifested more clearly than in the case considered ahove. This linear model is based on the same Ekman idea a correspond- _ ence between dissipative forces and Coriolis accelerations with a stipul- ated pressure gradient 2�/d y. We will write the initial system of equa- tions in the form _ y di 2 w w- 0 Y ~~Z 2 w ro tt a~ = 0 . ~ 16 ~ ~ tly + � ~r' = 0 Here in the first equation we have omitted the term - 2~ y/rp v for two reasons. The basic reason is a peculiarity characteristic for any BL: the ~ basic wall (underlying surface) effect is the appearance of a mass flow in the direction perpendicular to the wall. Accordingly, in the BL the terms with w are more important than the terms with v-- velocity "across" is less than the velocity "along" the main external flow ~v ~u~. Sec- ond, allowance for the discarded term complicates the algebraic part of formulation of the model, which is scarcely justified in such a rough (linear) model, presented for the most part to illustrate the overall qualitative picture and the difficulties in formulating a good model. This ~ simplification should not have a strong effect on the characteristic of particular interest w. We can formally introduce a small parameter in front of the term _2 v rp 16 ~ FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY and seek a solution in the form of a series in powers of this small para- meter: in the zero ap~roximation thia term disappears. Excluding v and w from system (16), we obtain one equation for u � a;" 4 � ? ri = - 2 ~ dy . (17) dz~ w- ay ru In the BL the pressure gradient is naturally considered to be independent of altitude. Assuming that the solution (17) can b~ represented in the form of the product u= ui (J) u~ (z) . ~18~ we find 1_ y (l-I~ro)* u C ~lg') [It will be assumed that ul is a dimensionless function and u2 and there- ~ after v2 and w2 are functions having the dimensionality of velocity.] where ~ is the separation constant. We will determine it from the condi- ~ tion that outside the BL the velocity of the zonal flow should be depend- ent on y in the same way that this is determined in the problem of the circulation mechanism without taking dissipative forces into account. In the simplest case, with a symmetric distribution of pressure u~ y2. Hence ~ _ -3/rp. This means, _ u~ - L? . (19) where L corresponds to the latitude at which u can be determined from the geostrophic relationship. As is easily checked, in this case the pressure gradient will be a function of the type d ~ y3 . ~y For u2(z) w.e obtain the equation dQz~ ~~ro u== ~ . ~20) . The roots of the characteristic equation are S5 + 12 w 2 _ p~ y r0 Sk_~ Y,~o [cos 5 (2k-1)-}-isin j (2k-1),,k-l,?,3,4,5. ~21~ 17 FOR OFPICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY It can be seen from this ~ormula that S1 and S5 have a positive real pa~rt. Accordingly, the arbitrary constants with these so]:utions must be assumed equal to zero. Thus, the solution (20) can be written in the fonn u: ~z) = U~. -F- Cs es' "-I- Ca es3 Z-~- C~ .es` Z, where U~ corresponds to a partial solution of the. inhomogeneous equation. Hence we find : 12 _ � ~2 ~o; - v~ r Voo - ~-J r ~Ia.. 0 0 In order to determine the integration constants it is necessary to stip- ulate three conditions; 1) z=0 u=0~uz=0~ which gives C2'I'C3"I"C4=-U~,. 2~ z=0 w=0-> aZ'-0, ~~e obtain ~ - CZ Cy S3 C4 S4 = O. ' 3) From ttie continuity equation dv dy I:_o = 0 ~ a~ I ~ as - ~ ~ az~ I:=o= 0. :-u Thus, the third equation will be C= SZ C~ S3 -F C~ S4 = O. Solving a system of three equations relative to C2, C3, C4, taking into ac- count the values of the roots Sk, we find their powers S7~ and make use of the Euler formulas et's = cos x�i sin x sin ~ ' Cx=C~=-U~ _ 5 3 ~ ~ -.0,'l76 U� 2 sln 5 T sin 5 � ' sin 3 " C3 U�� 5 s,~ - 0~447 U... 2 sin 5-~ sin 5 ~ 18 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY ~ Combining the reaults, we write u2 in tlie form 1 ' s uz = Um 1- 2 sin cos lU e- t coS 10 ( 2> 2 + 2 e C atn 10 sin 5~os (t cos . where, fo r the sake of brevity, we introduce the notation _ _ . C - s I'l z v~ rp Assuming the value 'V = 10 m2�sec-1 to be usual for the turbulence coef- ficient, we obtain Z.: 10-3z ( z in m) . The working formula for computing u(y, z) can now be represented as u- U~, 1- 0,447 e- s- 0,553 e- = cos (0,951 z)} , where z, y are in km, L^'S00 km. From the first equa~ion in system (16) we find . w � d~ r~ 2 W ds~ ' This means that the dependence of w on altitude will be determined using the formula v - 12 w- ~ 5 ~ ^ _ _ _ 1 ( cos " � e- c - w - U~, - . 2~ ro 2 sin 5-~- cos i~ l 10 ( 23) ~ - 2 sIn 5 e 10 cos (~-I- C cos i~ . From the continuity equation we find the dependence of v on altitude: s 9s = U~ L 1 ( - cos" e- ~ ~ "ro 2 sin ~ cosiu ~ !0 (24) ~ 3 r, a 2 sin 5 e~ sin sin ('1U C cos la 19 FOR OFFICIAL USE ONLY ~ APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-00850R040240040057-3 FOR OFFICIAL USE ONLY Expanding (22), (23) and (24) into a series of powers of z it is easy to . show that with small z we will have uZ I~V zJ v~ N Zl iQ12 Z~I. _ 2 . S. 4 7. 3 _ 4 . - , I 2 3 w I ~ I ~ ~ 11 u 2 ; t t ~ i~ ~ i ,~y y ~ ~ii~ ~ U 0, 5 0,4 0,3 0, 2 y ' j~ I i � ~ ~ ~ 'Q2-Q1 �0,1 0,2 !{J w I i'~3 ~ _ Z. . x ' Fig. 2. Vertical profiles of wind Fig. 3. Hodograph of wind velocity velocity components in dimension- (1), its projection onto the hori- less form using model (16). In order zontal plane (2) and the Ekman spir- to facilitate reading of the dia- al (3). The small cross indicates gram the positive direction for u the limiting value (with z--s c~o) of is indicated to the left. the velocity h~dograph projection. It can be seen that in accordance ~~ith model (16) a tendency to the limiting value occurs far more slow- ly than according to the Ackerblom ~ model. With respect to the dependence on y, for w it evidently will be the same as for u, but for 3 Y ~ N 3L3 . 20 FOR OFFICIAL USE ONLY ; ~ ~ I- APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000240040057-3 FOR OFFICIAL USE ONLY In general, from the point of view of BL theo ry the function U oo(y) is considered stipulated and it can be selected, to a certain degree, arbi- trarily. If in (18') it is assumed that 1- ~ r0 = 0, that is, if it is _ assumed that Un, = const, then only the one function v will be dependen.t on y. (In this case d~/~ y ti y). The dependence of the solutions on z does not change. Such a model is ~ustified by the fact that when y= 0, that _ is, on the equator, the boundary layer does not disappear: the components u and w are not different from zero. In such models w-?0 when z-y ~~o, which cax?no t be the case for a plane boun- dary layer; in a three-dimensional case there is a possibility for compen- sating the vertical flow due to a change in the second component of hori- - zontal velocity. In these models u can be of any sign. Figure 2 shows the u2( v2( and w2(~ ) profiles in dimensionless values. Normalization was carried out for the corresponding factors 'r 12 ,or -,r ~ - 3-3 ~ Um, uco ~-ro ~I v~!'u ) ' ~ vr3 . ~ 0 Figure 3 is a velocity hodograph for y= 1 and its pro~ection onto the plane xoy. It can be seen clearly from the figure that w-r0 when z-+~o due to the fact that aloft v fluctuates near zero with an attenuating amplitude. It should be noted that the main maximum w falls on the value ~ N 2.0. In dimensional values this mr:ans that in the layer 1.5 to 3.0 km vertical movements are particularly appreciable. This agrees fairly well with the results obtained during GATE (Atlantic Tropical Experiment - an internation- al program carried out in the summer of 1974 [12]). One of the characteristics which is usually determined in models of the at- mospheric boundary layer is the angle between the isobar and the limiting - position of the velocity vector when z= 0, that is, at the earth's sur- = face. In this case i~ is correct to speak of the angle between the wind di- ~ rection outside the BL and the limiting value when z= 0, although actual- ~ ly these concepts coincide here: the isobars run parallel to the equator and the principal external flow is zonal. In the classical Ackerblom model [7J, as is well known (when y= const) this. angle is a= 45�. However, in general it is dependent on the turbulence regime (on stratification) in the - boundary layer [9], but it is less than 45�. It is easy to compute this angle for the considered model as well: 5 l - sin-" tga = lim v= lim v.~ (z) = L~ 54 ~ lu - x..o u z~o u: (Z) ~ r~~ 2 tg ~ sin Z'~ . lU 10 , Digressing from the dimensionless factor, which with the selected values of the parameters is the value 0(1), we find a~ 31�. To be sure, this is only an approximate value whose main sense is that it is less than 45� the model value for the~PBL characteristic outside the equatorial zone. 21 ~ FOR OFFICIAL USE ONLY - APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-00850R040240040057-3 FOR OFFICIAL USE ONLY llespite a number of shortcomings the considered model of a three-dimension- al boundary layer reflects some qualitative aspects of the processes trans- ' piring in.the lower layers of the troposphere in the equatorial zone: a maximum of vertical movements in the layer 1.5-3.0 km; a tendency to its limiting values of the ~aind velocity components with altitude slower in camparison with the PBL; a lesser an gle between the main flow outside the BL and flow at the earth's surface. I consider it my duty to express appreciation to I. B. Kazitskaya for ae- sistance in finalizing the article. BIBLIOGRAPHY 1. Barenblatt, G. I.~ PODOBIYE, AVTOMODEL'NOST', PROMEZHUTOCHNAYA ASIMPTO- TIKA (Similarity, Self-Similarity, Intermediate Asymptotic Behavior), Ltningrad, Gidrometeoizdat, 1978. ~ 2. Gradshteyn, I. S., Ryzhik, I. M., TABLITSY INTEGRALOV, SjJr4f, F,YADOV I PROIZVEDENIY (Tables of Integrals, Sums, Series and Products), Moscow, . Fizmatgiz, 1962. ~ _ 3. Dobryshman, Ye. M., "Some Peculiarities of the Pr~ssure and Wind Fields in the Equatorial Region," METEOROLOGICHESKIYE ISSLEDOVANIYA (Meteoro- logical Investigatians), No 16, Moscow, 1968. 4. Dobryshman, Ye. M., "Determination of the Width of the Equatorial Zone," ME'iEOROLOGIYA I GIDROLOGIYA (Meteorology and Hydrology), No 12, 1973. 5. Dobryshman, Ye. M., "Dynamics of Atmospheric Processes in the Tropical Zone 5-15� in Latitude," METEORO~OGIYA I GIDROLOGIYA, No 4, 1975. ' 6. Zaychikov, B. P., Romanov, Yu. A., "Peculiarities of the Wind Regime, Temperature and Air Humidity in the Equatorial Region of Central Atlan- tic," TROPEKS-74 (TROPEX-74), Vol I, Leningrad, Gidrometeoizdat, 1974. 7. Koshmider, G., DINAMICHESKAYA METEOROLOGIYA (Dynamic Meteorology), Mos- cow, Gosizdat, 1938. 8. Krivelevich, L. M., Laykhtman, D. L., "Meridional Structure of the Plan- ~ etary Boundary Layer of the AtmosphQre in the Low Latitudes," IZV. AN SSSR, FIZIKA ATMOSFERY I OKEANA (News of the USSR Academy of Sciences, - _ Physics of the Atmosphere and Ocean), Vol 13, No 7[year not given]. = 9. rionin, A. S., Yaglom, A. M., STATISTICHESKAYA GIDROMEKHANIKA (Statis- tical Hydromechanics), Moscow, Nauka., Part I, 1966, Part II, 1967. 10. Pushistov, P. Yu., "Results of Numerical Modeling of Stationary Circula- tion of the Atmosphere in the Equatorial Region," IZV. AN SSSR, FIZIKA ATMOSFERY I OKEANA, Vol 9, No 3, 1973. 22 - _ FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 ~ FOR OFFICIAL USE ONLY 11. Roze, Kibel', Kochin, GIDRODINAMIKA (Hydrodynamics), Vol II, Moscow, Gostekhizdat, 1954. 12. TROPEKS-72, T':OPEKS-74 (Tropex-72~ Tropex-74), edited by M. A. P~etro- syants, Lenin�rad, Gidrometeoizdat, 1974, 1976. 13. Kuo, H. L., "On the Planetary Boundary Layer of the Equator," J. ATMOS. SCI., Vol 30, No 1, 1973. 23 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY - UDC 551.511.12 SMALL OSCILLATIONS OF THE POLYTROPIC ATMOSPHERE AND THE FILTERING ROLE OF THE IiYDROSTATIC APPROXIMATION rioscow METEOROLOGIYA I GIDROLOGIYA in Russian No 11, Nov 79 pp 23-33 [Article by Candidate of Physical and Mathematical Sciences V. M. Kadysh- ' - nikov, USSR Hydrometeorological Scientific Research Center, submitted for publication 23 May 1979] Abstract: The article gives a comparison of ~ the frequency spectra of wave moven~ents of the polytropic atmosphere in nonhydrostatic and hydrostatic cases. In the first case the _ dependence of the solution on altitude.is descrilied by confluent hypergeometric func- tions, and in the second case by Bessel functions. It is shown.that the hydrostatic approxitnation, as in an isothermic atmosphere, filters out acoustic oscillations and consid- erably distorts the first modes (that is, those changing least with altitude) of gravitational oscillations only with a wavelength shorter than 100 km. (Text] In [6] the authors obtained an anal~~i,;~1 s~Iu~~~i~ of the problem of small oscillations of a nonhydrostatic atmosphere under the condition that the linearization of the equations of hydrothermodynamics is carried out relative to an isothermic state. An analysis of this solution indicated that the atmosphere is characterized by wave processes of different nature, to wit: there are high-frequency acoustic and low-frequency gravitational � oscillationS. It was also demonstrated that the hydrostaticity hypothesis filters out acoustic wave~ and the spectrum of gravitational oscillations is distorted the lesser ths longer the corresponding wave. In.[2] a sim- ilar problem was examined for a model of the atmosphere with a real te~- perature stratification, but the solution obtained to a considerable de- gree is qualitative. In our article [4] we gave an analytical solution for a neutrally stratified atmosphere. In such an atmosphere there are no grav- itational waves. The objectivE of this article is solution of the problem of small oscillations in a stably stratified, nonhydrostatic, polytropic 24 FOR OFFICIAL USE ONLY ~ APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY _ atmosphere and a quantitative evaluation of the role of a hydrostatic ap- proximation in this case. A~ d~~monRtr~itc~d in (fi~~ the :;yatem of equatione tn hydrothermodynamice J.in- e,irL~~~~ rel.utlvc lu :i :~L.tCc ~~1~ t'cal h~tv lltr 1'ur.m rc~ pX lv, v~ py - lu, ~~e P: - P g, Pe = - (ux + vy ~s), Pe = - 2ro - e" (us vy w:), (1) . where u, v, w are F u*, P v*, ~ w*:/~ is the density of the main state (its _ pressure and temperature are p and T; all these three functions are depend- ent only on z), and uk, v*, w* are velocity vector components; p and Y are - the deviations of pressure and density from the corresponding values of the main state; 1~ is the Coriolis parameter (constant), g is the accelera- tion of free falling, c2 p/p ,~.is the ratio of heat capacities, (r = i ?~R( ~'a - r) is the stratification parameter, R is the gas constant, Ya is - the dry adiabatic vertical temperature gradient, ~'=-dT/dz. We will solve the problem with initial conditions periodic relative to x, y. Now we will consider formulation of boundary conditions for z. We will a~sume that the underlying surface is impermeable, that is, at the earth w = 0. This gives ~=0 (z=0). (2) For formulation of the second condition we note that the following formula follows from system (1) ~6]: ~ - ~ ~r + ~u# P)X + Iv'~ P),, -'t- P): _ ~ where the quadratic form (energy) is ~x + v~ ~ L 1 , ( ) E 2 0 2-c: [P'" B~P - P)",. 3 P If, accordingly, from the initial data it is required that , ~ _ ~ Edx dy ds < oo (4) (D is the periodicity region), then with condition (2) and the condition ~~"~*P~ z-r a~ ~ that is ~p ~ (Z -s c~o ( 5 ) the energy of the svstem will be conserved. Moreover, below we will confirm that if we take (4), rather than (5) for all t~ 0 as the second condition, condition (S) will be automatically satisfied, that is, from the limitation of energy follows its conservation. _ Thus, we will seek a periodic solution of the system of equations (1) under conditions (2) and (5). As already mentioned, the y value falling between 0 and ya is considered constant, so that r"= const. The altitude of the 25 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY polytropic atmosphere Ii is finite and equal to Tp/r , where T~ = Tlz a 0� We will represent all the functions in the form f~x, y, z, t) =f ~z)expi (k jz+k~ y-~. t). It follows from the equations of horizontal m~tion that the amplitude of plane divergence is, t7lk': ~ - ,,1 _ P, , where k2 = 1cX + ky. We introduce three-dimensional divergence ~G= uX + vy + wZ. We have - . _ 11? k= ~ . Y z= a- l, P-F- ( 6) The third equation of motion (with the continuity equation taken into ac- count) and the heat influx equation give, respectively -i?.z~-{-pr- ~ ~=U, ~ ~ ~ (8) -ilp-~I'w-}-c'X-O. Lde have a system of three ordinary differential equations (6)-(8) for the functions 'C, p and w. We will write an equation for some one of them. Since the coefficients of the system are variable (c2 is a function of 2), for each function there will be a,specific equation. It is convenient to examine the equation for '~C:; its eigenfunctions are single-term functions (it goes without saying that the final solution of the problem is not de- pendent on the choice). We will express p from (8). Then we differentiate this expression and we will substitute the p and pZ values into (6) and (7). We will solve the derived system of two equations as an algebraic system relative to w and wz. In particular, we have .(c? = s),. ~ xRY(a=--1=) ,g+1' ` ~ a~ (a~-n~-k~r�- xR7~ - 1 - ~9~ xRY -l'-') S~ ~-s%~'., Differentiating this expression and comparing it with the already deter- _ mined wZ, we obtain g + ~ ~ k~ . ~k= I' S7,~s- ( Sl - 2~ Si ~a=-t=) S + s2 + Sl ~A:_~_~,"/. _ (10) : where ~ _ s-c2=so-s~z (so=xR3`o, s~=xRy). . In the derivation of this equation we had to exclude from consideration the roots of the equation 26 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY ~.(7~~-l') ~~2~~2-1~)-k~I'2] =0. (11) They must be investigated separately. It is easy to show that the frequency 0 is the solution of our problem. It corresponds to a geostrophic hy- drostatic case. Howe^er, the ~lvalues, bringing the terms in parentheses and brackets to 0 in (11), are not solutions: in the first case from (7) and (8) for w and } we obtain a homogeneous algebraic system of equations with a determinant different from zer.o only at the one point where the linear function c2 is g(-/,~ 2, whereas according to the equations of motion p t'0; in the second case, from (6)-(8) we obtain a similar system for Y and xs which can be different from 0 also only at one point where c2 is r 2/ it therefore follows that there are no continuou~ nontrivial solu- tions in the two cases. The solution of the fundamental equation (10) of the problem is x=CiJi~S)+CsJs~s), ~12) where C1 and C2 are integration constants and the linearly independent special solutions are [5] _ as ` Y~ = e ~(1 - 2- 3, 2- 6, as . aS ~ 13~ _ Y= = S�-~ e 1~( 2 b+ as), cind the function of the three arguments ~(r, m, x) is a confl�ient hyper- geometric function (8), determined by a Pochhammer series ~~=1 ~ ~ r(r-i-1). . .(r+n-1) x" . m (m-r 1). . . (nc+n-1) n!' (14~ � n-1 Here we have introduced the nota~ions a= 2 x k S~~ a= ~ a= r-1 +~k= r b_ g+ r (15 ) S Si 1 ~i:= - ' 'l A ks, ) ' s~ Using (9), and in particular, the fact that when s~0 wNs ,L S- (b - 1)~1', and also that p�~~w (this, incidentally, confirms the advantage of reduction of system (6)-(8) to an equation specifically for x.), we substitute the resulting solution into condition (5), in a polytropic atmosphere exist- ing, as already noted, when s= 0. Since P-~ sb, this gives Qi~i +a~zC~C~+asC2 =0, _ and al is a linear combination of the powers s2-b, sl-b, s-b, a12 the powers sl, s~, s'1; a2 tlie powers sb, sb'1, sb-2. Since x 7=-1= ~ Yrt _~]2~ b - x-1 ~ 7 27 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY C1 = 0 and the C2 coefficient can be arbitrary. Thus, with an accuracy to a constant,factor the solution of equation (10), satisfying condition (5), has the form a % =S�-1 e `s ~ ( 2 - P~ b, a:)� (16) IE we found the remaining amplitudes from the function x, determined by Eormulas (12)-(15), and formed energy E from formula (3), condition (4) Cor all t> 0 would yield this same solution. Accordingly, we see that its con~ervation will follow from the condition of restriction on the Lot~l energy. We wi11 make two coimnents. First: the solution (16) was obtained proceed- ing on the basis that (13) are special solutions. But the latter is cor- rect under the condition that the difference in the roots of the determin- ing equation for (10), but the roots in this ca~e are b- 1 and 0., is not a whole number. Otherwise the second argument of the ~ function for one of yi, in accordance with (13), is a whole negative number, and - in accordance with (14) the eigenfunction does not exist. But, for e~ ample, with Y close to t:he real value 5.83 K/km, this is a whole number. Nevertheless, i.n this case the solution of equation (10) under the condi- tion (S) is the function (16). In actuality, in the neighborhood of the regular singularity s= 0 in this exceptional case as one independent so- luti.on, as before it is necessary to take Y2 (this corresponds to the greater root of the dete~nining equation), and as the other, not yl, but the expression P(s) + ayl ln s, where a is some number and P(s) is a power _ series with a free term [7]. But only this, and specifically the presence of a free term in the fun~ction yl, was used in the proof that C1 = 0. The second remark is as follows: we will assume that 0< Y< Y. When Y= ya there is also a solution of the formulated problem [4). At the same time, oscillations of the isothermic atmosphere with restricted energy are not possible. In actuality, in this case c~ = const and the equation for p is the same as for ~C,: g-F~ 1' k~ g i' Pu P: [ (1- P= 0. We denote g+ f/2 c2 by r and introduce R in such a way that the character- istic equation assumes the form x2 + 2rx + r2 - R= 0. There can be three forms of the dependence of p on z: e-~= (C, stn z j~- R-E- C, cos z~) (R < 0), e-.:(C~+C:z) (R=~)~ e-?z I C, es C= e- : vR ~~R ~ p~, . 2s FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY Taking into account that w-~ pZ + g/c2 p, and J~ e'2rZ, we conclude that con- dition (5) in the first two cases is not satisfied by the two apecial solu- tions, and in the third case, although it is satistied, the lower boundary condition (2) requires from it that r(2-~~x12 _ Ig+rl~ _ k~a: gr - L 2c' ~ ~ l~' ~ +~'-C'(1-~,~,1. / If this condition is satisfied, then w= 0. The corresponding ~l values can be found. However, it is clear that it is impossible to solve the problem with any general initial conditions. Now we substitute solution (16) into the lower boundary condition (2). Us- ing (9), we obtain the dispersion expression in the form ~1-2ti)~r? -3~-1, b~-l,al-(1- St~ c~(2 b~ al, ~17~ ~ i 1 \ / where ~ St = Si V~S~ ~)2 - g r, a=' as ~ s=so � ~18~ . In this case we considered that, in accordance with (15), i' - a S . k - YR~ -1~ . (19) - -"2H ' -2H ~ Equation (17) determines in the plane of the variables a, ~ two families of curves corresponding to the values S+ and S_. The first family corres- ponds to acoustic waves because it owes its existence to atmospheric com- pressibility. In actuality, we will assume that the atmosphere is incom- pressible, that is, X--~ ~[6]. This means that the parameters r, sp, sl also tend to infinity. This is in no way reflected in determination of the ~~G para~qeter in (15) . At the same time gk ( 7~ 1~1 ~ ~ 2~7i~ e=-l~ Accordingly, in accordance with (18), S+ ->c~o. Thus, in accordance with (19), ~~oo, that is, the assumption made actually filters out oscillations of this family. In this case _ S_~ k The second family corresponds to gravitational waves because it owes its existence to the presence of stratification. In actuality, with 0 the parameter S= 0, so that a= 0. It follows from the conservation of energy that all ~1 are real. Below it will be dempnstrated that ~l 2~~,2, that is, oC and Q are also real. It can be seen from (15) that oe~9 > 0. Now we will first examine the case a~ 0, 0, and then we will show that the case a< 0, ,8 0 is a sufficiently large whole number. We note that the several first solutions of equations (Zl) and (22) can be parasitic because sufficiently small a values under the condition ~~i b can lead to ~ values which are too small. The hyperbolic equations approx- imately following from (21) and (22) can be ~oined into one: 30 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY r +c / . 20~ I 11~ - L 2`n 4 = 9t� (23) Here and in the text which follows the subscripts and will relate ~ to the acoustic and gravitational families respectively. In case 2 we make direct use of the series (14). There are two possibilit- ies: is limited and ~--;~~o . For the first for acoustic waves there is no solution, whereas for gravitational waves we obtain f3 = b/2. And in general, if S= S, the entire straight line j3 = b/2 forroally satisfies equation (17). But the corresponding frequency curve, reducing the bracket- ed term in (11) to zero, as already noted, is not a solution. With the sec- ond possibility for acoustic waves there is again no salution, whereas for gravitational waves we have a~b+t~r a~ - g+( b-I- 1)1" ( 24 ) = In this case it is assumed that ~p ~ 1, that is (~-s 0. In case 3, in accordance with [8], we have ~ b b~ Q ~ (b-111 + (b-11! e� . , Z ~ 2 +a+l~~ (_a~~ ( ~ -~1--~~~ sb +a ~ Therefore, from equation (17)bwe obtain (-1)~ 1 a2 q 8a (25) ( 2 +~-t~i s+r ~l -~-i~! In order for the riglit-hand side of quation (25) with oc.-~ro to remain, like ttie left-hand side, finite (the a 2~ factor on the left-hand side cannot ensure the corresponding increase), it is necessary that the denominator on the right-hand side increase without limit, compensating the increase in the numerator. Hence - P=Pn-` I~R - ~ -}-iyl~ ~26~ ~ 1 where n is any whole non-negative number and I�I ~ 1. In this case b _ 1 ~ ~ ~n \l - n!` ' Substituting (26) into equation (25), we find e-a Qo+1R s_r E- (b+n-1)! n! S+~' With ~n ~ ' St =g 2 s, n�~(g 2 ~ I~ S,n (g i') 7. (s,~y)', / Therefore, for any n, taking into account that g>(', we have S+, Since, moreover, the root is a value not less than g-(~/2 + sln, for n~ 1 we have, on the contrary, S . This means that when a-s~~3+ --y ~n - 0, but -s ~n + 0 � - 31 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 r~rA vrrl~,lnL UDr, VriLi Since the different modes of the long-wave part of the spectrum are charac- terized by the parameter q, it is convenient to find 'lfrom the equation obtained by multiplying the first two formulas (15). Since k is small, we obtain ' - s' q+ + - H ~ (27) = 1= + g ~ k3. S~ 9_ (28) y - r= . In the short-wave part of the spectrum the modes are characterized by a constancy of the parameter ~(~~~n); using this parameter, directly from ~ (15) we have - l + 1' 4 +S�kz ' If, in particular, n is sufficiently large, then ~ l/ ~ = 1-~ Y q~- sl ~rt 2 g I' j k=, - ~2 _ l + ~ 4 + 4 t 'l k., . t We assumed OL ~ 0, p? 0. It follows from (15) that a simultaneous change in � t;l~e signs on oc. and ,B is equivalent to a replacement of a by - a. Since ~ a11 our formulas for a(k) contain ~ 2(that is, the frequency was deter- mined with an accur,3cy to the sign), we evidently obtain no new solutions. We note that the eigenfunction (~.6) in this case.also remains unchanged because ~(r, m, x) = eX ~(m - r, m, -x) [8]. We also assumed that a2 ~,Q.2. We will confirm this. We will assume that this is not so. Then a=-i o~ 1, ~_-i ~1~ 0~1 and ~1 are real positive numbers. It can be shown that if the wavelength does not exceed several tens of thousands of kilometers, even for very small f~ 0 O~C that is, it is possible to use the asymptotic form (20).. We again obtain the equations (21) and (22). Their left-hand side is complex, whereas the ~ right-hand side in the first case is real, but in the second case pure- - ly fictitious, because S= 2i ?crll~. Thus, in actuality, ~ 2 7.~-2~ It was demonstrated in [4] that this inequality also is observed with r= 0. Now we will proceed to a clarification of the problem of how the frequency � spectrum changes if the hydrostaticity hypothesis is~adopted. Here it is possible to use the results in [3J, where the corresponding spectrum was constructed. However, we will proceed differently and construct it by pro- ceeding directly from equations (1) and the boundary conditions (2j, (5). 32 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY W~~ wlll NuPply th~~ I~~ft-i~rind 9ide oE the thlyd c~quatiun of mc~t[an in (1) witli tiie Eactur E�1 ;;o ll~at Lhe transitton to liydroetatice will cnrrespond to the limiting transition ~2 -s0. First, the frequency ~1 = 0 as before is a solution. Second, we will again obtain equation (10) and the dispersion expression (17), but formulas (15) assume the form 2~ca a E2 + gk~ f as = S~ ~ _ , (15' ) sl y~~~ - . 2 ~c~ ks~ - that is, seemingly not entering into the combination ~2 -,Q,2, is sup- plied with the factor � 2. This means that in any range of wavelengths with - E-2 ~ 1 we have the case 1, that is, we have formulas (21)-(23). Similarly to (27), (28) we obtain . 2 _ S~ 4+ _ Hc~ , ~27~) _ l-}- s~ 9 g r k~' ( 28' ) ~ - c~ With E2 t0 the frequencies (27') tend to infinity, that is, the correspond- ing oscillations are filtered out. Thus, in a hydrostatic polytropic atmo- sphere all the solutions are given by the formula ~ a - ~a j~ S 9_ k, ~ 28� ~ H _ However, for its derivation it is not necessary to assume a smallness of k. It is sufficient from (15') to form q=� 8, assume E2 = 0 and solve the derived formula for a2. The q values can be found from equation (22), - bearing in mind that q= S 2/4. However, it is better (at the limit a, f 0 and the asymptotic formula (20) even with ~~-~-oo, generally speaking, is noe true) , to discard i12, except the combination -~Q,2, directl;~ in equation (10) and obtain a corresponding solution. Instead of (10), we ob- tain the Bessel equation ~ _ gk, r ~ s~s-(b-2)~+ S~ ~~~-t:~ ~=0. Its solutioii, satisfying the upper boundary condition, is b-1 " r- X�S ~ ~�-12 s, y ~,g ~ S~, and the lower boundary condition gives a dispersion equation in the form / _Stz ( 2k ~gf'so~ ~6 lz~ - 2 ~ ?b-: ~Z~ lz = Sl _ . (17' ) 33 FOA OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007/02108: CIA-RDP82-00850R000200040057-3 rvcc urrtt,tr~t~ u~~ ULVLT If, in particular, z~l, using the corresponding asymptotic formula [8], from (17') we obtain equation (22) for z. And formula (28") follows f rom tlia der~rmi,,,lti~~~ of z. We see that it virtually coincides with formula (28) for the frequency of long gravitational waves because the parameter ~,2 in the denominator of the latter can b~ neglected: for example, if T~ = 273 K, r= 5.83 K/km, then slq /H = 0.5�10-4 q; c~ > 1, whereas ~t2 = 10-8. A comparison of for- mulas (28) and (28") shows that the influence of the hydrostatic approx- imation decreases not only with an increase in wavelength, but also with an increase in the number of the mode denoting tlie depree of "vertical" variability of the solution. All the waves with the frequencies (28"), except the zero mode, which is designated so because ~j~ corresponds to its short waves, whereas the long = waves of all other modes are described by formulas (23), beginning with n = 1, are gravitational. In actuality, according to (28") with ~'~0 a2~.~2, but since these frequencies are not solutions, the corresponding oscilla- tions are filtered out. The zero mode, which in a general case consists of two branches (the short waves are acoustic, the long waves are gravitation- a1) cannot be regarded as purely gravitational: with j~=~ 0 it is gravita- t1UI1~11~ since in the long-wave range it corresponds to gravitational waves i~i a rionhydrostatic atmosphere, but with (--i0 it is not filtered out, but E,~~sses into a solution corresponding to acoustic waves in a nonhydrostatic - atmosphere. Now we will consider how this occurs. If r--3 0, then it also tlas one value q-~ 0, and specifically (24). Therefore, it follows from (28") that Az _ ls + BTo (r. - l) kZ t29~ y' ~ ~ i t}iat is, this mode actually exists also when ~'-i ~'a, but in this case it already becomes acoustic. In actuality, if 0, then, in accordance with (15'), the hydrostatic approximation means a s~ 1�a ~ 1. It follows from (17), in accordance with (14), that . b 4~b+1~ a. Therefore ~ ~ _ _ b 2 b a~ Nx 9= 4~b~-1) a- b+ 1~2-[- k2E~' ~ Substituting this into formula (27') for acoustic frequencies; we obtain (29) with the replacement of y by Ya. Thus, this mode must be considered acoustic-gravitational. Thus, the hydrostatic approximation filters out the acoustic waves and the ~ frequencies of the gravitational aald acoustic-gravitational waves are dis- torted the less the longer the waue and the higher the number of the mode. 34 _ FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-00850R040240040057-3 FOR OFFICIAL USE ONLY . ~7/C io t _ ~ Figure 1 shows the frequency curves 600 i computed using the dispersion equa- - ~ tions: (17) for a nonhydrostatic 0 i ~ atmosphere (acoustic part of the i i spectrum is tlot represented) and ' (17') for a hydrostatic atmo- S00- ~ sphere. In complete accordance with ~ the asymptotic formulas (28) and i ~ (28") the difFerence in the frequen- i cy spectrum decreases with an in- - y00 ~ crease in the mode nwnber and an increase in wavelength. Only the zero mode in a considerable part of ~ the spectrum is raore hydrostatic than, for example, the first mode. 300 ` But with the selected Y(5.83 K/km) ~ ` it is acoustic to L= 187 km, so - ~1 - ~ that its behavior is not describ- ed Uy formula (28). We see that - 2~ ~pi~~ gravitational waves longer than ap- 1~~~~, Proximately 100 km with great ac- curacy are hydrostatic. _ ~ 9 ~ - 0 ` ` 900 ~ ~ ~ � ~ _ , 0 SO 1J^ ;50 G r,rf Fig. 1. Wave frequencies of gravi- tational (1) and (for the zero mode) - the acoustic (2) nonhydrostatic os- cillations and wave frequencies of fiydrostatic oscillaLiuns (3) in de- pendence on wavelength in a polytropic ( Y= 5.83 K/km) atmosphere. The fig- . ures on the curves denote the mode number. BIBLIOGRAPHY ~ 1. Gradshteyn, I. S., Ryzhik, I. M., TABLITSY INTEGRALOV, SUNihf, RYADOV I - PROIZVEDENIY (Tables of Integrals, Sums, Series and Derivatives), Mos- cow, Fizmatgiz, 1352. 2. Dikiy, L. A., TEORIYA KGLEBANIY ZEMNOY ATMOSFERY (Theory of Os~illa- tions of the ~arth's Atmo~phere), Leningrad, Gidrometeoizdat, 1969. 35 - FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-00850R040240040057-3 FOR OFFICIAL USE ONLY - 3. Kadyshnikov, V. M., "Reason for Closeness of the Wind to Geostrophic," MET~OROLOGIYA I GIUROLOG.T..YA (Meteo~ology and Hydrology), No 9, 1977. 4. Kadyshnikov, V. M., "Modeling of Hydrostaticity of Large-Scale Atmo- spheric Processes," METEOROLOGIYA I GIDROLOGIYA, No 1, 1979. 5. Kamke, E., SPRAVOCHNIK PO OBYKNOVENNYM DIFFERENTSIAL'NYM URAVNENIYAM - (Handbook on Ordinary Differential Equations), Translated from German, _ Moscow, Nauka, 1971. 6. Monin, A. S., Obukhov, A. M., "Small Atmospheric Oscil?ations and Adap- tation of Meteorological Fields," IZV. AN SSSR, SER. GEOFIZ. (News ot the USSR Academy of Sciences, Geophysical Series), No 11, 1958. _ 7. Smirnov, V. I., KUItS VYSSHEY MATEMATIKI (Course in Higher Mathematics), ~~ol 3, Part 2, Moscow, Gostoptekhizdat, 1957. 8. .t~hnke, F., Emde, F., Loesch, F., SPETSIAL'NYYE FUNKTSII (Special Func- tions), translated from German, Moscow, Nauka, 1968. 36 FOR O~FICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY UDC 551.(512.2+54+55) _ � CORRELATION BETWEEN MINIMUM PRESSURE AND MAXIMUM faIND VELOCITY IN TROPICAL CYCLONES Muscow METEOROLOGIYA I GIDROLOGIYA in Russian No 11, Nov 79 pp 34-41 [Article by Candidate of Physical and Mathematical Sciences V. M. Radike- cich and Candidate of Geographical Sc3.ences G. G. Tarakanov, Leningrad Hydrometeorological Institute, submitted for publication 2 April 1979] The at~,thors propose a model for describing the dependence of maximum wind velocity V~ � in a tropical cyclone on the radius of the region with maximum wind velocity rp, the Coriolis parameter cJ Z, pressure P~ at the center of a tropical cyclone and pressure P1 at a great distance from the center rd. In addition to a precise expression for deter- mi.ning V~, the authors also derived an ap- ~ proximate and simplified expression. The re- sults of the computations agree satisfactor- - ily with observations. The systematic discrep- ancy between computations and observations falls within the range of discrepancies be- ~ tween the wind in the free atmosphere and the surface wind and varies from 0.78 to 0.65. ~ [Text] Tropical cyclones develop in the tropical and equatorial zones be- tween 22�S and 35�N, except for a narrow equatorial zone (about 2�N-2�S). ' However, the main mass of tr~pical cyclones (87%) is associated with a narrower zone taking in the region between 3 and 20�N and S[2J. Observa- tions show that in each tropical cyclone (TC) there is a so-called ring of naximum wind (if it is sufficiently narrow it is called the maximum wind circle). The maximum wind velocities can attain 80-100 m/sec. It is evi- dent that the instrumental measurement of such vElocities involves great difficulties. Measurements.and the presently known theoretical models o~ a tropical cyclone show that pressure decreases toward the center of the cyclone in conformity to a parabolic law. Taking this circumstan~e into 37 FOR OFFICIAL USE ONLY ~ APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY accuunt, and also that pressure measurements are usually inadequate for a precise determination of the pressure gradient in ~he neighborhood of the maximum wind zone, there have been numerous attempts to relate the maximum wind velocity (Vm) to the pressure drop at the center (Pp) and on the peri- ~ phery (P1) of the TC. A physical validation of the correlation between Vm and (P1 - Pp) can b~ and T as follows: Z _ g ~ dz T~:t-~t~Z~ ' P-poe . a~,~ R z= Zoa -I- Tn.~ r( P 1 g 1. . a~,2ll P~,21 ~ ~ a, R zo_O~ x'_ Q~~~PI ~ g _ 54 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY _ Computations sho~~ ti~at the relative ciianges in pressure in the mesosphere ~~nd u~~n~~r atr~~t~~:~~h~~re remuin fipproxim~~te].y conHtrint ~ which, r~y inclie�f~te~d ubuvey i;~ riLtrihutable to the fact that we are solving the Cauchy problem. The changes in temperature in p-coordinates in the presence of ascending and descending currents (3-5 cm/sec) are given in Fig. 1; it can be seen that: 1) considerable changes in the temperature of the mesopause (20 K) (c, d) cause a weak temperature change 1.5�K at z~75 km and about 0.6�K at other levels; the geopotential changes insignificantly (0.1-0.4 km); 2) a temperature increase at the stratopause by 10�K (b,c) causes a tem- perature change in the entire layer by approximately 0.5�K and virtually does not change tlie geopotential; 3) a temperature increase of the lower stratosphere by 10 K(a, c) does not change geopotential and weakly disturbs stratospheric temperature (0.3-0.4 K) . PHb mb ~ ~ '3 ~ ~ 9 ~ 9 _ ~ ~ ~ J 1 50 a 2~~\ b 2 a) ~ 9 b) ~ ~ Z ~ ~ a 2 b y r, O,SS / ~ r~. , ~a same not- = ations in ~J ~ ' ' ' ' ' ~ 0,01 ~ Fig. 2 u) ~ 2) ti BJ ~ Z~ \ J d ~~y 3 d 1 2 ~~~2 ~ 2 60- ` ~ ~ 0,55 ~ ~ - ~o ~ ' yo 210 230 250 270 210 230 ?50 T!f zn 210 23D 2~ 0 ~~4 230 250 T If Fig. 1. Temperature stratification of the atmosphere in p-coordinates with different initial conditions. 1) initial, 2) disturbed by ascending currents, 3) disturbed by descending currents. Fig. 2. Temperature stratification of the atmosphere in z-coordinates with different initial conditions. Thus, the maximum variations of temperature and~geopotential, caused by ver- ~ tical currents, are observed in the case of a very warm mesosphere (d), but they are quite small and it can be asserted that the role of the temperature gradients in the process of downward propa~ation of disturbances is also small. From (III) (2, 3): t (y-~) u~ = CaYe~' -au~, 2 1 ~ (y-~) d'z = 1 g~ e' yz ~-(a w)=1. 1 55 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY . r Here _ = y In p-F d-1' Cd - a u? - ut0 ~1 ~ Po ' _ It was demonstrated in I2] that the solution of system (III)-(VI) has the following peculiarity: all the even terms in the expansion u(u2, u4,...) and odd t~rms ~ ( ~t~, tY3 , . . . ) and ~ ( cAl, ~ 3, . . . ) become equal to zero . Therefore from (IV) (2, 3) u2 =~3 = 0. In the case of a linear stratificat:ion it is impossible to obtain an ana- - lytical solution for u3 and ~4, but estimates show that the temperature gradients do not change the behavior of u3, obtained in [~2] (~u3~~~ul ~ and of the same sign), but decrease ~ u 31 by .'.0-50% with differ~nt initial con- . ditions, and this means that ~ changes insignificantly. Since the contrib- ution of the terms A 3u3 and 8~ S~4 is small (less than 15~) and . they can only intensify the ef~ect, it can be assumed that the vortex for the mnst - part is determined by the terms ul, ~p and ~2. Thus, the transformation of the circumpolar vortex described in [2] is correct also for a real tem- perature stratification. Now we will cite variations of the temperature profile computed in z-coor- dinates (Fig. 2).~It can be seen that in this case the role of the t-emper- ature gradients is quite large. For example, at z~70 km the variations - 4T~2.5 K-- b(1, 2) and [jT=8 K-- d(1, 2); at z=40 km ~ T='15 K a(1, 2) and Q T= 20 K-- b(1, 2). Moreover, from the figure (b, d) it can be seen that an increase in the temperature gradient in the meso- sphere causes a decrease in L~T in it, and as indicated by computations, with great gradients (a ~ 2 I:/km) the ascending currents cause a.cooling of the mesosphere and heating of the stratosphere, which agrees well with ex- perimental data [3, 4). . BIBLIOGRAPHY 1. Bekaryukov, V. I., Purganskiy, V. S., "Precise Solution of a One-Dimen- sional System of Equations in Hydrodynamics and Some Possibilities for its Use," TRUDY TsAO (Transactions of the Central Aerological Observa- tory), No 115, 1973. 2. Zadvernyuk, V. M., "Some Results of Computation of Changes in the Ther- modynamic Parameters of tTie Strato-Mesosphere Caused by Disturbance of the Thermosphere," TRUDY GGO (Transactions of the rfain Geophysical Observatory), No 429, 1979. 3. Zadvernyuk, V. M., Mikhnevich, V. V., "Some Problems in Solar-Atmo- spheric Relationships," METEOROLOGIYA I GIDROLOGIYA (Meteorolo~y and Hydrology), No 9, 1973. 56 FOR OFFYCIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY 4. Ivanov-Kholodnyy, G. S., et al., "Discussion of the Results of the Con- ducted Experiment," ISSLEDOVANIYE ATMOSFERY I IQNOSFERY V PERIOD POVYSH- ENNOY SOLNECHNOY AK'�IVNOSTI (Investigation of the Atmosphere and Iono- sphere During a Period of Increased Solar Activity), Leningrad, Gidro- meteoizdat, 1970. 57 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 rUtt UFFLCiAL USE ONLY UDC 551.557(47-15)(571.5) CHARACTERISTIC DIURNAL VARIATIONS OF WINDS IN THE UPp~R MESOPAUSE REGION OVER CENTRAL EUROPE AND EASTERN SIBERIA ' , Moscow METEOROLOGIYA I GIDROLOGIYA in Russian No 11, Nov 79 pp 50-54 [Article by Doctor R. Schminder, Doctor of Physical and Mathematical Sci- ences E. S. Kazimirovskiy, Doctor D. Kurschner and.Candidate of Physical and Mathematical Sciences V. D. Kokourov and V. F. Petrukhin, :~ollm Geo- physical Observatory, Leipzig University, and Siberian Institute of Terres- trial Magnetism, Ionosphere and Radio Wave Propagation, submitted for pub- lication 26 February 1979] Abstract: On the basis of an analysis of simul- taneous wind observations in the upper mesopause regiQn over Central Europe and Eastern Siberia in the winter of 1977/1978 and in the spring of 1978 it is demonstrated that together with the coincid- ing nature of the diurnal variation there are differ- ences indicating the existence of ma~or regional structures in circulation and in the systems of tidal winds. A synoptic a~alysis requires an ade- quately dense network of observation stations. [Text] During the winter of 1974/1975 joint work was undertaken by the Kollm Geophysical Observatory at Leipzig University imeni K. Marks and the Badary Observatory of the Siberian Institute of Terrestrial Magnetism, Ionosphere and Radio Wave Propagation Siberian Department USSR,Academy of Sciences in the upper mesosphere region by the radiophysical :nethod by measuring iono- spheric drift with the spaced reception of signals of long-wave radio trans- mitters (D1 method). Continuing this program [1, 2, 4], it was possible to obtain results of simultaneous measurements for the periods 4-5 - 16-17 De- cember 1977 and 1-2 - 19-20 March 1978 (the measurement method was such that the observations were made during the nighttime hours local longitude time between sunset and sunrise). In addition to the results of ineasurements of long-wave drifts, the analysis was supplemented by wind measurement data for this same region obtained by the radiometeor method (D2) by specialists of the Kuhlungsborn Observatory of the Central Institute of Solar-Terrestrial Physics Academy of Sciences German Democratic Republic (Heinrich Hertz 58 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY ~ Institute). The characteristics of the measurement path~ are given in Table 1. � Table 1 Wind Measurement Paths in the Upper Mesopause Region at the Observatories in Kollm and Ku}ilun~sborn (German Democratic Regublic) and Badary (Soviet Union) - Iis~tcpiirenbxa~ PaccroA- MeTOA ~Iacrora 06o3xa- KoopAH}~ar~t rovKx Tpacca I Hxe, x.s ~ 3I ~ KexNe ~ orpaxtexaa 6 - 7 llc.i ~ex,~op~-konnM 170 D~ 185 xzy K 185 5.2� c, m, I3� s, A, $ 3ipKyres-&aAap~t 150 D~ 200 2 fi 200 b2 13 14 9 Bapmasa-KonnM 460 D~ 227 K 227 52 17 10 hY�~yHrc6opH-KonnM 500 D~ 245 K 245 53 12 11 paaHOHeTeopHme us~epe- xuA eeTpa 150 Dz 32,5 M~y K-PAM 54 ]5 5 56 12 KEY: 1. Measurement path 15. MHz 2. Distance, km 3. Method K = Kuhlungsborn 4. Frequency 6 = Badary 5. Notation PAM =.Radiometeor observations 6. Coordinates of reflection point 7. Zellendorf-Kollm . 3. Irkutsk-Badary - 9. Warsaw-Kollm 10. Kuhlungsbo rn-Kollm 11. Radiometeor wind measurements - 12. KHz > 13. N 14. E The selected interval of synchronous measurements did not make it possible to investigate disruptions of circulation associated with stratospheric heating, which according to measurements at Kollm reached the upper meso- pause region over Central Europe only on 1-2 January 1978. And here the March measurement period coincided with the onset of the.spring restructur- � ing of circulation, which despite the considerable distance between Kollm and Badary (about 5,000 km) began virtually simultaneously on 2-3 March 1978 and was manifested as a weakening of the westerly wind, a frequent change in the directions of zonal movement and an intensification of the easterly � wind. With respect to the semidiurnal tide there was a decrease in amplitude, a strong phase instability, and at the beginnin~ of the third 10-day period in March a rapid change in phase characteristic for the transition from a winter to a summer type of circulation [5j. 59 FOR OFFICT_AL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY � . Ol /0211! 1978 Od / 09d~ 9977 ~ - v H~cEelr m~se 1 ~ ~ I ir,~ 1 V"', r. I ~ l } wr ~ 40 - f Ij i.~ ,o, 1 ~ r 1 ,':ti a , t I B ~1: i 20 ~ _~~'~~~f ^ r' ~ r' ~b 1 l . - t ~ ~ r ..i�� l ; ~ y... ~ ~ ; ^ ~ ? W~- ~~'�!r's C 1 1 . ~ ' a ~ . . ~�l ~-'i-- ~ ~ A ~ Pj r 1 -1U ~ ~ ~ ~ + ~ r. -40 ~d ~ 3 I ~ - ti W ~ J? 40 � CN 1+.. 1! C~l � - ~ 6~ t ~'O,'~RI /ry !iM~~''v1 ~ 2~ ~ ' ; ,:u - Pl~~r ;~~e,~ ~ ~4~,�r,~;v~ d r~ =i _ ~ ~I�+~~, i I, I ; ~ f O ; � w ` - ~ VI , ' ~ , ' -~:l ` '�J ~ r' J ~ ~ ~ ~~r: i ~ - . ~ ~ , t ~ �'ti�~~~ ~ ~ - 9 - ~%o j y ~ - - ~o S ~ ~ ----~r s _ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~ ~6 18 ~ 20 ?2 24 02 09 OB � OB ~7 19 21 23 01 OJ OS 07 Oy MecmNOe BpeMA Local Time - Fig. 1. Diurnal variations of wind velocity vector in upner mesopause region o~~er Central Europe and Eastern Siberia at middle (a, b) and at end (c, d) of winter. a, b) zonal component, b, d) meridional compoiient. 1) K 185; 2) B 200; 3) K 227; 4) K 245; 5) K-PAM [K = Kollm; ~j = Badary; PAM = radio- meteor wind measurementsJ An analysis of wind measurement data on all paths in December 1977 and in March 1978 indicated a coincidence of the nature of the diurnal variations for all stations having an approximately identical geographical latitude, ~ regardless of longitude. However, in details it is possible to see an ex- cellent correlation and very great difference in both the prevailing and in the tidal o~inds. With respect to the short-per~od variations, which are interpreted as the effects of internal gravitational waves, here it is im- possible to detect any correlation, since the coherence scale for these variations is less than 200 km [6). - Figure 1 presents data from simultaneous measurements for one night in Decem- ber and one night in March, when the agreement of the results was rather good. The values are given in local longitude time. It must be remembered 60 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY that the difference in time between Central Europe and Eastern Siberia is . 7 hours. Figure la,h, which shows the wind values for all five measure- - ment paths (including measurements by the D2 method at Kuhlungsborn), gives some idea concerning the real differences in the results between the paths w3th reflection points situated very close together. For all the European paths they are all in a small volume with dimensions in longi- tude and latitude of the order of several degrees and with respect to al- titude not more than 10 km. It can be seen that the measured values occupy some interval and in general it would be very useful in many respects on the basis of ineasurements in Central Europe to ascertain the center of this "zone of values" and provide users not only the mean wind value, but also the width of this interval. Ztaelve years of experience with measurements at Kollm, where at first the measurements were made only on one path, en- ables us to note that in the interpretation of such measurements it is necessary to approach very carefully the evaluation of the accuracy in the measurements themselves and the representativeness of the results for one reflecti~n point relative to the mean chax�acteristic of the wind re- gime on a regional scale. " Figure l,c,d shows the diurnal variations on the last day of stable win- ter circulation for the paths K 185, 6 200 and K-PAM. The data from radiometeor measurements were furnished through the courtesy of the Kuhlungsborn ionospheric observatory (see Table 1). Figure lc,d clearly shows the excellent agreement of the results of ineas- urements by the Dl and D2 methods for Europe. The agreement is not always - so good, but already available observational data, making it possible to carry out their comparison, indicated that a high correlation exists with adherence to the following conditions. First, when the measurement condi- tions ensure statistical reliability of the result; second, when a stable circulatiori ensures a spatial uniformity of the wind regime; third, when there are sufficiently great amplitudes of the tidal winds, which facil- itates the relative analysis. If at least one of these conditions is vio- - lated, the correlation rapidly decreases. This means that if we observe differences in the measurement results by both methods, they can be at- tributed either to a certain uncertainty in the method.itself (for ex- ample, due to the lack'of data on the precise reflection altitude) or to some still not finally clarified pec~iliarities of these methods. Table 2 gives a sample of results for other days. It also gives some idea concerning both the coincidences and discrepancies in measurements at two points. After comparing the results, we must conclude that in addition to the coinciding nature of the diurnal variation for both the prevailing and the tidal wind in the upper mesopause region over Central Europe and East- ern Siberia, t~ere are differences indicating the presence of regional structures in circulation and in the systems of tidal winds. This means 61 FOR OFFICIAL USE ONLY , APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY Table 2 Result:~ of Analysis of Wind Measurements in Upper Mesopause Region at the - Observatories Kollm (GDR) and Badary (USSR) on Paths With Close Frequenc- ies (185 and 200 K~Iz) and Approximately Equal Distances from Transmitters 3oxaabx~ii xo~noxexz MepHltFCOeanbxd~ 06oaxa- 3 KowrtoHei+r q HOqb qeaxe 1 2 vo , v: I 1's Yo I vs I T~ 5 05/~ ,~exa6pA ] K 185 -05 38 20.45 -03 17 18.00 1977 r. 8 S 200 + 18 11 20.00 + 09 33 15.45 . OS/09 K 185 + 09 29 21.45 i-00 12 19.i5 S 200 +28 31 2L30 +04 1? �19.15 12/13 K 185 +28 33 21.15 +06 35 17.45 B 200 +32 12 17.30 -17 07 14.45 14/15 K 185 +24 32 20.00 +Ol 23 18.15 � S 200 + 17 20 17.00 -18 21 ! 3.15 . - Ol/02 uapTa K 185 +21 42 21.15 -11 23 16.30 6 1978 r. & 200 +07 35 21.15 -13 27 18.00 . : ~ O6/O7 K 185 +03 30 20.30 -04 18 15.45 B 200 -17 24 22.15 -15 ~ 34 17.15 C~EY : 1. Night 5. December 2. Notation 6. March 3. Zonal component 7. Kolmm - 4. Meridional component 8. Badary Note: V~ prevailing wind, positive to the north or to the east, V2 amplitude of the semidiurnal tidal wind, T2 (local time) phase of sersi- diurnal tidal wind, moment of maximum wind, directed to the north or east. the detailed coincidence of the curves of temporal wind variations can be either random (Fig. la,b) or associated with a marked seasonal change in the wind fields (Fig. lc,d). A study of changes in the prevailing wind from day to day could afford a possibility for investigating planetary _ 62 FOR OFFICIAL USE ONLY . , APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000240040057-3 I FOR OFFICIAL USE ONLY waves, but this would requira intermediate measurement points [2]. In the future, evidently, it is necessary to make a detailed study of the dis- covered regional large-scale structures separately in each regi~n and only then~ in the second stage, combine the results for constructing a general model of circulation. Real synoptic investigations require the organiza- tion of additional measurement points; the density of the network which exists at the present time for these purposes ts inadequate. BIBLICJGRAPHY _ 1. Grayziger, K. M., Kazimirovskiy, E. S., Kokourov, V. D., Petrukhin, " V. F., Shminder, R., Stirenher, K., "First Results of a Comparison of ~ Measurements of Ionospheric Drifts in the Long-Wave Range Over East- ern Siberia and Central Europe,".ISSLEDOVANIYA PO GEOMAGNETIZMU, AERO- NOMII I FIZIKE SOLNTSA (Investigations of Geomagnetism, Aeronomy and Solar Physics), Vol 38, 1976. 2. Shminder, R., Kyurshner, D., Kazimirovskiy, E. S., Kokourov, V. D., "Mean Monthly Diurnal Variations of Zonal and Meridional Components - of the tidind in the Mesopause Region Over Central Europe and Eastern Siberia in the Winter of 1976 and the Spring of 1977," GEOMAGNETIZM - I AERONOMIYA Geoma ~ ( gnetism and Aeronomy), No 18, 1978. 3. Mu11er~ I~. G., Nelson, L., "A Travelling Quasi-2-Day Wave in the Meteor Region," J. ATMOS. TERR. PHYS., Vol 40, 1978. 4. Schminder, R., Kazimirovskiy, E. S., Kurschner, I)., Petruchiny W. F., "Ein Vergleich der Ergebnisse von Ionospharendriftmessungen im Lang- wellenbereich uber Mitteleuropa und Ostsibirien in Winter 1975/76," _ Z. MET., b. 28, 1978. _ 5. Schminder, R., Kurschner, D., "On the Behavior of Wind Systems in the Upper Mesopause Region in Winter and During the Transition from Winter to Summer Conditions," J. ATMOS. TERR. PHYS., Vol 40, 1978. 6. Vidal-Madjar, D., "Gravity Wave Detection in the Lower Thermosphere ~dith the French Incoherent Scatter Facility," J. ATMOS. TERR. PHYS., Vol 40, 1978, - 63 FOR OFFICIAI, USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE OIJLY UDC 551.(521.32:524.7) EVALUATION OF ERRORS IN COMPUT"ING EFFECTIVE RADIATION - Moscow METEOROLOGIYA I GIDROLOGIYA in Russian tJo 11, Nov 79 pp 55-61 [Article by Doctor of Geographical Sciences A. I. Budagovskiy and L. Ya. Dzhogan, Institute of Water Yrohlems, submitted for publication 20 March 1979] Abstract: The paper describes a method for the indirect evaluation of errors in comput- ing effective radiation. It is based on a - comparison of the differences between the temperature of the soil surface and the air, determined by three independent methods. The merhod is applied to observational data from Takiatash actinometric station. The authors give a brief analysis of the results. [Text] Regular observations of the radiation regime have been made in the relatively dense network of actinometric stations in the Soviet Union - since the mid-1950's. The accumulated dat4 are necessary for the use of modern methods for climatological and hydrological computations, in par- ticular, for computations of evaporation and irrigation norms. However, the use of the mentioned materials involves considerable difficulties. They are caused, in particular, by the fact that the measured values of the radiation balance to a considerable degree are dependent on the albedo and temperature of the underlying surface at the measurement site. The = difficulties mentioned above are usually overcome by the use of the total - radiation values measured at actinometric stations and accumulated data on the characteristics of albedo for different underlying surfaces; in - this case effective radiation is determined by computations. ' It is evident that errors in such computations are in need of evaluation. It is desirable that this be done on the basis of quite extensive mater- ial so as to form somp idea concerning the statistical stability of the results. This problem can be solved indirectly. It involves essentially the following. ' On the basis of data from observations made at actinometric stations it is Possible to determine the effective radiation value. It is equal to the difference between the measured values of absorbed short-wave radiation 64 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY and the radiation balance. The latter is dependent, in particular, on the temperature of the underlying surface in the measurement sector beneath Lhe balancemeter. The effective radiation value obtained in this way wi11 t?~nceforth be ca:Lled the "measured" radiation and in case of necea- sity it can be assumed conditionally that it does not contain errors. It is convenient to represent the effective radiation value found in this way in the form of the sum of two terms. The first of these is equal to the effective radiation for a case when the temperature of the underlying surface beneath the balancemeter Tu is equal to the air temperature meas- ured in the meteorological booth. The second term takes into account the influence exerted on effective radiation by the difference between the _ iaentioned temperature values. III=1~+4a S (273--~ T2) 3 (T~-T:) . ~1) [ i7 = u ] Here Iu is the effec.tive radiation found by the method indicated above on the basis of the measurement results, I* is the effective radiation when Tu = T2. W* will denote the computed effective radiation by Icomp and assuming that _ I� Icomp~ on the basis of (1) we write etp=1~-ly =46S(273+Ts)3(T~~-rz)� (2) - [ p = comp; 1T = u] Since the computed values of effective radiation in the most general case can contain both random and systematic errors, they will automatically enter into Q I~omp . These errors can be evaluated by computing the Q I values by some other two independent methods. The first of these methods involves the use of ineasurements of soil surface temperature. At meteorological stations they are carried out on bare sec- tors of the soil. Therefore, their use for the purposes mentioned above in principle is possible only when the measurements of the radiation bal- ance are carried out over a sector with bare soil. On the basis of these " measurements the Q Iu value can be computed using the expression ~/i,=4vS(273+Tz)3(T�-T2). _ (3) The second possible method for indirect but independent determination of the difference Tu - T2 involves use of the formula P=a p ~pD(Tn-T2), ~4) 65 FOR OFFICIAL US~ ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY _ in which P is turbulent heat exchange, oC is a conversion factor, is air density, cp is the specific heat capacity of air at constant pressure, D is a coefficient having the dimensionality of velocity. It characterizes the turbulent conductivity of the air layer bet~een the soil surface and any fixed altitude (in this case 2 m). Henceforth we will call it the ex- change coefficient. The turbulent heat exchange value enters into the heat balance equation R=P+B~-LE, (5) in which R is the radiation balance, B is heat exchange in the soil, L is the latent heat of evaporation, E is the turbulent flux of moisture or evaporation. The day-to-day changes in heat exchange in the soil are usually small and rarely exceed 10% of the radiation balance. Therefore, in an approximate estimate of turbulent heat exchange this heat balance component can be omitted or estimated approximately using data from the literature [3, 5]. - Under extremely dry conditions, especially in the arid zone during the - summer months, evaporation from the surface of the bare soil can be assum- ~d equal to the precipitation (H). Accordingly, assuming E~ H and limiting ourselves to the above-mentioned anproximate estimate of heat exchange in the soil, it is possible to obtain turbulent exchange values which are employable for practical purposes. In their subsequent use for computing Che differences Tu - T2, and then L~1 it is necessary to know the values of the exchange coefficient. The literature gives estimates of this coef- ficient used in climatological computations. However, a high percentage of these coefficients were obtained on the basis of observations over the plant cover. Therefore, we used additional observational data given in [4]. They include information on all the components of the heat balance and meteorological elements, including on the temperature of the soil surface, measured using scattered temp2rature sensors [6], giving the spa-~ tially averaged temperature and having a small radiation error. The ob- ' servations were made in the Fergana valley in cotton.fields in 1956 be- _ ginning on 7 May. In the study observational data were used from early in - May through the first 10-day period of 3une inclusive, when the cotton " plants are extremely small and do not Exert a significant influence on for- mation of soil surface temperature. Computations of the exchange coefficient D were made using expression (4) a~d on their basis it was possible to construct the dependence D= f(U), where U is wind velocity at a height of 2 m(see Fig. 1). We note in pass- ing that an attempt to construct a curve of the dependence of the exchange coefficient on ~ind velocity, taking into account the influence of temper- ature stratification, did not give significant advantages. ~ 66 FOR OFFICIAL USE OPdLY . 1 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY Dca/ceK cm/sec - 1,3 0 0 0 ~ 0 0 1,0 0 0 0 o Q 0 0 J 4s o o ~ i i 0 Z U~r mf sec Fig. 1. Curve of dependence of exchange coefficient on wind velocity. Along y-axis: D cm/sec; along x-axis: U m/sec ~I? dl, ~ d J~ . 0 O ~ G~Q ~ uC0 O C ~O ~ ~A G~~ ~ ~ ^ ~ ~~O O 0 0 CA O ~ J c` `~j9 C O Op~ eS ~ O . C~ yCJ ~ . O O~ g C, ~ O O ~ O J 2 O ~ 00 O C � O n � O p � comp a u i J 2 !t ~.D ~ Z 9 C _ din r'ig. 2. Curves of correlation between Q Ii values computed by different methods. - In order to estimate the error in the dependence D= f(U) on the basis of data on the diurnal swns of turbulent heat exchange,by means of inverse computations we computed the temperature difference Tu - T2 and compared this with the corresponding measured values. Then after scaling the com- puted and measured Tu - T2 values into radiation units we determined the standard deviation characterizing the error in determining Q I, asso- ciated with the use oi the mentioned dependence. According to computation data, c7 = 0.46 Cal/(cm2�month). Taking into account what has been said above, on the basis of (4) and (5) we can write ~ ~!o = 4 ~S (273 T,)3 R-B-LH �pcoU ~6) Since the errors in computing effective radiation, acr_ording to the com- ment made above, are completely included in ~ Icomp. their evaluation in- volves the paired comparison of the determined d Ii values. It is possible to obtain some idea concerning the systematic differences of the ~ Ii values computed by the three above-mentioned metliods by using the expressions . 67 - - FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY ~lP=Q,~jo, ~11P = a_ ~ /n, ~Ip =Q~~/~~ \ ~~/P `'~lP ~1 = ~�S jD r - c, Ip ' n~ - E~1 . IP = comp; Tf = u] Graphs of the dependence (7) are shown in Fig. 2. It is evident that the measure of the relative value of the systematic dif- ~ ferences between the L~Ii values, determined by the two independent meth- ods, is the deviation of the corresponding ai values from 1. The random deviations between the paired comparable L1I values, computed by different methods, can be characterized by the values of the corres- ponding dispersipns ~2= -a ~8~ ~ ( ~ ~ o 02 = (t~ Ip - aZ /n)'-~ Q3 = (~fp- R3~Ia~=. Here the horizontal line is the averaging symbol. Since the L1I~omp ,[~ID and p Iu values were computed by independent meth- - ods, the errors in these computations are also independent. In actuality, the random errors in romputing Q I~o~P are related primarily to the devia- tion of the real vertical distribution of. air temperature and humidity, the deviation of clouds from their, "typical" values adopted in validation of the parameters in the computation formulas,.and also their approximate character. The error in computing ~ID is determined by the c~rrespond- ing error in ~he dependence D= f(U), incomplete adherence to the hydro- dynamic conditions for its use, and also the conditions necessary for ad- equately reliable determination of the monthly sums of turbulent heat ex- - change. Finally, the random error in computing Q Iu is related to the cor- responding errors in measuring the soil surface temperature. Accordingly, the va~ues of the dispersions of each of the two paired comparable values - must be the sum of the two corresponding dispersions, to wit 68 ~ - FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY y"~o Tab 1 qpe~~ae Q~ Q_ ~y Qi ~1 Q3 Qp I QD Jn a~ = aP + Jp~ Nu ber o te in s~r s a - - ie p cbmp n~ u 29 1.Oi 1.~~0 ' 1.49 1,10 1.26 0,9~ 0,84 O,G3 U,7-1 0! -'P + on+ - ~ n,`~~'~ I.~.; ! 1,43 1,03 t,33 0.88 0.36 0,54 U,77 ~9~ b8 I,i1'' I.:,_' ; 1,-~9 1,06 1,30 0,91 U,:,~ 0,58 O.iG �3-~Q+�Q' (T~= u; p= rad] Here O"rad is used to denote the dispersion Q Ip, Ct D denotes the dispersion Q ID and O~'u denotes the dispersion L~Iu. 1'I~e three equations of the very simple system (9) contain three unknown parameters O"rad, v"D and ~u, which can be determined using the already computed values cYi, O'2 and O'3. In order to apply the described method to an evaluation of the errors in computing effective radiation we used observational data for the actino- metric station Takhiatash (lower course of the Amudar'ya River), in whose description, given in [1] and in handbooks on climate, there is a clear indication that the observations were made in a sector with a bare soil surface, that is, there was satisfaction of the requirements necessary for computing the L~ID ~and ~ Iu values. In addition, this region of Central t\sia is characterized by very little precipitation (the annual precipita- tion norm is 98 mm). Therefore, here, more frequently than in other re- gions, for monthly time intervals th~re is satisfaction of the condition E~ H necessary for computing Q ID. In the computations we used the Brent formula, the values of whose coeffic- ients were refined by M. Ye. Berlyand [2] on the basis of theoretical com- putations. A linear dependence was used in taking into account the influ- ence of cloud cover. The value of the parameter c entering into it was also taken from M. Ye. Berlyand. The values ~ I~o~P, L1Iu and L~1D were determined using expressions (2), (3) and (6). In the computations we used observational data obtained during tt~e warm half-year (April-September) during 1955-1959 and 1961-1968, given in j2], in actinometric handbooks and in handbooks on USSR climate. Additional data on temperature of the soil surface, temperature and air humidity, wind velocity, total cloud cover and precipitation were taken from handbooks on climate and archival data. In computati,ons we excluded cases~'when satisfac- tion of the condition E~II caused well-founded doubt or when there were gaps in the observations of one of the elements. The results of observa- tions for a total of 58 months were used in the computations. In order to form some idea concerning the statistical stability of the computed charac- teristics the series was broken down into two equal parts. The results of paired comparison of the computed values are given in Fig. 2. - The computed ai and ~i values (in Cal/(cm2�month)) are given in Table 1. We must note that the differences between the ai and Cl i values obtained for the.first and second halves of the investigated geriod and also for the period as a whole do not exceed the limits of the admissible error 69 ~ FOR OFFICIAL USE ONLY ~ APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOK OFFICIAL USE ONLY in estimating the sample values of the corresponding parameters, that is, the ohtained ai and o'~ values have an adequate statistical stability. The second important conclusion following from an examination of the computed values is that the values of the coefficient al (the systematic discrepan- cies between the compared L~1~o~ and ~ID values) do not differ aignif- ~ icantly from unity not only for the investigated period as a whole, but al- so for it~ parts, that is, tfie two values virtually ~oincide. On the other hand, the values of tfie coefficients a2 and a3 differ substantially from unity, evidence of considerable systematic differences between the comput- ed d Icomp andLlID values, on the one hand, and L~1u on the other. The rea- son for this, as is clearly noted, is the systematic error in measuring temperature of the soil surface, data for which are considerably too lowo The latter for tfie conditions in Central Asia during the summer months is entirely natural since the surface of even dry soil has a considerably less- er reflectivity in comparison with the reflectivity of the thermometer res- ervoir. As a result the mean value of the characteristic temperature of the latter is considerably lielow the temperature of the dry soil surface. Thus, the results, even with the most cautious approach, give basis for drawing the conclusion that there are no significant systematic differences be- _ tween the effective radiation parameters obtained by computations and their measured values (under the conditton Tu = T2). _ In evaluating tfie determined values characteri~ing tlie rand4m error ~.I in computations attention must be given tQ the fact that for all three co~ putation methods they differ little from one another and in general are not great. A mare complete idea concerning the magnitude of the computation error ~ Icomp can be obtained by comparing it with the computed values cf effective radiation and the radiation balance. The mean value of the computed values of the effective radiation for April- September during the aFiove-mentioned time interval (1955-1959, 1961-1968) is 3.2 Cal/(cm2�month). Accordingly, the relative probable error in such computations is 18%, and the error with an 80% guaranteed probability - 34%. The mean value of the radiatioi: ~~lance, obtained on the basis of the re- sults of ineasurements,of absorbed short-wave radiation and the computed values of effective radiation, is 10.8 Cal/(cm2�month). In such a case the probable error and the error with 80% guaranteed probability, related to the use of the computed values of effective radiation fox the warm half- year, on the average will be 5.3 and 10.1% respectively, which does not ex- ceed the errors in measuring the monthly sums of the radiation balance, - which, in accordance with [7], are estimated at 10%. Next we note that the Oro value, cited atiove in the text, and ~D, cited in Table 1, determining the error in computing L11D, differ relatively little from one another. However, insignificant differences between them are en- tirely to be expected because cTO was ohtained on the basis of use of heat balance ohservations and d D on the basis of mass materials. It is more ~ 70 FOR OFFICIAi, USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY important that ~D is. appreciahly less than ~~o~ and u. This conclu- sion is important fo'r evaluating the error associated with the introduc- tion of a correction to the radiation halance in computations of evapora- - tion from the surface of the moistened soil. ~ l~in~i.lly, it should he noted that despite the wel..l.-known incorrectnese of standard observations of temperature of the soil surface carried out in the _ network of inetEOrological stations the 0'u value, dependent in the last analysis on the random errors of these observations, is relatively small. 1t ~s even somewhat less than the random error in computations of effec- ~ tive radiation characterized ~v the ~comp value. The described method for estimating the errors in computations of effective _ - radiation can also be used for other points if there is satisfaction of the above-mentioned conditions necessary for computing the Q Iu and ~ ID ~alues. BIBLIOGRAPHY 1. Barashkova, Ye. P., Gayevskiy, V. L., D'yachenko, L. N., Luchina, K. M., - Pivovarova~ Z. I., RADIATSIONNYY REZHIM TERRITOFII SSSR (Radiation Re- gime of the USSR), Leningrad, Gidrometeoizdat, 1961. 2. Berlyand, M. Ye., Berlyand, T. G., "Determination of the Earth's Effec- tive Radiation With Allowance for the Influence of Cloud Cover," IZV. AN SSSR, SER. GEOFIZ. (News of. the USSR Academy of Sciences, Geophys- ical Series), No 1, 1952. 3. Budagovskiy, A. I., ISPARENIYE POCHVENNOY VLAGI (Evaporation of Soil Moisture), Moscow, Nauka, 1964. 4. Budagovskiy, A. I., Mii;~yeva, Ye. N., "Results of Investigation of Evap- oration from Irrigated Fields in Central Asia," TRUDY GGI (Transactions of the State Hydrological Institute), No 151, 1968. 5. Budyko, M. I., KLIMAT I ZHIZN' (Climate and Life), Leningrad, Gidro- meteoizdat, 1971. 6. Kaganov, M. A., Ch~idnovskiy, A. F., "Instruments for Measuring Tempera- ture of the Soil Surface," SBORNIK TRUDOV PO AGRONOMICHESKOY FIZIKE (Collection of Yapers on Agronomic Physics), No 5, Moscow-Leningrad, 1952. 7. Lebedeva, K. D., Sivkov, S. I., "Accuracy in Measuring the Radiation Balance by Thermoelectric Balance Meters," TRUDY GGO (Transactions of the Main Geophysical Observatory), No 129, 1962. 71 FOR OFFICIAL USE ONLY - APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 ' FOR OFFICIAL USE ONLY UDC 551.593.5(263) OPTICAL CHARACTERISTICS OF THE ATMOSPHERE IN THE TROPICAL ZONE OF THE ATLANTIC OCEAN Moscow METEOROLOGIYA I GIDROLOGIYA in Russian No 11, Nov 79 pp 62-69 [Article by V. N. Adnashkin, L. K. Veselova and Candidates of Physical and Mathematical Sciences 0. D. Barteneva, A. G. Laktionov and N. I. Nikitin- skaya, Main Geophysical Observatory, Leningrad State University and Insti- tute of Applied Geophysics, submitted for publication 22 May 1979] Abstract: This paper presents the results of investigations of the spatial-temporal variability of the aerosol-optical character- istics of the near-water layer and atmospher- ic layer of the tropical zone in the Atlantic Ocean according to data from TROPEKS-72 and GATE-74. Also considered is the latitude var- iation of integral and spectral values in the region 0.35-1.OOfcm of atmospheric trans- parency, the selectivity index of aerosol at- tenuation of solar radiation and moisture content - of the atmospheric layer, and also a number of , characteristics of the near-water layer: meteor- ological ranRe of visibility, concentrations of large (d> 0.63�m) and giant (d > 10 �m) aerosol particles. It is shown that the principal fac- tor responsible for the inconstancy of optical weather ir~ the tropical region of the North At- lantic is the transport of dust from the des- erts of the African continent. [Text] Experimental investigations of the parameters of microstructure and optical characteristics~of aerosol in the near-water layer and atmo- spheric layer in the tropical latitudes of the Atlantic Ocean indicate their considerable temporal variability and spatial nonuniformity [1, 3, 4, 6, 7, 12-14, 16]. It is shown in this study that the principal factor responsible for the inconstancy of optical weather in the tropical region of the North Atlantic is the transport of great quantities of dust into 72 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY the ocean from the deserts of the African continent. During the period of GATE-74, from 27 June through 10 September, on the basis of data from Lhe geostationary satellite GMS-1, "dust cyclones" nine times crossed the Atlantic Ocean from east to west in the latitude zone from 10 to 20�N [2J. The eastern part of this region of the Atlantic Ocean, the so-called "sea of gloom," in which the transport of dust from the African continent by the NE Trades is systematically observed [10], - can also be seen clearly in the averaged latitude variation of optical and aerosol characteristics of the atmosphere from 50�N to the equator, represented in Fig. 1 on the basis of data from TROPEKS-72 [1, 6]. Spe- cifically in this region of the ocean all the aerosol-optical charactexL istics considered below attain their extremal values. For example, in the "sea of gloom" the following are noted: the lowest integral transparency PZ [9] and the highest value of the aerosol component of the optical lay- - er of the atmosphere 'G j~ for 1.OO�m; the maximum N concentrations of large (d~0.63~m) and giant (d>l0�m) particles and the minimtun val- ues of rhe meteorological range of visibility S in the near-water layer of the atmosphere. The n parameter, characterizing the degree of selec- _ tivity of aerosol attenuation of solar radiation in the well-known Ang- strom farmula 'Cj~= /3~-n has a value on the order of 0.2, that is, the spectral variation of the aerosol component of optical thickness of the atuasphere in the spectral range 0.35-1.00 ~ m is extremely close to neutral. Figure 1 shows that the air masses in the temperate latitudes are charac- terized by appreciably higher values of the transparency characteristics - and the selectivity of aerosol attenuation, and also a relatively small moisture content W of the atmospheric layer, which with the degree of ad- vance into the tropics increases and attains maximum values in the ICZ intertropical convergence zone [8]. _ The equatorial region from 10�N to the equator, within which the ICZ is situated, is also subject to the influence of the transport of contin- - ental dust from the central and southwestern regions of Africa. However, the effect of the "cloud filter" in the ICZ region and the arrival of pure oceanic air masses in the lower layers of the atmosnhere with south- erly and Anti-Trades westerly winds favor this picture: with movement from 10�N toward the equator there is a tendency to an increase in trans- parency of the entire layer of the atmosphere and especially its near- water layer [6]. Thus, Fig. 1 shows how nonuniform this region of the Atlantic Oceara is with respect to its optical properties. In Fig. 2, where the days of observations have ~een plotted along the x- axis, we show the temporal variability of the parameters of microstruc- ture and optical characteristics of aerosol considered above for the period of the three phases of GATE-74 from 28 June throuQh 19 September durin� which the "Passat" scientific research weather ship was situated 73 - FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY at a point on the equator 10�W at a distance of about 500 km from Africa. WcH S . ' J� > ~ r'_.~. x `h ~B ~ ~6' s n , , ~ O,G~ * ~ ta ~2 x R ~r~ % x ~`t ' x A 3 0,0 " Q `-x'x # x . I QS ~ ~ r SKM 6 x x f0 ~ ~~S 20 a \ x Q 'A~ , Ncn' a-J i,~~~~ s 1~ x~~ - ~oZ 1 ~ '�~l'\.,~7 10' SO y0 JO ?0 90 G y GrrG N Fig. l. Latitude variation of aerosol-optical characteristics of atmospher- ic layer and near-water layer over the Atlantic Ocean (Northern Hemisphere). Figure 2 shows that the amplitude of variation of aerosol-optica~. charac- teristics of the atmosphere during this period is extremely significant. The integral coefficient of atmospheric transparency P2 varied from 0.59 to 0.75. The concentration of large (d ~ 0.63 ~tm) aerosol particles in the near-water layer varied from 2 to 30 cm 3 and the concentration of par- ticles d~ 2 N.m was from C~.1 to 1.3 cm 3. The meteorological range of vis- ibility S varied in the range 10-100 km. The values of the aerosol component of the optical layer of the atmosphere '~~j, for a= 1.O�.m were in the ran~e 0.10-0.40. A fact of importance is that the nature of the spectral variation of aerogol attPnuation varied _ from close to neutral n= 0.1 to extremely selective n= 1.3. The water vapor content in the atmospheric layer was 2.5-4.5 cm. The concentration of giant particles with d> ZO ~.m during the entire period was very low and did not exceed 10-3 cm 3. 74 - - FOR OFFICIAL USE dNLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200040057-3 FOR OFFICIAL USE ONLY rrc~~ r pa3u d q~a~a ID~a3a r ~1~ ; --;-T r-r-r~~~.,,~T_ -`1a_r~n-rt-r-~--rT-r~T-- -~-T--, . T _ 4 ~ . , _ r _ F~ , . , TTrr~-i r ( . aJ 1~~~ ~ i ~ ~ ~ TTl- .Tr. ~ ~ ~ !2 I~ II~~ I~~~ ' i ~ ~ i i' I ~ ~ . 0;1 ' ae l T; ; ~ - , ! ;-~T~ ~1 r - ~ - - ~ t - - ~ _ i ~ Qo - - 0,7 ~ , i rr1 ~I ~ i . sr,r~ 6 - T--- ~ ~ ~ , ~ i- ^ ~ :IG1 ~ ~ I i I ~1I`~~_`-J . ' ` l.: ~'T~ 1 ~ i I 90' . _ ~ ~ 2~ ~ _ ~ . 10 _ - _ _ , � - - ~ ~ S _ . - - ' ' - - ' - ~ ~,_..-.-...:y_:-: L J . � L.1.,. _ ! ~ ' 18 .JOI 3 S".'' 15 ll T? ~3 J1' 3 S 7 9 1t 13 95 97~9 d1 ; 3 S 7 9 11 13 73 fl 19 f~yn 2 3 aery'rm 4 CeHma6pe 5 P Fig. 2. Temporal variatic,n of aerosol-optical characteristics of atmosphere . - during GATE-74 period ("Passat" scientific research weather ship, 0� lati- tude, 10�W). 1) moisture content of atmospheric layer; 2) aerosol compon- ent af atmospheric layer for 1.0 �m; 3) index of selectivity of aero- - sol attenuation of solar radiation; 4) integral coefficient of atmospheric attenuation; S) meteorological range of visibility in near-water layer; 6) concentration of giant particles (d> 2~m); 7) concentration of large par- ticles (d~0.63~m). KEY: 1. Phase ~ 2. June 3. July 4. Augus t 5. September We should note the presence of a hi;ii correlation between the values char- acterizing atmospheric turbidity P2 and 'G;~, and also between the concen- tration of large (d> 0. G3 �m) particles and the meteorological range of visibility; in ttie latter case the correlation coefficient was 0.94. In the considered region of the ocean the concentration of large particles in the near water layer is essentially dependent on the direction of the flows in the free atmosphere. The increase in the concentration of aerosol particles is associated, in particular, with advection of northeasterly ~ flows in the layer 2-3 km, which carry dust-laden air from Africa, wher~as 75 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-00850R040240040057-3 FOR OFFICIAL USE ONLY its decrease i~ associated with an intensification of monsoonal cirrula- tion in the layer 1-1.5 km [7]. These results agree with the data obtain- ed in TROPEKS-i2, on the basis of which the conclusion was drawn that the basic mass of large aeresol particles is of continental origin [6]. Table 1 Aerosol-Optical Characteristics of Atmosphere During Du::t Intrusion Periods . ~ _ _ ~ = U Q Q MeCio ~ ~ ~ x x _ I~aTa P.: S Ax ~ a_ o. s., xaxepex~ta ~ ~ 3 7~~ ~ F~ 5 1 'Z ~.a pS~ o=~ ii ~ e c~ x c. t- r 6 HNCII RIlaccaT~, 15 e~ons 1974 r. 11 (~.61 17 0,9 0,2J 16 0,5 ~ 0� m., 10' s. Cpe~aee 3a 1 c~asy 12 0,71 43 0,2 ~~,1~ 3,1 0,;~ 29 xwax 1974 r. 11 u,~9 15 1,3 O,Y~ 9 ~).4 30 a 0,58 iG 0,7 U.dU 2:i 1,~~ 11 aerycra 1974 r. 13 U,:'iU 8 - - 30 0,7 13 s G,59 1,0 (~,23 16 U,6 14 a 6.:~9 0.9 U.1ti. 14 U,8 Cpe~xee aa II ~asy 14 0,63 4~~ 0,5 O,Iti 8,5 0,6 g H~iC cAr.aaeMqK 1 aeryrra 1972 r. O,:i8 12 0,2 u,3�1 KypvaTOe 13 9 18,8� c. r . is,s 3. k. io KEY: 1. Place of ineasurement 11....July 2~. Date ~2. tiean for phase I 3. �m 13. ..August ~ Concentration of large particles 14. Mean for phase II 5. Concentration of giant particles 6. "Passat" scientific research weaLher ship 7. 0� latitude, 10�W - 8. "Akademik Kurchatov" scientific research ship A comparison of curves 2-7 in Fig. 2 shows that the greatest~variations in aerosol-optical characteristics of the atmosphere occurred on 14-15 and 29-30 July and 9-14 August, when dus*-laden air masses arrived in this regi~n of the Atlantic, some of which, moving from the arid zones of North and Southwest Africa in a southwesterly direction, under def in- ite synoptic conditions reached the equator. A~oint analysis of TV and IR images for this region of the earth from the American meteorological satellite GMS-1 indicated, for example, that on 30 July 1974 the centers of the dust storms were situated in regions with the coordinates SP = 18� N, ~l= 8�W and 31�N, ~l= 0�W at a distance of aBout 2000 km or more from the measurement site j15]. 76 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200040057-3 APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-00850R040240040057-3 FOR OFFICIAL USE ONLY - . _ N CH'J 102 90 1 9d' ~d2 3 ~ - . To~ ~ _ ~ 2 , 4 ~ i0 S ~ � p~y 1 2 4 10 20rMHM N` - Fig. 3. Spectra of size of aernsol particles in near-water layer obtainel - in different regions of tropical zone of Atlantic Ocean. 1) "sea of gloom," - ~14�N, 19�W); 2) "sea of gloom" (18.8�N, 16.5�W), 1 August 1972; 3) Equa- toriaJ. Atlantic (0� latitude, 10�W), 30 July 1974; 4) 1ow-dust region in ocea >30 > 10 ~7 2-2 17 Hoeo�A.~e~:cauipouch 2 It 19i;t > >(0 >10 `2 J-~ ~ 18 MNxafinoei:a 3 111 19iti ~ >30 >10 10:1 26-77 19 ~epThouo l X 197; c:?aGoe2 ~10 AHeM 2~ 2-2 I no~Nee ~ 20 ~4so6irnbxo~ X 197G ro xce 2 Tana~ 7- ~7 i-7 21 Kpacnoraap~eNCKOC 2 `CI 19F~~ ~ To Hce - 6i 12 -!3 22 ~'Iaoli~:?b}ioe 2 ~CI I~~a cyxas2 0,5 4 ~ -0,2 77,0 E Srp e VI-VII 75 40 14 16 5 a -5,4 -77,0 ~ Cp. t Q$ I-III 75 5'J 1 16 6 s -5,2 I,0 1 Cp. tQ sa npeAwecxeyw- 31 9 16 u~a~ roA 7 asrycr ~-6,9 98,0 . ~tas VI-VIII 73 :~3 1 l 15 8 To kce 6-7.8 72,0. Cp. t� a ~'I11 80 f7 13 15 9 , -1,23 118,0 Cp. S~P s X-IX 100 6i 14 15 f' 12 15 10 ~ -0,33 G5,9 S S~P a VI-VII fii ~:i 14 . l5 11 a -I,0 115,0 Qp, S~P s X-VI � - 12 s -5,0 31,0 Cp. taH I-VI 30 :,3 11 l5 13 ~ 0,5 63~U A~ p aall~ VI &.-A. 16 87 50 7 15 SAAt V11~7O 0,6 46,0 SA~~ ~ 0 � � 2,1 80,0 ~ np~t ~A p a I-VI K: A.