Search

Lower tropospheric thermal structure greatly affects atmospheric boundary-layer (ABL) stability and mixing processes with the free troposphere. In particular, in polluted urban zones, ABL stratification becomes a key variable in air quality research. This study focuses on generating a climatology (1990–2017) of the seasonal variability of ABL thermal structure in Mexico City by way of radiosonde analysis. Thermal inversion intensity and frequency are shown to be greater during winter and spring, a behaviour which coincides with greater pollutant concentrations. Higher concentrations are found during the dry season (November to May) than during the rainy months. In addition, significantly higher than normal surface pollutant concentrations are found on days with simple thermal inversion layers as well as during multilayer inversion days. Furthermore, stable layers, determined by potential temperature, are found throughout the year but more frequently during winter, whereas stable layers based on the virtual potential temperature prevail all year. In regions of complex terrain, such multiple stable layers have also been identified by previous authors. Additionally, the most unstable surface layers (in which the bulk Richardson number ( $${Ri}_{\mathrm {B}}$$ ) is small) develop during the rainy season, whereas during winter there are more levels in the vertical column with higher $${Ri}_{\mathrm {B}}$$ values. Although the Mexico City ABL and pollution episodes have been widely studied, this represents the first long-term investigation to consider the thermal stability of the ABL. Therefore, the present study provides a baseline for further research employing different observational techniques and high-resolution numerical models.

An open-pool type research reactor is designed and operated considering the accessibility around the pool top area to enhance the reactor utilization. The reactor structure assembly is placed at the bottom of the pool and filled with water as a primary coolant for the core cooling and radiation shielding. Most radioactive materials are generated from the fuel assemblies in the reactor core and circulated with the primary coolant. If the primary coolant goes up to the pool surface, the radiation level increases around the working area near the top of the pool. Hence, the hot water layer is designed and formed at the upper part of the pool to suppress the rising of the primary coolant to the pool surface. The temperature gradient is established from the hot water layer to the primary coolant. As this temperature gradient suppresses the circulation of the primary coolant at the upper region of the pool, the radioactive primary coolant rising up directly to the pool surface is minimized. Water mixing between these layers is reduced because the hot water layer is formed above the primary coolant with a higher temperature. The radiation level above the pool surface area is maintained as low as reasonably achievable since the radioactive materials in the primary coolant are trapped under the hot water layer. The key to maintaining the stable hot water layer and keeping the radiation level low on the pool surface is to have a stable flow of the primary coolant. In the research reactor with a downward core flow, the primary coolant is dumped into the reactor pool and goes to the reactor core through the flow guide structure. Flow fields of the primary coolant at the lower region of the reactor pool are largely affected by the dumped primary coolant. Simple, circular, and duct type discharge headers are designed to control the flow fields and make the primary coolant flow stable in the reactor pool. In this research, flow fields of the primary coolant and hot water layer are numerically simulated in the reactor pool. The heat transfer rate, temperature, and velocity fields are taken into consideration to determine the formation of the stable hot water layer and primary coolant flow. The bulk Richardson number is used to evaluate the stability of the flow field. A duct type discharge header is finally chosen to dump the primary coolant into the reactor pool. The bulk Richardson number should be higher than 2.7 and the temperature of the hot water layer should be 1 °C higher than the temperature of the primary coolant to maintain the stability of the stratified thermal layer.

Climate models still have deficits in reproducing the surface energy and momentum budgets in Arctic regions. One of the reasons is that currently used transfer coefficients occurring in parameterizations of the turbulent fluxes are based on stability functions derived from measurements over land and not over sea ice. An improved parameterization is developed using the Monin–Obukhov similarity theory (MOST) and corresponding stability functions that reproduce measurements over sea ice obtained during the Surface Heat Budget of the Arctic Ocean (SHEBA) campaign. The new stability functions for the stable surface layer represent a modification of earlier ones also based on SHEBA measurements. It is shown that the new functions are superior to the former ones with respect to the representation of the measured relationship between the MOST stability parameter and the bulk Richardson number. Nevertheless, the functions fulfill the same criteria of applicability as the earlier functions and contain, as an extension, a dependence on the neutral-limit turbulent Prandtl number. Applying the new functions we develop an efficient noniterative parameterization of the near-surface turbulent fluxes of momentum and heat with transfer coefficients as a function of the bulk Richardson number (Rib) and roughness parameters. A hierarchy of transfer coefficients is recommended for weather and climate models. They agree better with SHEBA data for strong stability (Rib > 0.1) than previous parameterizations and they agree well with those based on the Businger–Dyer functions in the range Rib ≤ 0.1.

Four cases of an overlying inversion imposed on a stable boundary layer are investigated, extending the earlier work of Hancock and Hayden (Boundary-Layer Meteorol 168:29–57, 2018), where no inversion was imposed. The inversion is imposed to one or other of two depths within the layer: midway or deep. Four cases of changed surface condition are also investigated, and it is seen that the surface and imposed conditions behave independently. A change of imposed inversion condition leaves the bottom 1/3 of the layer almost completely unaffected; a change of the surface condition leaves the top 2/3 unaffected. Comparisons are made against two sets of local-scaling systems over the full height of the boundary layer. Both show some influence of the inversion condition. The surface heat flux and the reduction in surface shear stress, and hence the ratio of the boundary-layer height to surface Obukhov length, are determined by the temperature difference across the surface layer (not the whole layer), bringing all cases together in single correlations as functions of a surface-layer bulk Richardson number.

Diurnal variations in the planetary boundary layer height (PBLH) at different latitudes over different surface characteristics are described, based on 45 years (1973−2017) of radiosonde observations. The PBLH is determined from the radiosonde data by the bulk Richardson number (BRN) method and verified by the parcel method and the potential temperature gradient method. In general, the BRN method is able to represent the height of the convective boundary layer (BL) and neutral residual layer cases but has relatively large uncertainty in the stable BL cases. The diurnal cycle of the PBLH over land is quite different from the cycle over ocean, as are their seasonal variations. For stations over land, the PBLH shows an apparent diurnal cycle, with a distinct maximum around 15:00 LT, and seasonal variation, with higher values in summer. Compared with the PBLH over land, over oceans the PBLH diurnal cycles are quite mild, the PBLHs are much lower, and the seasonal changes are less pronounced. The seasonal variations in the median PBLH diurnal cycle are positively correlated with the near-surface temperature and negatively correlated with the near-surface relative humidity. Finally, although at most latitudes the daytime PBLH exhibits, over these 45 years, a statistically significant increasing trend at most hours between 12:00 LT and 18:00 LT over both land and ocean, there is no significant trend over either land or ocean in the nighttime PBLH for almost all the studied latitudes.

Ilustraciones, mapas Esta investigación caracterizó la estructura térmica y dinámica de la capa límite nocturna para el año 2017 en el Valle de Aburrá–Colombia, un área urbana tropical con topografía compleja. Se utilizaron registros provenientes de un radiómetro de microondas, un radar perfilador de vientos, un ceilómetro, estaciones meteorológicas y el modelo WRF-ARW acoplado al modelo de parametrización urbana SLUCM. Este último fue ejecutado con seis esquemas de parametrización distintos de la capa límite atmosférica, para 33 noches distribuidas en el periodo de estudio. Un análisis exploratorio fue ejecutado para identificar procesos espacio-temporales usando variables de estado como los vientos, el número de Richardson aproximado, la temperatura potencial virtual y la intensidad de retrodispersión. Mediante un análisis de sensibilidad de los registros se encontró que el espesor de la capa límite nocturna corresponde a un número de Richardson crítico de 0,5. Además se evaluó el modelo para las horas de la noche y se encontró un desempeño aceptable del esquema de parametrización MYNN. También se identificaron patrones de circulación asociados a un jet de bajo nivel, inversiones térmicas, vientos catabáticos y acoplamiento de los vientos alisios con los vientos orográficos. Se observó que los trimestres junio-julio-agosto y septiembre-octubre-noviembre son más estables dinámicamente, mientras que los trimestres de diciembre-enero-febrero y marzo-abril-mayo lo son más desde el punto de vista estático. Finalmente, se concluye que los espesores de la capa límite nocturna en el Valle de Aburrá son relativamente bajos, con condiciones de velocidades significantes al principio de la noche, pero al final de la noche con velocidades muy cercanas a cero y con una estabilidad atmosférica cada vez fortaleciéndose más por el enfriamiento radiativo. (texto tomado de la fuente) This research characterized the thermal structure and dynamics of the nocturnal boundary layer for the year 2017 in the Aburrá Valley-Colombia, a tropical urban area with complex topography. Records from a microwave radiometer, a wind profiler radar, a ceilometer, meteorological stations, and the WRF-ARW model coupled to the SLUCM urban parameterization model were used. This last one was run with six different atmospheric boundary layer parameterization schemes, for 33 nights distributed in the study period. An exploratory analysis was performed to detect spatio-temporal processes using state variables such as winds, bulk Richardson number, virtual potential temperature and backscatter. Through a sensitivity analysis of the records, it was found that the thickness of the nocturnal boundary layer corresponds to a critical Richardson number of 0,5. In addition, the model was evaluated during night hours and an acceptable performance of the MYNN parameterization scheme was found. Circulation patterns associated with a low-level jet, thermal inversions, katabatic winds and coupling of trade winds with orographic winds were also identified. It was observed that the quarters of june-july-august and september-october-november are more dynamically stable, and whereas those of the december-january-february and march-april-may are more statically stable. Finally, it is concluded that the thicknesses of the nocturnal boundary layer in the Aburrá Valley are relatively low, with significant velocities at the beginning of the night, but at the end of the night with velocities very close to zero and with atmospheric stability becoming increasingly stronger due to radiative cooling. Maestría Magíster en Ingeniería - Recursos Hidráulicos Meteorología urbana y de montañas Área Curricular de Medio Ambiente

Atmospheric stability has a significant effect on wind shear and turbulence intensity, and these variables, in turn, have a direct impact on wind power production and loads on wind turbines. It is therefore important to know how to characterise atmospheric stability in order to make better energy yield estimation in a wind farm. Based on the research-grade meteorological mast at Alaiz (CENER's test site in Navarre, Spain) named MP5, this work compares and evaluates different instrument set-ups and methodologies for stability characterisation, namely the Obukhov parameter, measured with a sonic anemometer, and the bulk Richardson number based on two temperature and one wind speed measurement. The methods are examined considering their theoretical background, implementation complexity, instrumentation requirements, and practical use in connection to wind energy applications. The sonic method provides a more precise local measurement of stability while the bulk Richardson is a simpler, robust and cost-effective technique to implement in wind assessment campaigns. Using the sonic method as a benchmark, it is shown that to obtain reliable bulk Richardson measurements in onshore sites it is necessary to install one of the temperature sensors close to the ground where the temperature gradient is stronger.

Ceilometer measurements of aerosol backscatter profiles have been widely used to provide continuous PBLHT estimations. To investigate the robustness of ceilometer-estimated PBLHT under different atmospheric conditions, we compared ceilometer- and radiosonde-estimated PBLHTs using long term U.S. Department of Energy (DOE) Atmospheric Radiation Measurement (ARM) ceilometer and balloon-borne sounding data at three ARM fixed-location atmospheric observatories and from three ARM mobile observatories deployed around the world for various field campaigns, which cover from Tropics to Polar regions and over both ocean and land surfaces. Statistical comparisons of ceilometer-estimated PBLHTs from the Vaisala CL31 ceilometer data with radiosonde-estimated PBLHTs from the ARM PBLHT-SONDE Value-added Product (VAP) are performed under different atmospheric conditions including stable and unstable atmospheric boundary layer, low-level cloud-free, and cloudy conditions at these ARM observatories. Under unstable atmospheric boundary layer conditions, good comparisons are found between ceilometer- and radiosonde-estimated PBLHTs at ARM low- and mid-latitude land observatories. However, it is still challenging to obtain reliable PBLHT estimations over ocean surfaces even using radiosonde data. Under stable atmospheric boundary layer conditions, ceilometer- and radiosonde-estimated PBLHTs have weak correlations. Among different PBLHT estimations utilizing the Heffter, the Liu-Liang, and the bulk Richardson number methods in the ARM PBLHT-SONDE VAP, ceilometer-estimated PBLHTs have better comparisons with the Liu-Liang method under unstable and with the bulk Richardson number method under stable atmospheric boundary layer conditions.

Publisher Copyright: © Copyright © 2021 Huo, Wang, Paasonen, Liu, Tang, Ma, Petaja, Kerminen and Kulmala. Boundary layer height (BLH) is an important parameter in climatology and air pollution research, especially in urban city. We calculated the BLH with a bulk Richardson number (Ri) method over urban cities of China during 1980–2018 using European Centre for Medium-Range Weather Forecasts (ECMWF) ERA-interim data after carefully validation with sounding data obtained from two meteorology stations in eastern China during 2010–2018. The values of BLH between these two types of data have correlation coefficients in the range of 0.65–0.87, which indicates that it is reasonable to analyze long-term trends of the BLH from ERA data sets. Using ERA-interim calculated BLH, we found that there is an increasing trend in the daytime BLH in most cities of eastern China, particularly during the spring season. A correlation analysis between the BLH and temperature, wind speed, relative humidity and visibility revealed that the variability in meteorological parameters, as well as in aerosol concentrations over highly polluted eastern China, play important roles in the development of the BLH.

We currently have a limited understanding of the spatial and temporal variability in vertically stratified atmospheric layers over Namibia and the southeast Atlantic (SEA) Ocean. Stratified layers are relevant to the transport and dilution of local and long-range transported atmospheric constituents. This study used eleven years of global positioning system radio occultation (GPS-RO) signal refractivity data (2007–2017) over Namibia and the adjacent ocean surfaces, and three years of radiosonde data from Walvis Bay, Namibia, to study the character and variability in stratified layers. From the GPS-RO data and up to a height of 10 km, we studied the spatial and temporal variability in the point of minimum gradient in refractivity, and the temperature inversion height, depth and strength. We also present the temporal variability of temperature inversions and the boundary layer height (BLH) from radiosondes. The BLH was estimated by the parcel method, the top of a surface-based inversion, the top of a stable layer identified by the bulk Richardson number (RN), and the point of minimum gradient in the refractivity (for comparison with GPS-RO data). A comparison between co-located GPS-RO to radiosonde temperature profiles found good agreement between the two, and an average underestimation of GPS-RO to radiosonde temperatures of −0.45 ± 1.25 °C, with smaller differences further from the surface and with decreasing atmospheric moisture content. The minimum gradient (MG) of refractivity, calculated from these two datasets were generally in good agreement (230 ± 180 m), with an exeption of a few cases when differences exceeded 1000 m. The surface of MG across the region of interest was largely affected by macroscale circulation and changes in atmospheric moisture and cloud, and was not consistent with BLH(RN). We found correlations in the character of low-level inversions with macroscale circulation, radiation interactions with the surface, cloud cover over the ocean and the seasonal maximum in biomass burning over southern Africa. Radiative cooling on diurnal scales also affected elevated inversions between 2.5 and 10 km, with more co-occurring inversions observed at night and in the morning. Elevated inversions formed most frequently over the subcontinent and under subsidence by high-pressure systems in the colder months. Despite this macroscale influence peaking in the winter, the springtime inversions, like those at low levels, were strongest.

Recent research has demonstrated that hidden Markov model (HMM) analysis is an effective tool to classify atmospheric observations of the stably stratified nocturnal boundary layer (SBL) into weakly stable (wSBL) and very stable (vSBL) regimes. Here we consider the development of explicitly stochastic representations of SBL regime dynamics. First, we analyze whether HMM-based SBL regime statistics (the occurrence of regime transitions, subsequent transitions after the first, and very persistent nights) can be accurately represented by “freely running” stationary Markov chains (FSMCs). Our results show that despite the HMM-estimated regime statistics being relatively insensitive to the HMM transition probabilities, these statistics cannot all simultaneously be captured by a FSMC. Furthermore, by construction a FSMC cannot capture the observed non-Markov regime duration distributions. Using the HMM classification of data into wSBL and vSBL regimes, state-dependent transition probabilities conditioned on the bulk Richardson number (RiB) or the stratification are investigated. We find that conditioning on stratification produces more robust results than conditioning on RiB. A prototype explicitly stochastic parameterization is developed based on stratification-dependent transition probabilities, in which turbulence pulses (representing intermittent turbulence events) are added during vSBL conditions. Experiments using an idealized single-column model demonstrate that such an approach can simulate realistic-looking SBL regime dynamics.

An inflow generation method for large-eddy simulations of wall-bounded turbulent flows is presented. The method builds upon the work of Munoz-Esparza et al. (“A Stochastic Perturbation Method to Ge...

To investigate the evolution of the nocturnal boundary layer (NBL) and its impacts on the vertical distributions of pollutant particulates, a combination of in situ observations from a large tethered balloon, remote sensing instruments (aerosol lidar and Doppler wind lidar) and an atmospheric environment-monitoring vehicle were utilized. The observation site was approximately 100 km southwest of Beijing, the capital of China. Results show that a considerable proportion of pollutant particulates were still suspended in the residual layer (RL) (e.g., the nitrate concentration reached 30 μg m−3) after sunset. The NBL height calculated by the aerosol lidar was closer to the top of the RL before midnight because of the pollutants stored aloft in the RL and the shallow surface inversion layer, after midnight, the NBL height was more consistent with the top of the surface inversion layer. As the convective mixing layer gradually became established after sunrise the following day, the pollutants stored in the nocturnal RL of the preceding night were entrained downward into the mixing layer. The early morning PM2.5 concentration near 700 m in the RL on 20 December decreased by 83% compared with the concentration at 13:34 on 20 December at the same height. The nitrate concentration also decreased significantly in the RL, and the mixing down of nitrate from the RL could contribute about 37% to the nitrate in the mixing layer. Turbulence activities still existed in the RL with the bulk Richardson number (Rb) below the threshold value. The corresponding increase in PM2.5 was likely to be correlated with the weak turbulence in the RL in the early morning.

The behaviour of the turbulent Prandtl number ( Pr t ) for buoyancy-affected flows near a vertical surface is investigated as an extension study of Gibson & Leslie, Int. Comm. Heat Mass Transfer, Vol. 11, pp. 73–84 (1984). By analysing the location of mean velocity maxima in a differentially heated vertical planar channel, we identify an infinity anomaly for the eddy viscosity ν t and the turbulent Prandtl number Pr t , as both terms are divided by the mean velocity gradient according to the standard definition, in vertical buoyant flow. To predict the quantities of interest, e.g. the Nusselt number, a machine learning framework via symbolic regression is used with various cost functions, e.g. the mean velocity gradient, with the aid of the latest direct numerical simulation (DNS) dataset for vertical natural and mixed convection. The study has yielded two key outcomes: ( i ) the new machine learnt algebraic models, as the reciprocal of Pr t , successfully handle the infinity issue for both vertical natural and mixed convection; and ( i i ) the proposed models with embedded coordinate frame invariance can be conveniently implemented in the Reynolds-averaged scalar equation and are proven to be robust and accurate in the current parameter space, where the Rayleigh number spans from 10 5 to 10 9 for vertical natural convection and the bulk Richardson number R i b is in the range of 0 and 0.1 for vertical mixed convection.

Gustiness is examined for the wind speed, fluctuations, turbulence intensities and fluxes for a real urban topography. Using large-eddy simulation (LES) and an ensemble sampling approach, which allows a more comprehensive characterisation of the urban morphological features, a wide range of boundary-layer stabilities is considered: the bulk Richardson number, Rb ​∈ ​[−0.41, 0.82]. Ratios of the proposed gustiness statistics, G , over the conventional time-averaged flow and turbulence statistics are maximised for z/Have ​≲ ​1 (where Have is the mean building height). The strong linear scaling of G with the plan-area index (λp) for neutral stratification is found to persist for stably- and unstably-stratified flows (R2 ~ 0.8). By contrast, the non-dimensionalised building-height variability, σH/Have, and the effective frontal-area index, λ f ^ ≡ λ f H a v e / σ H , are argued to be of more appropriateness as scaling parameters for G compared to their original forms, σH and λf. While the sensitivity of G to Rb is well defined at greater heights, the influence of surface inhomogeneity may be strong enough to oppose the effect of thermal stratification in the lower surface layers. Qualitative differences in the sensitivities to the boundary-layer stability are narrowly distinguishable amongst the zeroth-, first- and second-order gustiness statistics. The results are relevant to the understanding of urban wind hazards in thermally-stratified boundary layers.

In this study, a climatology of two key boundary layer features, the Planetary Boundary Layer Height (PBLH) and the wind field over Greece is derived. The climatology is based on daily soundings collected in Athens, Thessaloniki and Heraklion and spanning a 32-year period. The PBLH is estimated using a method based on the gradient of potential temperature and a method based on the bulk Richardson number. The wind field is analyzed by calculating the wind shear and the turning angle of the wind vector between the surface and the top of the boundary layer. The PBLH of the daytime boundary layer over Athens and Thessaloniki is found to exhibit seasonal variability with summer maxima and winter minima and has annual median values in the range of 1.4&ndash, 1.7 km estimated using the gradient method. The PBLH over Heraklion is found to exhibit weak seasonal variability with a lower median value of 1.2 km. The nighttime boundary layer over all three sites is found to be much shallower with PBLH values in the range of 150&ndash, 200 m with no seasonal variations. In addition, the bulk Richardson number method is found to systematically underestimate the PBLH compared to the gradient method. The wind field in the daytime boundary layer at all three sites is found to have small shear of the order of 1 ms&minus, 1 and wind turning angles that are lower than 15 degrees, while in the nocturnal boundary layer it has larger shear of the order of 5&ndash, 10 ms&minus, 1 with turning angles lower than 20 degrees. In addition, for both the daytime and the nighttime boundary layer there is no general preference for veering or backing.

The accurate prediction of wildfire behavior and spread is possible only when fire and atmosphere simulations are coupled. In this work, we present a mechanism that causes a small fire to intensify by altering the atmosphere. These alterations are caused by fire-related fluxes at the surface. The fire plume and fluxes increase the convective available potential energy (CAPE) and the chance of the development of a strong pyroconvection system. To study this possible mechanism, we used WRF-Fire to capture fire line propagation as the result of interactions between heat and moisture fluxes, pressure perturbations, wind shear development and dry air downdraft. The wind patterns and dynamics of the pyroconvection system are simulated for the Horse River wildfire at Fort McMurray, Canada. The results revealed that the updraft speed reached up to 12 m/s. The entrainment mixed the mid and upper-level dry air and lowered the atmospheric moisture. The mid-level and upper-level dew point temperature changed by 5&ndash, 10 ∘ C in a short period of time. The buoyant air strengthened the ascent as soon as the nocturnal inversion was eliminated by daytime heating. The 887 J/kg total increase of CAPE in less than 5 h and the high bulk Richardson number (BRN) of 93 were indicators of the growing pyro-cumulus cell. The presented simulation has not improved the original model or supported leading-edge numerical weather prediction (NWP) achievements, except for adapting WRF-Fire for Canadian biomass fuel. However, we were able to present a great deal of improvements in wildfire nowcasting and short-term forecasting to save lives and costs associated with wildfires. The simulation is sufficiently fast and efficient to be considered for a real-time operational model. While the project was designed and succeeded as an NWP application, we are still searching for a solution for the intractable problems associated with political borders and the current liable authorities for the further development of a new generation of national atmosphere&ndash, wildfire forecasting systems.

Atmospheric stability has been studied for decades. There are several methodologies that evolved over the years. In this study, a special experimental meteorological mast that has been erected to a complex site has been used to calculate dimensionless Obukhov length ( &zeta, = z L ) , dimensionless momentum ( &phi, m ), and heat coefficients ( &phi, h ). The results are compared with the ones from average value approaches: Richardson number, flux-profile (F-P) relations, and wind shear exponent methods. The results show that the estimated &zeta, values, using the bulk Richardson number, get along well with the reference &zeta, within the neutral and stable regimes. F-P relations and wind shear exponent methods result in the best agreement for stable and neutral regimes. Nevertheless, average oriented methods are not reliable for the other regimes.

Micrometeorological observations from a tower, an eddy-covariance (EC) station and an unmanned aircraft system (UAS) at the WINSENT test-site are used to validate a computational fluid dynamics (CFD) model, driven by a mesoscale model. The observation site is characterised by a forested escarpment in a complex terrain. A two-day measurement campaign with a flow almost perpendicular to the escarpment is analysed. The first day is dominated by high wind speeds, while, on the second one, calm wind conditions are present. Despite some minor differences, the flow structure, analysed in terms of horizontal wind speeds, wind direction and inclination angles shows similarities for both days. A real-time strategy is used for the CFD validation with the UAS measurement, where the model follows spatially and temporally the aircraft. This strategy has proved to be successful. Stability indices such as the potential temperature and the bulk Richardson number are calculated to diagnose atmospheric boundary layer (ABL) characteristics up to the highest flight level. The calculated bulk Richardson values indicate a dynamically unstable region behind the escarpment and near the ground for both days. At higher altitudes, the ABL is returning to a near neutral state. The same characteristics are found in the model but only for the first day. The second day, where shear instabilities are more dominant, is not well simulated. UAS proves its great value for sensing the flow over complex terrains at high altitudes and we demonstrate the usefulness of UAS for validating and improving models.

s The release of airborne hazardous substances in the urban roughness sublayer (RSL) is a risk to human health. The atmospheric dispersion of these materials should be contained to mitigate their adverse consequences. This study investigates the effects of atmospheric stability on the air, heat, and pollutant transport within RSL by a validated 3-D RANS code. A building morphology with height variance is proposed for better representing the real urban environment and characterizing the complex interactions between the roughness elements and the flow regimes. The thermal stratification stability test subject to the ground-inflow temperature difference is extended to span from strongly unstable to moderately stable. The quantitative indicators, including air exchange rate, heat removal rate, pollutant removal rate, heat transfer coefficient, and pollutant transfer coefficient, are introduced and analyzed. The correlations between indicators and thermal stabilities are established, which provides explicit expressions to describe the influence of the changing bulk Richardson number (Rb). Results show that the design of height-variance elements enables a stronger lateral momentum towards the target street canyon, which suppresses the spanwise dispersion of pollutants. The discussions upon heat and pollutant transport based on stability categories show a high Rb-dependence of heat/mass transfer efficiency. The stability threshold is demonstrated to be Rb ≈ 0.7, where the flow speed near the ground approaches zero. As a result, the high temperature gradient is formed and acts as positive feedback to facilitate a more stratified condition. The accurate and rational correlations obtained in this study will save the computational cost considerably for further research. Also, the available results will be a reference for environmentally friendly designs.

The availability of 5-year time series of velocity and temperature data from two sonic anemometers installed at Jang Bogo Station, Antarctica, allowed a systematic investigation of the turbulence features in a stable layer affected by submeso motions and characterized by the vertical divergence of some second-order moments for a large fraction of time (quite a non-ideal surface layer). The investigation of the effect of the averaging time interval on the statistics of the second-order moments showed that this is greater for the variances of the velocity components with respect to that for the vertical fluxes. This corresponds to a greater contribution from low-frequency motions. The turbulence statistics were investigated and compared with current literature results in terms of vertical structure, share of energy between horizontal and vertical components, skewness of the vertical velocity and turbulent velocities. As a general result, all the normalized second-order moments show a clear change passing from neutral to stable conditions, passing through the range of bulk Richardson number equal to 0.1–1.

The seasonal thermocline in shallow shelf seas acts as a natural barrier for boundary-generated turbulence, damping scalar transport to the upper regions of the water column, and controlling primary production to a certain extent. To better understand turbulence and mixing conditions within the thermocline, two unique 12- and 17-day datasets with continuous measurements of the dissipation rate of turbulent kinetic energy (e) collected by autonomous underwater gliders under stratified to well-mixed conditions are presented. A highly intermittent e signal was observed in the stratified thermocline region, which was mainly characterized by quiescent flow (turbulent activity index below 7). The rate of diapycnal mixing remained relatively constant for the majority of the time with peaks of higher fluxes that were responsible for much of the increase in bottom mixed layer temperature. The water column stayed predominantly strongly stratified, with a bulk Richardson number across the thermocline well above 2. A positive relationship between the intensity of turbulence, shear and stratification was found. The trend between turbulence levels and the bulk Richardson number was relatively weak, but suggests that e increases as the bulk Richardson number approaches 1. The results also highlight the interpretation difficulties in both quantifying turbulent thermocline fluxes, as well as the responsible mechanisms.

The paper investigates buoyancy impact on the vertical flow over a backward-facing step at low Prandtl number by Direct Numerical Simulation. In particular, the very low Prandtl number of liquid sodium, 0.0088, is considered in the regime of mixed convection, i.e. for Richardson numbers below unity. The effects of buoyancy on mean flow, heat transfer and turbulence are assessed. Buoyancy is found to attenuate recirculation and, consequently, increase heat transfer. Turbulence is decreased in the attached boundary layer for moderate buoyancy impact but surpasses the levels found in forced convection at the largest Richardson number investigated. Beyond the mean flow and second moments, the budgets of turbulent kinetic energy, Reynolds shear stress, temperature variance, and turbulent heat flux components are studied and related to the alterations in the mean field quantities. Due to scale separation, production and dissipation nearly balance for temperature variance while this is not the case for turbulent kinetic energy. Similar findings for the turbulent heat fluxes show that the correlation between temperature and pressure gradient is the most important contribution to the budget aside from production and dissipation. In addition to the physical insight into this flow, the data presented may be used for the validation and improvement of turbulence models for liquid metal flows.

In climate and weather prediction models the near-surface turbulent fluxes of heat and momentum and related transfer coefficients are usually parametrized on the basis of Monin–Obukhov similarity theory (MOST). To avoid iteration, required for the numerical solution of the MOST equations, many models apply parametrizations of the transfer coefficients based on an approach relating these coefficients to the bulk Richardson number $$Ri_{b}$$ . However, the parametrizations that are presently used in most climate models are valid only for weaker stability and larger surface roughnesses than those documented during the Surface Heat Budget of the Arctic Ocean campaign (SHEBA). The latter delivered a well-accepted set of turbulence data in the stable surface layer over polar sea-ice. Using stability functions based on the SHEBA data, we solve the MOST equations applying a new semi-analytic approach that results in transfer coefficients as a function of $$Ri_{b}$$ and roughness lengths for momentum and heat. It is shown that the new coefficients reproduce the coefficients obtained by the numerical iterative method with a good accuracy in the most relevant range of stability and roughness lengths. For small $$Ri_{b}$$ , the new bulk transfer coefficients are similar to the traditional coefficients, but for large $$Ri_{b}$$ they are much smaller than currently used coefficients. Finally, a possible adjustment of the latter and the implementation of the new proposed parametrizations in models are discussed.

The Miles–Howard theorem states that a necessary condition for normal-mode instability in parallel, inviscid, steady stratified shear flows is that the minimum gradient Richardson number, $Ri_{g,min}$, is less than $1/4$ somewhere in the flow. However, the non-normality of the Navier–Stokes and buoyancy equations may allow for substantial perturbation energy growth at finite times. We calculate numerically the linear optimal perturbations which maximize the perturbation energy gain for a stably stratified shear layer consisting of a hyperbolic tangent velocity distribution with characteristic velocity $U_{0}^{\ast }$ and a uniform stratification with constant buoyancy frequency $N_{0}^{\ast }$. We vary the bulk Richardson number $Ri_{b}=N_{0}^{\ast 2}h^{\ast 2}/U_{0}^{\ast 2}$ (corresponding to $Ri_{g,min}$) between 0.20 and 0.50 and the Reynolds numbers $\mathit{Re}=U_{0}^{\ast }h^{\ast }/\unicode[STIX]{x1D708}^{\ast }$ between 1000 and 8000, with the Prandtl number held fixed at $\mathit{Pr}=1$. We find the transient growth of non-normal perturbations may be sufficient to trigger strongly nonlinear effects and breakdown into small-scale structures, thereby leading to enhanced dissipation and non-trivial modification of the background flow even in flows where $Ri_{g,min}>1/4$. We show that the effects of nonlinearity are more significant for flows with higher $\mathit{Re}$, lower $Ri_{b}$ and higher initial perturbation amplitude $E_{0}$. Enhanced kinetic energy dissipation is observed for higher-$Re$ and lower-$Ri_{b}$ flows, and the mixing efficiency, quantified here by $\unicode[STIX]{x1D700}_{p}/(\unicode[STIX]{x1D700}_{p}+\unicode[STIX]{x1D700}_{k})$ where $\unicode[STIX]{x1D700}_{p}$ is the dissipation rate of density variance and $\unicode[STIX]{x1D700}_{k}$ is the dissipation rate of kinetic energy, is found to be approximately 0.35 for the most strongly nonlinear cases.

A non-iterative analytical scheme is developed for unstable stratification that parametrizes the Monin–Obukhov stability parameter $$\zeta $$ ( $${=}z{/}L$$ , where z is the height above the ground and L is the Obukhov length) in terms of bulk Richardson number ( $$Ri_B$$ ) within the framework of Businger–Dyer type similarity functions. The proposed scheme is valid for a wide range of roughness lengths of heat and momentum. The absolute relative error in the transfer coefficients of heat and momentum is found to be less than 1.5% as compared to those obtained from an iterative scheme for Businger–Dyer type similarity functions. An attempt has been made to extend this scheme to incorporate the similarity functions having a theoretically consistent free convection limit. Further, the performance of the scheme is evaluated using observational data from two different sites. The proposed scheme is simple, non-iterative and relatively more accurate compared to the schemes reported in the literature and can be used as a potential alternative to iterative schemes used in numerical models of the atmosphere.

The mixing properties of statically stable density interfaces subject to imposed vertical shear are studied using direct numerical simulations of stratified plane Couette flow. The simulations are designed to investigate possible self-maintaining mechanisms of sharp density interfaces motivated by Phillips’ argument (Deep-Sea Res., vol. 19, 1972, pp. 79–81) by which layers and interfaces can spontaneously form due to vertical variations of diapycnal flux. At the start of each simulation, a sharp density interface with the same initial thickness is introduced at the midplane between two flat, horizontal walls counter-moving at velocities$\pm U_{w}$. Particular attention is paid to the effects of varying Prandtl number$\mathit{Pr}\equiv \unicode[STIX]{x1D708}/\unicode[STIX]{x1D705}$, where$\unicode[STIX]{x1D708}$and$\unicode[STIX]{x1D705}$are the molecular kinematic viscosity and diffusivity respectively, over two orders of magnitude from 0.7, 7 and 70. Varying$\mathit{Pr}$enables the system to access a considerable range of characteristic turbulent Péclet numbers$\mathit{Pe}_{\ast }\equiv {\mathcal{U}}_{\ast }{\mathcal{L}}_{\ast }/\unicode[STIX]{x1D705}$, where${\mathcal{U}}_{\ast }$and${\mathcal{L}}_{\ast }$are characteristic velocity and length scales, respectively, of the motion which acts to ‘scour’ the density interface. The dynamics of the interface varies with the stability of the interface which is characterised by a bulk Richardson number$\mathit{Ri}\,\equiv \,b_{0}h/U_{w}^{2}$, where$b_{0}$is half the initial buoyancy difference across the interface and$h$is the half-height of the channel. Shear-induced turbulence occurs at small$\mathit{Ri}$, whereas internal waves propagating on the interface dominate at large$\mathit{Ri}$. For a highly stable (i.e. large$\mathit{Ri}$) interface at sufficiently large$\mathit{Pe}_{\ast }$, the complex interfacial dynamics allows the interface to remain sharp. This ‘self-sharpening’ is due to the combined effects of the ‘scouring’ induced by the turbulence external to the interface and comparatively weak molecular diffusion across the core region of the interface. The effective diapycnal diffusivity and irreversible buoyancy flux are quantified in the tracer-based reference coordinate proposed by Winters & D’Asaro (J. Fluid Mech., vol. 317, 1996, pp. 179–193) and Nakamura (J. Atmos. Sci., vol. 53, 1996, pp. 1524–1537), which enables a detailed investigation of the self-sharpening process by analysing the local budget of buoyancy gradient in the reference coordinate. We further discuss the dependence of the effective diffusivity and overall mixing efficiency on the characteristic parameters of the flow, such as the buoyancy Reynolds number and the local gradient Richardson number, and highlight the possible role of the molecular properties of fluids on diapycnal mixing.

Observations of the heat flux over the open sea and in the coastal zone are analyzed to reexamine the relation of the heat flux to the air–sea temperature difference and wind speed. The study begins by examining problems with different methods for estimating the air–sea temperature difference. The difference between the air and within-water temperature is found to be most suitable for one dataset, while the difference between the air and the radiometrically measured sea surface temperature must be used for the second dataset. On average, the heat flux is linearly proportional to the product of the air–sea temperature difference and wind speed corresponding to approximately constant transfer coefficient for heat CH. Deviations of the heat flux from this simple relationship were generally weakly related or unrelated to the surface bulk Richardson number Rb, even for the coastal zone site. Similar results are also found for the moisture flux. In contrast to the general success of constant transfer coefficient for heat, CH plotted directly as a function of Rb is systematically related to the square root of Rb. The role of the wind speed as a shared variable between CH and Rb also predicts such a square root dependence, suggesting that the relationship could be largely due to self-correlation, which is also supported by a nominal study of self-correlation. However, confident isolation of the influences of stability, self-correlation, uncertainties in the air–sea temperature difference, and other physics requires more extensive data.

We present the linear stability properties and nonlinear evolution of two-dimensional plane Couette flow for a statically stable Boussinesq three-layer fluid of total depth$2h$between two horizontal plates driven at constant velocity$\pm \unicode[STIX]{x0394}U$. Initially the three layers have equal depth$2h/3$and densities$\unicode[STIX]{x1D70C}_{0}+\unicode[STIX]{x0394}\unicode[STIX]{x1D70C}$,$\unicode[STIX]{x1D70C}_{0}$and$\unicode[STIX]{x1D70C}_{0}-\unicode[STIX]{x0394}\unicode[STIX]{x1D70C}$, such that$\unicode[STIX]{x1D70C}_{0}\gg \unicode[STIX]{x0394}\unicode[STIX]{x1D70C}$. At finite Reynolds and Prandtl number, we demonstrate that this flow is susceptible to two distinct primary linear instabilities for sufficiently large bulk Richardson number$Ri_{B}=g\unicode[STIX]{x0394}\unicode[STIX]{x1D70C}h/(\unicode[STIX]{x1D70C}_{0}\unicode[STIX]{x0394}U^{2})$. For a given bulk Richardson number$Ri_{B}$, the zero phase speed ‘Taylor’ instability is always predicted to have the largest growth rate and to be an inherently two-dimensional instability. An inherently viscous instability, reminiscent of the ‘Holmboe’ instability, is also predicted to have a non-zero growth rate. For flows with Prandtl number$Pr=\unicode[STIX]{x1D708}/\unicode[STIX]{x1D705}=1$, where$\unicode[STIX]{x1D708}$is the kinematic viscosity, and$\unicode[STIX]{x1D705}$is the diffusivity of the density distribution, we find that the most unstable Taylor instability, maximized across wavenumber and$Ri_{B}$, has a (linear) growth rate which is a non-monotonic function of Reynolds number$Re=\unicode[STIX]{x0394}Uh/\unicode[STIX]{x1D708}$, with a global maximum at$Re=700$over 50 % larger than the growth rate as$Re\rightarrow \infty$. In a fully nonlinear evolution of the flows with$Re=700$and$Pr=1$, the two interfaces between the three density layers diffuse more rapidly than the underlying instabilities can grow from small amplitude. Therefore, we investigate numerically the nonlinear evolution of the flow at$Re=600$and$Pr=300$, and at$Re=5000$and$Pr=70$in two-dimensional domains with streamwise extent equal to two wavelengths of the Taylor instability with the largest growth rate. At both sets of parameter values, the primary Taylor instability undergoes a period of identifiable exponential ‘linear’ growth. However, we demonstrate that, unlike the so-called ‘Kelvin–Helmholtz’ instability that it superficially resembles, the Taylor instability’s finite-amplitude state of an elliptical vortex in the middle layer appears not to saturate into a quasiequilibrium state, but is rapidly destroyed by the background shear. The decay process reveals$Re$-dependent secondary processes. For the$Re=600$simulation, this decay allows the development to finite amplitude of the co-existing primary ‘viscous Holmboe wave instability’, which has a substantially smaller linear growth rate. For the$Re=5000$simulation, the Taylor instability decay induces a non-trivial modification of the mean velocity and density distributions, which nonlinearly develops into more classical finite-amplitude Holmboe waves. In both cases, the saturated nonlinear Holmboe waves are robust and long-lived in two-dimensional flow.