A STUDY OF NUMERICAL MODELLING ON THE DEVELOPMENT OF DEPRESSIONS IN THE SOUTH CHINA SEA

The comparison analyses between two tropical depressions in the South China Sea are completed by numerical ex-periments with a limited area model developed in Guangzhou Institute of Tropical and Oceanic Meteorology. One of thedepressions develops and finally becomes a typhoon within 24 hours of the analyzing period (defined as DVD hereafter),and the other not (defined as UNDVD) The analyses show that their initial structures of temperature, humidity, pres-sure. wind and stability are obviously different from each other. There are a very deep wet layer and a clear warm-coredstructure in the mid-lower troposphere in the depression area in the DVD case, but with the warm core in the upperrather than lower troposphere arid a very shallow wet layer in the lower troposphere in the depression area in theUNDVD case. The developing and non-developing processes are successfully simulated by the model, showing that theCISK mechanism plays the most important role in the development of SCSDs (Tropical Depressions in the South ChinaSea). Several numerical experiments show that the structures of humidity and temperature in the depression area haveimportant effect on the development of SCSDs. It is favourable to the development of SCSDs when a very deep wet layerexists in the mid-lower troposphere or a warm-cored structure exists in the mid-lower troposphere instead of in uppertroposphere, and conversely, it is unfavourable to the development of SCSDs when the wet layer is very shallow in thelower troposphere or the warm-cored structure is in the upper troposphere rather than in mid-lower troposphere. Thestructures of stability for each case are also analysed, which shows evident differences between the two cases, with adeeper instability layer in the DVD case and a shallower one in the UNDVD case. Finally, the sensitivity of the devel-opment of the SCSDs to the vertical structures of humidity and temperature in the depression area is discussed.     

AN ANALYSIS OF THE STRUCTURE OF THUNDERSTORM IN THE ATMOSPHERIC BOUNDARY LAYER

The structure of thunderstorm in the atmospheric boundary layer in Beijing area is analysed by using three-year data of tower. It is indicated that the outflow current of the thunderstorm in the lower layer is a sort of density current. An area of evident wind direction shear is found at about half an hour to one hour before the arrival of the gust front.The maximum intensity of the shear can reach 0.35sec-1.The inner structure within the density current is also very complicated. At the nocturnal stable boundary layer in summertime, the development of the convective motions is often triggered due to the instability of the Kelvin-Helmholtz wave.     

NUMERICAL SIMULATION OF THE INFLUENCE OF VEGETATION COVER FACTOR ON BOUNDARY LAYER CLIMATE IN SEMI-ARID REGION

A numerical model has been developed for simulating land-surface processes and atmosphericboundary layer climate of vegetation and desert in semi-arid region.Dynamically,thermal andhydrological processes take place in the atmospheric boundary layer.Vegetation and surface layerof soil are included in the soil-vegetation-atmosphere coupled system,in which,vegetation isconsidered as a horizontally uniform layer,soil is divided into 13 layers and the horizontaldifferences of variables in the system are neglected.The influence of local boundary layer climateby vegetation cover factor is simulated with the coupled model in the semi-arid region of NorthwestChina (around 38°N,105°E).Results indicate that due to significant differences of water andenergy budgets in vegetation and desert region,the air is colder and wetter over the vegetation andcorrespondingly an obvious local circulation in the lower atmosphere is formed.Simulating results also show that maximum updraft and downdraft occur around thevegetation-desert marginal area,where the dynamical and thermodynamical properties of PBL(Planetary Boundary Layer) are uncontinuous.It is stronger at daytime,weaker and reverse atnighttime.In the simulation,the moisture inversion phenomena are analyzed.Finally.theinfluences of vegetation cover factor exchange on local boundary layer climate are simulated.Thesimulating results bring to light that water may be conserved and improved by developing treeplanting and afforestation,and improving cover factor of vegetation in local ecoenvironment,andthis is an important way of transforming local climate in arid and semi-arid area.Results indicatethat the coupled model can be used to study the soil-vegetation-atmosphere interaction and localboundary layer climate.     

NUMERICAL SIMULATION FOR THE EFFECTS OF PBL AND THE SURFACE ON POLLUTANT CONCENTRATIONS

On the basis of theoretical and experimental results of study of planetary boundary layer (PBL), the physical parameters describing the structure of PBL are calculated by using the data obtained from a meteorological tower and the effects of PBL and the surface on pollutant concentrations are numerically simulated with a time-dependent two-dimensional advection and diffusion equation.It is shown that the diurnal variation of PBL results in that of concentration. The height of mixing layer is an important factor to determine the ground-level concentration. As for an elevated point source, the height of mixing layer, growing from lower to higher than the releasing height is a necessary condition for the phenomenon of fumigation. It is also shown that the surface may be considered as a boundary with perfect reflection when Vd ≤0.001 m s-1, but has an important effect on concentration and must be carefully dealt with when Vd≥0.01 m s-1.     

THE STUDY OF THE CHARACTERISTICS OF BOUNDARY LAYER OVER HILL TERRAIN

Taking advantage of the actual continuous materials of hilly area, we analyse the structures of the wind and temperature fields of planetary bounday layer over hill terrain under different weather conditions. As the results show that being similar to the situation in flat area, the planetary boundary layer has striking characteristics of daily variation and distinct regularity. We also find that in certain weather situation or under certain weather condition, these characteristics are principally determined by complex effect of such factors as the features of specific radiation field and intensity of turbulent exchange, and the advection of wind, and the thermodynamic effect of valley or lake-land, and the extending direction of valley, as well as difference between states of sunlessness and sunshine of the surface. Therein, the local environmental situation, where the reference station stands, is of considerable importance to the formation of the characteristics of planetary boundary over hill terrain (shortened as HTBL). These results not only help us to disclose the characteristics of planetary boundary layer over hill terrain in theory, but also are of great value in practice, like developing hilly area, and implementing environment conservation etc.     

ADJUSTMENT OF WIND AND PRESSURE FIELDS IN THE BOUNDARY LAYER

The vertical transport of mass at the top of the boundary layer is considered as a link between theboundary layer and free atmosphere.The adjustment of the wind and pressure fields in the boundary layeris studied under the consideration of the interaction between the boundary layer and free atmosphere.Thevertical motion at the top of the boundary layer is evaluated.The results show that the distinguished differ-ences of the present results from classical Ekman layer do exist and they are discussed in the paper.     

THE INFLUENCES OF OROGRAPHY UPON THE FLOW WITHIN EKMAN BOUNDARY LAYER UNDER THE APPROXIMATION OF GEOSTROPHIC MOMENTUM

In this paper, the influences of orography on the boundary layer flow with the approximation of geo strophic momentum are studied. The wind velocity at the lower boundary will not always be zero when the orography exists. So the structure of the boundary layer flow, as well as the vertical velocity at the top of the boundary layer, is affected. There are three factors affecting the vertical motion at the top of the boundary layer: lifting due to orography; divergence due to Ekman flow, and advection of the geostrophic momentum. These effects and the features of the flow within the boundary layer are discussed in detail.     

THE WIND IN THE BAROCLINIC BOUNDARY LAYER WITH THREE SUBLAYERS INCORPORATING THE WEAK NONLINEAR EFFECT

In considering the weak non-linear effect, and using the small parameter expansion method, the analytical expressions of the wind distribution within PBL (planetary boundary layer) and the vertical velocity at the top of the PBL are obtained when the PBL is divided into three layers and different eddy transfer coefficients K are adopted for the three layers. The conditions of barotropy and neutrality for the PBL are extended to that of baroclinity and non-neutral stratification. An example of a steady circular vortex is used to display the characteristics of the horizontal wind within the PBL and the vertical velocity at the top of the PBL. Some new results have been obtained, indicating that the magnitude of the speed in the lower height calculated by the present model is larger than that by the model in which k is a constant within the whole boundary layer, for example, in the classical Ekman boundary layer model and the model by Wu (1984). The angle between the wind at the top of the PBL and the wind     

A NUMERICAL EXPERIMENT OF THE PBL WITH GEOSTROPHIC MOMENTUM APPROXIMATION

In this paper, a numerical experiment of the motion in the PBL (planetary boundary layer) is performed with geostrophic momentum approximation, in which a nonlinear eddy transfer coefficient is used. Some results are obtained for the boundary layer winds in cyclone-anticyclone and trough-ridge systems. This treatment improves W-B's work. The effects of geostrophic wind tendency and the advection of the geostrophic wind on the winds in the PBL are also discussed.     

OBSERVATIONS OF TEMPERATURE MICROSTRUCTURE IN THE ATMOSPHERE

In order to provide data of atmospheric temperature microstructure for the investigation of light propagation we measured fluctuations of atmospheric temperature below the height of 300 m with a platinum wire thermometer in Tianjin in May and September, 1980. The results measured in daytime revealed some properties of the temperature structure parameter and spectrum. It has been confirmed that there is a maximum in the profile of the structure parameter produced probably by the entrainment in the interfacial layer at the top of convective boundary layer. The average of CT2 in the interfacial layer and its Wyngaard calculating method are discussed, and the thickness of the interfacial layer is obtained.It is shown by spectrum analysis that a wide inertial subrange exists in the convective boundary layer and the strong turbulent zone in the free atmosphere. The spectral law with the power of -2.5 was measured within the upper half of boundary layer over the sea in vicinity of Tanggu.     

NUMERICAL SIMULATIONS OF THE EFFECTS OF AEROSOLS IN THE ATMOSPHERIC BOUNDARY LAYER ON THE CLIMATE

The climatic effects of the atmospheric boundary aerosols are studied by the use of a three-dimensional climatemodel.Simulated results show that the climate states both at the surface and in the atmosphere change remarkably whenthe aerosols with different optical thicknesses and properties are introduced into the atmospheric boundary layer of themodel.The aerosols absorb and scatter the solar shortwave radiation,therefore,they reduce the solar energy reachingthe ground surface and decrease the surface and the soil temperatures.The temperature in the boundary layer increasesbecause of the supplementary absorption of radiation by the boundary aerosols.In the atmosphere,the temperatures atall isobaric surfaces rise up except for the 100 hPa level.The atmospheric temperatures below the 500 hPa level aredirectly influenced by the boundary aerosols,while the atmospheric temperatures above the 500 hPa level are influencedby the heating due to convective condensation and the changes in the vertical motion field.Cyclonic differential circula-tions appear over the desert areas at the low levels,and anticyclonic differential circulations exist at the upper levels inthe horizontal flow fields.The vertical motions change in correspondence with the differential circulations.The changesin precipitation are directly related to that of vertical motions.The mechanisms of climate effects of the boundaryaerosols are also discussed in this paper.     

A MODEL STUDY OF THE NOCTURNAL BOUNDARY LAYER

In this paper, a second-order model is proposed for the study of the evolution of the nocturnal boundary layer (NBL). The model is tested against the Wangara data on atmospheric boundary layer. The computer results show that the model can simulate some important characters observed in the NBL, and that a kind of sudden change may occur in the developing process of NBL.     

The influence of convergence movement on turbulent transportation in the atmospheric boundary layer

Classical turbulent K closure theory of the atmospheric boundary layer assumes that the vertical turbulent transport flux of any macroscopic quantity is equivalent to that quantity‘s vertical gradient transport flux. But a cross coupling between the thermodynamic processes and the dynamic processes in the atmospheric system is demonstrated based on the Curier-Prigogine principle of cross coupling of linear thermodynamics. The vertical turbulent transportation of energy and substance in the atmospheric boundary layer is related not only to their macroscopic gradient but also to the convergence and the di-vergence movement. The transportation of the convergence or divergence movement is important for the atmospheric boundary layer of the heterogeneous underlying surface and the convection boundary layer.Based on this, the turbulent transportation in the atmospheric boundary layer, the energy budget of the heterogeneous underlying surface and the convection boundary layer, and the boundary layer parameteri-zation of land surface processes over the heterogeneous underlying surface are studied. This research offers clues not only for establishing the atmospheric boundary layer theory about the heterogeneous underlying surface, but also for overcoming the difficulties encountered recently in the application of the atmospheric boundary layer theory.     

THE NUMERICAL EXPERIMENTS ON APPLICATING GEOSTROPHIC MOMENTUM APPROXIMATION TO THE BAROCLINIC AND NON-NEUTRAL PBL

In this paper, Wu and Blumen's boundary layer geostrophic momentum approximation model (Wu and Blu-men, 1982) is applied to baroclinic and non-neutral PBL, the motion equations for the PBL under the geostrophic momentum approximation are solved, in which the eddy transfer coefficient is a function of the distributions of the wind and temperature. The results are compared with those in barotropic and neutral conditions with the geostrophic momentum approximation. It is found that in the baroclinic condition, the wind distribution has both the characteristics of a steady, homogeneous and baroclinic PBL and those caused by the geostrophic momentum approximation. Those in non-neutral conditions show that they retain the intrinsic characteristics for the wind in non-neutral PBL, at the same time, the effects of the large-scale advection and local variation are also included. We can predict the wind in the non-neutral and baroclinic PBL by use of the geostrophic momentum approximation when the temporal and sp     

THERMODYNAMIC PROPERTIES AND ROLE OF THERMALS IN THE STRATUS-TOPPED BOUNDARY LAYER

From results of two large-eddy simulation of stratus-topped boundary layer,the structure,thermodynamic prop-erties and role of thermals are investigated by using conditional sampling methods,which divided the thermals into thewarm/moist,cool/dry,warm/dry and cool/moist events.The results show that the main turbulent circulation in thestratus-topped boundary layer is composed of the warm/moist updraft and cool/dry downdraft.Below entrainmentregion,the warm/moist updrafts and cool/dry downdrafts are,respectively,positively and negatively buoyant andcontribute most to total fluxes and variances.Evaporative cooling has important effect on the structure and thermodynamic properties of thermals instratus-topped boundary layer.     

A Study of the Atmospheric Boundary Layer Structure During a Clear Day in the Arid Region of Northwest China

The local climate and atmospheric circulation pattern exert a clear influence on the atmospheric boundary layer （ABL） formation and development in Northwest China. In this paper, we use field observational data to analyze the distribution and characteristics of the ABL in the extremely arid desert in Dunhuang, Northwest China. These data show that the daytime convective boundary layer and night time stable boundary layer in this area extend to higher altitudes than in other areas. In the night time, the stable boundary layer exceeds 900 m in altitude and can sometimes peak at 1750 m, above which the residual layer may reach up to about 4000 m. The daytime convective boundary layer develops rapidly after entering the residual layer, and exceeds 4000 m in thickness. The results show that the deep convective boundary layer in the daytime is a pre-requisite for maintaining the deep residual mixed layer in the night time. Meanwhile, the deep residual mixed layer in the night time provides favorable thermal conditions for the development of the convective boundary layer in the daytime. The prolonged periods of clear weather that often occurs in this area allow the cumulative effect of the atmospheric residual layer to develop fully, which creates thermal conditions beneficial for the growth of the daytime convective boundary layer. At the same time, the land surface process and atmospheric motion within the surface layer in this area also provide helpful support for forming the particular structure of the thermal ABL. High surface temperature is clearly the powerful external thermal forcing for the deep convective boundary layer. Strong sensible heat flux in the surface layer provides the required energy. Highly convective atmosphere and strong turbulence provide the necessary dynamic conditions, and the accumulative effect of the residual layer provides a favorable thermal environment.     

DYNAMICAL CHARACTERISTICS OF THE STREAM FIELD OF COLD FRONT IN BOUNDARY LAYER

The motion equation of atmospheric boundary layer with cold front surface under geostrophicmomentum approximation is solved and some characteristics of the stream field for the cold frontsurface in the boundary layer are derived,for example,the slope of the cold front surface increaseswith the increases of geostrophic vorticity and the temporal tendency of geostrophic wind speed,and also increases with the decrease of the component of thermal wind speed along the motiondirection of the front:the stream field above and below the cold front surface in the boundary layerdepends on the slope of the front surface,geostrophic wind speed and its temporal and spatialdistributions.A common characteristic is that there exist updraft motion above the cold frontsurface and downdraft motion below it.     

EFFECT OF WIND VERTICAL SHEAR ON DIFFUSION CHARACTERISTICS IN THE MESOSCALE RANGE

Using the turbulent statistical form of the vertical vortex diffusion coefficient K2, in the planetary boundary layer (PBL) and Ekman spiral wind profile, the three-dimensional diffusion equation is solved by the numerical method. The influences of vertical shear of both wind direction and wind speed on pollution trajectory and horizontal diffusion parameters σy are numerically analysed. The expressions of both pollution trajectory and σy, including the factor of wind shear, are obtained. The results show that the vertical shear of wind is important among all factors affecting the mesoscale dispersion. Specifically, from neutral to stable atmospheric conditions, vertical shear of wind makes greater contribution to σy than turbulence, thus it is the most important factor. In this paper, we have compared horizontal dispersion pattern with both Pasquill's dispersion pattern considering wind direction shear, and experimental data collected at 9 different sites ranging from 10 to 100 km, and the results show that our dispersion pattern is closer to the experimental values than Pasquill's results, and his correction to shear of wind direction is too large under the stable conditions.     

Evaluation of the Ocean Feedback on Height Characteristics of the Tropical Cyclone Boundary Layer

In this study, the interaction between the tropical cyclone（TC） and the underlying ocean is reproduced by using a coupled atmosphere-ocean model. Based on the simulation results, characteristics of the TC boundary layer depth are investigated in terms of three commonly used definitions, i.e., the height of the mixed layer depth（HVTH）, the height of the maximum tangential winds（HTAN）, and the inflow layer depth（HRAD）. The symmetric height of the boundary layer is shown to be cut down by the ocean response, with the decrease of HVTH slightly smaller than that of HTAN and HRAD. The ocean feedback also leads to evident changes in asymmetric features of the boundary layer depth. The HVTH in the right rear of the TC is significantly diminished due to presence of the cold wake, while the changes of HVTH in other regions are rather small. The decreased surface virtual potential temperature by the cold wake is identified to be dominant in the asymmetric changes in HVTH. The impacts of ocean response on the asymmetric distributions of HTAN are nonetheless not distinct, which is attributed to the highly axisymmetric property of tangential winds. The HRAD possesses remarkable asymmetric features and the inflow layer does not exist in all regions, an indication of the inadequacy of the definition based on symmetric inflow layer depth. Under influences of the cold wake, the peak inflow area rotates counterclockwise distinctly. As a consequence, the HRAD becomes deeper in the east while shallower in the west of the TC.     

THE INVERSE PROBLEM OF KLEIN-GORDON EQUATION BOUNDARY VALUE PROBLEM AND ITS APPLICATION IN DATA ASSIMILATION

Since the solution of elliptic partial differential equations continuously depends on the boundary condition, the Euler equation derived from variational method cannot be solved without boundary condition. It is often difficult to provide the exact boundary condition in the practical use of variational method. However, in some application problems such as the remote sensing data assimilation, the values can be easily obtained in the inner region of the domain. In this paper, the boundary condition is tried to be retrieved by using part solutions in the inner area. Firstly, the variational problem of remote sensing data assimilation within a circular area is established. The Klein-Gordon elliptic equation is derived from the Euler method of variational problems with assumed boundary condition. Secondly, a computer-friendly Green function is constructed for the Dirichlet problem of two-dimensional Klein-Gordon equation, with the formal solution according to Green formula. Thirdly, boundary values are retrieved by solving the optimal problem which is constructed according to the best approximation between formal solutions and high-accuracy measurements in the interior of the domain. Finally, the assimilation problem is solved on substituting the retrieved boundary values into the Klein-Gordon equation. It is a type of inverse problem in mathematics. The advantage of this method lies in that it overcomes the inherent instability of the inverse problem of Fredholm integral equation and alleviates the error introduced by artificial boundary condition in data fusion using variational method in the past.