Different aspects of the stable boundary-layer structure are contrasted between the very stable and the weakly stable regimes from a new point of view. This study finds a limit wind speed, referred to as the crossover threshold, when the average vertical gradient of the turbulent kinetic energy switches sign at all observational levels. When the wind speed exceeds this transition, the entire stable boundary layer becomes vertically fully coupled. Consequently the very stable boundary layer in this study is considered as a decoupled regime, while the weakly stable state is referred to as a coupled regime. It is shown that the vertical profiles of other quantities, such as friction velocity, heat flux and thermal gradients are strikingly different between the two coupling states.Decomposition of turbulent kinetic energy and heat flux into temporal scales indicates overlapping of non-turbulent sub-mesoscale flow with turbulence in the decoupled case, while there is a clearer scale distinction between the two types of motions when coupling takes place. The turbulent kinetic energy budget is dominated by dissipation and shear production in both coupling states. However, the relative importance of the buoyant destruction term is shown to be appreciably larger in the decoupled regime. In the heat flux budget equation, buoyant destruction is larger in magnitude than production by the thermal gradient in the decoupled case, but not when there is full coupling. These results indicate that the surface heat flux plays a major role in controlling the stable boundary-layer state, as previously proposed. For the entire dataset, the frequency distributions of turbulence quantities near the surface are shown to be bimodal. The two modes are associated with the two coupling states, each well described by independent log-normal distributions.

Previous observational studies in the stable boundary layer diverge appreciably on the values of dimensionless ratios between turbulence-related quantities and on their stability dependence. In the present study, the hypothesis that such variability is caused by the influence of locally dependent nonturbulent processes, referred to as submeso, is tested and confirmed. This is done using six datasets collected at sites with different surface coverage. The time-scale dependence of wind components and temperature fluctuations is presented using the multiresolution decomposition, which allows the identification of the turbulence and submeso contributions to spectra and cospectra. In the submeso range, the spectra of turbulence kinetic energy range increases exponentially with time scale. The exponent decreases with the magnitude of the turbulent fluctuations at a similar manner at all sites. This fact is used to determine the smaller time scale with relevant influence of submeso processes and a ratio that quantifies the relative importance of such nonturbulent processes with respect to turbulence. Based on that, values for the local stability parameter that are unaffected by nonturbulent processes are found. It is shown that the dimensionless ratios do not usually converge to a given value as the time scale increases and that it is as a consequence of the locally dependent submeso influence. The ratios and their stability dependence are determined at the time scales with least influence of nonturbulent processes, but significant site-to-site variability persists. Combining all datasets, expressions for the dependence of the dimensionless ratios on the local stability parameter that minimize the role of the submeso contribution are proposed.

Observations of the vertical and temporal structure of the nocturnal boundary layer before and after a transition from the weakly to the very stable regime have been presented in Part I. Here, similar transitions are investigated using a one‐dimensional second‐order closure numerical model, with an energy budget solved at the surface. The transition is driven by a decreasing mean wind at the top of the domain, and simulations with different cloud covers and surface thermal properties are considered. The time of the transition depends on the wind speed at the top of the domain and on the “coupling strength” between the surface and the atmosphere, which is affected by the cloud cover and surface thermal properties. The vertical profiles and temporal evolutions of the terms of the budgets of turbulent kinetic energy (TKE), heat flux and temperature variance are presented. Of these, only TKE budget presents the same dominant terms in both regimes. Absolute heat flux in the model is proportional to the cube of the wind speed in the very stable regime.

The nocturnal boundary‐layer regime transition from weakly stable (strong wind) to very stable (weak wind) is analyzed using 10 levels of turbulence observations made at a 140 m micrometeorological mast near the southeastern Brazilian coast. The combination of synoptic and local flow favors the systematic occurrence of such a transition, typically 5 to 7 h after sunset. The regime transition is marked by decreases in temperature, wind speed, turbulent kinetic energy (TKE) and absolute heat flux. The decrease in temperature is often abrupt and the inflection point in the temperature series marks the regime transition. Absolute heat flux peaks before the transition during the weakly stable period, while temperature variance peaks near the transition. Composites from 36 cases when the cooling rate exceeded 2 °C/h are used to describe the vertical structure of the stable boundary layer (SBL) in both regimes. For these abrupt transitions, dimensionless variables that relate thermal and mechanical properties of the flow are compared as indicators of the SBL regime, and the gradient Richardson number is found to be better for that purpose. The absolute heat flux is shown to be proportional to the cube of the wind speed only in the strong wind limit of the very stable regime. Simulations of similar transitions using a second‐order model are described in Part II.

A 140-m micrometeorological tower provides detailed observations of the vertical structure of the mean and turbulent fields of meteorological variables of a coastal region in southeastern Brazil, and reveals the extent to which a nearby power plant affects the local atmospheric boundary layer.

Even though mesoscale convective systems (MCS) are of great importance in precipitation regimes besides being related to severe weather events, they are still not easy to predict. This study builds an objective index for South American MCS based on synoptic features present before the initiation of the MCS between 2005 and 2010. The National Centers for Environmental Prediction (NCEP) Climate Forecast System Reanalysis (CFSR) product is used to access environmental conditions near the MCS-initiation centroid and in unorganized convective environments to obtain a South American MCS index (SA-MCS index). Upper and middle atmospheric levels presented the main patterns used to identify the MCS environments and unorganized convection, with the jet streak downstream of a short trough at medium levels frequently present in the convectively organized environment. The SA-MCS index is shown to depend on a vertical wind shear between 0 and 6 km high, temperature advection at 775 hPa, lifted index (LI) and vertical velocity omega at 800 hPa. When compared with the MCS index developed through the United States MCS climatology, the SA-MCS index was able to predict more than twice as many MCS at a distance of less than 1 ∘ from the point of maximum intensity of the index.

Two contrasting flow regimes exist in the stable boundary layer (SBL), as evidenced from both observational and modeling studies. In general, numerical schemes such as those used in numerical weather prediction and climate models (NWPCs) reproduce a transition between SBL regimes. However, the characteristics of such a transition depend on the turbulence parameterizations and stability functions used to represent the eddy diffusivity in the models. The main goal of the present study is to detail how the two SBL regimes occur in single-column models (SCMs) by analyzing the SBL structure and its dependence on external parameters. Two different turbulence closure orders (first order and an E–l model) and two types of stability functions (short and long tail) are considered. The control exerted by the geostrophic wind and the surface cooling rate on the model SBL regimes is addressed. The model flow presents a three-layer structure: a fully turbulent, weakly stable layer (WSL) next to the surface; a very stable layer (VSL) above that; and a laminar layer above the other two and toward the domain top. It is shown that the WSL and VSL are related to both SBL regimes, respectively. Furthermore, the numerically simulated SBL presents the two-layer structure regardless of the turbulence parameterization order and stability function used. The models also reproduce other features reported in recent observational studies: an S-shaped dependence of the thermal gradient on the mean wind speed and an independence of the vertical gradient of friction velocity δu* on the mean wind speed.

Occasionally, storms that share many features with tropical cyclones, including the presence of a quasi-circular “eye” a warm core and strong winds, are observed in the Mediterranean. Generally, they are known as Medicanes, or tropical-like cyclones (TLC). Due to the intense wind forcings and the consequent development of high wind waves, a large number of sea spray droplets—both from bubble bursting and spume tearing processes—are likely to be produced at the sea surface. In order to take into account this process, we implemented an additional Sea Spray Source Function (SSSF) in WRF-Chem, model version 3.6.1, using the GOCART (Goddard Chemistry Aerosol Radiation and Transport) aerosol sectional module. Traditionally, air-sea momentum fluxes are computed through the classical Charnock relation that does not consider the wave-state and sea spray effects on the sea surface roughness explicitly. In order to take into account these forcing, we implemented a more recent parameterization of the sea surface aerodynamic roughness within the WRF surface layer model, which may be applicable to both moderate and high wind conditions. The implemented SSSF and sea surface roughness parameterization have been tested using an operative model sequence based on COAWST (Coupled Ocean Atmosphere Wave Sediment Transport) and WRF-Chem. The third-generation wave model SWAN (Simulating Waves Nearshore), two-way coupled with the WRF atmospheric model in the COAWST framework, provided wave field parameters. Numerical simulations have been integrated with the WRF-Chem chemistry package, with the aim of calculating the sea spray generated by the waves and to include its effect in the Charnock roughness parametrization together with the sea state effect. A single case study is performed, considering the Medicane that affected south-eastern Italy on 26 September 2006. Since this Medicane is one of the most deeply analysed in literature, its investigation can easily shed some light on the feedbacks between sea spray and drag coefficients.