JPRS ID: 9121 USSR REPORT METEORLOGY AND HYDROLOGY NO.3, MARCH 1980

Document Type: 
Collection: 
Document Number (FOIA) /ESDN (CREST): 
CIA-RDP82-00850R000200090007-3
Release Decision: 
RIF
Original Classification: 
U
Document Page Count: 
187
Document Creation Date: 
November 1, 2016
Sequence Number: 
7
Case Number: 
Content Type: 
REPORTS
File: 
AttachmentSize
PDF icon CIA-RDP82-00850R000200090007-3.pdf9.45 MB
Body: 
APPROVEE3 FOR RELEASE: 2007/02/08: CIA-ROP82-00850R000200090007-3 3 JUNE 1980 ME NO. 3f MAi4Ci-i 1980 I OF 2 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-44850R000200094407-3 FOR OFFIC[AL USE ONI.Y JPRS L/9121 3 June 1980 USSR Report METEOROLOGY AND HYDRUIOGY No. 3, March 1980 ~BIS FOREICN BROADCAST INFORMATION SERVICE FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 NOTE JPRS publications contain information primarily from foreign ' - newspapers, periodicals and books, but also from news agency transmissions and broadcasts. Materials from foreign-language 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 IL are supplied by JPRS. Processiag indicators such as [Text] _ or [ExcerptJ in the 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- 4 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 appropriate in context. Other unattributed parenthetical notes within 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 farther information on report content ca11 (703) 351-2938 (economic); 3468 (political, sociological, military); 2726 (life sciences); 2725 (physical sciences). COPYRIGHT LAWS AND REGULATIONS GOVERNING OWNERSHIP OF MATERIALS REPRODUCED HEREIN REQUIRE THAT DISSEMINATION OF THIS PUBLICATION BE RESTRICTED FOR OFFICIAL USE ONLY. APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY I USSR REPORT METEOROLOGY AND HYDROLOGY No. 3, March 19$0 JPRS L/9121 3 June 1980 Selected articles from the Russian-language journal METEOROLOGIYA I GIDROLOGIYA, Moscow. F CONTENTS Dynamic-Statistical Parameterization of the Ocean's Thermal Memory ~ (Sh. A. Musayelyan) 1 Zonal Dietribution of Cloud Cover Over the Earth (T. G. Berlyand, et al.) 13 Determinatione of the Mass Concentration of Aerosols in the Plumes of Industrial Plants Using a Lidar (I. M. Nazarov, et al.) 25 Change in Microstructure of Stratiform Clouds Under Influence of Condensation Nuclei (K. B. Yudin) 37 - Probabiliatic Approach to the Objective Classificatiou Problem (A. A. Burtsev) 45 = Structure of Tropical Cyclone "Carmen" Determined from Aerological Sounding Measurements Over the Ocean - (V. V. Galuahko, et al.) 53 Accuracy in Determining Atmospheric Ozone Content Using Data from . Measurements of Outgoing Radiation (Yu. M. Timofeyev, et a1.) 61 Model for Computing Ttnickness of the Quasihomogeneous Layer in the Ocean C (T. R. Kil'matov and S. N. Protasov) 7) - a - [III - USSR - 33 5& T FOUO) FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 , Computation of Characteristics of the Quasi-Isothermal Layer in the ' Equatorial Zone of the Ocean (A. B. Polonskiy) 80 Opposition in Ice Re3istribution in the Waters of the Foreign Arctic (V. I. Smirnov) .............e................................, 91 Methods for Regime and Operational Determination of River Runoff (I. F. Karasev) 98 Objective Method for Evaluating the Probable Success af an Individual Hydrological Forecast (G. N. Ugreninov) 113 Distribution of Current Velocity and Turbulent Friction Near the Water- Air Interface (N. K. Shelkovnikov, et al.) 118 . Method for Measuring the Thickness of a Petroleum Film on the Surface of a liater Basin (T. Yu. Sheveleva, et al.).................................... 124 Computation of Mean Oblast Yield of Winter Wheat (A. R. Konstantinov and D. N. Peradze) 130 Computation of Atmospheric Counterradiation (N. I. Rudnzv) 137 Nephelometric Method for Dctermining Light Attenuation in the UV and Visib].e Spectral Regions (T. P. Toropova) 141 Transmission of Intpgral Solar Radiation by a Grass Cover (T. K. Tammets) 146 Optical Properties of Crystalline Clouds . (0. A. Volkovitskiy) 150 - Sixtieth Birthday of Solomon Moiseyevich Shmeter 164 _ At the USSR State Committee on Hydrometeorology and Environmental Monitoring (V. N. D�rozdov and V. M. Voloshchuk) 166 Conferences, Meetings and Seminars � (L. S. Speranskiy, et al.) 169 _ News from. Abroad (B. I. Silkin) .............o................................. 178 - b - FOR OFFICIAL USE ONLY I APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY c PUBLICATION BATA English title : METEOROLOGY AND HYDROLOGY - Russian title ; METEOROLOGIYA I GIDROLOGIY A Author (s) . - Editor (s) : Ye. I. Tolstikov Publishing House : Gidrometeoizdat Place of Publication ; Moscow Date of Publication : March 1980 Signed to press 18 Feb 80 Copies ; 3780 COPYRIGHT . "Me'teorologiya i gidrologiya", 1980 � c - FOR OFFICIAL USE ONLY � i APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 F'i; : OFFICIAL USE ONLY UDC 551.509.33 DYNAMIC-STATISTICAL PARAMETERIZATION OF THE OCEAN'S THERMAL MEMORY Mb scow METEOROLOGIYA I GIDROLOGIYA in Russian No 3, Mar 80 pp 5-14 (Article by Doctor of Phyaical and Mathamatical Sciences Sh. A. Musayelyan, USSR Hydrometeorological Scientific Research Center, submitted for publica- tion 5 July 1979] Abstract: The problems of predictability of atmospheric processes and long-range weather forecasting are discusaed. Particular atten- tion is devoted to the role of external ener- gy sources. The author emphasizes the role - of cloud cover over the ocean as a regulator of the heat influx in the process of absorp- tion of solar radiant energy by the ocean. A dynamic-atatistical method for parameter- - ization of the ocean's thermal memory is out- 1 ined . [Text] In the short-range forecasting of ineteorological fields the atmo- sphere in the first approximation can be regarded as an isolated me33um and - the processes transpiring in it can be regarded as adiabatic. In other words, in the first approximation the short-range forecasting of ineteor- ological fields is the Cauchy problem a problem with initial data. Ac- cordingly, in numerical shcrt-range weather forecasting the principal fac- tor is the quality of analysis of the fields of initial data. [T11e pro'..lem of bnundary conditions will not be discussed here (Author's note).] In ac- tuality, however, the processes transpiring in the atmosphere are strictly nonadiabatic, and when one speaks of study of atmospheric behavior over ~ long and very long time periods the external energp conditions must be re- $arded as determining factors. And this means that the problem of modeling of short-period climatic varigtions or long-range predictiofl of ineteorolog- ical fields is the Cauchy problem, but with a source. In other words, in . this case it is insufficient to have only qualitatively analyzed initial - fields; it is also neceasary to take into account, and rhis is partic-lnr_- ly important, the principal macroscale nonadiabatic factors. _ 1 tvin /1L'0TnrA1 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 L.,A.. - The contributions of initial data and external energy sources to numerical weather forecasting for different times in advance have been examined by a number of authors (for example, see [1, 4]). These contributions can alsu be inveatigated directly on the basis of the firat law of thermodynamica. Qualitative Analysis of Contributions of Initial Data and Nonadiabatic - Factors As is well known, the mathematical expression for the first law in thermo- dynamica, applicable to atmospheric conditions, is the energy equation, - which can be written in the following form: _ Or dT di k = E. - UF t u d H + a Sin N d A a'= A Here t is time, /k is geographic longitude, e= 7T12 - so, Pis latitude, ve ( 9, 21 , t) and va (B, A , t) are the velocity components al.ong the ~ B and A axes respectively, T( e, X , t) is temperature, a is the earth's ~ mean radius, k is the coefficient of horizontal macroturbulent exchange, A is the Laplace operator in spherical coordinates and the function ~ E( B, i1 , t) integrally describes all the heat influxes. We will assume that it is necessary to obtain a solution of equation (1), periodic relative to A , symmetric relative to the equator and limited at ' the pole, satisfying the initial condition T(d, X, t) jr=o= To (g,. i.). (2) For convenience in the subsequent presentation we introduce the notation - vp dT vX dT - a d 9 a sin d d?. - I(9, i., t) and equation (1) is rewritten in the following form: ar L A T~E-}-I. (3) dt a= Aesuming that the functions E(8, a, tj, I( e,'1., t) and Tp( e,~) can be repreaented in the �orm of series in spherical functions E-_ Re r~ EA ~t~ E-'m a Pn (COS 6), ~.r n m ~4 ~ I:= Re ,~,~J: I� (t) C- 'm I Pn (COS e)r nm T, = Re J:G, Totn e' r"' a Pn (COS fj), /t fi we will seek a solution ot equation (3) in the form T= IZe 11 Tm (t) C- iin a E7- (COS H)� (5) nm 2 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY In formulas (4) and (5) E1�(t), In(t) and Tn(t) are complex functions, TO n are complex numbers, and ~(cos e) is the adjoint Legendre p,~lynom~al. It ie easy to show that the considered problem is reduced to integration of the equation dTn"l k _L dl ~ tt= 11 ~J1 -j- ] ) 7''" _ ~ni -t j~ n (6) with the stipulated condition jm (t) I t-0 = Tomn� (7) _ The solution of problem (6)-(7) can be written in the form � /.n (t) - akn c Xp I a, 1l (lt 1) ll t ~ J) ~8) tl Ton ex p In Q, (I[ 1) t~. where ETn and In are values of the functions En(t) and In(t), averaged in . the considered time interval. Formula (8) is prognostic. The second term on the right-hand side of this formula shows that in the predictlon of temperature the role of initial data actually decreases.with time, and this decrease occurs exponentially. However, the first term on the right-hand side of (8) shows that the con- tribution of the energy sources increases with time in conformity to the law 1- eXp n, !t (lt -t- 1) t~. L It is entirely obvious that both these processes (decrease in the role of initial data with time and increase in the contribution of energy sources) occura differently for different wave numbers n. In addition, the intensity of both these processes ia essentially dependent on the value of the caef- ficient of horizontal macroturbulent exchange k, which, as is well known, variea in a rather broad range. For waves of the scale of middle-latitude cyclones nZ 5 and for the value k= 106 m2/sec the contributions of the above-mentioned processes are rep- resented in Fig. 1. In this figure we have plotted time in days along the horizontal axis, and along the vertical axis some conditional rp value characterizing the contribution of some factor to the precomputed meteorological f ield (in easence, this is the parameter characterizing the success of the predic- - tion, for example, the correlation coefficient). 3 F(1R (1FFTf;TAT. ttCF (luT v APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 o,e' \ - 2 0.61 - 04 - - A' ZC JO 40 JD t r - Fig. 1. Contributions of initial data (1) and hest influxes (2). The following cor.clusiona can be drawn on the basis of an analysis of Fig. 1: 1) In the time interval 0< tG 5 days the success of the forecast is deter- mined pr-Imarily by the quality of analysis of the initial field and the contribution of the energy sources is secondary. In this case the fore- casting problem in the first approu;.mation can be regarded as the classic- al problem with initial data. According to the terminology which is now - in use forecasts for up to five days are called short-range. 2) In the time interval approximately from one to two weeks the role of the above-mentioned factors is approximately identically important, that is, the success of the forecast is determined to an approximately equal degree both by the quality of analysis of the initial field and the con- trib ution of energy sources. In this case the forecasting problem can be cons idered as the Cauchy problem with a source; both these factors are identically important. This is a case of intermPdiate forecasts. 3) In the time interval uing decrease in the ro tion of energy sourcea. the Cauchy problem with by energy sources. This and a season. from two weeka to one-two months there is a contin- Le of initial data and an increase in the contribu- In this case the forecasting problem is evidently a source, but the main role is nevertheless played is the case of long-range forecasts for a month 4) Finally, in the time interval frnm one-two months to infinity the role - of initial data becomes negligib].e and everything (or almost everything) is determined by the external energy conditions. Th{s is the case of a superlong-range or climatic forecast. Othe r classificstions also exiat. We will return to formula (8). It is easy to show that the time 20, at whoae end the role of initial data becomes equal to the contribution of the energy sources, is determined using the formula 0,7 a2 - To kn (n + 1) ' 4 . FOR OFFICIAL USE ONLY ~ APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY Below we give the 0 valuea for different wave nt;mbera with k- 106 m2/ _ sec. n 1 2 3 4 5 6 7 8 9 10 -6 o days 162 54 27 16 11 8 6 5 4 3 The cited data show that the ultralong components of the initial fields _ (n - 1, 2, 3) are conserved for a very long time in the atmospheric mem- ory and only after one or more mpnths their role becomes less than the contribution of energy so+srces. Waves of a synoptic acale (n s,*-5) in the atmoepheric memory, however, are conserved only for approximately two weeks, after which the role of nonadiabatic factors increases. Finally, for the high-frequency components the role of initial data de- creases very rapidly and literally after several days the external energy sources become the determining factor. Another conclusion which follows from formula (8) relates to the predict- ability problem. As is well known, the traditional interpretation of the problem of predic- tability of atmospheric processes is based on the conclusion that the err- ors contained in the initial fields are doubled each two to four days, as a reault of which predictability is a definite limit equal to five-aeven _ days (for example, see [6]). This, to a certain degree, agrees with the just mentioned formula. Moreover, curve 1 in Fig. 1 shows that the con- tribution of initial data decreases with time independently of their ac- curacy such is the structure of the equations in hydrothermodynamics. But, according to this same formula, with a decrease in the role of ini- tial data the contribution of energy sources increases. Therefore, to apeak of atmoapheric predictability, taking into account only the role of initial data and laying aside the contribution of energy sources,means - from the very beginning to set a limit to predictability of 5-7 days. ' Z'his follows from formula (8) and can be seen clearly in Fig. I, The na- ture of atmospheric processes is evidently such that allowance for only one factor the contribution of initial data leads to a sharp limita- tion on the predictability limit, whereas the addition of another factor energy sources makes this limit longer. In any case, such is the nature of the equationa of hydrothermodynamics, and this rigorously fol- lows from formula (8) and can be seem clearly in Fig. 1. The structure of the prognostic formula (8) is such that when having even precise initial data, atrictly speaking, it is impossible to precompute the predicted field precisely even for one day. However, if in addition precise data are available concerning the sources, uaing this formula it is possible to precompute the temperature ffeld for any ti.me with the ac- curacy with which equation (3) is integrated. However, the initial fie1d ia known very approximately and only for a part of the earth constituting 5 F(1R (1FFTr7AT. iTCx Our v APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 about 20% of its entire surface. With respect to the energy sources, ,judg- ing from everything, our ideas concerning them are extremely vague. There- fore, when reference is to long-range weather forecasting the principal = problem is that of discriminating and parameterizing the most important ~ large-scale energy sources controlling the behavior of the global atmo- , sphere over long time intervals. Everything set forth above indicates that when refereace is to short-range weather forecasting the term "atmo- ~ spheric predictability" in general is an apt expression. However, if ref- erence is to long-range forecasting of ineteorological fields the tradi- tional problem of predictability must evidently be dealt with in a new way and the *_erm "atmospheric predictabil ity" itself must be replaced by some other, more suitable term, such as the term "long-range predictabil- , ity of the eai�*_h-atmosphere system." In equation (3) and formula (8) the temperature T( e,X , t) and its ex- panaion coefficients n(t) are functions highly dependent~on time. Now - we introduce into consideration the climatic temperature T( e,,A ) and we will ass ume that it is slightly dependent on timp (as it actually is), so that a T/ a tN 0. Then, if^we denote the expansion coefficienta T( e ,~k) in spherical functions by 1n it is easy to show that 2 FE" + IR ~ (g) n kn (n + 1) . dn t�z other har.d, on the basis of (8) it is easy to find (10) lim Tn (t) r-oo kit (+t t I) ' - Comparing (9) and (10), we have - , T~'= lim Trt (t). (-U) Thus, according to formula (8), modern climate is some asymptomatic state of the real atmosphere with real boundary conditions and external energy - eources. 'Lt should be noted that the conclusions drawn above do not have a categor- ical nature due to a number of simplifications of the considered problem (in particular, no allowance ia made for the process of energy exc?-�--:ge between waves of different scales). However, it appears that these C-on.- - clugiQns, at least qualitativelq, agree with modern concepts. The enumerated classifications are based on a very important fact: the irr itial data on the state of the atmosphere sooner or later "are �orgotten" and the factors controllino atmospheric behavior over long time intervals will be the external energy sources. 6 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY Thus, for developing a long-range weather forecasting method it is neces- - sary that we first discriminate and parameterize those factors which con- trol the processes transpiring in the "earth-atmosphere" system over the course of Iong time intervals. ! L7 9 P Y! FIr Pd! LT d dI dIl - S,1 S.2 S,? 5,1 3,7 5,1 5,2 5,1 �,7 S,? 5,2 5,2 t :Jw oa� oi�� c;~ oa. or� -oa�,c res7~o 'erti;e;ro Qz. ,e,m or~ q:ur~,se.~y 1 _ ;GU' J== � ;Ji 0~ D :OJ~~�~~ j.:w ~ ~0� %-0' ' " Ob 0~=~' :i~ ia 1�`D ' - - - - ~ xti Jb0 S;U 5[U n di ~v 700 6 lA,~ 1s 0 I DmG'FPIIPIDILddIdII1 Fig. 2. Mean monthly vertical water temperature profiles and Qi values. av-t c 012 -012 -0,2 0 D -42 _ 02 0,30 -0,4 'T -H -0,2 -0,41 IX -.T -0, 2 i7D-!X 01f 0.2 p OZ �OES '02 0,,'tlf Y YI ~J � 0, 6 -0, J4 02 0 p jv-P - os -~z D, 56 o,~ L149 �D,? ~ 59 \\0,6 1 1 1 , ~nz 1 ~ 4 6 B 10 z Fig. 3. Cross-correlation matrix of moving averages of the S'(tm) and T'(tnl - +t) values for two months. ~ Comments on Thermal Memory of Ocean In an investigation of large-scale atmospheric processes it can be assuined that the sun is the only energy source. 7 FOR OFFT(:TAT, TTSF nNr v 0 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/48: CIA-RDP82-00850R000200094407-3 rvn vrVi%,.LnL .,..L. The nature of atmospheric-oceanic processes is such that the ocean in the warm half of the year absorbs solar radiant energy, transforms it into thermal energy and "remembers" it, that is, accumulates it in its depths. In the cold half of the year, using different heat exchange processes, the ocean transmits this heat to the atmosphere. In this case the cloud _ cover is the principal regulator with a feedback [3-5, 7]. In order to trace this process we will turn to Fig. 2, where on the basis of data for weather ahip M(66�00'N, 02�00'E) we have shown the mean month- ly vertical water temperature profiles for all the months of 1953. Along the vertical axie we have plotted depth in meters, and along the horizontal axie time t(the months are denoted by Roman numerals) and temperature T in degrees Celsius. "Qp" denotes the heat content of an ocean layer with a thickness 150 m with a temperature conatant with height, equal to 5.2�C. "Qi" (i = 1, 2,.... 12) denotes deviations of the heat content of the men- tioned layer from Qo (heavy shading). An analysis of this figure indicates that an increase in the heat content of the active layer of the ocean occur_a in the warm half of the year, where- _ as in the winter the accumulated heat ia transmitted to the atmoaphere. This can be seen particularly well in the inset to Fig. 2, where the solid curve represents the approxi-mate depth of the homogeneous layer of the ocean, whereas the dashed curve represents the ratio Qi/Q4� For greater clarity we will cite the following well-known example [4]. If the temperature of the upper 100-m layer of the world ocean is reduced by 0.10C ' and it is assumed that all the released heat is expended on heating of the atmospheric air, the temperature of the entire atmosphere to an altitude 30 km ie raised by approximately 5-6�C. - Thus, in developing a method for long-range weather forecasts an acute need arisea for parameterization oE the thermal memory of the ocean. Role of the Cloud Cover of the Oceans The process of absorption of solar radiant energy by the ocean, its trans- - formation into internal energy stored in the warm half of the year and the release of this accumulated heat to the atmosphere in the cold half-year cannot for the time being be parameterized in general form with the accur- acy necessary for the development of long-range weather forecasting schemes - on thie basis. The author of [5, 7] gave a certain phenomenological interpretation of this problem. There, instead of the temperature of sea water, he proposes the use of information on the cloud cover of the oceans received from meteorological satellites. In particular, it was demonstrated that there is a quite close asynchronous correl3tion between anomaliea of summer cloud _ cover over the ocean, averaged for three-month time intervals, and devia- tiona of winter air temperature on the continent from the norm. 8 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL U5E Oi1LL Further investigations indicated that such a correlation also exists with ' two-month averaging of anomalies of the above-mentioned meteorological el- ements [2]. As an example, we will examine Fig. 3, where by means of isc,- carre.lates we have represented the correlation matrix for Leningrad. in computing this matrix on the basis of data from meteorological satellites for the 10 years from 1965 through 1974 we formed the mean monthly values for total cloud cover at the points of a geographic grid 5� x 5� for the part of the North Atlantic bounded by latitudes 40 and 65�N. Using these data we determined the norms for each month in the year and then computed the corresponding anomalies and averaged them for the considered part of the ocean. The monthly temperature anomalies used were those employed in the operational work of the USSR Hydrometeorological Center. iJsing these data we computed the moving mean two-month anomalies of the considered meteorological elements. Along the horizontal axis in Fig. 3.we have plotted the time shift 'C , and along the vertical axis the months of the year. The nature of the field of isocorrelates in this figure is approximately the s.ame as for the case of three-month averaging [S, 71, to wit: the field of isocorrelates is rep- resented by alternating positive and negative correlation regions extend- ing along the mairi diagonal of the matrix. It is easy to discriminate a region of negative correlation between the characteristics of summer cloud cover and winter temperature. Such cross-correZation matrices were computed for each of 31 stations uni- formly distributed over the European USSR. An lnalysis of these matrices made it possible to discriminate the most informative indicators for for- mulating a statistical scheme for predic.ting air temperature anomalies for the European USSR on the basis of cloud data registeYed over the North At- lantic. These prediEtors are listed in the table. One of the important peculiarities of this table is that in the prediction of the temperature anomaly for the cold half of the year the number of predictors invariably includes the cloud cover anomaly for the warm half-year. A prognostic re- gression equation was written on the basis.of this table. The regression coefficients were determined using data from a teaching sample for 10-12 years by the ieast squares method. Then, using data from rin examination sample we computed a series of forecasts of two-month air temperature anomalies for the European USSR. These predictions are usually Punrly suc- cessful with respect to sign, especially in those cases when tllere is a predominance of zonal circulation in the lower half of the troposphere. At the present time the prediction method is in the stage of operatio:al tests. The preliminary results of these tests are encouraging. s Thus, like the earlier published studies [5, 71, the results presented in this article convincingly demonstrate that there is a quite close negative asynchronous correlation between summer cloud cover anomalies over thc, North Atlantic and deviations of winter air temperature from khe norm ii 9 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 the European USSR. And this means that the heat absorbed and rlccumulaled in the active layer of the ocean during the warm half of the year iti rv- tained for a very long time in its depths and is released to the atmo- sphere during the winter months. Predictor Months Using Cloud Cover Anomaly for Temperature Anomalies Prediction months Predictor months .Tanuary-February' April-May ' February-March August-September March-April September-October . April-May September-October rtay-June February-March June-July February-March Ju].y-August January-February August-September June-July September-October . April-May October-November June-July November-Decemter May-June December-January July-August Dynamic-Statistical Approach to the Problem of Parameterization of the Oceanic Thermal Memory On the basis of the results obtained in the preceding sections we will at- tempt to formulate some phenomenolog ical approach to the problem of dynam- ic-statistical parameterization of t he process of the thermal effect of the ocean on the.atmosphere [5]. Bearing in mind the. above-r^ntioned conclusions as a first approximation, we will assume that the quantity of heat Qr( 0,2l, t) accumulated in the active layer of the ocean during some summer month is proportional to the function - N, (P, t) = 1-0,1 S, (H, t), ~ where r is the numbering of the months during the warm half-yeaX, Sr(0, P1, t) is the corresponding total quantity of clouds. It then can be written that C' Qr (H, ~'r t) _ ~r ~vr L) fY, (0, Xe t), (11) are some still unknown functions to be determined and which are where , r assumed to be quasiuniversal. It is,obvious that y., when Sr = 0/10 Q~ - ^ ' i o when Sr = 10/10 , , 10 FOR OFFICIAL USE ONLY \ APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY It evidently can be postulated that the quEntity of heat Qq released by the ocean into the atmosphere in the q-th winter month consists of the quantity of heat accumulated in the active layer of the hydrosphere dur- ing each of all the summer months. The contribution of each month of the warm half-year must enter with its weight. Now we will assume that as a teaching sample we have an archivea containing observational data for yeara. For each pair of succeseive warm and cold halves of the year belong- ing to this sample we have, on the basis of (11) Q�%C/ N' (12) Accordingly, the heat influx equation for predicting air temperature (with stipulated initial and boundary conditions and known horizontal velocity components) of the q-th winter month, with (12) taken into account, can be written in the form dT vd dT uk dT k dt + a d b + n sin 8 d?~ - nT ~ T= X, N,. Q (13) r For convenience in exposition, we will call the functions Y r( e, the asynchronous influence functions. ~ Our objective is a determination of elements of the matrix 1f xr(~, A)II functions for each pair of numbers r c y~ and q c yl possible for real at- mospheric-oceanic conditions. Equation (13) describes the change in air temperature caused by the thermal effect of the ocean on the atmosphere, advection and horizontal macrotur- bulent exchange. A distinguishing characteristic of this equation is that its right-hand side was written in a form with a"lagging" argument and its computation requires the availability of data on cloud cover only for the times which have elapsed (that is, cloud cover need not be predicted). - Since our objective is a determination of the asynchronous influence func- - tion )Cr, therefore, assuming that in equation (13) the functions T, vg , vIN _ and Nr are known, and x r are to be found, we denote its right-hand side by Fq( Q, ,X) and rewrite it in the following form: ~ N, F`i. r Hence the asynchronous influence functions )tr can be determined by differ- ~ ent methods. However, it seems that in this case the spectral method is the most natural. In this case the problem essentially involves solution of a system of linear algebraic equations for the expansion coefficients of the sought-for functions [5]. The degree of practical applicability of the asynchronous influence func- tions determined in this way can be Judged only after carrying out exten- sive diagnostic and prognostic experiments during which there should also be a atudy of the nature of the temporal variability of these functions. 11 RnR nl?FTrTeT rTeV rnvrv APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 CUK Vrrl%..LnL u.._ Uaing the proposed method for parameterization of the oceanic thermal mem- ory it ia evidently impossihle to study atmospheric-oceanic proceBSes with any detail. However, we hope that on the basis of such parameterization it will be possible, in a rough approximation, to describe the principal large-acale procesaes of the thermal effect of the ocean on the atmosphere during the cold half of the ,year. - - BIBLIOGRAPHY 1. Blinova, Ye. N., "Hydrodynamic Prediction of Mean Monthly Temperature Anomalies for the Earth's Northern Hemiaphere Using IGY Data, DOKLADY AN SSSR (Repo-ts of the USSR Academy of Sciences), Vol 131, No 2, 1960. 2. Zadorozhnaya, T. N., "Evaluations of the Correlations Between Cloud Anomaly Fields Over the Ocean and Temperature Anomalies Over the Con- tinent With a Two-Month Averaging Interval," TRUDY GIDROMETTSENTRA SSSR (Transactiona of the USSR Hydrometeorological Center), No 192, 1977. 3. Marchuk, G. I., Musayelyan, Sh. A., "Methods for Computing Variations of the Total Flux of Radiant Energy for the Purposes of Long-Range Prediction of Large-Scale Meteorological Fields," METEOROLOGIYA I GIDROLOGIYA (Meteorology and Hydrology), No 8, 1974. - 4. Monin, A. S., "Physical Mechanism of Weather Changes," METEOROLOGIYA I GIDROLOGIYA, No 8, 1963. 5. Musayelyan, Sh. A., 0 PRIRODE NEKOTORYKH,SVERKFIDLITEL'NYKH ATMOSFER- NYKH PROTSESSOV (Nature of Some Super-Prolonged Atmospheric Processes), - Leningrad, Gidrometeoizdat, 1978. 6. Lorenz, E. N., "The Predictability of Hydrodynamic Flow," TRANS. NEW YORK ACAD. SCI., Ser. 2, Vol 25, 1963. 7. Muasaelyan, Sh. A., "The Use of Cloud Cover Satellite Information for Quantitative Long-Term Forecasting," PROCEEDINGS OF THE SYMPOSIUM ON i ME'rEOROLOGICAL OBSERVATIONS FROM SPACE, Philadelphia, Pennsylvania, ` USA, 8-10 June 1976 (19th COSPAR). ; 12 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 s FOR OFFICIAL USE ONLY UDC 551.(58+576.2)(100) ZONAL DISTRIBUTION OF CLOUD COVER OVER THE EARTH Moscow METEOROLOGIYA I GIL.^.^LOGIYA in Russian No 3, Mar 80 pp 15-23 [Article by Doctor of Geographical Sciences T. G. Berlyand, Candidate of Geographical Sciences L. A. Strokina and L. Ye. Greshnikova, Main Geophys- ical Observatory, aubmitted for publication 11 September 19791 Abstract: A atudy was made of the latitudinal diatribution.of the quantity of clouds on the basis of new global monthly maps of cloud cover compiled on the basfs of data from the world network of ineteorological stations. It is shown that in both hemispheres the mean quantity of clouds over the oceans is greater than over the land (by approximately one-tenth). The aonthly quantity of clouds for the earth as a whole is virtually constant, regardless of season of the yeaf. The results are compared with similar data from earlier studies made for climatological generalization of cloud cover distribution. [Text] Status of problem and data used. The study of characteristics of _ the global cloud cover regime is now attracting considerable attention in connection with work on estimating the earth's energy balance, model- ing of climate and its changes and also solution of many other timely _ problems of a scientific and practical nature. - One of the first climatological generalizations of cloud cover data was made by Brooka, who in 1927, on the basis of a relatively limited volume of obaervational data obtained the zonal distribution of the quantity of clouds over the oceans, land and the earth's surface as a whole [9]. Af- ter this study the zonal characteristics of cloud cover were virtually not investigated for a lonq time. Only after three decades did K. Tele- gadas and J. London [10] compute the zonal values of the quantity of clouds over the northern hemiaphere for two months (April, October) and two aeasons (winter, summer). In characterizing cloud cover over the oceans they employed data from the atlas of climatic maps of the oceans 13 Fnu nFFTrrer TrcL' n*rrv APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 L'VL\ VL'aiva..u+ prepared by McDonald [12]. The results of these studies have found use in a number of climatological investigations, in particular, devoted to the modeling of climate. During the last 10-15 Veare the volume of cloud cover information has in- creased considerably as a result of the increasing series of observations, the implementation of international geophysical programs, determining the expansion of the ground network of stations in poorly studied regions, and also the launching of ineteorological earth satellites. This facilitated ~ the further atudy of the patterns of spatial-temporal distribution o� cloud cover. Using data from two years of observations from a satellite (1965-1966), J. Sadler [14] obtained the zonal distribution of the quantity of clouds for the tropical zone (30�N-30�S). - H. van Loon [18] used both ground and satellite information on cloud cover. On the basis of a generalization of [9, 14 and others] he investigated the characterietica of the latituclinal distribution of the quantity of cloada in the southern hemiephere for Jaizuary and July. In the fundamental publication ATLAS OF THE OCEANS [2], recently published in the USSR, there is a series of monthly cloud cover maps of the Pacific, Atlantic and Indian Oceans. They were compiled using data from regular and expeditionary observationa on ships over a long period of time (from the end of the last century through 1965-1970). Z. M. Makhover and L. A. Nudel'man [1] prepared an atlas containing a con- siderable number of different climatic characteristica of cloud cover for the southern hemisphere. In the course of preparation of the atlas maps they generalized for all montha of the year an extenaive volume of daily data from ground observations for a five-year period (1967-1972). Data on the distributior, of the quantity of clouds over the earth's surface, except for the equatorial zone (13�N-5�S), are cited in studiea by C. Shutz and W. Gates [15, 161. They were prepared for the purpose of obtaining init- ial material for investigations for the modeling of climate. Aa the charac- teristics of the cloud cover regime in the northern hemisphere the authors used the results of surface and satellite observations. The latter were made only at two times (midday and midnight). For the southern hemisphere the data were taken from the mentioned study by van Loon [18j. They also drew upon the reaults of generalization of a four-year series of observa- tiona (1967-1970) from a satellite, taken from the atlas published by D. Miller [13], relating only to the interval 1400-1600 hours local solar time and characterized the cloud cover conditions for a limited part of the day. The use of inadequately uniform observational data on cloud cover from the point of view of inethods for collecting information concerning it, the ob- servations made once or twice during the day and the relatively short (not 14 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL IrSE ONLY - alway$ corresponding to climatological requirements) obaervation periods all this could not but exert an influence on the accuracy of the valuea and completeness of the characteristica in the above-mentioned atudies. World monthly maps of the dietribution of the total quantity of clouda based on more uniform inaterial, covering a long period of observationa, were com- piled by T. G. Berlyand and L. A. Strokina. Some of these maps were publish- ed in 1975 [3]. The initial data used in their compilation was the informa- tion included in numerous handbooks, monographs, atlases, monthly summaries and yearbooks published by the meteorological aervices o.f different coun- tries. The results of modern investigations of the cloud cover regime in individual regions of the earth were also taken into account [4, 5, 7, 8, 17, 191 in this investigation. In the course of analysis of the various - initial data, in addition to surface data, use was made of data rrom sat- ellite obaervations of cloud cover [6, 13, 14]. In compiling monthly maps - of the quantity of clouds on the contine*ts the authors of [3] used data for about 3,500 points, of which 950 were for Europe, 1,150 for Asia, 450 for Africa, 300 for North America, 460 for Central and South America, 140 for Australia and Oceania and 30 for Antarctica (we note that in the men- _ tioned study by Brooks data for a total of 1,000 stations were uaed for _ compiling such maps [9]). At most stations the observation perioda exceed- _ ed 30 yeara, whereas for etations situated in the low latitudes from 15 to 20 yeare. Z'he distribution of cloud cover over the oceans was determin- ed for the most part uaing data from the MARINE CLIMP,TIC ATLAS [17], which contains data on the probability of the total quantity of clouda for indi- vidual regions characterizing the principal climatic regions of the oceans. For the area of the Atlantic Ocean the initial values were given for 93 regions, for the Pacific Ocean for 111 regions, for the Indian Ocean for 48 regions. This article is devoted to a further development of [3]. In that paper, on the basis of data taken from new world monthly maps, the authors computed the zonal values of the quantity of clouds for the continents, oceans and the earth's surface as a whole. Computation of Zonal Characteristics of Cloud Cover The mean latitudinal cloud quantities were determined by averaging the val- ues for 35 circles of latitude 30i = i-0, 5, 10, 15,...,85� with a longitude interval 5�. A total of more than 30,000 values were taken from all the monthly maps; these characterize the quantity of clouds at the intersec- tions in a regular 5� grid. - The mean value for a 10� i-th latitude zone, taking the oceans and contin- ents separately into account, was determined using the formula Ni _ 5 Nt ,vl+5 0.5 cos Tt_ 5~ ny_S) I + cos yr ~ n jj + 0, 5 cos v;~; ~ n0+s> ! > > I n! 0,.5 Ni-cos x_ N cos Y' 0 ~ s t s i i + ,5 Ni+5 co~ ft+5 15 FOR OFFICIAL USE ONLX APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 v..,.. l P Ult ur r i%. icil. U.._ where nij is the quantity of clouds at the j-th point of intersection of a regular grid at the circle of latitude ~Pi; Ni is the number of points of grid intersection at the circle 5pi, the subscript i= i0, 5, 10,..., - 85. The zonal value of the quantity of clouds over the earth's surface was com- puted as the mean weighted value of the corresponding data for the oceans and continenza, taking into account the areas of the latter. The results are given in tables accompanying the paper. The data in the table make it poeaible, more completely than up to this time, to ascertain the peculiar- ities of the latitudinal and annual variation of the quantity of clouds. over the earth's : urface. Tables 1 and 2 give the ar,nual variation of the zonal values of cloud cover in tenths over the land, oceans and for the earth as a whole respectively. Table 3 contains data on the mean monthly quantity of clouds for individ- ual continents and individual oceans. Analytical Reaults Table 1 indicates that during the course of the year both for the land and for the oceans in both hemispheres an identical latitudinal variation in the quantity of clouds with maxima in the temperate and equatorial lati- tudes and minima in the tropics is characteristic. The difference between the extremal values sometimes exceeds 3/10. The position of the maximum and minimum values does not remain constant during the course of the year. In each hemisphere during the warm half-year they shift in direction toward the higher latitudes; during the cold season the shift is toward the low latitudes. Thie ehift ia quite great, being 10-15� in latitude. The great- Pst withdrawal of the maxima of the quantity of clouds in the direction from the equator to the north is observed in the period from May through September; in the opposite direction from November through March. Ztaice a year (April, October) the extrema are situated in both hemispheres al- most symmetrically relative to the equator. The results of this analysis of the mean monthly cloud quantities for each of the 10� zones for the most part show that during the course of the year there is one maximum and one minimum. A maximum in the warm season and a minimum in the cold seaeon is the most commonly observed type of annual variation. It is characteristic of latitude zones situated above 50� in both hemispheres; with advance toward the polea there is an increase in the annual amplitude. The amplitude attaina its maximum values (4/10) over the north pole; this amplitude is twice as great as ovEr the sour', pole. The noted type of annual variation ia also characteriatic of the equatorial zone (20�N-20�S), where the annual amplitude is about 3/10 over the land and 1/10 over the oceans. 16 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/48: CIA-RDP82-44850R000200094407-3 FOR OFFICIAL USE ONLY 5ap~�~ Tenths . ~ � 8 a) . , y\~ . . . , . 7 / v r r 6 `J r \ I /x.11 ~l 1 / , S 1 ~ ~ w 9- ~ ~ ~ 6) B - ^ . ~ ' : Y/ : i�''.' . i / �~r , . S - ,�f } i ~i � 1 - 2 r j e-a9 ---5 Y .~6 7 ~ ~ ~ ~ ~ ~ ~_-1-_----�_�~.-=~ ~ ~ ' 90i 70 SO 747 0 10 70 g Fig. 1. Zonal distribution of cloud quantities according to data from dif- ferentauthora, a) January; b) July; 1) C. Brooks; 2) J. London; 3) J. Sad- ler; 4) H. ve.: Loon; 5) C. Shutz and W. Gates; 6) Z. M. Makhover and L. A. Nudel'man; 7) T. G. Berlyand, L. A. Strokina and L. Ye. Greshnikova 17 _ FOR OFFICIAL iTSR ONi.Y APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/48: CIA-RDP82-44850R000200094407-3 r Vn vr r i%,..na, w.,.. Table 1 Annual Variation of Cloud Quantities LUtipotal ( [ I II I !II I !V I V I VI I V11 I Vllll i`C I 1 I XI I XII I fo;t G 7 Han cyureA 80-70�c2 5,1; 6 3 ~i,fi 6 1 5,1 5 9 I 5,7 6 2 7,0 7 1 7,2 7 1 7,3 7 1 7,2 7 3 7,7 7 8 7,~4 7 8 6,3 9 6 5,3 6 5 G,;~ 8 6 70-60 . , , , , , , , , , , , , 60- b0 6,3 G,l 6,1 6,5 6,7 6,7 6,6 6,6 6,8 7,1 7,0 6,7 6, 6 " .50-40 5,7 5,8 5,9 6,1 5,9 5,5 5,0 4,7 4,7 5,1 5,8 6,1 5.5 40-30 5.1 5.0 5,0 4,9 4,5 3,9 3.3 3.7 3,4 3,5 4,0 4,7 4,3 30-20 3,7 3,6 3,7 3,6 3,6 3,7 4,0 4,2 3,7 3,0 3,4 3,6 3,7 20-10 3, z-, 3,6 3,8 4,0 4,6 5,2 6.0 6,2 5,6 4,1 3,9 3,6 4,1, .10-0 5,4 5,6 6,0 6,6 6,7 h,% 4,9 7,1 6,7 6,0 5,8 5,1 6,1 0-10� 103 7,2 7,3 7,3 6,8 6,4 5,6 .5,0 5,4 fi,0 f;,2 6,8 7,0 6.4 10-20 6,4 6 ,3 6,0 5,4 4,0 3,~ 3,1 3,0 4,2 :~.2 5,7 5,3 4,~+ 20-30 4,5 -a,i, 4,3 3,3 3,5 3,3 3,0 2.7 3,2 3,9 4,1 4,2 3,8 30-40 11 3 3 4 7 3 3 8 4 7 4,8 4.7 3,9 3,9 4,2 3,8 3,8 4,1 40-50 , 4,!) , 3,3 , 4,9 , 5,4 , 6,6 6,1 5,8 5,8 5,9 5,3 5,9 5,8 5,6 50-GO 60--70 7,1 6,7 6,7 6,5 5,8 5.7 5,9 6,4 6,4 6,4 6,2 6,3 6,3 70--80 6,2 6,3 5,4 5,7 5,0 4,9 �`t 0 5,7 5,8 6.2 5,8 6,2 5,7 60-90 5,4 5,7 5,9 4,5 4,2 4,1 3,9 5,0 5,0 5.U 5,4 5,7 5,0 'Ceecpnor 4 nonywapac 5,2 .5 .:i 5,2 5,4 5,6 5,5 5,6 3,6 5,5 5,3 5,3 5,2 5,4 .IOHaioc 5 no7yiuapF+e 5,9 5,9 5,8 5,3 4,8 4,4 4,2 4,4 4,9 5,3 5,4 5,6 5,2 Cywa n .ucnom 5,4 5,4 5,4 5,4 5,3 5,2 5,1 5,2 5,3 5,3 5,4 5,3 5,3 6 g Ha1c oxeauaMjs 90-80�C 5,3 5,2 5,3 5,7 7,7 8,5 3.7 9,0 I US ~ 8 5 8,0 1 8 6,2 5,4 ' 6 8 6 0 7,1I 1 7 80--70 6,0 5,9 5,7 6,3 7,7 8,0 8,1 8,4 , , , , , 70-60 7,0 6,3 6,3 7,0 7,6 7,7 7,7 7,6 3,U 7,8 7,4 7,1 7.4 -60-50 7,6 7.5 7,5 7,7 8,1 8.2 8,3 8,1 7,8 7,6 7,7 7,7 7,8 50-90 7,5 7,4 7,3 7,5 7,7 7,9 7.7 7,4 7,1 7,0 7,2 7,4 7,4 40-30 6,5 fi.r, 6,5 6,5 6,6 6,6 6.1 5,7 5,0 ~ G,I 6,2 6,4 6,3 30-20 J, r) l1,3 5,3 5,3 5,3 f,,[; 5,7 5,3 :i,'Z ~ 5,4 5,3 5,5 20-10 .5,2 5.0 :,U 5,3 5,4 6,1 6,3 6,1 (i,Q I 5,8 5,4 5,4 :),G 10�-0 5,8 Fi,6 53 5,9 6,1) 6,3 G,? 6.1 6,0 (>,O 5,$ ri.9 5,51 0-10� w 5,9 5,7 5.6 n,6 5,3 5,2 5,3 ~i,3 5,6 F,,fi 5,7 5,8 ' 5,& 10-20 5,7 i,5 5.6 5,3 M 5,0 ;5,3 :>,4 5,5 I 5,7 5,8 .`i,i ;r,.`, 20-30 5,G 5,5 5,7 5.6 :i, 5 5,4 5,6 5,6 5,8 6,0 .5,> (P, ~t.) = 3 K~; i:~ j)Z (1 2 n'� nt I''R ri P 1 ' (4) K (r, !t.) - � G Pn Kis (P, ~t) = j nn~K -}-2(1 - f ~n)~-~- (Yt - 1) g ta-n~ l where K(P) is the Van de Hulst function. ' > nrt' h~P~ 2 sin c sin-' ~ 1 - P ' ~ o= 4 r is a dimensionleas parameter) 30 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY (n)~ 1,05; n 1,5: C~ Thie function is correct ior any real values of the refractive index of perticlee. The mean appruximation error is 4% for a refractive index 1.5. Failure to take abeorption into account in K(p , n) gives errors oF a few percent for polydieperse media. The particle size function was atipulated in the form of a three-parameter gamma function f (r) = Rr'` e-bi (;L b (5) r which with a broad change in the b and i-L parameters describes well the diverae apectra of particles observed in the atmosphere. With (4) and (5) taken into account, the expression of the scattering function for the adopted model has the form [6] Nt~�~= 47tN + 1)(FL -F2) y�.ta b2 2 4' Y" (Et + ~ + y=) r) (6) -f-1/1 t y` cos (;.L + 1) ~1 1 + 2 b2N (IL+ 1)(lL -f- 2) X r n l~n - i X L 1 f(n) c�ty 1' ~ b n ~ p Y = 2=; C=1- 8nb x; x= r An analytical model of the volume scattering coefficient (6) can be used in determining the mass concentration and parameters of an aerosol in the interpretation of data from laser sounding. It makes it possible to reduce the problem of determining the micrnstructure of an aerosol to the optimization problem in parameter space by means of minimizing the functional k F r, N) _ ~ [~i Z) - P; (~�)1', (7) constituting the sum of the squares of nonclosures of the differences in the measured values of the scattering coefficients and the corresponding valuea in the model. In principle, for determining the four parameters of the diatributiion it is necessary to know the measured scattering co- efficients at four wavelengths. However, taking into account that the 31 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 variations of thP refractive index in the range from 1.3 to 1.7 lead to a change in the scattering coeff icients A) by not more than 10-15Y [10], with an accuracy in measuring )8(/A, z) of 20% the value of the refractive index can be selected taking into account the specific condi- tions of the investigated object, additional measurements of compoai- tion or taken as a mean. In other words, with a f ixed value of the re- fractive index in determining the aerosol parameters it is possible to limit ouraelves to three wavelengths. Having apecific experimental data on fi(,;1, z), uaing an electronic computer we obtain the numerical values of the selected P (a) model. The end result of the calculations is evaluations of the model parameters with a atipulated accuracy, deter- mined by the sum af the squarea of the difference in the values of the parameters i:, two successive iterations. With a fixed value of the re- fractive index n, the accuracy of the method gives an error in determin- ing the parameters of 3-5X [lO]. An assumption concerning the particle- size diatribution function does not limit the possibilities of the model employed because using the variations of r and � it is possible to obtain the most diverse particle distribution spectra. In order to obtain multi- modal spectra it is necessary to have measurements at a great number of wavelengths. Analyais of Inversion Results Tiue to the fact that two wavelengths were used in our experiment in sound- ing the plume, our method makes it possible to determine two parameters of the model. In this case the most interesting parameters of the plume should be the mean particle size and the numerical concentration. Precise- ly these parameters were determined in this experiment. Two other para- meters distribution width � and the refractive index of particles ft were selected from the accompanying samplings. The results of the samplings cited in Table 1, on the basis of the dispersion value, made it poseible to assume that the investigated aerosols have a quite broad size spectrum and therefore the � value should be selected equal to 0 or 2. The value of the real part of thp refractive index for particles of zinc, iron and copper, according to the Ivlev model [4] and cited in the book by Van de Hulst, can vary in the range 1.40-1.47. As mentioned above, such changes in n cannot introduce significant errors into the determination of mass concentration. - In additiony for work with the optimum parameterization method using an electronic computer it ie necessary to select two initial points contain- ing the surmised sought-for parameters n, the mean particle size r, � numerical N or mass M= N b 4/3110 concentration, where b is aerosol den- aity. In our experiment these data were also taken from microphysical _ measurementa. We note that in principle we can select any reasonable val- - uea of the parameters as the initial points. Their choice influencea only the convergence of the problem when seeking the minimum of the F function and the expenditure of computer time. In order to stipulate these points it is not at all mandatory to carry out simultaneous microphysical meas- urements, but with the availability of the latter the choice of initial points is simplified. 32 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 ; FOR OFFICIAL USE ONLY Thus, as a result of inversion of equation (3) by the optimum parameteriza- - tion method we obtained evaluations of model parameters as indicated in Table 2. Table 2 Evaluation of Aerosol Parameters by Optimum Parameterization Method � n r.41 KM _ I N c.u-3 I I eixa/M~~ ~,if % �m 2,0 1,40 0,20 5,0� 10' 3,2 S~0 0 1,47 0,50 3,7�107 74 40 2,0 1,44 0,58 5,7�107 985 50 Table 3 _ Values of Volume Scattering Coefficienta and Aerosol Parameters for Pure Atmosphere . - [IapaMer- I Afn.+t P(X. z) ceH Pb. �pac- i _ rtpe2e.7e- 1 HHA 2 Z , ~ � 0,.532 4�10-6 7-0,22 790 10 ' 0,694 3�10-6 �=2.0 1,064 1,7 � 10-6 n=1,33 0,532 2,1 � 10-e i=0.20 530 20 0,694 1,5�10-6 li-2,0 - 1,064 7�10-7 rr=1,33 KEY: - 1. �-m 2. Distribution parameters The magnitude of the error in determining the mass concentration !S M(with a density S= 2 g/cm3) was determined from L�he minimum nonclosure L1F, which characterizes the minimum deviation of our model p T(a) from the set of experimental data. - The mean radius of the particles falls in the region 0.3-0.6 � m and varies in dependence on the spectral variation of the scattering coefficient [2]. The concentration was equal to 5�107 cm 3. The mass concentration with a 50% error was 75 � g/m3, which corresponds to data from the sampling. TiLle 2 shows that the parametera most sensitive to variations of experimentai 33 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 data are N and the mean size of particles, and hence, to a high degree, also the mass concentration. Such high errors in determining M are at - tributable to this circumstance. For example, the inversion of the volume acattering coeff icients (Fig. 3) for a case when there were ac- companying,aircraft sa2plings with stipulated values 0, n= 1.47 made it possible to obtain r= 0.5~xm and a value of the mass concentration M = 74 � g/m3. The minimum value of the mass concentration, corresponding to the background value A (;1 , z) = 5�10-4 cm-1, is 2.6 �10'2 j.L. g/m3, whereas the concentration N is 1.3�104 cm 3. The maximum M value is 400 � g/m3/m3 and was observed with A(z) = 9.0 cm 1. The computed value of the scattering coefficient in the model for the assumed particle distrib- ution with N= 1 cm 3 is P(9l) = 3.77�10-8 cm 1f1.5�10-8 cm 1. The accur- acy in determininq the mass concentration in this case is 40%. ,a "r -",).o.u2 Ji't0 JIBO J6: U 3Bo: [ M Fig. 4. Dependence of the volume scattering coefficient on depth of pene- tration. Similar computations were made for other series of laser sounding of the plume; these are shown in Fig. 4. The mean values of the distribution and mass concentration parameters for these series are given in Table 2. As an illustration of the method described above we wi1Z give an example of inversion of the results of three-frequency laser sounding of the frAe atmosphere at Leningrad carried out in 1978. In this experiment we used an Al-Na garnet laser ( T= 1.06�.m and 0.532�-m) and a ruby laser (A _ 0.69 �.m) with a common receiving antenna designed in a Cassegrainian sys- tem with a diameter of the large mirror 0.26 m. Table 3 gives data on the composition of aerosol obtained with values of the meteorological range of visibility 10 and 20 km. 34 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY - 1'he refraction coefficient for the pure atmosphere was assumed equai to 1.33. The use of three sounding wavelengths made it possible to evaluate three aerosol parameters. The n value was stipulated, taking the experi- mental conditions into account. The determined values of tlie concentration N are in agreement with data in the literature for the free atmosphere [10]. Conclusion The results of this investigation of aerosol in smoke plumes using two-fre- quency laser sounding and the employed method fcr inversion of the volume scattering coefficients indicated that the mass concentration of particles obtained by this method agrees well with the results of sampZing. This confirms the considerable possibilities of remote sensing and monitoring of amoke plumes using laser techniques. In the case of three-frequency sounding the mass concentration can be obtained in the absence of ac- companying microphysical measurements. T~ao-frequency sounding also can be ~ used without accompanying microphysical measiirements if the element com- - position of the effluent is known. With an increase in the number of sounding frequencies the method describ- ed here makes it possible to obtain multimodal aerosol distribution spec- tra. Taking into account that enterprises of the same class give close spectra ' of aerosol particles in effluent, for solution of ~ypical problems in mon- itoring effluent it is sufficient to obtain statistically representative characteristics of effluent for individual branches of industry. This will make it possible ta organize remote routine monitoring of effluent by means of laser systems. BIBLIOGRAPHY 1. Badayev, V. V., Georgiyevskiy, Yu. S., Pirogov, S. M., "Aerosol At- tenuation in the Spectral Region 0.25-2.2 � m," IZV. AN SSSR, FIZIKA ATMOSFERY I OKEANA (News of the USSR Academy of Sciences, Physics of the Atmosphere and Ocean), Vol II, No 5, 1975. 2. Zhitkov, L. V., Kaul', B. V., et al., "Use of a Lidar for Investi.gat- ing the Dynamics of a Smoke Plume,"TRUDY 5-go VSESOYUZNOGO SIMPOZIUMA PO LAZERNOMU T AKIJSTICHESKOMU ZONDIROVANIYU ATMOSFERY (Transactions of the Fifth All-Union Symposium on Laser and Acoustic Sound:i.ng oi = the Atmosphere), Tnmsk, 1978. - 3. Zuyev, V. Ye., et al., "Optical Experiment and Results of ::nversion of Data from Multifrequency Laser Sounding of the M4 crostructur2 of 'Aerosol in the Surface Layer," PROBLEMY DISTANTSIG`NNOGO 7.ONDIROVANIYA _ ATMOSFERY (Problema in Remote Sensing of the Atmoapheie), Tomsk, SO AN SSSR, In-t Optiki Azmosfery, 1976. 35 Vnu nVVTrrnT rrcV nlknv APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 rvn vrri%.iti1j u.... ....L~ 4. Ivlev, L. S., Popova, S. I., "Complex Refractive Index of Matter of - the Disperse Phase of Atmospheric Aerosol," IZV. AN SSSR, FIZIKA ATMO- SFERY I OKEANA, Vol 9, No 10, 1973. 5. Kovalev, V. A., "One Method for Pracessing a Laser Signal," TRUDY GGO ('I'ransactions of the Main Geophyeical Observatory), No 312, 1973. 6. Martynova, V. I., "Polydisperse Scattering Coefficient for the Mie Approximation Function," TRUDY IPG (Transactions of the Institute of Applfed Geophysics), No 17, 1973. 7. Nazarov, I. M., Nikolayev, A. N., Fridman, Sh. D., DISTANTSIONNYYE I EKSPRESSIVNYYF METODY OPREDELENIYA ZAGRYAZNENIYA OKRUZHAYUSHCHEY SREDY (Remote and Express Methods for Determining Environmental Con- tamination), Mnscow, Gidrometeoizdat, 1977. 8. Pakhomov, P. V., Martynova, V. I., "Mie Approximation Function for the Factor of Effectiveness of Scattering," TRUDY IPG, No 17, 1973. 9. Rozhdestvenskaya, V. I., "Determinatiun of the Concentration and Para- meters of Aerosols from the Results oi Tailight Measurements," TRUDY IPG, No 32, 1977. 10. Rozhdestvenskaya, V. I., RASEYANIYE SVETA NA POLI7I:iPERSNYKH SISTEMAKH _ AEROZOL'NYKH CHASTITS (Light Scattaring on Polydfsperse Systems of Aero- sol Particles), Moscow, IPG, 1977. 1 I I ~ ' I ~ 36 ~ FOR OFFICIAL USE ONLY ~ ~ APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY UDC 551.509.616 CHANGE IN MICROSTRUCTURE OF STRATIFORM CLOUDS UNDER INFLUENCE OF CONDENSATION NUCLEI Moscow METEOROLOGIYA I GIDROLOGIYA in Ruasian No 3, Mar 80 pp 33-38 (Article by K. B. Yudin, Institute of Experimental Meteorology, submitted for publication 19 Suly 19791 - Abstract: The article gives the results of f ield experiments for the modification of warm strati- - form cloud cover by artificial hygroacopic eon- densation nuclei. The author gives some charac- teristics of a modified cloud medium and describes the nature of restructuring of the spectrum of - cloud element sizes ae a result of modification. [Text] The problem of artificial modification of warm clouds and fogs is of great scientific and practical interest. One of the objectivea of - artificial modification is an improvement in visibility in clouds and foga. Their total dissipation ie not always necessary. Zt follows from the theory of light attenuation of dispersive media [10] that it is pos- aible to improve viaibility by a restructuring of the microstructure of a cloud or fog. After Yu. S. Sedunov [6] proposed a method for the mod- ification of microstructure by the introduction of additional finely dis- - , perse condensation nuclei directly prior to the onset of cloud formation, aseries of theoretical studiea [1, 2, 5] was carried out in which a study was made of the effect of nuclei with different characteristics. The , experimental studies carried out in an aerosol chamber [3, 8] indicated that the considered method affords a possibility for modification of the _ process of forming of the droplet spectrum. A merit of the method is a high degree of use of the reagent, attributable to the fact that in the - proceas of formation of the microstructure the total mass, which can be extremely small, participates, and to a certain degree determines its - large number of condensation nuclei. The ob,jective of this study was a field investigation of the regularities in changes in the microstructure of the cloud medium with introductian of a reagent. Artificial hygroacopic particlea were introduced into the cloud near the condensation level. They must be entrained into a cloud by ascending flowa, and competing with natural nuclei, change the course 37 LV1D A'VrTl,T AT T7Vn n*rr V APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 L'VIl VL'l lVlau~ v.. . v.�..~~ r-1 d ~ ~ E~ Of N Q~ ar tD u 1+ G! ~ cd rl r-I a w 44 V o~o ~ ~ 7 G! U Gl O tl1 -rl y p 0 ~0 ,~C G b J M Gl ,W U P 0 i4 MH N y H ~O h o0 ON O r"i ~ ~ ~ ~ a~ ~d ro ~ n a o w -W Q4 ~5 q ~n O ~ Ir ri � cNd -W O p ~z � � � rlNM~t V1 ~ APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY the modification zone can be represented as an elliptical cylinder with a horizontal axis of 7-8 km. The lengtha of the axes of the ellipse, lying at the base of the cylinder, 5 minutes after modification were 600 and 200 m, and after 10 minutes 1,200 and 400 m. The mean concentration of reagent after 5 and 10 minutes was 20 � g/m3 and 5 N.g/m3 reapectively. However, the reagent, diffuaing in the cloud under the influence of dynamic turbulence, doee not fill the entire volume uniformly and at individual points of ita concentration can differ coneiderably from the mean values. For the mon- itoring and regietry of changes in microstructure the aircraft repeatedly returned to the modification track at an altitude exceeding by 20-30 m the altitude at which the reagent is introduced. The flight trajectory, the moments of entry into the modification zone and emergence from it were regiatered by a ground radar. The time of presence of the aircraft in the modification zone was about two minutes. The principal instruments for investigating microstructure and its changes were a photoelectric instrument for measuring the sizes and concentration of cloud droplets in the range from 0.5 to 401km [4] and an instrument for - measuring aerosol scattering in clouds. During presence in the modification zone 100-120 microstructure samples were taken for analysis. The number of droplets in the sample varied from 500 to 6,000, depending on their concen- - tration in the cloud. The background data on the state of clouds were col- lected in regions adjaCent to the modification zone during aircraft man- euvers for return to the modification track. - Table 1 gives data on the synoptic conditions under which the experiments were made and some integral characteristics of microstructure computed from the measured apectra. A preliminary analyaie indicated that in the modification zone there is an increase in the concentration of hydro- meteora and aerosol scattering. However, there is no direct proportional- ity between theae changes, which indicates a restructuring of the spectra. In order to explain the nature of this restructuring each of the measured spectra was divided into two parts water-enveloped nuclei and dropleta [7]. The part of the spectrum containing water-enveloped nuclei was de- scribed by a power-law diatribution in the form Y'' f (r)=A ,~where A is a normalizing factor numerically equal to the value of the above power function with r= rp = 1 � m. The A and -V parameters were computed by the least squares method. The-parameters of the distribution function, de- - scribing the droplet part of the spectrum, were computed from the modal values of the radius and distribution function and the liquid water content � value, computed from the measured spectrum. Table 2gives the parameters of the distribution function (and not the parameters of the mean spectrum), averaged from individual spectra, and their variability in background meas- urements and in the modif ication zone 5-7 minutes after introduction of ttie reagent. The last four lines in the table generalize data on change in the parameters as a result of modification and on the significance levels f,)r theae changes, evaluated using the Bartlett and Fisher tests. 41 - FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 An analysis of the table indicates that as a result of introduction of the reagent there was an increase in the concentration n of water-encased nuclei, that is, most of the additional nuclei were encased in water, but could not form droplete. The droplet concentration m and the aerosol acat- tering coefficient o4increase somewhat. The changes of the exponent in the inverse power-law diatribution -V, the mean radius r of drop.lets and the relative dispersion of the droplet part of the apectrum p,/r are in- significant at the level of the natural inhomogeneity of clouds. A sub- stantial increase in fluctuatione of some parameters ia alao observed in the wake of the effect. f!r) Nxrf I� cn"i J000 - a) 1000 � \ J00 ~ 100 � JO � 10 ~ 45 1 2 S 10 0,5 b) Z \ - - J \ \ \ l . ~ J 1 r M,Y,y t ~ Fig. 1. Spectra of size of hydrometeors before modification and 5-7 minutes after cloud modification. a) background, b) wake of modification; 2) in etratocumulus cloud, 3) in haze under clouds, 5) in stratiform cloud The additional information on restructuring of the spectra is shown in the figure. The numbers on the curves in this figure correspond to the numbers of the experiments. The figure shows that in the modification zone the value of the distribution fu:,::tion in the region 0.5-0. 7 � m increasea by a factor of 2-3, there is a deepening of the "dip" between the parts of the spectrum, and droplets which are larger than those observed in background measurements appear. However, droplets capable.of causing coagulation (r ) 20pL,a) have not been regiatered. Our work demonstrated that the introduction of additional artif icial con- densation nuclei somewhat modifies the microstructure of clouds. Under natural conditions the process is traced for 13-15 minutes after intro- duction of the reagent. At the end of this period the concentrations of 42 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY the diffueing reagent and modified droplets become so small that the mod- ification effect ia insignificant and all the measured parameters of the cloud medium approach the backgrouad levels. On the other hand, thia circumstance indicates that repeated flights of an a3rcraft through the _ modification zone do not result in eubstantial changes in etratiform cloude and exert no influence on "purity" of the eaperiment. It ie evi- dently poasible to trace the further course of the process only with a - _ marked increase in the seeding area. A significant effect from modifica- tion which could be used for practical purposes has not yet been estab- lished. The author expresses appreciation to A. V. Lachikhin, V. P. Bedrenko and B. S. Yurchak, apecialiats at the Institute of Experimental Meteorology, for assiatance in carrying out field experiments and measurements. BIBLIOGRAPHY 1. Aleksandrov, E. L., Klepikova, N. V., "Influence of Artificial Conden- sation Nuclei on the Kinetics of Formation of Cloud Droplets," TRUDY IEM (Transactions of the Institute of Experimental Meteorology), No 8(46), 1974. 2. Aleksandrov, E. L., Klepikova, N. V., "Effect of Artificial Condensa- tion 2Juclei on Development of the Cloud Spectrum," TRUDY IEM, No 9(52), 1975. 3. Aleksandrov, E. L., Yasevich, N. P., "Results of Preliminary Experi- tnents for Fog Modification in a Condensation Nuclei Chamber," TRUDY IEM, No 20, 1971. 4. Aleksandrov, E. L., Lachikhin, A. V., Posadskiy, V. I., Yudin, K. B., "Aircraft Photoelectric Instrument for Measuring Cloud Droplets," TRUDY IIIM, No 19(72), 1978. 5. Klepikova, N. V., Sedunov, Yu. S., "Kinetics of the Droplet Spectrum in the Initial Stage of Condensation," TRUDY IEM, No 3(37), 1973. 6. Sedunov, Yu. S., "Numerical Experiment for Modifying the Kinetics of Formation of the Cloud Spectrum by the Introduction of Additional Con- densation Nuclei," TRUDY IEM, No 6, 1969. - 7. Sedunov, Yu. S., FIZIKA OBRAZOVANZYA ZHIDKOKAPEL'NOY FAZY V A7.`MOSFERE (Physics of Formation of the Liquid Droplet Phase in the Atmosphere), Leningrad, Gidrometeoizdat, 1972. 8. Smirnov, V. V., "Restructuring the Microstructure of Fogs Under the Iri- fluence of Hygroscopic Particles," TRUDY IEM, No 25(93), in press. 43 FOR OFFICIAT iTSF f1NT.Y APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 - 9. Skhirtladze, G. I., Yurchak, B. S., "Measurement of the Coefficient of Horizontal Diffusion in Cumulus Clouda by the Radar Method," IZV. AN SSSR, FIZIKA ATMOSFERY I OKEANA (News of the USSR Acadeary of Sciences, Phyaics of the Atmoaphere and Ocean), Vol 15, No 2, 1979. _ 10. Shifrin, K. S., RASSEYANIYE SVETA V MUTNOY SREDE (Light Scattering in a Turbid Medium), Moscow, Goatekhizdat, 1951. 44 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY UDC 551.509.314 PROBABILISTIC APPROACH TO THE OBJECTIVE CLASSIFICATION PROBLEM Moscow METEOROLOGIYA I GIDROLOGIYA in Russian No 3, Mar 80 pp 39-44 [Article by A. A. Burtsev, USSR Hydrometeorological Scientific Research Center, submitted for publication 13 July 1979] Abstract: In the probabilistic approach to objective classification the meteorological fields are not characterized as belonging un- ambiguously to some class, but by degrees of reliability of belonging to all classes. It ia demonstrated ttiat the proposed probabilistic approach ia more universal and has a number of advantaaes in comparison with the traditional determir.istic approach. In order to ascertain the degrees of reliability of such a determin- ation the author proposes the synthesis of a Golovkin algorithm with a potential functions method algorithm and its realization using ac- tual data. The results indicate an objective character of the restored degrees of.reliability. [Text] Recently in meteorology ever-increasing attention is being devoted to the problem of objective classif ication of the fields of ineteorolog- _ - ical elements. The timeliness of this problem is attributable not only - to theoretical and cognitive considerations, but al.so to the fact that the results in one form or another can be used in weather forecasting. Many authors, concerned with classification, indicate the difficulty [2], and some [7] even suggest the impossibility of creating a universal ob- jective classification in meteorology. The author of [6] makes an attempt _ to create such a classification, but by virtue of the factors which will be considered below, this-approach can be considered universal only tc a limited degree. In this article we propose a new approach making it pos-- sible to come close to solution of this problem. _ Among the objectiva classification algorithms special attention shotll.c' be given to the B. A. Golovkin algorithm [5]. For the first time in meteor- ology it was applied in [4] and then with some modifications in [3]. 45 L'/1D !NL'VT/~TAT Ti[+T ALTrtI APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02108: CIA-RDP82-00850R000200090007-3 This algorithm, constituting an approximate solution of the problem of linear programming with Boolean variables, makes it possible to discrim- inate in the set of objects X the centers of clustering of objects xi E X, and specifically, obtain a set of the most typical, or as they are uaual- ly called, "standard" objecta. Then the entire ensemble of objecta xi C X, i= 1,...,N (N is the number of objects) is broken down into classes by as- signing each xi to that class the distance to whose standard object is minimum. A merit of this algorithm is that it is completely free of any a priori assumptions concerning the spatial structure of the set of ob- jects to be classified and in the best way possible evaluates the number of classes and unambiguously distributes all the objects by classes. In this algorithm there is no need to stipulate any threshold distances; the only input data are the matrix (measuring N x N) of distances between each pair of objects. If the spatial structure of the set X is such that all xi can be broken - down into classes which are adequately compact and distant from one anoth- er, the Golovkin algorithm rather qualitatively solves the classification problem. The situation is far more complex in a case when the character- istic regions of the classes are situated close to one another and partic- ularly when the characteristic regions intersect. In this case objects _ which are extremely close to one another, situated in a zone of intersec- tion of the characteristic regions, can be assigned to different classes - and in th is case the probability of a classification error will be very - high. And precisely such a pecu7.iarity is characteristic for some large- scale fields of ineteorological elements. For example, it was established - in [3] that the fields of summer mean inonthly anomalies of H-SOQ (regis- - tered in the latitude zone 75-40�N during the period 1948-1977) have the property of intersection of characteriatic regions. This circumstance followed from the fact that after constructing a classification on the _ basis of the Golovkin algo:-ithm the zero hypothesis of equality of intra- - class means could not be refuted on the ba.sis of the Fisher test with a 5% significance level. The classification became significant orly after exclusion of a number of objects falling in the zone of intersection of the characteristic regions. But such an artificial procedure evident- ly cannot be considered justified because we obtain a number of objects i_ ' unrelated to any one of the classes. Returning to [6], it should be noted that the theory of a universal objec- tive classification constructed there in the considered case is inapplic- able since at its basis is the hypothesis of nonintersection of the char- acteristic regfons of classes. Thus, the classification in [6J can be re- garded as universal only in a limited class of fields of ineteorological elements. - In this article the author does not attempt to dispute the position that - in the atmoaphere there are some stable states around which the entire ensemble of temporal realizations of atmoapheric processes tends to be 46 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02108: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY grouped. However, the objectively established fact of the existence of some individual meteorological f ields not similar to one of the standard fields of determined classes forces us to seek a more flexible approach ta the objective classification problem than used now and before. Now we caill examine still another aspect of objective classification theory. In developing different meteorological classifications it is cus- - tomary to limit ourselve to the construction of "standard" fields of classes, assuming in this case that these standards are an exhaustive char- acteristic for all fields combined with them into a single class. Now we will consider the problem of the informational properties of the stan- dards, for clarity in description examining the case of two classes. We will denote these classes by A and B, and their standarda by xA and xB. The quantity of information I(y) concerning some field y EX consists of two ~arameters, to wit: the quantity of information abou* y included in xA and the quantity of information on y included in xB. I (y) _ / (v/xA) + f (.vlxR). The Iky/x) parameter is determined in the following way: 1 (J/�r) = H (x) +J/ (J) ---fl (X, Y), where c, H(x) f p(x) IogP (x) ttx is the statistical entropy of the x field and p(x) is probability density. = If y belongs to the region where p(x*A , y) is close to P(x*B, y), it is ob- vious that I(y/x ~)/I(y/.x*B) do not differ so much thar one of them can be neglected. Should we formally assign y to that class for which the p(x, y) value is maximum, we lose a considerable part of the information, which is scarcely justifiable. Thus, in the considered case of inter- sections of the characteristic regions the informational properties of the standards do not enable ue to solve the classification problem suc- cessfully if we employ traditional methods for breaking down the objects by classes. Another problem which must be dealt with in solving classification problems is the availability of a teaching sequence of sufficient length. V. N. Gli- venko has demonstrated the theorem of uniform convergence of the empirical distribution curve to the distribution function with an increase in the volume of the sample. It follows from this theorem that the length of the teaching sequence I ad adequate so that after ending the teaching process it will be possible to recognize the newly appearing objects is proportion- al to the dimensionality of criterial space and inversely proportional to the probability of classification error. Applicable to our data, in partic- ular, for the above-mentioned mean monthly 6H500 fields, it was establish- ed in [3] that lad N 6000, which is an order of magnitude greater than the archives of field data at our disposal. In order to reduce lad it is 47 FnR nFFTrTer TiCF nrrrv APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02108: CIA-RDP82-00850R000200090007-3 ~ . _ - _ _ _ necessary to reduce considerably the dimensionality ot criterial space. This can be achieved by the approach of expansion in natural orthogonal functions. However, since in order to obtain lad in our data we were forced to limit ourselves to a very small number of the first expansion terms, some perceiltage of the variabili[y (dispersion) of the fields re- mains not taken into account. In addition to what has been stated above, this circumstance indicates that salution of L�he classification problem, which in essence is a com- pression of the useful information, in the considered case a transforma- tion of this lnformation into a loss. It seems that this difficulty also cannot be overcome by any effective method if we remain within the framework of the traditional, that is, de- terministic theory of breakdown by classes. Accor3ingly, we will proceed to fo:mulation of a fundamentally new, probabilistic classification prob- lem: we will characterize each field xi E X not by an unambiguous assign- ment to the class A or B, but instead we will introduce the concept of the - degree of reliab ility of belonging DA(ai) and DB(xi), that is, we will say - that the xi field belongs to the class A with the degree of reliability DA(xi) and belongs to the class B with the degree of reliability Dg(xi). In the case of M classes each xi will be characterized by the vector D(xi) Dil,..., DiM 3, i= 1,..., M. Such an approach is far more universal a than the deterministic approach, and in addition, with such a representa- tion there is no informa.tion loss. In actuality, if we f ind all the degrees of reliability, the determined vector D(xi) will unambiguously characterize xi, that is, rigorously fix its position in the space X, not leavfng room for uncertainty. It is con- venient to uae a probabilistic characteristic as the degree of reliabil- ity. Thus, we will assume that in all X there are functions DA(a) and DE(x) which are the conditional probabilities that x belongs to class A or class B respectively. We note that such a probabilistic formulation of the prob- lem is a special case uf the deterministic approach, characterized by the fact that DA(x) and DB(x) only assume values equal to or close to 0 or 1 in objects of A and B, as occurs in the case of adequate compactness and a considerable distance between A and B. Now our problem is, using the objects xi E X appearing in the teaching pro- cess, to restore DAW and Dg(x) = 1- DA(x}. In order to approximate DA(x) we will use the adaptive algorithms approach, and specifically the poten- - tial functions method [1]. This method is attractive, first of all, because it is free of any a priori assumptions concerning the nature of the distrib- ution of objects, Lhat is, we need only assume that there is some non-zero probability of the :ppearance of objects xiE X; second, the method is at- Cractive due to its simplicity and the economy of computer time; third, the potential functions method allows successful synthesis with a Golovkin al- gorithm. Applying the Golovkin algorithm, we will seek the optimum set of standard objects and then we will "fix" them as the most typ3.ca1 represent- ativea of the classes, and after this, applying the potential functions 48 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02108: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY approach, we will ascertain the degree of reliability of whether each xi belongs to these standar.d objects. . TFius, aseuming that DA(x) is representable by a series in the form DA (x) cj?I (t), . ~ _ where (Pi(x) is some camplQte system of functions, we will use the recur- _ rent procedure fA 1 (x) =Jq (x) -f ~11 K ('CA+1, X), where fA(x) is the n-th approximation to DAx), K(xn+l, x) is a potential function. In general form K has the following form: _ C6 K )�i (_X) 'c(Y), (2) = where the followi..E, conditions are imposed on T i ~ -n : k,T--- L i t.ti'S; -10 n FE~;?'�-:~^,~�:~ �cnI 6 Fig. 2. Hydrodynamic and thermodynamic characteristics of atmosphere in _ zone of influence of typhoon "Carmen." 1) cloudless Trades zone, 2) north- ern periphery, 3) central cloud mass, 4) tail convective zone. KEY: 1. mb 2. m/sec 3. cal/g A11 the experimental data obtained cordance with the synoptic situatii sible to discriminate seveYul such 4. cm/sec 5. sec 6. cal/(cm2�sec) during this period were grouped in ac- 3ns in the polygon region. It was pos- situations: 1) Northern periphery of typhoon "Carmen"; 2) CentrRl cloud mass; 3) Tail convective zone. It ahould be noted that the proposed classification of synoptic situations reflects quite well the peculiarities of the state of the troposphere in the region of the aerological polygon. The distribution of the thermohydro- dynamic parameters of the atmosphere is qualitatively similar within each group, which makes it possible to consider the reaulta of the computations as objective characteristics of the corresponding synoptic eituations. The results of computations of the thermo- and hydrodynamic characteristics of the troposphere for the above-mentioned situations are presented in Fig. 2. The horizontal components of wind velocity u and v are determined in a Cartesian coordinate system with a reading origin situated at the center of Che polygon. The coordinate syatem is oriented in such a way that the OX- 56 FOR OFFICIAL USE ONLY f' . APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY axis is directed from the center of the typhoon. Here u and v are the radial and tangential components of wind velocity, directed along the OR and OY axes respectively. The positive direction of the tangential velocity com- ponent is counterclockwise. The positive sign for the vertical velocity component w corresponda to an upward direction, whereas relative vorticity ROT corresponds to counterclockwise rotation. Figure 2 also shows the thermodynamic characteristics of the atmosphere: latent heat LQ, static energy EL, and also their influxes (FLQ and FEL) per unit time in an air column with a baee 1 cm2 and an altitude 50 mb. _ Now we will examine the most significant characteristics of the troposphere atructure over the waters traver.sed by typhoon "Carmen" in its different zones. The aerological polygon of 10-11 August was to all intents and pur- poses aituated outside the zone of active influence of the typhoon (Table 1, Figure 1), on its northern pexiphery, at a distance of about 800 km from the center in a cloudless zone. As a result of typhoon movement on 12 Auguet the aerological polygon was within the cloud mass at a distance of about 400 km from the center. Then the typhoon moved and the aerological polygon occupied a position in its tail cloud system at distances from the _ � center of 800 km on 13 August and from 1,200 to 1,400 km on 14 August. - The greatest changes in the characteristics of the troposphere occurred _ with movement of the polygon from the northern periphery of the typhoon into its central cloud masa. The northern periphery of the typhoon, in all probability, was a slightly disturbed zone, subjected only to a limited in- fluence from the typhoon. This can be seen most clearly from the profile of the radial component of wind velocity, which does not have traits char- acteristic for zones of strong influence of typhoons (small influx in low- er troposphere, no upper tropospheric outflow). Other distinguishing char- acteristics of the considered zone are posi*_ive relative vorticity, the = presence of weak but noticeable ascending air flows, small heat influx ` and virtually zero moisture influx, on the average, in the entire thick- _ neae of the atmosphere. = ' The central cloud mass differs appreciably in its characteristics from the ~ considered situation. This difference was manifested most clearly in the following facts: there is a considerable change in the profile of wind radial velocity there is a great influx in the lower troposphere and an outflow in the upper troposphere; . there is a change in the stage of vorticity from positive on the north- ern periphery to negative (anticyclonic) within the cloud mass, especially at great altitudes; the descending flows in the entire thickness of the troposphere in- crease by a factor of 4-5; in comparison with the preceding situation the total tropospheric mois- ture content increases bq a factor of 1.4; 57 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 1'Vl\ there was a considerable warming of the troposphere a minimum in the profile of static energy (in the case of the northern periphery sit- uated in the layer 600-700 mb) was smoothed and the atatic energy in the entire air column increased. There was a subatantial change in the energy flux in a vertical column of . the tropoaphere. There was an increase in the influx of static energy FEL. It was distributed for thz most part in two layers� 950-700 and 700-300 mb. In the first layer the influx averaged about 2�10-~ cal/(cm2�sec), in the second approximately 5�10'3 cal/(cm2�sec), However, the influx of latent heat FLQ became negative and in the middle layers of the troposphere at- tained values 8�10-3 cal/(cm2�sec) (the integral moieture outflow from a vertical unit air column in this case attained 0.85 g/day). It should be noted that tt�e energy influx of 10'3 cal/(cm2�sec) in an air column with a height of 50 mb and with a base of 1 cm2 can heat (under the condition that energy is expended only on heating) air by 7�C per day. Thus, in the central cloud zone of the typhoan the tropesphere was intensively heated. Tha tail convective zone of the typhoon, in which the aerological polygon was present on 13 August, differed little with respect to hydrodynamic characteristics from the central cloud mass. As on the preceding day, there were descending vertical flows, but their velocitiea decreased by a factor of.4-5. Vorticity remained negative. The profile of radial wind velocity retains features characteristic of the central cloud mass, but the outflow in the upper troposphere substantially decreased. There wae a considerable change iz Lhe energy balance of the tropoephere. . Although the vertical distribution of EL and LQ virtually did not differ from the central cloud zone, the energy balance (that is, the influxea of heat and moisture into an air column) decrea.sed by a factor of 6-7. The total influx of all types of heat (static energy) was close to zero. - Table 2 Cobrdinates of Polygons and Observation Times (Cloudless Situation) - - ~ - KUQpitilnlbi $ C\1A Koor~,u~~r~ ~ - BPC\IA i10811C01(2, U ' 4itcno Ttuc.iol n ~rnHi� �ontiroiia, zpa~ 3 - ceFiTH6p� ~rpuneFi, zpao 3 CC}ITAUPA ( P I choe), ir cKoc), v - ~1 Z c. ut. 1 B. :1� 5 _ c. In. B. a. ~ 9 12 15,1 143.5 29 12 12,2 146,0 9 00 1511 143.5 29 18 12,2 146,7 pp 13,3 138,8 30 00 12,2 146,9 12 00 13,2 141,6 30 OG 12,1 146,?, 12 pp 12,8 141,3 30 12 11,5 141,4 12 12 15,8 143.5 2 12 12,3 145,5 29 OG 12,1 196,8 12 00 12,G 145.5 KEY: 1. Day in September 3. Polygon coordinates 5. E _ _ 2. Time (GMT), houra 4. N 58 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY Thus, an analysis of the thermodynamic characteristir_s of typhoon "Car- men" makea it possible to draw the following conciusiona: 1. The tail convective zone of a typhoon ie characterized by a coneiderable heat content, differing little from the heat content of the central cloud mass, but differing appreciably from the Iittle-diaturbed northern peri- phery of the typhoon. 2. The tail convective zone is relatively conservative in the sense that the influx of latent heat FLQ and eepecially the total energy balance of the influxes FEL into a vertical air column are not great. 3. Both in the central cloud mass and in the tail convective zone of the ty- phoon there is a considerable anticyclonic vorticity in the entire thick- ness of the troposphere. It is particularly great in the layer 400-200 mb. 4. The vertical velocities both on the northern periphery and in the tail convective zone of the typhoon its eastern periphery are great. In this way the periphery differs appreciably from the central cloud mass of the typhoon, in which the vertical velocities are considerable. Some thermodynamic characteristics are given in Figures 2 and 3 for a com- parison of the thernadynamic characteristics of typhoon "Carmen" with the characteristics of the troposphere under undisturbed conditions. These thermodynamic characteristics are for the cloudless Trades zone, whose frequenc,y of recurrence during the time of the expedition was rather great. In order to compute these characteristics we averaged the experimental data obtained at different times during the course of the expedition and _ characterizing a relatively extensive cloudless zone. This fact was es- tabliehed ueing satellite photographs, nephanalysis maps and meteorolog- ical evaluationa of cloud cover over the polygon. Table 2 gives the coor- dinates of the polygons and the observation times for such cases. It should be noted that in this case (in the abaence of typhoons near the polygon) the coordinate system is oriented in such a way that the OX-axis is directed to the east and the OY-axia is directed to the north, that is, u is the ~ zonal component of wind velocity, v is the meridtonal component. We note the moat characteristic features of the thermohydrodyriamic parameters rep- resented on the graphs. These include: 1) Weak positive vorticity in the boundary layer of the troposphere and a considerable negative vorticity at great altitudes. 2) Weak, virtually zero ascending flows in the entire thickness of the tro- posphere. 3) A virtually complete absence (on the average) of exchange of heat and ua isture. The influx of static energy into an air column is virtually equal to zero. Thus, as follows from the computations, the cloudless Trades zone is char- acterized by a considerable conaervatiam, that is, a constancy of the ver- tical profiles of the principal thermodynamic characteristics and small (on the average) heat influxes into an air column. 59 F(1R l1FFTCTAT TTCF l1TTT V APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 BIBLIOGRAPHY 1. Minina, L. S., Arabey, Ye. N., "Structure of the Troposphere Under Tropical Cyclogeneais Conditions," METEOROLOGIYA I GIDROLOGIYA (Meteor- ology and Hydrology), No 12, 1979. 2. Petrova, L. i., Nesterova, A. V., "Dynamic and Energy Characteristics ~ of the Tropoephere on the Periphery of Typhoons and in Zones of Tteir Maximum Frequency of Recurrence," TAYFUN-75 (Typhoon-75), Vol I, Len- - ingrad, Gidrometeoizdat, 1977. 60 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY UDC 551.521.3 ACCYJRACY IN DETERMINING ATMOSPHERIC OZONE CONTENT USING DATA FROM - MEASUREMENTS OF OUTGOING RADIATION = Moscow METEOROLOGIYA I GIDROLOGIYA in Russian No 3, Mar 80 pp 51-58 - [Article by Candidates of Physical and Mathematical Sciences Yu. M. Timo- = feyev and M. S. Biryulina and V. P. Kozlov, Leningrad State University, _ aubmitted for publication 15 January 19791 - Abstract: Computations of the informativeness and - = optimum conditions for measurements of outgoing thermal radiation in the 03 absorption band 9.6 m are given in the problem of indirect determination _ of the vertical profile and total content of ozone - in the earth's atmosphere. The article gives an analysis of the physical reasons for the low in- formativeness and accuracy of the indirect method relative to the q(p) profile. Different ways are - proposed for increasing the accuracy in determin- ing the characteristics of ozone content using - measurements from meteornlogical artificial earth - _ satellites. It is shown that heterodyne measurements of outgoing radiation (with a resolution 10'2-10-3 = cm 1) make possible a substantial increase in the accuracy of the indirect method in the upper sttato- _ sphere. - [TextJ Measurements of the'total content of ozone U03 and the vertical pro- ` file of its concentration in the atmosphere, already carried out for a - quite long time, have become particularly timely during recent years, in = particular in connection with the problem which is arising of possible changes in the thickness of the ozone layer as a result of change in the chemical composition of the atmosphere [10, 18]. Present-day requirements for obtaining global information on atmospreric characteriatics stimulated the development of special satellite methods for determining ozone content based on interpretation of data on outgo-ing _ radiation in different spectral ranges. In one of these methods use is made 61 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY of ineasurements of outgoing thermal radiation in the region of the 03 ab- sorption band 9.6 � m. Sources [7, 12, 16] give the results of a theoret- ical analysis of the accuracy of this indirect method. The meteorological artificial earth sate113.tes "Nimbus-3" and "Nimbus-4" were used in measur- ing outgoing radiation, making it possible to evaluate its real accuracy, and also to obtain a great volume of valuable information on the global field of ozane content [13]. After analyzing the principal results of these studies it is possible to come to the frllowing conclusions: _ a) Theoretical investigations and satellite experiments indicated that the - informativeness of the considered indirect method when using an average epectral resolutian (Q y= 2-5 cm 1) is small. In particular, with the pres- ent-day accuracy of ineasurements of outgoing radiation it is possible to determine only 1-2 independent parameters of the vertical structure of the ozone content profile (for example, 1-2 coefficients in the expansion of the profile in the r_haracteristic empirical base). As a result, it is not possible to obtain sufficiently reliable data on the vertical structure - of the 03 content profile in the atmosphere. b) The indirect method makes it possible to determine the total ozone con- tent with a relatively hi;h accuracy. In particular, both theoretical in- vestigationa and the itice,:pretation of satellite data indicate the pos- sibility of determining U03 with an accuracy of 6-10% relative to the true or mean U03 value. c) In not a aingle one of the studies with which we are acquainted was any attempt made to analyze the physical reasores for a low informativeness of the indirect method or to find wa7s to increase it. In particuiar, until - now there have been no computations of the optimum canditions for measure- ments of outgoing radiation in the absorption band 9.6 �-m. Experience has ' indicated that a correct solution of this problem in a number of cases will considerably inerease tiie accuracy of indirect methods [3]. In this study an attempt is made to fill the existing gaps. The mathematical basis for solution of inverse problems in the theory of transfer of thermal radiation is the well-known integral form of the trans- - fer equation [4]. After linearization of the equation the considered prob- lem is reduced to solution of the following operator equation of the first kind: 8l=AXZq, (1) _ [It is assumed that the atmospheric temperature profile is known on the basis of solution of the problem of theroal sounding of the atmosphere [4l]� ~ _ where b I is the vector of deviations of the measured (with different fre- quenciea yi, i= 1, 2,...,m) radiation Z( v i) from the radiation I( V i), correaponding to the mean profile of of the mixture ratio 03 - q(p), 62 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY Sq is the corresponding vector of deviation of the mixture ratio (at dif- ferent atmospheric levels p1,j - 1, 2,...,n), A is a matrix m x n, whose elements have a simple physical sense they characterize the sensitiv- ity of outgoing radiation in the i-th spectral range to variationa of the q(p) pref ile at the level with the preseure pi. n uF Fig. 1. Vertical temperature profile T and mean ozone profile q for which - Computationa were made (a), the kernels A(p) for the absolute variations of ozone concentration in three registry channels (b) and the kernels A'(p) for relative variations of the ozone concentration (c). 1) y= 1055 cm 1, ' 2) V =995cm1, 3)~) =1125cm1. KEY: 1) mb 2) 8/8 ' 3) erg/(cm2�sr�cm 1)/g/g Figure lb shows examples of computation of Ai(p) for three spectral ranges at the center and in the wings of the absorption band with a spectral reso- lution of the measurements Q,y = 5 cm'l. The computations were based on use of a atatiatical absorption ;nodel [9] and for the mean profiles q(p) and T(p) corresponding to the period May-September [17]. An analysis of the behavior of the Ai(p) curves shows that for all the con- sidered apectral intervals (we examined 20 intervals in the range from 965 to 1155 cm 1) the ma.ximum influence on outgoing radiation is exerted by absolute variations of ozone content at approximately one and the same level, situated in the troposphere, that is, svbstantially below thc, level where the maximum 03 content is situated (see Fig. la, which gives mean q(p) profile. A change in the optical density of the spectral :~lte_r- vals (with movement from the center of the band to its edges) leads only 63 - FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 ~ v y�iu c~c 2> 0 -A�~d~ e 3pt/(CM1CeKGO�Gti'l �AO0'30tAcHt�ctKtP-f4 41 ?OD ' ?SO T K t/t 3 ) APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-00850R040240090007-3 I to changes in the absolute values Ai(p), but not to a displacement of the curvea in the vertical coordinate. This leads to the conclusion that for the considered indirect method the information available for layers above and near the ozone content maximum is very low. A similar behavior of the _ Ai(p) curvea is also obaerved for other aeaeons. The Ai(p) curves shown in Fig. 1 characterize the sensitivity of I(v) to absolute variations of q(p). However, due to the strong variability of q(p) within the limits of the ozone laper (q(p) varies by two orders of magnitude) it is useful to analyze the behavioi. the Ai(p) curves cor- responding to an ~erational equation relative to the ~q/q values (it is easy to see that Ai(p) = q�Ai(p). Figure lc shows Ai(p) curves for these same spectrai intervals. It can be seen that although the Ai(p) curves are situated higher in the atmosphere, they are characterized by the above- mentioned peculiarities. Under conditions of a low information yield of the indirect method it is _ desirable to make computations of the optimum conditions for measuring outgoing radiation for the purpose of maximum use of the information con- tained in the I(v) spectrum. In this study we used the optimization meth- od proposed by one of the authors [13 and involving the optimum combining of the radiation c+f different spectral intervals in different measurement - channels. The choice of the optimum number of ineasurements (N) is accom- _ plished on the basis of the known criterion [2] ~k~/(y I~ (2) where /'\k is.the k-th eigenvalue of the covariation matrix of radiation KII = AKqqAT, O'I is the dispersion of random error, AT is a transposed matrix. Table 1 Eigenvalues ~k(erg/(cm2�sec�sr�cm 1))2 for Different Seasons k V-IX X-XII I-IV 1 0.4135�102 0.4375�102 0.1558�102 2 0.3694 0.2382 0.1565 3 0.3167�10'2 0.1480�10-2 0.1700�10-2 Table 1 gives the first three eigenvalues for three seasons when using co- . variation matrices of ozone content ICq.q from [17]. - It follows from the data in Table 1 that with the preseat-day accuracy in measurements ( O" =0.1-1.0 erg/(cm2�sec x sr�cm-1) with a spectral resolu- tion 6v= 5 cm"l) the crfter.ion (2) is satisfied by 1-2 eigenvaLues. 64 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY of, .�1 , 71 /65 1000 � IOSS 1010 11J5 o 961 1010 lPi010SJ 1,0 6) 0" , 96S 1D10 1055 1010 1135 1 ^ ~ S6S ~ 1010 f0~y0 fOSS n. . 96' 1010 1055 1070 1135 Fig. 2. Configuration of optimum plan for combining spectral intervals. a) for two-channel instrument (N = 2); b) for three-channel instr.ument (N = 3) (with a value of the ordinate equal to unity a channel is cut in t_o }he combining operation). Table 2 KEY: o, a0,1 3pe1(c.x..cen�cp.c.u-1) 2 c, =i,u~p2l (c.,,=.cen - �cp.c.s--1) . - . n ~ (1-7Qj7Q)�1M% - n (1--:Q17q)�100ti p ,04.~ 9, o ! 1009~ I nQ q�I00% ~ N M Nc3 ~ M N=3 N N II II M II t `1 N N II M II M II i 21 3 1a 1 .9 s 7 1 s 9 1 io 0.045 o 0 o 16,3 77 o 0 o 11,7 77 0,425 u u o 37,4 27 o u u 21i,6 27 s 7,0 a,c> 13,7 �10,5 16,5 3,8 ~1,0 4,3 17,4 18,7 10 10,2 17 5 6,3 28 2 19,8 4 28 61;,0 :31 5 10,1 9 1,8 J S 5,1 } 2 2 .5,9 25 6 25,8 25 6 11,9 i 9 20 , , , , ,1 , : , : , , 100 33,2 13,3 33,3 33.8 21,0 29.4 31,9 33,0 33,�1 21,1 200 Sfi,:i -A.3 57,2 59,2 26,2 ;52.5 53,9 :,i ,:3 27,0 400 24,5 18,5 31,6 40,6 19,4 12,1 25,5 2G,7 :3,1.7 21,7 800 4,7 2.6 9,1 17,6 32,3 l,'l IO,i 11,4 11.�S 31,6 A n A 1 -aL o 92,0 94,90 Jf,,16 9G,Eil =0,47% U 63,3 9~1,75 9:,,d~l or~.i,:! =0,55% G 7 U I p:~ I I t 0 1. mb 2. erg/(cm2�sec�sr�cm 1) It can be aeen in Fig. 2, which shows examples of the optimum plans with N= 2 and 3 for the season May-September~that the optimum measurement channels occupy rather broad and overlapping spectral regions and also 65 FOR OFFICIAL USE ONT.Y APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 include sectors in different parts of the band. Computations for other seasons indicated an absence of significant variations in the optimum plans. Now we will proceed to an analysis of the results of evaluations of errors in the indirect method for different measurement schemes. It should be _ noted that in a large number of studies (for example, see [12]) as the parameters characterizing the error in the method use is made of the relative errors in restoration, determined relative to the true or mean q(p) or U03 values. However, taking into account that the relative natur- al variations of the characteristics of ozone content are relatively small (12-60%) f17j, such an approach can create the impression of a rather high uccuracy of the indirect method. For example, in evaluating accuracy it is more correct to use the parameter ~ (v) Yo ip1= �y (p) characterizing the relative decrease in a priori uncertainty in knowledge _ of the profile q(p) after carrying out a satellite experiment ( O'q and ~ -ri ~ P4 04 D ~ b 0 ~ ~ P4 ~ 1. i r-~I ctl Cl O' .~G G a1 'z+ H D 00 O rI > a, ~ L ~ Cr' �rl "Cy 1J r+ 7 co ~ ~.I ~ ~ ~H r l C O U N p ~ ~ 0 ~ N ~ U T1 l.r Q cO i.+ cU . ~ ~ U u G1 �r l - p ~ ri ~ U w ~ co ~ G l + J ~ b0 Gl > ~f c0 H A; O V p rn 107 _ FOR (1FFT('. TAT. TTSF. nur.v APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 . _ 0 o C) ~ N N r--I ~O C w u o .,4 Ai ~ ~ ~ Q~N I I, I ,,..r U O as p O~7 r-i ~ O O~ Ctl r-i Q1 F-+ U i-+ O U O U G M 0 ~ O M u O r~ co II ~D ~ a0 ~1 N "O O 2~ ~ rn v ~ 0 00 CT% r- M O G �n O 0 SG ' \.O r-I ~ I C1. H C > r-~ c~~1 ri f~ rl O~ ON ~ C) 41 I r- 00 > rl MH r- 0 w ~4 O m r-i M O ~ 41 D, 00 II cn * � -rl r-I o Q) M � ~ --i ~ o~~ 41 co o c s o~~ ~ a H r-I oo � ~ i i ro ' N N v " ~ ~ > o 0 i c O u a a a~ y z o ~ b i i b 0 . > ~ > 9 " c~ � ri Cd 3 CY) ~ u 0 0 O F H O) d O U ' D+ d ~ O' m ~ O -H N 4.) cd cti 0 ~ U cC ~ ~ y O cd ~ ~ co V V 0 b0 41 H ~ ~ l U ~ bo W r w ~ 44 p P~G O A 10 8 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY extrapalation interval 'teX, so that Pj-),0 with 7 km-1 for Sc, As relate to clouds of the mixed phase. Measurements of spectral transmission in [18], made in the atmosphere in the presence of crystals In the form of needles, indicated that attenuation in the region of wavelengths X = 0.5-12 � m is close to neutral, although in some spectra there is a weak minimum at /^k= 9.2-10 � m. It is reported in - [38] that the attenuation of IR radiation with A = 10.6 � m in crystalline clouds is somewhat greater than in the case of visible radiation. On the basis of these data it can be concluded that in contrast to droplet clouds, attenuation of radiation in the range A = 0.5-12�.m by crystalline clouds is approximately identical. Therefore, Table 2 gives some idea concerning , the attenuation indices of crystalline clouds in the IR spectral region as well. 151 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY In order to compute the attenuation of radiation by crystalline clouds it is important to know whether it is possible to approximate crystals of dif- ferent configuration by ice spheres. At present there is no univerQal agree- _ ment on thia subject. Computations [15, 161 of the factora involved in the effectiveness of attenuation ICO(10.6) for radiation with a= 10.6~.cm by ice spheres, cylindera and platelets ahow that the influence of config- uration is manifeated for particlea with aizes ~o - 27Tr/,A < 6(where r is the radius of a aphere or cylinder and the half-thickness of the plate- let). The authors of [3, 36] compare the measured Kp(10.6) values with the . computed values for ice spheres and come to the conclusion that regardless of the configuration of the crystals, in computations of attenuation at 10.6 � m it ia possible to use a model of ice spheres of an equivalent sec- tion. In [15], on *he basis of a comparison of the results of their meas- urements of spectral transparency in crystalline fogs and computations for large water droplets, the conclusion is drawn that a model of spherical ; particles cannot be used in computing spectral attenuation. However, the authors of this study neglected the difference in the optical constants of water and ice. Accordingly, further investigation is required for an- swering the formulated question. ; Radiation Scattering by Ice Crystals One of the distinguishing characteristics of the scattering properties af j- _ crystalline clouds is manifested most clearly in observation of a number ' of optical phenomena accompanying the transmission of solar radiation _ through a veil of cirrostratus clouds (less frequently, cirrus, cirrocum- ulus and altocumulus). These include large and small halos, false suns, tangential and circumzenith ares, light pillars, etc. These phenomena have long attracted the attention of researchers and have been described in detail in [11, 30 and in a number of other studies]. All thesE pheno- mena are caused by the reflection of light rays from plane surfaces or by light refraction in a crystal, which occurs the same as in a prism [11]. An ice crystal, having the configuration of a regular hexagonal platelet or pillar, for the rays constitutes a prism with refracting angles 60 and 90�. A ray, entering a prism, emerges from it after double refraction, forming small and large halos in the direction of angles 22 and 46�. Due to light dispersinn during the transmission of rays in the prism the halos must be brightly colored (inner edge red, outer edge violet)., ` in contrast to the phenomena caused by reflection (for example, a hori- zontal circle). However, the color of a halo is frequently not pure or is totally absent. In this case there is a"white" halo. The reasons for this may be both oscillations or rotation of the crystals about their axes and the superposing of reflected, refracted and diffracted rays [11, 30]. All the enumerated phenomena arise only when clouds contain crystals and for observers serve as one of the methods for determining the phase state of cloud cover. Moreover, the appearance of most of them indicates that in clouds conditions sometimes arise for the predominant orientation of crystal axes [11]. The process of radiation scattering by particles of 152 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY any configuration is described by a transformation matrix (scattering matrix), by means of which it is possible to express the Stokes parameters - of the scattered wave through the Stokes parameters of an incident wave [19]. The Stokes parameters characterize the intensity and polarization of radiation. In a general case the scattering matrix has 16 components. Spherical particles of an optically inactive aubstance have a scattering matrix of the aimpleat type (the tranaformation matris containa only three independent parametera). The asphericity of particles exerts an influence on all the light scattering matrices. The predominant orientation of aspherical particles can cause linear and circular refraction effects in combination with dichroism, even if the particles are not optically active [19]. Due to the great mathematical difficulties no expression has been - derived for computing the components of the scattering matrix, but in [33] a study was made of the theoretical principles for their experimental de- termination. %7' iU I , P% ~ IN` ~ fA 20 0 II -20 -yo 0 60 9C 120 8� Fig. 2. Polarization of light scat- tered by crystals cylinders (1) and platelets (2) according to data in [35]. I1 F~t ~ 0 Fig. 1. Scattering indicatrices for radiation with a= 0.63 � m. 1) com- puted for prism and horizontal scat- tering plane, 2) same for vertical scattering plane; 3) measured in a crystalline medium (data from [1]). 153 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY Table 1 Phase Structure of Clouds in Different Regions of the Earth Cloud form Eastern zone of Temperate Arctic tropical Atlantic latitudea Cu droplet - St droplet, mix2d, mixed droplet Sc droplet droplet, mixed, mixed dropl et Ns mixed mixed Ac droplet droplet, mixed, mixed droplet As mixed, mixed, crystalline, droplet crystalline mixed Cs mixed, crystalline crystalline crystalline C1 crystalline crystalline crystalline Table 2 - Values of Attenuation Index in Clouds of Different Forms, km71 Cloud form Eastern zone of Temperate Arctic trop{cal Atlantic latitudes Sc 25 As 7-20 25 4 Cs, Ci sp 8.3 2.5 2.5 Crystalline fogs 0.5-1.2 2.3-32 Until recently the attention of researchers has been concentrated only on measurement of the components fll and f12, which characterize the intensity - and degree of polarization of scattered light. Tte angular dependence of fll (scattering indicatrix) in crystalline fogs [27], measured for the firat time, did not exhibit intensity maxima in the range of halo angles, which in [28] is attributed to an inadequate angular resolution of the in- etrument. Later maxima of different intensity were noted at scattering angles of about 22� and 46� [4, 8, 9, 34, 35], corresponding to the amall and large halo circles observed in nature, and also with e= 142� [35] 154 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY Bouguer halos [ll]. A halo with 6a;46�, also under natural conditions, was obaerved considerably less frequently. In additton, in these investi- gations greater scattering in the region of angles close to 90� was observed in comparieon with droplet clouds (for the first time in [12, 27]). However, the scattering indicatricea meaeured by different authors reveal a substantial scatter of values, especially in the region of lateral scatter- ing angles. An analysis of the data indicated that thia scatter cannot be attributed to a difference in the shape and size of particles. The reseon for the discrepancy in the eaperimental data was expleined in [13, 141, where it was demonetrated that the presence of droplets or ice spheres in the medium, together with crystals, exerts an influence on the angular dis- tribution of scattered radiation. With a clear identification of the phase state and discrimination of cases of a purely crystalline medium the scatter of scattering indicatrices for crystals of different shape and size in the range of angles 10-180� was 1-30X. In [1, 14] the authors also demonstrated that an increase in scattering by crystals in the region of lateral angles occurs due to a decrease in the region of angles 2< 6< 40�. These data were obtained when measuring scattering in a single plane. However, as dem- _ onetrated in [2, 8, 9, 36, 44], the predominant orientation of crystals - leads to an asymmetric distributi4n of the scattered radiation in different scattering planes, although in the range of angles 10-180� the difference _ in fll( B) for different planea does not exceed the discrepancy in the in- tenaity values measured in a single plane. Figure 1 shows the averaged - scattering indicatrix from [1, 13], and Figure 2 shows the degree of plane polarization of scattered light [35]. Xn) Hpucmonnu11r.cKUe _ . - - - - - ~ " 1 u6nann � V, o z d,o ~ ~q~ o s - O's ` 0 7 _ . ~ NcKyccmDerrNOn � 2� cp~da croetuaNyou , ~ ` ~'anena~reie p ~~-~-r3 06nara -ZO -10 0 ~ 10 t 'C Fig. 3. D('rf) values for clouds of different phase state according to data from [41]. 1) snowflakes, 2) snowflakes covered with hoarfrost, 3) fall- ing anow, 4) graupel, 5)-6) thawing particles, 7) freezing droplets of aspherical configuration. KEY: 1. Crystalline clouda 3. Drnplet clouds _ 2. Artificial medium of mixed phase 155 L�!1D AL'L`TrT AT ttOV nwn v APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 Measurements of the angular variation of other matrix components have been made only in [7], whose basic conclusion is that the matrix of light scat- tering by crystals is characteristic for asymmetric particles with a ran- dom orientation. In addition, it is noted that measurement of ite compon- ents can become a means for the remote atudy o� clouda, for example, study of the depolarization of acattered light. The conclusion that it is possible to identify the phase state of clouds on the basis of the depolarization ratio D(31) was drawn earlier in [20]. D(31) is determined as the ratio of the perpendicular and parallel compon- ents in the backscattering signal in the casp of irradiation of particles by plane-polaiized radiation. 1'he experiments revealed that for water droplets with a diameter up to 5.5 mm the depolarization ratio D(J7)< 0.03 [40], whereas in crystalline clouds and snowfalls from 0.5 to 1.0 [20, 41]. However, it is difficult to characterize clouds of differ- ent composition by any definite D(71) values. It is only possible to es- tablish some limits of these values in dependence on the phase state of cloucis at different temperatures, as is demonstrated by Fig. 3 from [41]. Thus, the singularity of the light scattering matrix for a crystalline me- - dium is caused for the most part by the configuration of the crystalline particles. Therefore, the use of the scattering indicatrices in computa- _ tions as the model of a crystalline cloud, as models of ice spheres [24, _ 39] or circular cylinders [31] does not give a satisfactory result. In [1] there was a study of light scattering by hexagonal ice prisms of finite length in the approximation of geometrical optics and diffraction and analytical expressions were derived for computing the scattering in- dicatrix. Figure 1 showa the computed scattering indicatrices in the hor- - izontal and verticaZ planes for prisms oriented in the horizontal plane - (the radiation beam is propagated in this same plane) [1]. The consider- ed case ma.kes iti possible to evaluate the maximum possible asymmetry of light scattering becauGe all other cases of orientation should lead to its decYeas2. The a.greement with the experimental data in [13] can be consfdered satisfactory. The discrepancy in the computed and experimental ~ [1, 421 indicatrices in the region of backscattering angles can be attrib- uted to the fact that in the computations no allowance was made for a aumber of effects, such as surface waves [19]. Thus, the scattering indicatrices shown in Fig. 1 can evidently character- ize the angular distribution of apparent radiation scattered by a crys- talline cloud with e,> 10�. Interrelationship Between Backscattering and Optical Thickness In interpreting the results of lidar sounding of clouds at different wave- lengths it is of interest to study the interrelationship between backscat- tering Crand the optical thickness 'G~ . Experimental investigations [3, 156 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02108: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY 13, 14] have shown that this interrelationship for ~k= 0.63~t.m, inde- pendently of the phase composition of the medium, is described by a re- gression equation with identical coefficients. For Z 0.63 < 0.4 the in- terrelatioriship between 0'it and optical thickneas can be considered linear. The deviation from linearity in the case of large 2 63is attributable to the increasing contribution of multiple scatterin$: The equality of the regrebsion coefficients means that the value of the backscattering in- dicatrix for ^A = 0.63 in a cryatalline medium is the same as in a droplet medium: fll(51 ='m 0.05 (f30Y). Measurementa [6, 42] for wave- lengths 0.57, 0.63 and 10.6 � m, the same as in [3, 141, did not re- veal C~rdifferences in droplet and crystalline fogs with identical 2~. Investigations were not made at other waveleagths. Thus, from the backscattering signal for A = 0.63, 0.57 and 10.6 � m it is evidently impossible to determine the phase composition of cluud cover. This requires that measurements of polarization, such as D(Tt), be made. Transmission and Scattering of Radiation by a Crystalline Cloud Crystalline clouds, the same as droplet clouds, constitute an optical medi- um which is quite dense in which single scattering is accompanied by mul- tiple scattering processea. The problem of the propagation of radiation in cloud media in which multiple ecattering occura and with allowance for ab- sorption by atmoapheric gases and aerosol has not been solved at all. In theoretical investigations it has been special methods for solving this problem which for the most part have been developed and the results of _ computations of the fluxes of solar radiation reflected, absorbed and - transmitted by layers of droplet clouds are given. Individual character- - istics of radiation transfer were also examined for crystalline clouds, a layer of ice spheres or circular cylinders being used as the model [23, 24, 31, 32, 391. It was noted earlier that these models cannot be re- garded as satisfactory for computations of the indicatrix of light scat- - tering by ice crystals. However, the results of some theoretical and ex- perimental studies give basis for assuming that a model of ice spheres can be employed in computing the characteristics of radiation transfer. For example, it was demonstrated in [32] that the computed values for radiation and transmission in the window 8.3-12.5 � m for layers of different thickness containing ice cylinders and spheres of equivalent cross section differ insignificantly. Measurements of the reflectivity of artificial crystalline fogs in the near-IR spectral region were made in [45]. The author of [24] demonstrated that the results of these measure- ments agree well with the camputed reflectivity values for spherical ice particles. The transmission and reflection of solar radiation by clouds is dependent both on the optical properties of droplets and crystals and on the char- acteristics of radiation absorption in and outside a cloud [21, 22, 261. The influence of absorption by atmospheric gases in the range a= 2.4-3.0 �m is illustrated by Fig. 4, which shows the solar radiation spectra 157 Vnv nvVTrrAT ttcV rnnv APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 for radiation transmitted by crystalline clouds Ci [21]. The solar spec- trum was measured in a clear sky at the same altitude (6.2 km) at which in another flight the base of Ci clouds with a thickness from 0.5 to 3 - - km was si[uated. 1 lDO U 1, 4 v /r ~ ~ i i3 2, 9 J 1,4 MKM ~A m Fig. 4. Transmission of solar radiation by the pure atmosphere at an aYti- tude of 6.2 km (1) and by Ci clouds thin (about 0.5 km), (2) intermedi- ate thickness (3) and dense (about 3 km) (4). Ie11.a'yDm/(CNr� c,) � rrn.y} 1 M 1000I1 100 � 1CG ',`l ~ ip . f I ~ i,0 i,S 1,4 ?,.1,?AMIr+I 2 Fig. 5. Reflectivity of solar radiation by dense Ci (1) and dense Cu (2). KEY: 1. Irefl , 1A W/ (cm2 � sr �tkm) 2. m The inf luence of absorption by atmospheric gases is also manifested in the spectra of reflectivity of solar radiation. There are reflectivity minima corresponding to the absorption bands of o.rygen A = 0.76 and 1.26~i.m, the strong bands of water vapor ~1, = 0.95, 1.13, 1.38, 1.47, 1.86, 1.9 and 2.58 �,m, as well as -the H20 and C02 bands with = 2.01, 2.06 and 2.68~t,m. In addition, there are minima corresponding to the water or icE absorption bands [21, 22]. - 158 FOR OFFICIAL USE ONLY LA APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY The latter is attributable to the fact that the absorption of radiation by particles decreasea scattering. Figure 5 ehows the averaged reflectiv- ity spec tra for cryatalline Ci and droplet Cu clouds from [21]. In aome intervals of wavPlengths the reflectivity of these clouds is different. The difference is maximum where the complex refractive indices of water and ice are most different. Thus, the measurement of reflectivity in the IR spectral region can be used in identifying the phase state of clouds. This conclusion was later confirmed by computations for plane-parallel layers of the medium con- taining spherical particles of water and ice [24, 29, 39]. _ Summary - Thus, the optical properties of crystalline clouds at present have not been = fully studied. Some components of the scattering matrix of visible radi- - ation, especially the scattering indicatrix with 0> 10� and the polariz- = ation of scattered radiation for angles 15-150; can be considered known. - An analyais of published investigations shows that the singularity of the - scattering properties of crystalline clouds, distinguishing them from droplet clouds, includea: ~ 1) manifestation of such optical phenomena as halos, light pillars, false suns, different ares, etc., caused by refraction, reflection or diffrac- - tion of Iight on ice crystals in the case of a random or predaminant ori- - entation; _ 2) an increase in the fraction of radiation scattered in lateral direc- tions in compar3son with drops; 3) the partial depolatization of backscattered radiation. ' All these effects are caused by a differeace iri the shape of the crystals from spherical and therefore cannot be precisely modeled by means of poly- disperse ice spheres ar cylinders. In computations of the scattering indi- - catrix for visible radiation it is possible to use the analytical expres- - sions proposed in [1]. These were derived for a model medium consisting of ice hexagonal prisms. In addition, the spectral crystall ine and droplet c able for the most gart to and ice. Aecordingly, for ity of crystalline clouds corresponding size. variation of reflectivity of IR radiation by louds is different. This is evidently attribut- the infZuence of the optical constants of water computing the spectral variation of reflectiv- it is possi'ble to use a uodel of ice spheres of The characterist3cs of the optical progerties of crystalline clouds given above give basis for remote determi.nation of the presence of the crystall- ine phase in elouds and the relative cnntent of cryst.als and droplets. 153 rnID AVL�T/�7A� rrOa n*nv APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-44850R000200094407-3 rvn vrri...Lru., v.... BIBLIOGRAPHY 1. Volkovitskiy, 0. A., Pavlova, L. N., Petrushin, A. G., "Light Scat- tering by Ice Crystals," IZV. AN SSSR, FIZIKA ATN:OSFERY I OKEANA (News of the USSR Academy of Sciencea, Physica of the Atmosphere and Oceun), Vol 16, No 2, 1980. 2. Volkovitskiy, 0. A., Pavlova, L. N., Snykov, V. P., "Asymmetry of the Scattering Properties of a Crystall ine Cloud Medium," IZV. AN SSSR, FIZIKA ATMOSFERY I OKEANA, Vol 11, No 7, 1975. 3. Volkovitskiv, 0. A., Nikiforova, N. K., Pavlova, L. N., Petrushin, - A. G., Snykov, V. P., "Model Investigations of the Optical and Aero- dynamic Properties of Ice Crystals," VOPROSY FIZIKI OBLAKOV (Problems in Cloud Physics), Gidrometeoizdat, Leningrad, 1978. 4. Dugin, V. P., Golubitskiy, B. M., Mirumyants, S. 0., Paramonov, P. I., Tantashev, M. V., "Experimental Investigat ions of the Optical Charac- - teristics of Artificial Ice Clouds," IZV. AN SSSR, FIZIKA ATMOSFERY I OKEANA, Vol 7, No 8, 1911. ~ - 5. Dugin, V. P., Volkovitskiy, 0. A., Maksimyuk, V. S., Mirumyants, S. - 0., Snykov, V. P., "Spectral Transmission of Artificial Crystalline Cloud Formations," IZV. AN SSSR, FIZIKA ATMOSFERY I OKEANA, Vo1 12, No 4, 1976. ' 6. Dugin, V. P., Mirumyants, S. 0., Pavlova, L. N., "Experimental Inves- tigations of Backscattering at Wavelengths 10.6 and 0.57 �.m by Artif- icial Cloud Formations," TRUDY IEM (Transactions of the Institute of Experimental Meteorology), No 13(58), 1976. 7. Dugin, V. P., Mirumyants, S. 0., "Scattering Matrices of the Light of Artificial Crystalline Clouds," IZV. AN SSSR, FIZIKA ATMOSFERY I OKEANA, Vol 12, No 9, 1976. 8. Dugin, V. P., Volkovitskiy, 0. A., Mirumyants, S. 0., Nikiforova, N. K., "Anisotropy of Light Scattering by Artificial Crystalline Cloud Formations." IZV. AN SSSR, FIZIKA ATMdSFERY I OKEANA, Vol 13, rro i, 1977. 9. Dugin, V. P., Maksimyuk, V. S., Mirumyants, S. 0., Nikiforova, N. K., "Anisotropy of Light Scattering by Artificial Cryatalline Cloud For- mations (Vertical Illumination of Medium)," TRUDY IEM, No 13(58), 1976. 10. Kosarev, A. D., Mazin, I. P., Nevzorov, A. N., Potemkin, V. G., Shug- ayev, V. F., "Comparison of Some Microphysical Characteristics of Clouds in Different Geographical Regions," VOPROSY FIZIKI OBLAKOV, Leningrad, 1978. 160 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2047102108: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY 11. Mindart, M., SVET I TSVET V PRIRODE (Light and Color in Nature), Mos- cow, Nauka, 1969. 12. Nikiforova, N. K., "Experimental Study of Light Scattering by Droplet and Crystalline Fogs," TRUDY IEM, No 10, 1970. 13. Pavlova, L. N., "Investigation of the Scattering of Visible Radiation in Cloud Media Containing Crystals," IV VSESOYUZNYY SIMPOZIUM PO RAS- PROSTRANENIYU LAZERNOGO IZLUCHENIYA V ATMOSFERE (Fourth All-Union Symposium on the Propagation of Laser Radiation in the Atmosphere), Summaries, Tomsk, 1977. 14. Pavlova, L. N., "Measurement of the Coefficients of Directed Light - Scattering in Droplet and Crystalline Cloud Media," IZV. AN SSSR, FIZIKA ATMOSFERY I OKEANA, Vol 14, No 9, 1978. 15. Petrushin, A. G., "Attenua.tion, Scattering and Absorption of Radia- tion at-10.6 m by a Model Cloud Medium Containing Ice Circular Cyl- inders of Infinite Length," TRUDY IEM, No 11(54), 1975. 16. Petrushin, A. G., "Attenuation of IR Radiation in Crystalline Cloud Media," TEZISY DOKLADOV NA IV VSESOYUZNOM SIMPOZIUME PO RASPROSTRAN- ENIYU LAZERNOGO IZLUCHENIYA V ATMOSFERE (8ummaries of Reports at 4th All-Union Symposium on the Propagation of Laser Radiation in the At- mosphere), Tomsk, 1977. 17. FIZIKA OBLAKOV (Cloud Physics), edited by A. Kh. Khrgian, Leningrad, Gidrometeoizdat, 1961. 18. Filippov, V. L., Ivanov, V. P., Makarov, A. S., "Variations of Aero- sol Attenuation of Radiation Under Conditions of an Ice Fog," TEZISY DOKLADOV NA IV VSESOYUZNOM SIMPOZIUME PO RASPROSTRANENIYU LAZERNOGO IZI,UCHENIYA V ATMOSFERE, Tomsk, 1977. 19. Van de Khyulst, G., RASSEYANIYE SVETA MALYMI CHASTITSAMI (Light Scat- tering by Small Particles), Moscow, IL, 1961. 20. Shupyatskiy, A. B., Shlyakhov, V. I., Kravets, V. V., Tyabotov, A. Ye., "Use of Laser Technology in Polarization Investigations of Meteorolog- ical Formations," METEOROi.OGIYA I GIDROLOGIYA (Meteorology and Hydre,l- ogy), No 2, 1967. _ 21. Blau, H. H., Espinola, R. P., Reifenstein, E., "Near-IR Scattering by Sunlit Terrestrial Clouds," APPL. OPTICS, Vol 5, No 4, 1966. 22. Blau, H. H., Espinola, R. P., "Spectral Property of Cloude from 2.5 to 3.5 Microns," APPL. OPTICS, Vol 7, No 10, 1968. 161 xnu nFFTrTer TrcV n7%nv APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-00850R040240090007-3 23. Dave, J. V., "Intensity and Polarization of the Radiation Emerging from a Plane-Parallel Atmosphere Containjng Monodispersed Aerosol.y," APPL. OPTTCS, VoJ 9, Na 12, 1970. 24. Hansen, J. E., Cheynay, H., "Theoretical Spectral Scattering of Ice Clouds in the Near Infrared," JGR, Vol 74, No 13, 1969. _ 25. Heymsf ield, A. J., Knollenberg, R., "Properties of Cirrus Generating Cell," J. ATMOS. SCI., Vol 29, No 7, 1972. 26. Hovis, W. A., Blaine, L. R., Forman, M. L., "Infrared Reflectance of High-Altitude Clouds," J. APPL. OPTICS, Vol 9, No 3, 1970. 27. Huffman, P. J., Thursby, W. R., "Light Scattering by Ice Crystals," J. ATMOS. SCi., Vol 26, No 5, 1969. 28. Huffman, P. J., "Polarization of Light Scattering by Ice Crystals," J. ATMOS. SCI., Vol 27, No 8, 1970. 29. Irvine, W. W., Follack, J. B., "Infrared Optical Properties of Water and Ice Spheres," ICA1tUS, Vol 8, 1968. - 30. Lenggenhager, K., "Eine Meteorologische Erscheinung. Z. Erklarung der - "unchten" Spektralfarben des Sonnenhalos von 22� Radius," UNIVERSUM, " Vol 30, No 3, 1975. 31. Liou, K. N., "Transfer of Solar Irradiance Through Cirrus Cloud Lay- ers," JGR, Vol 78, No 9, 1973. 32. Liou, K. N., "On the Radiative Properties of Cirrus in the Window Re- . , gion and Their Influence on Remote Sensing of the Atmosphere," J. ATMOS. SCI., Vol. 31, No 3, 1974. 33. Liou, K. N., "lfieory of the Scattering Phase-Matrix Determination ; for Ice Crystals," JOSA, Vol 65, No 2, 1975. 34. Liou, K. N., Baldwin, R., Kaser, T., "Preliminary Exgeriments on the Scattering of Polarized Laser Light by Ice Crystals," J. ATMOS. SCI., - Vol 33, No 3, 1976. 35. Morita, Y., "Scattering Cross Section of Freely Suspended Ice Crys- ; tals for Visible Light," TENKI, Vol 20, No 3, 1973. - 36. Nikiforova, N. K., Pavlova, L. N., Petrushin, A. G., Snykov, V. P., Volkovitskiy, 0. A., "Aerodynamic and Optical Properties of Ice Crystals," J. AEROSOL. SCI., Vol 8, No 3, 1977. 37. Ohtake, T., Huffman, P. J., "Visual Range in 7ce Fog," J. APPL. METEOROL., Vol 8, 1969. 162 FOF OFFICIAL USE OtdLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY 38. Pembrook, J. D., Gryvnak, D. A., Burch, D. B., "Attenuation of Solar Radiation by Natural Clouds at Several Wavelengths in the Visible and Infrared," JOSA, Vol 61, No 4, 1971. 39. Plass, G. N., Kattawar, G. W., "Radiative Transfer in Water and Ice Clouds in the Visible and Infrared Region," J. APPL. OPT., Vol 10, No 4, 1971. 40. Sassen, K., "Depolarization of Laser Light Backscattered by Artif- icial Clouds," J. APPL. METEOROL., Vol 13, No 8, 1974. 41. Sassen, K., "Laser Depolarization 'Bright Hand' from Melting Snow- flakes," NATURE, Vol 255, 1975. 42. Sassen, K., "Backscattering Cross Sections for Hydrometeors: Measure- ments at 6328 A," J. APPL. OPT., Vol 17, No 5, 1978. 43. Schaaf, I. W., Williams, D., "Optical Constants of Ice in the Infra- red," JOSA, Vol 63, I?o 6, 1973. 44. Thucnan, W. C., Brown, A. G., "Preliminary Studies of Intensity af Light Scattering by Water Fogs and Ice Fogs," SCIENCE, Vol 120, No 3128, 1954. - 45. Zander, R., "Additional Details on the Near-Infrared Reflectivity of Laboratory Ice Clouds," JGR, Vol 73, No 20, 1968. 163 Fnu nFFTrrer TiCV l1TTT V APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 SzXTIETH BIRTHDAY OF SOLOMON MOISEYEVICH SHMETER Moacow METEOROi,OGIYA I GIDROLOGIYA in Russian No 3, Mar 80 p 122 [Article by personnel of the Central Aerological Observatory] [Text] Solomon Moiseyevich Shmeter, a well-known scientist in the field of atmospheric physics, a Doctor of Physical and Mathematical Sciences, _ head of the Division of Cloud Physics and Atmospheric Dynamics of the Central Aerological Observatory, marked his 60th birthday on 10 March _ 1980. Solomon Moiseyevich came to meteorology af�er graduating from the Mechanics- Mathematics Faculty of Moscow State University and then Khar'kov Hydro- ~ meteorological Institute. During the period 1944-1947 he worked at the Yak- utsk Geophysical Observatory, first as an aerological engineer and then as head of the aerology division. In 1947 he undertook graduate studies at the 164 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY ' Central Aerological Observatory. S. M. Shmeter was the first graduate stu- dent at the obaervatory and after then all his scientific activity over a period of 32 yeare hae been closely aseociated with the Central Aerological - Observatory. S. M. Shmeter dealt with a very broad range of problems during these years. Original investigations of the chemical composition of clouds were the sub- ject of his Candidate's dissertation, which he successfully defended in 1951. During 1951-1957 he carried out a series of important theoretical and experimental studiea devoted to an improvement in aerological measure- mente. Since 1957 the scientific activity of Solomon Moiaeyevich for the most part hae been concentrated on experimental investigations of atmo- ~ spheric turbulence and cloud physics. His !nany years of investigations of the dynamics and mesostructure of convective clouds later served as a ba- sis for his doctoral disaertation, successfully defended in 1968. The book by S. M. Shmeter, entitled FIZIKA KONVERTIVNYYE OBLAKOV (Physics of Con- - vective Clouds), received.broad international recognition. S. M. Shmeter repeatedly headed flight expeditions and participated in flights carried out under highly complex meteorological conditions. Dur- ing theae flight expeditions he created a new method for aircraft inves- tigationa of the atmosphere and obtained unique materials substantially enriching our concepts concerning turbulence, convection, and especially - the physics of convective clouda. The results of the studies of S. M. Shmeter have been published in more than 60 scientific articles. He is the author and coauthor of six mono- graphs on cloud physics and atmospheric turbulence. Three of his mono- graphs have been translated and published abroad. The rich experience of thia working aerologiat and At the same time hia excellent masL�ery of - theory assisted him in creating a text for uae in the aerology course (1965), he being the coauthor. Solomon Moiseyevich devotes great attention to teaching work. Under his , direction several Candidate's dissertations were prepared on different aepects of atmospheric physics. S. M. Shmeter is engaged in much public and political work. A member of the CPSU aince 1945, he devotes much time and energy to the indoctrina- tion of young scientists. Now S. M. Shmeter is at the height of his creative forces and we wish him good health and retention of the same vigorous activity and further creative successes. 165 FOR OFFTrrAT, T1GF nNr.v APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR ur'r'lulAL uoL unLY A'1' THE USSR STATE COP4fITTEE ON HYDROMETEOROLOGY AND ENVIRONMENTAL MONITORING Moscow METEOROLOGIYA I GIDROLOGIYA in Russian No 3, Mar 80 pp 122-123 [Article by V. N. Drozdov and V. M. Voloshchuk] _ [Text] A regular session of the Scientific Council of the State Committee on Hydrometeorology and Environmental Monitoring on the problem "Study of - the Oceans and Seas" was held on 28 November 1979. A report entitled "Preliminary Results of an Interdepartmental Expedition ~ - _ in the Baltic Sea and Prospects for its Hydrometeorological Investigation" was presented by I. N. Davidan. In his report he noted that the Leningrad Division of the State Oceanographic Institute and the Northwestern Admin- istration of the Hydrometeorological Service, within the framework of the Commission "Study of Hydrology and Water Contamination of the Baltic Sea," - during 1978-1979 carried out much work for bringing together the efforts of different departments in the field of organization and implementation of comprehensive investigations of the Baltic Sea, a generalization of the reaults of two interdepartmental eapeditions. The expeditions were carried out to study the most timely problems of the Baltic Sea (processes on a synoptic scale, variability and spatial nonuniformity of contamination fields) in accordance with programs developed at the Leningrad Division _ of the State Oceanographic Institute. ; Due to the close contacts between the Leningrad Division of the State Oceanographic Institute, the Northwestern Administration of the Hydrometeor- ological Service and the Arctic and Antarctic Scientific Research Institute, the following organizations became involved in the expeditions: Baltic Sci- entific Research Institute of Fisheries, Central Scientific Research Inst- tute of Geodesy, Aerial Mapping and Cartography, All-Union Scientific Re- search Institute o� Ma.rine Geology and the TEF and GS DKBF [expansions un- - known]. The expeditions were carried out successfully and made it possible to obtain new, original results, important for study of the Baltic Sea. Then the speaker discussed in detail the results obtained from the observa- tions and those preliminary conclusions which have been drawn. In conclu- sion he analyzed the principal directions in hydrometeorological investiga- tions in the Baltic Sea. 166 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY The Council noted that the results obtained from the interdepartmental ex- peditions in the Baltic Sea and also the results of intradepartmental eaped- itions in Tallinn Gulf, in which several shipa and aircraft participated eimultaneously, make it poseible to coneider such an organization of field investigationa to be the most rational in light of new tasks for preserva- tion of the environment. In its resolution the Council approved the organization and implementation of auch work and recommended that the Leningrad viviaion of the State Ocean- ographic Institute prepare materials on the organization and implementation of such expeditions for extending the experience to other basins. The coun- cil also decided to consider an interdepartmental program for investigations _ of the Baltic Sea as a unified geographical feature over the extended period - of 15-20 yeara (Project "Baltika"). A representative of the State Oceanographic Institute, A. N. Ovsyannikov and a representative of the Arctic and Antarctic Scientific Research,Institute, A. T. Bozhkov reported on the second problem "Status and Prospects of De- - velopment of the Sea Coastal Network." After hearing and discussing the reports, the Council noted that the marine administrations of the Hydrometeorological Service are devoting great at- - tention to supporting the normal work of the marine hydrometeorological net- work. The quality of observations in the marine network on the whole is good. The State Oceanographic Institute and the Arctic and Ar.tarctic Scientific Reaearch Institute are carrying out a systehiatic checking of operation of the marine network and are rendering the necessary assistance. During 1980 the State Oceanographic Institute, in collaboration with the Arctic and Antarctic Scientific Research Institute, Far Eastern Scientific Research Hydrometeorological Institute and hydrometeorological observator- ies of marine administrations of the Hydrometeorological Service will carry out scientific research work on the theme "Development of a Scheme for the Distribution of the Marine Coastal Hydrumeteorological Network in the Seas of the USSR," which will make possible a scientific validation of the prospects for development of the marine network. V. N. Dro zdov An unofficial conference of experts on statistical planning of the Precip- itation Enhancement Project was held at the State Committee on Hydrometeor- ology during the period 29 October-2 November 1979. It was devoted to the International Precipitation Enhancement Project being carried out in the basin of the Douro River (Spain). The scientific objectives of this project were defined by the 28th Seseion of the Executive Committee WMO in 1976. There have already been a series of unofficial conferences of experts at 167 T/~T I~TTTnT JT n~~ t� APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 which there was a discussian of different aspects of the Precipitation En- hancement Project. The materials of these conferences were given in dif- ferent WMO bulletins and reparts (Weather Modification Programme. Precip- itation Gnhancement Project). - At the conference of experts on statistical planning of the Precipitation ' Enhancement Project there was discussion of the present status of work on the project, the specialists analyzed the aspects of statistical planning and evaluation of fielci experiments with artificial modification for the purpose of obtaining additional precipitation and formulated recommenda- tions for further preliminary investigations which must be made for the purposea of the Precipitation Enhancement Project in the Douro basin. The Soviet delegation familiarized the foreign specialists with the prob- lems involved in the planning of field experiments for enhancing precipi- tation in the steppe region of the Ukraine and in the basin of Lake Sevan, with the results af experiments already carried out, and also with plans for further improvement of statistical evaluations of the results of field experiments being carried out at the present time in the USSR under the direction of Academician A. N. Knlmogorov. Foreign specialists presented concise information on studies discussed at the Sixth Conference on the Use of Probabilistic and Statistical Methods in the Atmospheric Sciences held by the American Meteorological Society during the period 9-12 October 1979 in Canada. The conference materials will be published in the next WMO bu'Lletin (Weather Modification Programme. Precipitation Enhancement Pro- ject). V. M. Voloshchuk 168 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY CONFERENCES, MEETINGS AND SEMINARS Moacow METEOROLOGIYA I GIDROLOGIYA in Russian No 3, Mar 80 pp 124-128 [Article by L. S. Speranskiy, K. Ya. Vinnikov, N. F. Dement'yev and V. D. - Komarovj _ [Text] An All-Union Conference on "Hydrodynamic and Statistical Methods for Local Weather Forecasting and Problems in Mesometeorology" was held in Novosibirsk during t%e period 15-19 October 1979. '1'he conference was attended by representatives of 15 acientific institutes and colleges in the country. A total of 36 reports were presented. A large group of repoxts was devoted to hydrodynamic methods for local weather forecaeting. Here, in particular, we should mention a report by Academician G. I. Marchuk, V. V. Penenko, A. Ye. Aloyan and G. L. Lazri- yev (Computation Center Siberian Department USSR Academy of Sciences) on nwnerical modeling of the microclimate of cities, with anthropogenic fac- turs taken into account. Individual aspecte of thie problem were examined in reports by A. Ye. Aloyan, D. L. Iordanov and V. V. Penenko, entitled - "Numerical Modeling of Tranaport of Passive Admixtures," A. Ye. Aloyan, G. L. Lazriyev and V. V. Abramenko, "Methods for Taking the Turbulent En- ergy Equation into Account," and A. Ye. Aloyan and G. V. Isayev, "Propaga- tion of an Admixture in the Surface Layer from Fixed and Moving Sources." A series of reports presented by specialists of the USSR Hydrometeorolog- ical Center was devoted to studies for creating a local forecasting model on the basis of deep convection equations. A numerical algorithm for in- tegrating a system of deep convection equations and investigating the properties of this system was diacuased by V. Z. Kisel'nikov, Ye. M. Pekelis and D. Ya. PreBSman. A report by M. Alautdinov and N. P. Vel'tishch- ev presented a method for computing radiation heat influxes and gave the = results of modeling of the diurnal variation cf ineteorological elements under moist deep convection conditions. A. A. Zhelnin and A. A. Bregman told of an algorithm for the parameterization of microphysical processes related to the formation of clouds and precipitation. V. M. Losev reported on computation of stationary mesoscale temperature and wind fields on the baeis of synoptic information using linearized deep coavection equations. - 169 FnD nL`L`rnTAT ifc`c nIkrrv APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 . ~ . _ _ P. Yu. Pushistov, L. S. Speranskiy, T. A. Toloknova and Ye. G. Sementsova (West Siberian Regional Scientific Research Hydrometeorological Institute) reported on a model of local forecxsting for a lowland territory and pre- sented some results of forecasts made on the basis of this model. V. K. Arguchintsev (Irkutsk State University) told about a two-dimensional model `l for investigating breeze effects, constructed without the simplification of quasistatics. B. G. Vager (Leningrad Institute of Construction Engin- eers) reported on allowance for horizontal diffusion in a model of the eurface layer intended for calculation of the heat and moisture transfer reglme in the preaence of inhomogeneities of a small ].inefir scale (canal in the desert, flight atrip). His second report was devoted to experience in the use of cplines in boundary layer problems. M. S. Akhmetov (Pskov State University) told of some aspects of interaction between open pit mines and the atmosphere. A number of reports were devoted to an investigation of convection and tur- bulence. R. S. Pastushkov and S. A. Vladimirov told about numerical model- ing of the convective cloud cover accompanying a squall. G. G. Goral' and 0. I. Chapovskaya (High-Mountain Geophysical Institute) reported on exper- imental investigations of convective pracesses accompanied by squalls. The synoptic and thermodynamic conditions determining the localization and intensity of convection in the Northern Caucasus were eaamined by L. M. Fedchenko and V. A. Belentsova (High-Mountain Geophysical Institute). The report of A. N. Koval`chuk and T. N. Terskovoy (High-Mountain Geophysical Institute) was devoted to an investigation of the dynamics of cumulonimbus clouds and the conditions for the formation of well-developed hail cells. An investigation of a model of a cloud ensemble on the basis of GATE data was preaented by A. I. ral'kovich (USSR Hydrometeorological Center). The subject of two reports by P. Yu. Pushistov, V. M. Mal'bakhov and S. M. Kononenko was the numerical modeling of penetrating convectton in the boundary layer dnd an inves~igation of two-level convective ensembles and the mechanism of formation of cumulus clouds. A repart by A. S. Gavrilov (Leningrad Hydrometeorological Institute) was devoted to an investigation of the convective boundary layer using a dif- ferent tool: systems of Fridman-Keller moment equations. In this report, as in the communication of the above-mentioned authors, great attention was devoted to the problem of antigradient heat transfer. A. G. Tarnopol'- skiy (State Oceanographic Institute) and V. A. Shnaydman (Odessa Hydro- meteorological Institute) made a comparative investigation of the charac- teristics of turbulence obtained employing two closing methods: on the basis of an equation for the rate of dissipation of turbulent energy and using the generalized Karman hypothesis. An investigation of the struc- ture of vertical movements in the boundary layer on the basis of experi- mental data was the subject of a report by Yu. P. Perevedentsev and Yu. G. Khabutdinov (Kiev State University). The reaults of experimental in- vestigationa and methods for computing turbulence in the layer 50-500 m were reported by V. N. Barakhtin and E. A. Morozov (West Siberian Regional Scientific Research Hydrometeorological Institute). 170 FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02108: CIA-RDP82-00850R000200090007-3 FGR OFFICIAL USE ONLY A report by B. N. Sesgeyev (Central Aerological Observatory), tell ing about a two-dimensional model of an atmospheric front, devoted great at- tention to the parameterization of phase transitions and the formation of precipitation. A model of formation of a frontal zone, taking into ac- count nonadiabatic factors, was examined by B. Ya. Kutsenko (Central Aero- logical Observatory). A report by Z. N. Kogan and N. Z. Pinus (Central Aerological Observatory) was devoted to an evaluation of the budget of the principal types of energy in a middle-latitude cyclone. The subject of a aecond communication by Kogan was a atudy of the kinetic character- istics of a cyclone on the basis of radiosonde data. Inveatigations of the interaction between the boundary layer and the free atmosphere for processes on a subsynoptic scale were reported by V. A. Shlychkov and P. Yu. Pushistov. T~ao reports were devoted to problems involved in flow around barriers: "Lee Wave Resistance Exerted by a Real Mnuntain System" V. N. Kozhevaikov, N. N. Zidlev and N. N. Pertsev (Moscow State University), "Numerical Mod- eling of the Problem of Flaw Around Obstacles by an Incompressible Fluid" I. G. Granberg (Institute of Physics of the Atmosphere). A number of reports examined models of prediction of local phenomena on the basis of statistical methods. The methods for ordering and sampling of parameters in local alternative models were investigated in a report by T. A. Anikina, G. M. Vinogradova and L. N. Romanov (West Siberian Regional Scientific Research Hydrometeor- ological Institute). One atatistical forecasting method, based on the ex- prinHlon of the FieldH of ineteoroloAical elemente into serieA in u syAtem of linezirly independent Eunctiona, was discusaed by V. V. Kostyukov (West Siberian Regional Scientific Research Hydrometeorological Institute). A _ model of prediction of thP frequency of recurrence of breezes was examin- _ ed in a report by E. A. Burman and F. Ya. Stupina (Odessa Hydrometeorolog- _ ical Institute). The following reports of specialists of the West Siberian Regional Scientific Research Hydrometeorological Institute were also de- voted to the prediction of individual weather phenomena: M. Ya. Kogan and L. N. Romanov "Quasilinear Models for Prediction of Diurnal Temperature Variation," G. G. Polyakov and Z. V. Torbina "Prediction of Precipita- tion Using a Piecewise-Linear Model," D. I. Zenkevich "Prediction of Visibility in Snowfal.ls," I. P. Prokop'yev "Prediction of Signif icant Snowfalls in the Southeastern Part of Western Siberia." " L. S. Speranskiy A bilateral Soviet-American symposium on the modeling of climate, climatic changes and the statistical processing of climatic data was held at Tbil- isi during the period 15-22 Octaber 1979. - It was the sixth of a series of symposia held during 1976-1979 within the framework of project 02.08.11 "Effect of Changes in the Heat Balance ,-)n Climate" of Working Group VIII ("Effect of Environmental Change on Climate"; 171 nnn nnnrnrAT y.r.,. n+.... APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-44850R000200094407-3 rVn vrri%..Lnu - of the Soviet-American Commiasion on Scientific Cooperation in the Field - of Preservation of the Environment. The project is headed on the Soviet side by Corresponding Member USSR Acad- emy of Sciences M. I. Budyko, and on the American side Professor W. L. Gates. Well-known American scientists participated in the symposium: W. L. Gates, E. N. Lorenz, R. D. Sess, R. M. Chervin, S. Esbensen, M. Halem, I. M. Held, R. D. Jenne, R. Madden, G. R. North, A. Robock, M. J. Suarez, J. M. Wal- lace, E. R. ;.eiter and E. W. Bierley. On the Soviet side there were about 30 specialists, for the most part from the inatitutes of the State Com- mittee on Hydrometeorology and the USSR Academy of Sciences. The symposium program included 37 reports, including 20 reports of Soviet - - scientists and 16 reports of scientists of the United States delegation. = One report was presented jointly by the scientists of tne USSR and the United States. The most general problems were considered in the reports of M. I. Budyko and E. N. Lorenz. The report of M. I. Budyko contained a synthesis of modern concepts con- cerning impending changes in global climate and the biosphere under the _ influence of anthropogenic change in the content of carbon dioxide in the _ atmosphere. The report expressed the assurance that already at the pres- ent time it is possible to prepare forecasts of climatic change necessary for working up long-range economic development plans. The problem of the predictability of behavior of a climatic system was discussed in a report by E N. Lorenz. Touching upon the problems involv- ed in prediction of climatic changes in relation to the increase in con- tent of atmospheric C02, E. N. Lorenz surmises that this is precisely the type of climatic problem whose solution will be most feasible by ap- plication of the numerical modeling method. Empirical evaluations of recent changes in the thermal and ice regimes of the northern hemisphere were examined in ajoint report of specialists of four institutes of the StaCe Committee on Hydrometeorology: K. Ya. Vin- nikov, G.V. Gruza, V. F. Zakharov, N. P. IGovyneva, K. M. Lugina and E. Ya. Ran'kova. It was demonstrated in the report that beginning from the middle 1960's there has been a warming of the northern hemisphere with a mean intensity 0.1-0.2� C/10 years. This warming is accompanied by a sub- stantial reduction in the total ice content in the Arctic Ocean. In a report by A. Robock, on the basis of an analysis af empirical data characterizing the change in the mean temperature of the northern hemi- - sphere during the last 400 years, carried out using a numerical zonal model of a climatic system with parameterized dynamics, the author - 172 - FOR OFFICIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2047102108: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY refuted the hypothesis that "solar activity" is one of the imporctant fac- tors in change in global climate. Empirical data on recent changea in climate and the results of a statis- tical analysis and generalization were also presented in the reporta of E. R. Reiter, K. Ya. Vinnikov and P. Ya. Groysman, G. V. Gruza and E. Ya.Ran'kova, R. D. Madden, 0. A. Drozdov, J. M. Wallace and D. S. Gutsler, Ye. M. Dobryshman, M. I. Fortus and Ya. M. Kheyfets, N. V. Koby- sheva and L. P. Naumova, L. S. Petrov, et al., V. N. Adamenko, Yu. L. Rauner. The climate modeling problem was discussed in more than ten reports. Eval- uations of the influence of an increase in the content o� C02 in the at- mosphere on the climatic regime, obtained using the Oregon University model, were presented in a report by W. L. Gates. Although the model in- cludes a realistic description of the topography of the earth's surface and the seasonal variation cf solar radiation, these evaluations cannot be used for predicting c'limatic changes, since in this investigation the temperature of the ocean surface is considered stipulated and does not change with an effect on the climatic system. 'Che necessity for allowance for the influence of thermal inertia of the ocean on the change in air temperature accompanying an anthropogenic in- crease in the atmospheric content of C02 was demonstrated in a report by R. D. Sess. In another report by this author it was stated that the pro- cess of global warming can be accelerated as a result of the increasing anthropogenic discharge of the compounds C0, NOx, CH4, forming from fuel combustion. In a report read by V. P. Dymnikov, Academician G. I. Marchuk and his col- leagues gave detailed information cuncerning the status of studies on con- struction of models in the theory of climate at the Computation Center Siberian Department USSR Academy of Sciences. I. V. Trosnikov told about a model of general circulation of the atmosphere developed at the USSR Hydrometeo rological Center. It was demonstrated in a report by V. P. Meleshko and R. T. Vezerold that allowance for the real geographical distribution of cloud cover exerts a substantial influence on the local characteristics of the climatic and circulatory regime of the atmosphere. Important problems in study of the laws of behavior and sensitivity of a global climatic system by means of models were also examined in reports by Corresponding Member USSR Academy of Sciences G. S. Golitsyn, G. North, M. Suarez, I. L. ICarol' and I. M. Held. The reports of V. A. Aleksandrov and V. Ya. Sergin, R. M. Chervin and M. Halem were devoted to a study of the statistical characteristica of the meteorological regime, reproducible using models of general circulation 173 'VAD /1VL�TIIT AT iinn /%%n lr APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007/02148: CIA-RDP82-44850R000200094407-3 of the atmosphere. The studies of these authors give additional informa- tion on how realistic climatic models are and make possible a more cor- rect formulation of numerical experiments for evaluating the sensitivity of cl.imate to eaternAl effects. Without dwelling on the other reports (the materials contained in them have been partially published), we note a paper by I. M. Held, which was discussed from the floor, telling about a new study by Manabe and Weth- erhold which gave an evaluation of the influence of doubling and quadrupl- ing of atmospheric C02 content on the mean annual climatic regime over the continents and oceans separa*_ely. According to these evaluations, the warm- ing of climate can be accompanied by a substantial change of moistening conditions i^ agriculturally important regions of the continents. The participants in the symposium expressed the opinion that it is neces- sary to summarize the results of the series of Soviet-American symposia held during the period 1976-1979 at a conference of Soviet and American experts on the problem of the influence of C02 on climate. Upon completion of the symposium this proposal was adopted at a session of the Eighth Working Group of the Soviet-American Commission. R. Ya. Vinnikov The Tenth Conference of the Danubian Countries on Hydrological Forecasts was held in Vienna during the period 11-14 September 1979. It was attend- - ed by delegations from Austria, Bulgaria, Hungary, USSR, West Germany, Czechoslovakia and Yugoslavia, as weZl as two hydrologists from Switzerland and one each from East Germany and the Netherlands. Three international or- ganizations were represented at the conference: Danubian Commission, World Meteorol.ogical Organizatio.L and International Institute of Systems Analy- sis. The total number of participants was 134. The working languages were Russian and German. - In evaluating the organization of the conference work we should note the high level of its preparation and implementation. Before the conference began two volumes of mater.ials were published which contained national _ reports on the methods used in each country for predicting the water - and ice regimes of the Danube (one volume) and the reports presented at the conference (second volume). Also included were author's simmaries of reports and detailed programs of the scientific sessions and field trips. Each conferee received these materials. We should also note the exceptionally high quality of simultaneous translation. All this favored the complete success of the conference. Thirty-four reports were presented and discussed at six scientific ses- sions. Most of these were devoted to the following problems and aspecte of hydrological forecasts for the Danube and its tributaries: short-range 174 FOR OFFIGIAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2047102108: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY forecasting of levels and discharges, long-range forecasting of levels and discharges, forecasting of ice phenomena, forecasting of channel processea and runoff of sediments, evaluation of the accuracy and econom- ic effectiveneas of forecasts. The greatest number of reports (15) dealt with methods for the short-range forecasting of water discharges and levels in the Danube and its tributar- y ies. The basis for the methods is mathematical models of the channel phase of runoff, diffe:ing in complexity, and models known under the name "pre- cipitation-runoff." Among the first of these models attention was centered on an adaptive lin- ear model of channel runoff whose parameters vary continuously in de- pendence on the observed error in forecasts on the preceding days. Pre- computations of river discharges and levels using an electronic computer on the basis of such a"self-adjusting" forecasting method are probably of considerable practical interest. The model was developed by the Hungar- ian hydrologist Sh. Ambrush for the reach of the Danube between Dunafoldvar and Baja, but it naturally has quite general significance. In the fore- casting methods which were discussed in the reports of Kresser, Gutknecht and Dreher (Austria), Engel and Mendel, Schiller and Tissa (West Germany) and Stefanova (Bulgaria) the emphasis was on computations of the travel- time curve of water on the basis of models of a cascade of linear and nonlinear capacitances. In general, these models are quite well known. Nevertheless, the reports were of considerable interest because they fam- iliarized the conferees with some modifications of these models and pro- - cedures for determining their parameters. It is also important that river reaches with a cascade of reservoirs were among the investigated objects. - "Precipitation-runoff" models were discussed in reports by Lutz, Koch and ~ Rosemann (West Germany), Lukachova (Czechosluvakia), Georgevic (Yugoslav- ia) and Nachtnebel (Austria). In scientific respects the most interesting - reports were an investigation of the influence exerted by the intensity of precipitation, runoff losses and some other factors on the form of the - unit hydrograph (report by Lutz) and a nonlinear model of the "runoff- precipitation" type (reports by Lukachova and Nachtnebel). In the reports presented at the preceding conference in Budapest it could be seen that in the "precipitation-runoff" model there is still a part ~ which has not been adequately developed; reference is to computations of runoff losses with time. Unfortunately, this situation still persists, as follows from the conference reports in Vienna. The reason for this is the _ absence of extensive investigations af the problem of infiltration of rain water, especially investigations on the basis of experimental and representative basins. The use of a mathematical approach alone does not solve this complex problem. Attention was given to this circumstance _ at the final session of the conference. The reports on the considered problem also contained important informa- tion relating to remote measurements of hydrometeorological elements and the automated collec*ion of the results of these measurements at the 175 Ff1R !1VFTf�TAT TTCF r1NTV APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102108: CIA-RDP82-00850R000240090007-3 rviN. ul i forecast preparation center. Since the maximum possible advance time for predicting rain-induced high waters in many cases is less than 24 and even 12 hours, an acceleration of the collectior of initial data for preparation of a forecast is of very great practical importance. A pre- diction of precipitation, to be sure, makes it possible to increase the advance time for predicting high water. In this connection informational reports on precipitation forecasts were presented at the conference. - Ir should be noted that in connection with the construction of a large nur~lber of reservoirs on the Danube tributaries the problem of their con- Lrol ts becoming increasingli, ucute. For this reason, and Also in rela- tlon to the prob~.ems involved in monitoring water quality (two reporta from the USSR and Bulgaria were devoted to water quality problems), short-range Forecasts are acquiring ever-greater importance. Naturally, the prediction of exceptionally high waters, constituting a danger for the population and enterprises along the shores, also retains its im- portance in the presence o'L reservoirs. At this conference there was not a single report on the prohlem of long- range forecasting of low levels of the Danube, having great importance for navigation. For that reason particular interest was shown in a re- port by V. M. Mukhin (USSR Hydrometeorological Center) on a method for predictinR the 1 evels of the Danube for 10-20 days in advance. The basis - of the method is a mathematical model which takes into acc:,unt the role of water distribution along the length of the river and the role of the dtstribution of precipitation over the area of the basin in the formation of c�iter discharges at the lowest-lying station on the river. Such fore- casts make it possible to know the depth in the main channel during the - mc,vemeiit of a ship in the Danube reach from Bezdan to Tulcea, with an ex- tent oC about 1,400 km. The report caused intereat as well due to the - fact that for solving the nroblem V. M. Mukhin used a modern approach (component analysis, theory of solution of incorrectly formulated prob- lems). - Four reports were devoted to a superlong-range forecasting of the water voltune in the Danube and an analysis of its changes over a period of many years. 'i'hey did not contain new approaches to solution of the prob- lem oF superlong-range forecasting of river runoff. We note that among the hydrologists of the Soviet Union the opinion prevails that there is _ no firm scientif ic basis for the superlong-range forecasting at this - time and that investigat'ons of this problem have a strictly explora- tory nature. The Soviet delegation presented no reports on the super- long-range forecasting of river volutae, although at the preceding con- ference at Budapest we did present reports on this problem. _ Six reports were presented on the problem of channel processes and the runoff of sediments. However, only the reports from the Soviet Union = were devoted to a quantitative pred iction of river channel deformations (g. F. crishcherko, State Hydrological Institute and V. N. Mikhaylov, 176 FOR OFFIr,IAL USE ONLY APPROVED FOR RELEASE: 2007/02/08: CIA-RDP82-00850R000200090007-3 APPROVED FOR RELEASE: 2007102/08: CIA-RDP82-00850R000200090007-3 FOR OFFICIAL USE ONLY Moscow State University). The reports on sediments, in essence, contain- ed no methods for predicting the solid runoff of the Danube and its trib- utaries. Neverthelesa, they were interesting since they gave the resulta of investigations which without question will be used in developing meth- ods for predicting the runoff of sediments. For example, in a report pre- _ sented by Bulgaria there were data characterizing the distribution of the granuiarity of bottom deposits in a reach of the Danube with an ex- - tent of about 2,000 kil.ometers and casting light on the evolution of the c eompoHitton oC depoAits alonR the lenRth of thiR river. The acouRtLe _ method for meayuring bottom sedimenta is quite novel; it was deecribed in a report by Sehlagge (West Germany). It is well known that the groblems involved in predicting the ice regime of rivers in the Danube basin are being studied primarily in the Soviet - Union. It is therefore not surprising that only our reports were present- ed at the conference: one about thP hydrological conditions of winter navigation and another about the method for long-range forecasting of the duration of ice phenomena and the times for clearing of the Danube _ from ice. There were only two special reports on the problem of evaluating the econ- - omic effectiveness of forecasts. The matter of economic effectiveness was _ touched upon in a number of other reports in which data were cited quanti- - _ tatively characterizing the advantage gained from different specific fore- casts of the water volumes in the Danube and its tributaries. Upon complet.ion c>f the scientific program of the conference an excursion w;iq orKiinlzed for its piirtic-iPunts to the Danuhe Me]k and ihba-PersenhoyK llydroelectric I'uwer S