Evaluation of the impact of observations on blended sea surface winds in a two-dimensional variational scheme using degrees of freedom

This paper presents an evaluation of the observational impacts on blended sea surface winds from a two-dimensional variational data assimilation (2D-Var) scheme. We begin by briefly introducing the analysis sensitivity with respect to observations in variational data assimilation systems and its relationship with the degrees of freedom for signal (DFS), and then the DFS concept is applied to the 2D-Var sea surface wind blending scheme. Two methods, a priori and a posteriori, are used to estimate the DFS of the zonal (u) and meridional (v) components of winds in the 2D-Var blending scheme. The a posteriori method can obtain almost the same results as the a priori method. Because only by-products of the blending scheme are used for the a posteriori method, the computation time is reduced significantly. The magnitude of the DFS is critically related to the observational and background error statistics. Changing the observational and background error variances can affect the DFS value. Because the observation error variances are assumed to be uniform, the observational influence at each observational location is related to the background error variance, and the observations located at the place where there are larger background error variances have larger influences. The average observational influence of u and v with respect to the analysis is about 40%, implying that the background influence with respect to the analysis is about 60%.     

On the problem of velocity of long continuous waves in overland flows

The velocity of one-dimensional long continuous wave is computed from the condition of continuity of water both with its variable and constant discharge along the way passing through the control volume of water moving with the velocity V _W . The wave stability on the free surface of an overland flow is analyzed.     

Statistics of the <Emphasis Type="Italic">Z</Emphasis>–<Emphasis Type="Italic">R</Emphasis> Relationship for Strong Convective Weather over the Yangtze–Huaihe River Basin and Its Application to Radar Reflectivity Data Assimilation for a Heavy Rain Event

The relationship between the radar reflectivity factor (Z) and the rainfall rate (R) is recalculated based on radar observations from 10 Doppler radars and hourly rainfall measurements at 6529 automatic weather stations over the Yangtze–Huaihe River basin. The data were collected by the National 973 Project from June to July 2013 for severe convective weather events. The Z–R relationship is combined with an empirical q_r–R relationship to obtain a new Z–q_r relationship, which is then used to correct the observational operator for radar reflectivity in the three-dimensional variational (3DVar) data assimilation system of the Weather Research and Forecasting (WRF) model to improve the analysis and prediction of severe convective weather over the Yangtze–Huaihe River basin. The performance of the corrected reflectivity operator used in the WRF 3DVar data assimilation system is tested with a heavy rain event that occurred over Jiangsu and Anhui provinces and the surrounding regions on 23 June 2013. It is noted that the observations for this event are not included in the calculation of the Z–R relationship. Three experiments are conducted with the WRF model and its 3DVar system, including a control run without the assimilation of reflectivity data and two assimilation experiments with the original and corrected reflectivity operators. The experimental results show that the assimilation of radar reflectivity data has a positive impact on the rainfall forecast within a few hours with either the original or corrected reflectivity operators, but the corrected reflectivity operator achieves a better performance on the rainfall forecast than the original operator. The corrected reflectivity operator extends the effective time of radar data assimilation for the prediction of strong reflectivity. The physical variables analyzed with the corrected reflectivity operator present more reasonable mesoscale structures than those obtained with the original reflectivity operator. This suggests that the new statistical Z–R relationship is more suitable for predicting severe convective weather over the Yangtze–Huaihe River basin than the Z–R relationships currently in use.     

Formulation of stability-dependent empirical relations for turbulent intensities from surface layer turbulence measurements for dispersion parameterization in a lagrangian particle dispersion model

Season- and stability-dependent turbulence intensity (σ _u /u _*, σ _v /u _*, σ _w /u _*) relationships are derived from experimental turbulence measurements following surface layer scaling and local stability at the tropical coastal site Kalpakkam, India for atmospheric dispersion parameterization. Turbulence wind components (u′, v′, w′) measured with fast response UltraSonic Anemometers during an intense observation campaign for wind field modeling called Round Robin Exercise are used to formulate the flux–profile relationships using surface layer similarity theory and Fast Fourier Transform technique. The new relationships (modified Hanna scheme) are incorporated in a Lagrangian Particle Dispersion model FLEXPART-WRF and tested by conducting simulations for a field tracer dispersion experiment at Kalpakkam. Plume dispersion analysis of a ground level hypothetical release indicated that the new turbulent intensity formulations provide slightly higher diffusivity across the plume relative to the original Hanna scheme. The new formulations for σ _u , σ _v , σ _w are found to give better agreement with observed turbulent intensities during both stable and unstable conditions under various seasonal meteorological conditions. The simulated concentrations using the two methods are compared with those obtained from a classical Gaussian model and the observed SF₆ concentration. It has been found that the new relationships provide comparatively higher diffusion across the plume relative to the model default Hanna scheme and provide downwind concentration results in better agreement with observations.     

Characteristics of Turbulent Transfer during Episodes of Heavy Haze Pollution in Beijing in Winter 2016/17

We analyzed the structure and evolution of turbulent transfer and the wind profile in the atmospheric boundary layer in relation to aerosol concentrations during an episode of heavy haze pollution from 6 December 2016 to 9 January 2017. The turbulence data were recorded at Peking University’s atmospheric science and environment observation station. The results showed a negative correlation between the wind speed and the PM_2.5 concentration. The turbulence kinetic energy was large and showed obvious diurnal variations during unpolluted (clean) weather, but was small during episodes of heavy haze pollution. Under both clean and heavy haze conditions, the relation between the non-dimensional wind components and the stability parameter z/L followed a 1/3 power law, but the normalized standard deviations of the wind speed were smaller during heavy pollution events than during clean periods under near-neutral conditions. Under unstable conditions, the normalized standard deviation of the potential temperature σ _θ /|θ_*| was related to z/L, roughly following a –1/3 power law, and the ratio during pollution days was greater than that during clean days. The three-dimensional turbulence energy spectra satisfied a –2/3 power exponent rate in the high-frequency band. In the low-frequency band, the wind velocity spectrum curve was related to the stability parameters under clear conditions, but was not related to atmospheric stratification under polluted conditions. In the dissipation stage of the heavy pollution episode, the horizontal wind speed first started to increase at high altitudes and then gradually decreased at lower altitudes. The strong upward motion during this stage was an important dynamic factor in the dissipation of the heavy haze.     

Raindrop Size Distribution Parameters Retrieved from Guangzhou S-band Polarimetric Radar Observations

According to the statistical shape–slope (μ–Λ) relationship observed for the first time by several 2D-Video-Distrometers (2DVD) in southern China, a constrained gamma (C-G) model was proposed for the retrieval of rain drop size distributions (DSDs) from Guangzhou S-band polarimetric radar observations. Two typical precipitation processes were selected to verify the accuracy of the retrieval scheme. The μ–Λ relationship: Λ = 0.0241μ² + 0.867μ + 2.453 was obtained based on the 2DVD observation results from at Huizhou Longmen station, which is a very representative location in the area. Relying on the Guangzhou polarimetric radar measurements of radar reflectivity (Z_HH) and differential reflectivity (Z_DR), the gamma (Γ) size distribution parameters (N₀, μ, and Λ) can be retrieved by the C-G model retrieval scheme. The results show that the Guangzhou polarimetric radar retrievals of DSDs were close to the 2DVD observations at Guangzhou Maofengshan station. The rain rate, mass mean diameter, and normalized intercept parameter of radar retrievals were in good agreement with the 2DVD observations, and the relative errors were less than 10%. The overall accuracy of the retrieval scheme was high. The retrieval scheme has established the relationship between the polarimetric radar measurements and gamma size distribution parameters. It will be helpful to in-depth research and application of the dual-polarization radar data in microphysical precipitation processes analysis, as well as convection-resolved numerical model data assimilation and prediction effect evaluation.     

Wind-Speed Undulations Over Swell: Field Experiment and Interpretation

Results of field measurements of the swell-induced undulation of the wind speed taken from a Black Sea platform are presented. The wind speed and its fluctuations were measured at several heights between 1.3 and 21 m above the mean sea level under various wind and swell conditions. Parameters of the swell-induced undulations were derived from cross spectra of the wind-speed fluctuations and the sea-surface displacement. As found, the phase and the amplitude of the wind speed undulation in the layer from k _p z = 0.1 to k _p z = 3 (k _p is the swell wavenumber) are in good agreement with the theory of inviscid shear flow over a wavy surface. The main feature of the vertical profile of the swell-induced undulation is the exponential attenuation of its amplitude with height typical for the potential flow over the fast running waves. At the lowest levels the potential undulations are significantly distorted by the wind-speed variations caused by the vertical displacements of the shear airflow relative to a fixed sensor. No direct impact of swell on the mean properties of the turbulent boundary layer at k _p z > 0.1 is revealed. In particular, the mean wind-speed profile and spectra of the horizontal velocity in the inertial subrange obey Monin-Obukhov similarity theory.     

On the extension of the wind profile over homogeneous terrain beyond the surface boundary layer

Analysis of profiles of meteorological measurements from a 160 m high mast at the National Test Site for wind turbines at Høvsøre (Denmark) and at a 250 m high TV tower at Hamburg (Germany) shows that the wind profile based on surface-layer theory and Monin-Obukhov scaling is valid up to a height of 50–80 m. At higher levels deviations from the measurements progressively occur. For applied use an extension to the wind profile in the surface layer is formulated for the entire boundary layer, with emphasis on the lowest 200–300 m and considering only wind speeds above 3 m s^?1 at 10 m height. The friction velocity is taken to decrease linearly through the boundary layer. The wind profile length scale is composed of three component length scales. In the surface layer the first length scale is taken to increase linearly with height with a stability correction following Monin-Obukhov similarity. Above the surface layer the second length scale (L _MBL ) becomes independent of height but not of stability, and at the top of the boundary layer the third length scale is assumed to be negligible. A simple model for the combined length scale that controls the wind profile and its stability dependence is formulated by inverse summation. Based on these assumptions the wind profile for the entire boundary layer is derived. A parameterization of L _MBL is formulated using the geostrophic drag law, which relates friction velocity and geostrophic wind. The empirical parameterization of the resistance law functions A and B in the geostrophic drag law is uncertain, making it impractical. Therefore an expression for the length scale, L _MBL , for applied use is suggested, based on measurements from the two sites.     

Cold air outbreaks along a non-frozen sea channel: effects of wind on snow bands

Wintertime cold air outbreaks along a non-frozen sea channel or a long lake can become destructive if the related bands of heavy snowfall hit onto land. The forcing for such bands is studied with a 2D numerical model set across an east–west sea channel at 60^oN (‘Gulf of Finland’), varying the basic geostrophic wind V _g. Without any V _g opposite coastal land breezes emerge with convergence. This results in a quasi-steady rising motion w _max ~ 7.5 cm/s at 600 m in the middle of the gulf, which can force a snow band. During weak V _g, the rising motion is reduced but least so for winds from 60^o to 80^o (~ENE), when modest alongshore bands could exist near the downstream (Estonian) coast. During V _g of 4–6 m/s from any direction, the land breezes and rising motions are reduced more effectively, so snow bands are not expected during moderate basic flow. In contrast, during a strong V _g of 20–25 m/s from 110^o to 120^o (~ESE) the land breeze perturbations are intense with w _max up to 15–18 cm/s. The induced alongshore bands of heavy snowfall are located in these cases at the sea but quite close to the downstream (Finnish) coast. They can suddenly make a landfall if the basic wind turns clockwise.     

Comprehensive modeling of monthly mean soil temperature using multivariate adaptive regression splines and support vector machine

Soil temperature (T _s) and its thermal regime are the most important factors in plant growth, biological activities, and water movement in soil. Due to scarcity of the T _s data, estimation of soil temperature is an important issue in different fields of sciences. The main objective of the present study is to investigate the accuracy of multivariate adaptive regression splines (MARS) and support vector machine (SVM) methods for estimating the T _s. For this aim, the monthly mean data of the T _s (at depths of 5, 10, 50, and 100 cm) and meteorological parameters of 30 synoptic stations in Iran were utilized. To develop the MARS and SVM models, various combinations of minimum, maximum, and mean air temperatures (T _min, T _max, T); actual and maximum possible sunshine duration; sunshine duration ratio (n, N, n/N); actual, net, and extraterrestrial solar radiation data (R _s, R _n, R _a); precipitation (P); relative humidity (RH); wind speed at 2 m height (u ₂); and water vapor pressure (V_p) were used as input variables. Three error statistics including root-mean-square-error (RMSE), mean absolute error (MAE), and determination coefficient (R ²) were used to check the performance of MARS and SVM models. The results indicated that the MARS was superior to the SVM at different depths. In the test and validation phases, the most accurate estimations for the MARS were obtained at the depth of 10 cm for T _max, T _min, T inputs (RMSE = 0.71 °C, MAE = 0.54 °C, and R ² = 0.995) and for RH, V _p, P, and u ₂ inputs (RMSE = 0.80 °C, MAE = 0.61 °C, and R ² = 0.996), respectively.     

Multivariate drought frequency estimation using copula method in Southwest China

This study reveals the impacts of climatic variable trends on drought severity in Xinjiang, China. Four drought indices, including the self-calibrating Palmer drought severity index (sc-PDSI), Erinç’s index (I _m), Sahin’s index (I _sh), and UNEP aridity index (AI), were used to compare drought severity. The ensemble empirical mode decomposition and the modified Mann-Kendall trend test were applied to analyze the nonlinear components and trends of the climatic variable and drought indices. Four and six climatic scenarios were generated in sc-PDSI, I _m, I _sh, and AI with different combinations of the observed and detrended climatic variables, respectively. In Xinjiang, generally increasing trends in minimal, average, and maximal air temperature (T _min, T _ave, T _max) and precipitation (P) were found, whereas a decreasing trend in wind speed at 2 m height (U ₂) was observed. There were significantly increasing trends in all of the four studied drought indices. Drought relief was more obvious in northern Xinjiang than in southern Xinjiang. The strong influences of increased P on drought relief and the weak influences of increased T _min, T _ave, and T _max on drought aggravation were shown by comparing four drought indices under different climate scenarios. Decreased U ₂ had a weak influence on drought, as shown by the AI in different climate scenarios. The weak influences of T and U ₂ were considered to be masked by the strong influences of P on droughts. Droughts were expected to be more severe if P did not increase, but were likely milder without an increase in air temperature and with a decrease in U ₂.     

Classification of Vertical Wind Speed Profiles Observed Above a Sloping Forest at Nighttime Using the Bulk Richardson Number

Wind speed profiles above a forest canopy relate to scalar exchange between the forest canopy and the atmosphere. Many studies have reported that vertical wind speed profiles above a relatively flat forest can be classified by a stability index developed assuming wind flow above a flat plane. However, can such a stability index be used to classify vertical wind speed profiles observed above a sloping forest at nighttime, where drainage flow often occurs? This paper examines the use of the bulk Richardson number to classify wind speed profiles observed above a sloping forest at nighttime. Wind speed profiles above a sloping forest were classified by the bulk Richardson number Ri _B . Use of Ri _B distinguished between drainage flow, shear flow, and transitional flow from drainage flow to shear flow. These results suggest that Ri _B is useful to interpret nighttime CO₂ and energy fluxes above a sloping forest. Through clear observational evidence, we also show that the reference height should be high enough to avoid drainage-flow effects when calculating Ri _B .     

Simplification and Validation of a Spectral-Tensor Model for Turbulence Including Atmospheric Stability

A spectral-tensor model of non-neutral, atmospheric-boundary-layer turbulence is evaluated using Eulerian statistics from single-point measurements of the wind speed and temperature at heights up to 100 m, assuming constant vertical gradients of mean wind speed and temperature. The model has been previously described in terms of the dissipation rate \(\epsilon \), the length scale of energy-containing eddies \(\mathcal {L}\), a turbulence anisotropy parameter \(\varGamma \), the Richardson number Ri, and the normalized rate of destruction of temperature variance \(\eta _\theta \equiv \epsilon _\theta /\epsilon \). Here, the latter two parameters are collapsed into a single atmospheric stability parameter z / L using Monin–Obukhov similarity theory, where z is the height above the Earth’s surface, and L is the Obukhov length corresponding to \(\{Ri,\eta _\theta \}\). Model outputs of the one-dimensional velocity spectra, as well as cospectra of the streamwise and/or vertical velocity components, and/or temperature, and cross-spectra for the spatial separation of all three velocity components and temperature, are compared with measurements. As a function of the four model parameters, spectra and cospectra are reproduced quite well, but horizontal temperature fluxes are slightly underestimated in stable conditions. In moderately unstable stratification, our model reproduces spectra only up to a scale \(\sim \) 1 km. The model also overestimates coherences for vertical separations, but is less severe in unstable than in stable cases.     

Improving Wind-Ramp Forecasts in the Stable Boundary Layer

The viability of wind-energy generation is dependent on highly accurate numerical wind forecasts, which are impeded by inaccuracies in model representation of boundary-layer processes. This study revisits the basic theory of the Mellor, Yamada, Nakanishi, and Niino (MYNN) planetary boundary-layer parametrization scheme, focusing on the onset of wind-ramp events related to nocturnal low-level jets. Modifications to the MYNN scheme include: (1) calculation of new closure parameters that determine the relative effects of turbulent energy production, dissipation, and redistribution; (2) enhanced mixing in the stable boundary layer when the mean wind speed exceeds a specified threshold; (3) explicit accounting of turbulent potential energy in the energy budget. A mesoscale model is used to generate short-term (24 h) wind forecasts for a set of 15 cases from both the U.S.A. and Germany. Results show that the new set of closure parameters provides a marked forecast improvement only when used in conjunction with the new mixing length formulation and only for cases that are originally under- or over-forecast (10 of the 15 cases). For these cases, the mean absolute error (MAE) of wind forecasts at turbine-hub height is reduced on average by 17%. A reduction in MAE values on average by 26% is realized for these same cases when accounting for the turbulent potential energy together with the new mixing length. This last method results in an average reduction by at least 13% in MAE values across all 15 cases.     

Field Experiments and Numerical Modeling of Wind Speed and Surface Waves in Medium-size Inland Reservoirs

An attempt is made to apply the modern methods of surface wave simulation developed for oceanic conditions to the modeling of waves in medium-size inland reservoirs (10–100 km). The results of field measurements of wind speed and waves are described, and on their basis the parameterization C _D (U ₁₀) is proposed. WAVEWATCH III spectral wave model was adapted to the conditions of a medium-size inl and reservoir. The simulated data are compared with the field data. The use of the new parameterization C _D (U ₁₀) allowed reducing the values of the wind wave growth rate that improved consistency in data from the field experiment and numerical modeling concerning the height of significant waves. Further steps towards improving the quality of prediction of the adapted WAVEWATCH III model are discussed.     

The Impact of Landscape Fragmentation on Atmospheric Flow: A Wind-Tunnel Study

Landscape discontinuities such as forest edges play an important role in determining the characteristics of the atmospheric flow by generating increased turbulence and triggering the formation of coherent tree-scale structures. In a fragmented landscape, consisting of surfaces of different heights and roughness, the multiplicity of edges may lead to complex patterns of flow and turbulence that are potentially difficult to predict. Here, we investigate the effects of different levels of forest fragmentation on the airflow. Five gap spacings (of length approximately 5h, 10h, 15h, 20h, 30h, where h is the canopy height) between forest blocks of length 8.7h, as well as a reference case consisting of a continuous forest after a single edge, were investigated in a wind tunnel. The results reveal a consistent pattern downstream from the first edge of each simulated case, with the streamwise velocity component at tree top increasing and turbulent kinetic energy decreasing as gap size increases, but with overshoots in shear stress and turbulent kinetic energy observed at the forest edges. As the gap spacing increases, the flow appears to change monotonically from a flow over a single edge to a flow over isolated forest blocks. The apparent roughness of the different fragmented configurations also decreases with increasing gap size. No overall enhancement of turbulence is observed at any particular level of fragmentation.     

Turbulent Structure of a Simplified Urban Fluid Flow Studied Through Stereoscopic Particle Image Velocimetry

Stereoscopic particle image velocimetry was used to provide a three-dimensional characterization of the flow around a simplified urban model defined by a 5 by 7 array of blocks, forming four parallel streets, perpendicular to the incoming wind direction corresponding to a zero angle of incidence. Channeling of the flow through the array under consideration was observed, and its effect increased as the incoming wind direction, or angle of incidence (AOI), was changed from \(0^{\circ }\) to \(15^{\circ }\), \(30^{\circ }\), and \(45^{\circ }\). The flow between blocks can be divided into two regions: a region of low turbulence kinetic energy (TKE) levels close to the leeward side of the upstream block, and a high TKE area close to the downstream block. The centre of the arch vortex is located in the low TKE area, and two regions of large streamwise velocity fluctuation bound the vortex in the spanwise direction. Moreover, a region of large spanwise velocity fluctuation on the downstream block is found between the vortex legs. Our results indicate that the reorientation of the arch vortex at increasing AOI is produced by the displacement of the different TKE regions and their interaction with the shear layers on the sides and top of the upstream and downstream blocks, respectively. There is also a close connection between the turbulent structure between the blocks and the wind gusts. The correlations among gust components were also studied, and it was found that in the near-wall region of the street the correlations between the streamwise and spanwise gusts \(R_{uv}\) were dominant for all four AOI cases. At higher wall-normal positions in the array, the \(R_{uw}\) correlation decreased with increasing AOI, whereas the \(R_{uv}\) coefficient increased as AOI increased, and at \({\textit{AOI}}=45^{\circ }\) all three correlations exhibited relatively high values of around 0.4.     

Parameterization of sheared entrainment in a well-developed CBL. Part I: Evaluation of the scheme through large-eddy simulations

The entrainment flux ratio A _e and the inversion layer (IL) thickness are two key parameters in a mixed layer model. A _e is defined as the ratio of the entrainment heat flux at the mixed layer top to the surface heat flux. The IL is the layer between the mixed layer and the free atmosphere. In this study, a parameterization of A _e is derived from the TKE budget in the firstorder model for a well-developed CBL under the condition of linearly sheared geostrophic velocity with a zero value at the surface. It is also appropriate for a CBL under the condition of geostrophic velocity remaining constant with height. LESs are conducted under the above two conditions to determine the coefficients in the parameterization scheme. Results suggest that about 43% of the shear-produced TKE in the IL is available for entrainment, while the shear-produced TKE in the mixed layer and surface layer have little effect on entrainment. Based on this scheme, a new scale of convective turbulence velocity is proposed and applied to parameterize the IL thickness. The LES outputs for the CBLs under the condition of linearly sheared geostrophic velocity with a non-zero surface value are used to verify the performance of the parameterization scheme. It is found that the parameterized A _e and IL thickness agree well with the LES outputs.     

Long-term variation of rainfall erosivity in Calabria (Southern Italy)

The changes in rainfall erosivity have been investigated using the rainfall erosivity factor (R) proposed for USLE by Wischmeier and Smith (R _W-S ) and some simplified indexes (the Fournier index modified by Arnoldus, F, a regional index spatial independent, R _Fr , and a regional index spatial dependent, R _Fs ) estimated by indirect approaches. The analysis has been carried out over 48 rainfall stations located in Calabria (Southern Italy) using data collected in the period 1936–2012 and divided in three sub-periods. The series of the erosivity indexes and of some precipitation variables have been analyzed for evidence of trends using standard methods. The simplified indexes suggested a general underestimation of the rainfall erosivity with respect to R _W-S . The mean underestimation ranged between 23 and 54 % for R _Fr and from 10 to 15 % for R _Fs . Both the sign and the magnitude of the trends were different for the different stations depending on the variable and sub-period considered. In general, the erosivity increased during the period 1936–1955 (1^st sub-period) and during the more recent sub-period (1992–2012, 3^rd sub-period), whereas it decreased during 1958–1977 (2^nd sub-period). The evidence of trends was generally higher for R _W-S than for R _Fr and R _Fs . Focusing on the most recent sub-period (3^rd sub-period), all the variables analyzed showed mainly increasing trends but with different magnitude. More particularly, R _W-S showed a mean increment of 29 %; F, R _Fr and R _Fs increased by 11, 15 and 18 %, respectively; the maximum intensity of 0.5-h precipitation increased by 5 %; and the annual precipitation increased by 22 %. Consequently, it remains difficult to define which precipitation variable plays the dominant role in the temporal variation of rainfall erosivity in the region. However, the overall results suggest that the indexes estimated by indirect procedures (F, R _Fr , and R _Fs ) should be used with caution for climate change analysis, despite they are used for practical purposes considering they are based on easily available information.     

Retrieval of eddy thermal conductivity in the weakly nonlinear Prandtl model for katabatic flows

Because the nonlinearity of actual physical processes can be expressed more precisely by the introduction of a nonlinear term, the weakly nonlinear Prandtl model is one of the most effective ways to describe the pure katabatic flow (no background flow). Features of the weak nonlinearity are reflected by two factors: the small parameter ε and the gradually varying eddy thermal conductivity. This paper first shows how to apply the Wentzel–Kramers–Brillouin (WKB) method for the approximate solution of the weakly nonlinear Prandtl model, and then describes the retrieval of gradually varying eddy thermal conductivity from observed wind speed and perturbed potential temperature. Gradually varying eddy thermal conductivity is generally derived from an empirical parameterization scheme. We utilize wind speed and potential temperature measurements, along with the variational assimilation technique, to derive this parameter. The objective function is constructed by the square of the differences between the observation and model value. The new method is validated by numerical experiments with simulated measurements, revealing that the order of the root mean squre error is 10^–2 and thus confirming the method’s robustness. In addition, this method is capable of anti-interference, as it effectively reduces the influence of observation error.