Vertical Wind Profiles Research Articles

Modulating local winds and turbulence around a single building obstacle with the obstruction of tall vegetation

Strategic vegetation placement can significantly alter airflow patterns and turbulence, fostering desired wind environments. By comparing scenarios where vegetation is placed upstream, downstream or absent (treeless) relative to a single building using large-eddy simulation, this study provides detailed insights into the sensitivity of flow dynamics to the positioning of the vegetation. Upstream vegetation more significantly disrupts the flow patterns around the building obstacle, altering vertical wind profiles and modifying wake circulations, compared to downstream vegetation. A small shear layer developed at the plant top for upstream vegetation markedly influences turbulent kinetic energy (TKE) on both the leeward and windward sides of the building, shifting the inflection point in vertical TKE profiles by up to 0.13H. By contrast, smaller tree-building separations lead to an effective merging of their aerodynamic profiles, whereas larger separations confine the streamwise breadth of turbulent fluxes, amplifying flux exchanges in the spanwise direction. Spectral analyses reveal that upstream vegetation consistently results in higher power spectral densities of the streamwise turbulence in the residential area than downstream vegetation. While small-scale spanwise velocity fluctuations are found to be comparably energetic at the building's windward side for upstream vegetation, the power becomes substantially concentrated on large-scale eddies in the building wake region, providing specific insights into modulating turbulent eddy motions within the residential zone.

How horizontal transport and turbulent mixing impact aerosol particle and precursor concentrations at a background site in the UAE

Abstract. The optical, physical, and chemical properties of aerosol particles have been previously studied in the United Arab Emirates (UAE), but there is still a gap in the knowledge of particle sources and in the horizontal and vertical transport of aerosol particles and their precursors in the area. To investigate how aerosol particle and SO2 concentrations at the surface responded to changes in horizontal and vertical transport, we used data from a 1-year measurement campaign at a background site where local sources of SO2 were expected to be minimal. The measurement campaign provided a combination of in situ measurements at the surface and the boundary layer evolution from vertical and horizontal wind profiles measured by a Doppler lidar. The diurnal structure of the boundary layer in the UAE was very similar from day to day, with a deep, well-mixed boundary layer during the day transitioning to a shallow nocturnal layer, with the maximum boundary layer height usually being reached around 14:00 local time. Both SO2 and nucleation-mode aerosol particle concentrations were elevated for surface winds coming from the east or western sectors. We attribute this to oil refineries located on the eastern and western coasts of the UAE. The concentrations of larger cloud condensation nuclei (CCN)-sized particles and their activation fraction did not show any clear dependence on wind direction, but the CCN number concentration showed some dependence on wind speed, with higher concentrations coinciding with the weakest surface winds. Peaks in SO2 concentrations were also observed despite low surface wind speeds and wind directions unfavourable for transport. However, winds aloft were much stronger, with wind speeds of 10 m s−1 at 1 km common at night and wind directions favourable for transport; surface-measured concentrations increased rapidly once these particular layers started to be entrained into the growing boundary layer, even if the surface wind direction was from a clean sector. These conditions also displayed higher nucleation-mode aerosol particle concentrations, i.e. new particle formation events occurring due to the increase in the gaseous precursor.

Proof of concept demonstration of a metastable helium fluorescence lidar for thermospheric wind and temperature measurements.

We report on the first, to the best of our knowledge, spectral measurements of terrestrial thermospheric metastable helium using ground-based lidar. By stimulating fluorescence of He(23S) at four closely spaced wavelengths within the He line around 1083 nm and measuring the lidar returns, we measured the He(23S) spectrum at 600 km, providing coarse constraints on the He(23S) temperature and vertical wind speed. This work serves as a proof of concept and precursor experiment for future, more powerful helium lidar systems capable of measuring vertical profiles of neutral wind and temperature in the upper terrestrial thermosphere.

Acoustic Propagation in the Near‐Surface Martian Atmosphere

AbstractThis work introduces a comprehensive model of sound propagation on Mars, in light of the recent operation of several microphones on the Martian surface. The main outcome of this work is an operational acoustic model capable of simulating the sound field created by any source, at any location on the Martian surface, at any time. Expanding on the result of previous work (Gillier et al., 2024, https://doi.org/10.1029/2023je008257), we use the parabolic equation method for sound propagation in order to obtain the overall sound field produced by a source, in a given atmospheric composition and state, and accounting for ground properties. The resulting model enables the study of acoustics on Mars, and has the potential also to be used to probe the properties of the Martian environment using acoustic measurements with known sources. We investigate the effects of the Martian ground and the vertical profile of temperature and wind, on sound propagation. We find that the ground has a minor effect on sound propagation, and the wind profile strongly influences sound propagation as on Earth. However, the midday near surface temperature profiles on Mars are shown to cause refraction, which generates non‐negligible acoustic losses that are an order of magnitude stronger than typical refraction‐related acoustic losses on Earth. We show that the effect of the Martian atmospheric turbulence is to slightly reduce the acoustic losses due to refraction. Finally, we apply our model to show that refraction and atmospheric turbulence have a negligible effect on the propagation of sound from Ingenuity to the Perseverance rover.

Atmospheric flow simulation over stationary waves with various steepness using WRF-LES

We investigated the capability of the Weather Research and Forecasting (WRF) model to perform large-eddy simulations (LESs) of the structure of the marine atmospheric boundary layer and turbulence exchange between the sea surface and the atmosphere. The WRF-LES model results were compared to results from Sullivan (2008) [1] for two cases: flow over a flat surface and over an idealized stationary sinusoidal surface. We also investigated the sensitivity of the WRF-LES model in simulating turbulence through the marine atmosphere under a range of horizontal grid spacings and wave steepnesses. All simulations were performed in a conventionally neutral atmosphere with a geostrophic wind of 5 m s−1, aligned with the direction of wave propagation. The simulated wind vertical profiles from the WRF-LES model exhibit similar behavior to those of Sullivan (2008) over both wavy and flat surfaces, but the differences in Sullivan’s vertical wind profiles for the two cases are considerably larger in comparison to those derived from WRF-LES. Furthermore, spectral analysis clearly shows the effect of horizontal grid spacing. The impact of wave steepness on turbulence distribution is clearly noticeable in the vertical profiles of resolved velocity variances and covariances. The impact is largest close to the surface, but it is important throughout the bulk of the boundary layer.

Balloon drift estimation and improved position estimates for radiosondes

Abstract. When comparing model output with historical radiosonde observations, it is usually assumed that a radiosonde has risen exactly above its starting point and has not been displaced by wind. This changed only relatively recently with the availability of Global Navigation Satellite System (GNSS) receivers aboard radiosondes in the late 1990s, but even then the balloon trajectory data were often not transmitted, although this information was the basis for estimating the wind in the first place. Depending on the conditions and time of year, radiosondes can sometimes drift a few hundred kilometres, particularly at the middle latitudes during the winter months. The position errors can lead to non-negligible representation errors when the corresponding observations are assimilated. This paper presents a methodology to compute changes in the balloon position during its vertical ascent, using only limited information, such as the vertical profile of wind contained in the historical observation reports. The sensitivity of the method to various parameters is investigated, such as the vertical resolution of the input data, the assumption about the vertical ascent speed of the balloon, and the departure of the surface of Earth from a sphere. The paper considers modern GNSS sonde data reports for validation, for which the full trajectory of the balloon is available, alongside the reported wind. Evaluation is also conducted by comparison with ERA5 and by conducting low-resolution data assimilation experiments. Overall, the results indicate that the trajectory of the radiosondes can be accurately reconstructed from original data of varying vertical resolutions and that the more accurate balloon position reduces representation errors and, in some cases, systematic errors.

Investigating the impact of coupling HARMONIE-WINS50 (cy43) meteorology to LOTOS-EUROS (v2.2.002) on a simulation of NO2 concentrations over the Netherlands

Abstract. Meteorological fields calculated by numerical weather prediction (NWP) models drive offline chemical transport models (CTMs) to solve the transport, chemical reactions, and atmospheric interaction over the geographical domain of interest. HARMONIE (HIRLAM ALADIN Research on Mesoscale Operational NWP in Euromed) is a state-of-the-art non-hydrostatic NWP community model used at several European weather agencies to forecast weather at the local and/or regional scale. In this work, the HARMONIE WINS50 (cycle 43 cy43) reanalysis dataset at a resolution of 0.025° × 0.025° covering an area surrounding the North Sea for the years 2019–2021 was coupled offline to the LOTOS-EUROS (LOng-Term Ozone Simulation-EURopean Operational Smog model, v2.2.002) CTM. The impact of using either meteorological fields from HARMONIE or from ECMWF on LOTOS-EUROS simulations of NO2 has been evaluated against ground-level observations and TROPOMI tropospheric NO2 vertical columns. Furthermore, the difference between crucial meteorological input parameters such as the boundary layer height and the vertical diffusion coefficient between the hydrostatic ECMWF and non-hydrostatic HARMONIE data has been studied, and the vertical profiles of temperature, humidity, and wind are evaluated against meteorological observations at Cabauw in The Netherlands. The results of these first evaluations of the LOTOS-EUROS model performance in both configurations are used to investigate current uncertainties in air quality forecasting in relation to driving meteorological parameters and to assess the potential for improvements in forecasting pollution episodes at high resolutions based on the HARMONIE NWP model.

Towards urban wind utilization: The spatial characteristics of wind energy in urban areas

Urban wind energy is gaining increasing recognition for its environmentally friendly nature, particularly in the context of energy shortages. Predicting wind speed in urban environments is challenging due to the varying roughness and drag caused by obstacles on the ground. The complexity of urban morphology significantly disrupts the wind field and hampers the utilization of wind energy. This study aimed to investigate the diverse impacts of highly heterogeneous urban environment on the urban wind energy by innovatively establishing spatially varying power law wind profiles. Firstly, long-term LIDAR observations were conducted to measure the vertical wind profile in the study area. Then, a high-resolution LCZ map was created and integrated into the Weather Research and Forecasting (WRF) model to reproduce the urban wind field, which was validated using the observational data. Finally, the spatial characteristics of the urban wind field and wind energy were analyzed based on 3 representative points and the urban LCZ categories. The results show that the average power law exponents of each LCZ classification ranges from 0.29 to 0.75, with the average power law exponent of the high-volume building area being larger than 0.6. The highest wind power density in the study area can reach 343.3W/m2 at 200m height. Further, significant negative correlation coefficients were found between wind power density and building height, as well as building density, with values of approximately −0.6 and −0.35, respectively. Overall, the high-volume-ratio built-up areas with higher roughness had a notable impact on the wind speed by increasing the power law exponents. These spatially varying wind profiles are expected to provide technical support for the utilization of urban wind energy.

Investigating the Role of the Low-Level Jet in Two Winters Severe Dust Rising in Southwest Iran

The dust storms with local and non-local dust sources mostly affect Khuzestan province in southwest (SW) Iran. In this study, the role of the low-level jet in the activation of the internal dust events in SW Iran during two severe dust cases was investigated. For this purpose, the fifth-generation ECMWF reanalysis for the global climate and weather (ERA5) data was used to identify the synoptic patterns and the low-level jet (LLJ) characteristics in the study area. Furthermore, the images of the moderate resolution imaging spectroradiometer (MODIS) sensor, the outputs of the hybrid single-particle Lagrangian integrated trajectory (HYSPLIT) model, and a weather research and forecasting model coupled with chemistry (WRF-Chem) were used to investigate the propagation and transport of the dust particles. The results of the synoptic analysis in both dust cases show the simultaneous occurrence of the divergence zone associated with cyclonic curvature in the subtropical jet stream (STJ) at 300 hPa, causing convergence at 925 hPa, upward motion, and the development of low surface pressure in SW Iran. Examining the vertical wind profile shows the existence of the maximum horizontal wind speeds of 975 to 875 hPa, along with the positive and negative shear below and above it, respectively, which emphasizes the existence of the LLJ and its role in local dust emission. The results of the comparison between the satellite images, WRF-Chem, and HYSPLIT model outputs show the formation and transportation of dust particles from the inner regions of Khuzestan in SW Iran. The horizontal dust surface distribution, vertical raised dust mass, and kinetic energy transfers are well simulated by the WRF-Chem model when LLJ broke at 09:00 to 12:00 UTC. The most important finding of this research is that, for the first time, the role of low-level jet is investigated in the activation of internal dust events in SW Iran.

Vertical profiles of temperature, wind, and turbulent fluxes across a deciduous forest over a slope observed with a UAV

To contribute to a better understanding of the dynamics of the atmosphere inside and above a forest, vertical profiles are flown with a remotely-controlled multicopter in the Steinkrug forest. This area is located over a slope in the Solling natural area in Lower Saxony (Germany), composed mostly of deciduous trees about 30 m tall. Fifteen vertical flights made near sunset between summer 2019 and spring 2020 were inspected from the surface to 100 m above ground level. These measurements provide information on the vertical structures of wind and temperature within and above the canopy, including the effects of shallow slope flows near the ground. Contrasting measurements downhill outside the forest were also made. The gathered data allow estimated profiles of the turbulent fluxes of sensible heat and momentum to be obtained by computing averages and fluctuations for layers of 5 m depth. A leaf area density profile in both leafy and leafless conditions could also be produced. The presence of a slope flow is inspected at both sites, and the applicability of existing theories is explored.

Applicability of the generalized wind profile model over mountainous forests.

The analytical equation based on Monin-Obukhov (M-O) similarity theory (i.e., wind profile equation) has been adopted since 1970s for using in the prediction of wind vertical profile over flat terrains, which is mature and accurate. However, its applicability over complex terrains remains unknown. This applicability signifies the accuracy of the estimations of aerodynamic parameters for the boundary layer of non-flat terrain, such as zero-displacement height (d) and aerodynamic roughness length (z0), which will determine the accuracy of frequency correction and source area analysis in calculating carbon, water, and trace gas fluxes based on vorticity covariance method. Therefore, the validation of wind profile model in non-flat terrain is the first step to test whether the flux model needs improvement. We measured three-dimensional wind speed data by using the Ker Towers (three towers in a watershed) at Qingyuan Forest CERN in the Mountainous Region of east Liaoning Province, and compared them with data from Panjin Agricultural Station in the Liaohe Plain, to evaluate the applicability of a generalized wind profile model based on the Monin-Obukhov similarity theory on non-flat terrain. The results showed that the generalized wind profile model could not predict wind speeds accurately of three flux towers separately located in different sites, indicating that wind profile model was not suitable for predicting wind speeds in complex terrains. In the leaf-off and leaf-on periods, the coefficient of determination (R2) between observed and predicted wind speeds ranged from 0.12 to 0.30. Compared to measured values, the standard error of the predicted wind speeds was high up to 2 m·s-1. The predicted wind speeds were high as twice as field-measured wind speed, indicating substantial overestimation. Nevertheless, this model correctly predicted wind speeds in flat agricultural landscape in Panjin Agricultural Station. The R2 between observed wind speeds and predicted wind speed ranged from 0.90 to 0.93. The standard error between observed and predicted values was only 0.5 m·s-1. Results of the F-test showed that the root-mean-square error of the observed and predicted wind speeds in each secondary forest complex terrain was much greater than that in flat agricultural landscape. Terrain was the primary factor affecting the applicability of wind profile model, followed by seasonality (leaf or leafless canopy). The wind profile model was not applicable to the boundary-layer flows over forest canopies in complex terrains, because the d was underestimated or both the d and z0 were underestimated, resulting in inaccurate estimation of aerodynamic height.

Improvement of Stable Atmospheric Boundary Simulation with High-Spatiotemporal-Resolution Nudging over the North China Plain

The WRF model often struggles to accurately replicate specific characteristics of the atmospheric boundary layer, particularly under highly stable conditions. In this study, we reconstructed an OBS-nudging module using meteorological data with high spatiotemporal resolution, then coupled it in the WRF model (WRF-OBS) to improve stable boundary layer (SBL) simulation over the North China Plain (NCP). The results showed that WRF-OBS improved the simulation of SBL characteristics and reduced the deviation from observations significantly. The correlations (R2) between WRF-OBS simulations and observations of 2 m temperature, relative humidity, and 10 m wind speed at 460 stations across the NCP were 0.72, 0.56, and 0.75, respectively, which were much higher than the values for results from the unassimilated WRF model (WRF-BS). The simulated vertical profiles of temperature, relative humidity, and wind were generally consistent with observations at Pingyuan station. The meteorological factors which caused heavy air pollution was also investigated based on WRF-OBS simulation. The SBL characteristics obtained from WRF-OBS showed that light wind persisted over the NCP region during the period of heavy pollution, and Pingyuan was affected by warm and humid air. Vertically, the persistent temperature inversion at Pingyuan station was one of the main drivers of the heavy pollution. The WRF-OBS simulation captured the characteristics of the two temperature inversion layers very well. The two inversion layers covered the NCP, with a horizontal scale of approximately 200 km, and created very stable conditions, preventing the vertical diffusion of pollutants and maintaining high PM2.5 concentrations.

Profiling the molecular destruction rates of temperature and humidity as well as the turbulent kinetic energy dissipation in the convective boundary layer

Abstract. A simultaneous deployment of Doppler, temperature, and water-vapor lidars is able to provide profiles of molecular destruction rates and turbulent kinetic energy (TKE) dissipation in the convective boundary layer (CBL). Horizontal wind profiles and profiles of vertical wind, temperature, and moisture fluctuations are combined, and transversal temporal autocovariance functions (ACFs) are determined for deriving the dissipation and molecular destruction rates. These are fundamental loss terms in the TKE as well as the potential temperature and mixing ratio variance equations. These ACFs are fitted to their theoretical shapes and coefficients in the inertial subrange. Error bars are estimated by a propagation of noise errors. Sophisticated analyses of the ACFs are performed in order to choose the correct range of lags of the fits for fitting their theoretical shapes in the inertial subrange as well as for minimizing systematic errors due to temporal and spatial averaging and micro- and mesoscale circulations. We demonstrate that we achieve very consistent results of the derived profiles of turbulent variables regardless of whether 1 or 10 s time resolutions are used. We also show that the temporal and spatial length scales of the fluctuations in vertical wind, moisture, and potential temperature are similar with a spatial integral scale of ≈160 m at least in the mixed layer (ML). The profiles of the molecular destruction rates show a maximum in the interfacial layer (IL) and reach values of ϵm≃7×10-4 g2 kg−2 s−1 for mixing ratio and ϵθ≃1.6×10-3 K2 s−1 for potential temperature. In contrast, the maximum of the TKE dissipation is reached in the ML and amounts to ≃10-2 m2 s−3. We also demonstrate that the vertical wind ACF coefficient kw∝w′2‾ and the TKE dissipation ϵ∝w′2‾3/2. For the molecular destruction rates, we show that ϵm∝m′2‾w′2‾1/2 and ϵθ∝θ′2‾w′2‾1/2. These equations can be used for parameterizations of ϵ, ϵm, and ϵθ. All noise error bars are derived by error propagation and are small enough to compare the results with previous observations and large-eddy simulations. The results agree well with previous observations but show more detailed structures in the IL. Consequently, the synergy resulting from this new combination of active remote sensors enables the profiling of turbulent variables such as integral scales, variances, TKE dissipation, and the molecular destruction rates as well as deriving relationships between them. The results can be used for the parameterization of turbulent variables, TKE budget analyses, and the verification of large-eddy simulations.

Vertical Wind Profiles in the Mesosphere and Lower Thermosphere Driven by Meteor Radar and Ionospheric Connection Explorer Observations Over the Korean Peninsula

AbstractMeteor radar observations provide wind data ranging from 80 to 100 km altitude, while the Michaelson Interferometer for Global High‐resolution Thermospheric Imaging (MIGHTI) onboard the Ionospheric Connection Explorer satellite offers wind data above 90 km altitude. This study aims to generate wind profiles in the mesosphere and lower thermosphere by combining the winds derived from meteor radar and MIGHTI observations over the Korean Peninsula from January 2020 to December 2021. The wind profiles derived from the two instruments are continuous at night, but they show discrepancies during the day. The atomic oxygen 557.7 nm (green line) emission intensity measured by MIGHTI peaks at approximately 100 km during the day and 94 km at night. The vertical gradient of the airglow volume emission rate is more pronounced during the day. These differences can cause day‐night differences in the MIGHTI wind retrieval accuracy, potentially leading to discrepancies during the day.

Research on attitude correction algorithm for mobile wind lidars

The laser wind measurement technology is remarkable for detecting clear-sky wind fields. The Doppler beam swinging algorithm for wind lidars has been developed to obtain vertical wind profiles based on fixed observation methods. However, the Doppler frequencies are superposed due to the self-motions of lidars caused by carrier motions when lidars are used on motion carrier platforms. Meanwhile, the emission directions of laser beams are uncertain due to changes in carriers’ motion directions and tilts. Thus, a new wind measurement correction model must be studied with lidar attitudes. This study considers the influences of the motion velocities, the carrier’s tilt angles, and the laser beams’ yaw angles at the 0° azimuth angle on the measured results under lidar motions, a correction model of motion attitudes for mobile wind lidars was designed. Sensitivity simulation tests for motion attitude parameters were carried out, and the influences of different attitude parameters of the carrier on the measured results were investigated to evaluate and verify the effects of the correction model. Results indicated that the wind measurement correction model could correct data errors caused by the carrier’s motion and tilts. The motion velocities, carrier directions, and the yaw angles of the laser beams at the 0° azimuth angle had an essential influence on the wind velocity measurements. Besides, the carrier’s pitch angles and the roll angles, which did not influence the wind velocity measurements, only affected the altitudes of the wind field data. Furthermore, the pitch angles exerted more significant influences on the data altitudes than the roll angles.

A novel empirical model for vertical profiles of downburst horizontal wind speed

AbstractThis study proposes an empirical model for preliminary wind‐resist design of downburst flow. Existing empirical models were compared with field data and found to underpredict horizontal wind speed below the height corresponding to the maximum radial velocity, due to the neglect of viscous effects and the evolution of vertical wind profiles along radial direction. To address these deficiencies, semi‐empirical piecewise functions including wall shear effects in the local turbulent boundary layer and interpolation functions were proposed to improve the accuracy of existing models. The wind profile based on Coles' theory was found to agree well with field data, with the parabola interpolation function being the most desirable. Using the proposed method, the vertical profile of horizontal wind speed at different local radial locations can be predicted for wind resist design given the inlet wind speed of the downburst flow. Overall, this model improves upon existing empirical models and allows for more accurate wind‐resist design.

Comparisons and quality control of wind observations in a mountainous city using wind profile radar and the Aeolus satellite

Abstract. Observations of the vertical wind profile in Chongqing, a typical mountainous city in China, are important, but they are sparse and have low resolution. To obtain more wind profile data, this study matched the Aeolus track with ground-based wind observation sites in Chongqing in 2021. Based on the obtained results, verification and quality control studies were conducted on the wind observations of a wind profile radar (WPR) with radiosonde (RS) data, and a comparison of the Aeolus Mie-cloudy and Rayleigh-clear wind products (Aeolus winds measured in cloudy and aerosol-rich atmospheric conditions from Mie-channel-collected data and winds measured in clear-air conditions from Rayleigh-collected data) with WPR data was then performed. The conclusions can be summarized as follows: (1) a clear correlation between the wind observations of WPR and RS was found, with a correlation coefficient (R) of 0.71. Their root mean square deviation increased with height but decreased at heights between 3 and 4 km. (2) After quality control using Gaussian filtering (GF) and empirical orthogonal function construction (EOFc; G=87.23 %) of the WPR data, the R between the WPR and RS reached 0.83 and 0.95, respectively. The vertical distribution showed that GF could better retain the characteristics of WPR wind observations but with limited improvement in decreasing deviations, whereas EOFc performed better in decreasing deviations but considerably modified the original characteristics of the wind field, especially regarding intensive vertical wind shear in strong convective weather processes. (3) In terms of the differences between the Aeolus and WPR data, 56.0 % and 67.8 % deviations were observed within ±5 m s−1 for Rayleigh-clear and Mie-cloudy winds (Aeolus winds measured in cloudy and aerosol-rich atmospheric conditions from Mie-channel-collected data and winds measured in clear-air conditions from Rayleigh-collected data) vs WPR winds, respectively. Vertically, large mean differences of both Rayleigh-clean and Mie-cloudy winds versus WPR winds appeared below 1.5 km, which is attributed to the prevailing quiet and small winds within the boundary layer in Chongqing; in this case the movement of molecules and aerosols is mostly affected by irregular turbulence. Additionally, large mean differences at a height range between 4 and 8 km for Mie-cloudy versus WPR winds may be related to the high content of cloud liquid water in the middle troposphere of Chongqing. (4) The differences in both Rayleigh-clear and Mie-cloudy versus WPR winds had changed. Deviations of 58.9 % and 59.6 % were concentrated within ±5 m s−1 for Rayleigh-clear versus WPR winds with GF and EOFc quality control, respectively. In contrast, 69.1 % and 70.2 % of deviations appeared within ±5 m s−1 for Rayleigh-clear versus WPR and EOFc WPR winds, respectively. These results shed light on the comprehensive applications of multi-source wind profile data in mountainous cities or areas with sparse ground-based wind observations.

Overview of Wind Parameter Radar Applications Around the World

Relying on the development of radio electronics and computer technology, a direction has emerged in radar meteorology related to the development of methods and means for remote measurement of the vertical wind profile in the atmosphere by determining the Doppler frequency shift of the reflected signal from meteorological irregularities. The article provides a brief overview of the use of wind parameter radars (profilers) that operate in the decimeter wavelength range (400–1700 MHz), since it is in this range that wind parameters are most fully determined under absolutely clear sky conditions. The distribution of profilers across countries and continents is considered, their main technical characteristics are given, and the current situation in the Russian Federation is described.

Evaluation of the GRAMM/GRAL model for high-resolution wind fields in Heidelberg, Germany

Independent verification of mitigation efforts for climate and air quality action in cities relies on inferring emissions from atmospheric concentration measurements. As emissions are dispersed in the atmosphere before they reach an instrument, the quantitative estimation of emissions requires an understanding of the atmospheric transport and associated uncertainties. In this study, we analyse the catalogue of steady-state flow fields generated by the Graz Mesoscale Model (GRAMM) coupled to the Graz Lagrangian Model (GRAL) for an entire year in Heidelberg, Germany. We use a loss function for the wind field selection, which assigns a best-matching catalogue entry to any given hour by exploiting observation data. We introduce a new loss function which finds an optimal balance between differences in wind speed and wind direction. We evaluate the performance of the model based on 15 meteorological measurement sites, of which 14 are in the inner high-resolution and building-resolving GRAL domain (12.5 km × 12.5 km, 10 m resolution). Performance metrics include mean bias (MB) and root mean square errors (RMSEs) of simulated and observed wind speed and wind direction for all individual stations. On average, we find a mean underestimation of wind speed of 0.14 ms−1 corresponding to about 7 % of the mean wind speed and a mean RMSE of 1.03 ms−1. For wind direction, a mean overall bias smaller than 1° is achieved, but individual stations show larger biases (mean absolute bias: 37°), especially at stations where wind speeds are low on average. Evaluation benchmarks for mean biases of wind direction and wind speed of mesoscale models provided by the European Environmental Agency (EEA) are met at 11 and 14 out of 15 stations at low measurement heights, respectively. Recently suggested extended benchmarks for complex terrain are met at almost all stations. Additionally, for the first time, we analyse the model's ability to simulate the vertical wind profile and we analyse the benefit of implementing a wind profile measurement into the process. We find that the model does not fully capture the vertical profile in our setting. We further study the required measurement network size and find that a high number (> 6) of meteorological stations improves the selection of flow fields over the entire GRAL domain substantially. The conducted comprehensive analysis of the wind fields in the GRAL domain are the basis for detailed quantitative analysis of greenhouse gas and air pollutant emissions using the GRAMM/GRAL modelling framework.

Comparison of the LBM snowdrift model output with the observation results

In this work, snowdrift experiments which are equivalent to one drifting snow event are performed by the snowdrift model. The model consisted of the computational fluid dynamics part of the large-eddy simulation with the lattice Boltzmann method and the drifting snow part of the conventional advection algorithm for representative Lagrangian particles. The observed vertical wind profile of a 4 h drifting snow event in Teshikaga Town was used as the inflow boundary conditions in the model to compare the results of the snowdrift estimated by the model and the observed snowdrift distribution. Parallelization enabled us to simulate the snowdrift distribution in a realistic domain and on the time scale of a single drifting snow event. We demonstrated that the upgraded model could quantitatively reproduce the height and position of the observed snowdrift along the center of a three-dimensional fence.