New modified and extended stability functions for the stable boundary layer based on SHEBA data

2020 ◽  
Author(s):  
Vladimir M. Gryanik ◽  
Andrey Grachev ◽  
Christof Lüpkes ◽  
Dmitry Sidorenko
Keyword(s):  
Boundary Layer ◽  
Sea Ice ◽  
Stable Boundary Layer ◽  
Richardson Number ◽  
Prediction Models ◽  
Weather Prediction ◽  
Bulk Richardson Number ◽  
Transfer Coefficients ◽  
Near Surface ◽  
Stability Functions

<p>The calculation of the near-surface turbulent fluxes of energy and momentum in climate and weather prediction models requires transfer coefficients. Currently used parametrizations of these coefficients are based on stability functions derived from measurements over land and not over sea ice. However, recently, a non-iterative parametrization has been proposed by Gryanik and Lüpkes (2018), which can be applied to climate and weather prediction models as well but uses stability functions of Grachev et al. (2007). These functions had been obtained from measurements during the Surface Heat Budget over the Arctic Ocean campaign (SHEBA) and thus from measurements over sea ice. A drawback of the scheme of Gryanik and Lüpkes (2018) is that there is still some complexity due to the complexity of the SHEBA based functions.</p><p>Thus new stability functions are proposed for the stable boundary layer, which are also based on the SHEBA measurements but avoid the complexity. It is shown that the new functions are superior to the former ones with respect to the representation of the measured relationship between the Obukhov length and the bulk Richardson number. Moreover, the resulting transfer coefficients agree slightly better with the SHEBA observations in the very stable range. Nevertheless, the functions fulfill the same criteria of applicability as the earlier functions and contain furthermore as an extension a dependence on the neutral Prandtl number. Applying the new functions, an efficient non-iterative parametrization of the near-surface turbulent fluxes of momentum and heat is developed where transfer coefficients result as a function of the bulk Richardson number (Ri<sub>b</sub>) and roughness parameters. The new transfer coefficients, which are recommended for weather and climate models, agree well with the SHEBA data in a large range of stability (0< Ri<sub>b</sub><0.5) and with those based on the Dyer-Businger functions in the range Ri<sub>b</sub> <0.08.</p><p><strong>References</strong></p><p><span>Grachev A.A., Andreas E.L, Fairall C.W., Guest P.S., Persson POG (2007) Boundary-Layer Meteorol., </span><span><strong>124</strong></span><span>, 315–333</span><span>.</span></p><p><span>Gryanik, V.M. and Lüpkes C. (2018) An Efficient Non-iterative Bulk Parametrization of Surface Fluxes for Stable Atmospheric Conditions Over Polar Sea-Ice,Boundary-Layer Meteorol 166:301-325</span></p>

scholarly journals New Modified and Extended Stability Functions for the Stable Boundary Layer based on SHEBA and Parametrizations of Bulk Transfer Coefficients for Climate Models

2020 ◽  
Vol 77 (8) ◽  
pp. 2687-2716
Author(s):  
Vladimir M. Gryanik ◽  
Christof Lüpkes ◽  
Andrey Grachev ◽  
Dmitry Sidorenko
Keyword(s):  
Sea Ice ◽  
Richardson Number ◽  
Climate Models ◽  
Turbulent Fluxes ◽  
The Arctic ◽  
Bulk Richardson Number ◽  
Transfer Coefficients ◽  
Near Surface ◽  
Stability Functions ◽  
Stable Surface Layer

Abstract 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.

A Package of New Universal Non-Iterative Parametrizations for Stable Surface Layer Transfer Coefficients of Momentum, Heat and Moisture in Numerical Earth System Models

2021 ◽  
Author(s):  
Vladimir Gryanik ◽  
Christof Luepkes ◽  
Andrey Grachev ◽  
Dmitry Sidorenko
Keyword(s):  
Climate Models ◽  
Similarity Theory ◽  
Weather Prediction ◽  
Climate Projections ◽  
Transfer Coefficients ◽  
Near Surface ◽  
Stability Functions ◽  
Stable Surface Layer ◽  
Heat And Moisture ◽  
The Impact

<p><span>Results of weather forecast, present-day climate simulations and future climate projections depend among other factors on the interaction between the atmosphere and the underlying sea-ice, the land and the ocean. In numerical weather prediction and climate models some of these interactions are accounted for by transport coefficients describing turbulent exchange of momentum, heat and moisture. Currently used transfer coefficients have, however, large uncertainties in flow regimes being typical for cold nights and seasons, but especially in the polar regions. Furthermore, their determination is numerically complex. It is obvious that progress could be achieved when the transfer coefficients would be given by simple mathematical formulae in frames of an economic computational scheme. Such a new universal, so-called non-iterative parametrization scheme is derived for a package of transfer coefficients.</span></p><p><span>The derivation is based on the Monin-Obukhov similarity theory, which is over the years well accepted in the scientific community. The newly derived non-iterative scheme provides a basis for a cheap systematic study of the impact of near-surface turbulence and of the related transports of momentum, heat and moisture in NWP and climate models. </span></p><p><span>We show that often used transfer coefficients like those of Louis et al. (1982) or of Cheng and Brutsaert (2005) can be applied at large stability only with some caution, keeping in mind that at large stability they significantly overestimate the transfer coefficient compared with most comprehensive measurements. The latter are best reproduced by Gryanik et al. (2020) functions, which are part of the package. We show that the new scheme is flexible, thus, new stability functions can be added to the package, if required. </span></p><p> </p><p> <span>Gryanik, V.M., Lüpkes, C., Grachev, A., Sidorenko, D. (2020) New Modified and Extended Stability Functions for the Stable Boundary Layer based on SHEBA and Parametrizations of Bulk Transfer Coefficients for Climate Models, J. Atmos. Sci., 77, 2687-2716</span></p><p><br><br></p>

scholarly journals Impact of planetary boundary layer turbulence on model climate and tracer transport

2014 ◽  
Vol 14 (23) ◽  
pp. 31627-31674
Author(s):  
E. L. McGrath-Spangler ◽  
A. Molod ◽  
L. E. Ott ◽  
S. Pawson
Keyword(s):  
Boundary Layer ◽  
Planetary Boundary Layer ◽  
Richardson Number ◽  
Surface Wind ◽  
Specific Humidity ◽  
Boreal Summer ◽  
Tracer Transport ◽  
Bulk Richardson Number ◽  
Turbulent Length Scale ◽  
Near Surface

Abstract. Planetary boundary layer (PBL) processes are important for weather, climate, and tracer transport and concentration. One measure of the strength of these processes is the PBL depth. However, no single PBL depth definition exists and several studies have found that the estimated depth can vary substantially based on the definition used. In the Goddard Earth Observing System (GEOS-5) atmospheric general circulation model, the PBL depth is particularly important because it is used to calculate the turbulent length scale that is used in the estimation of turbulent mixing. This study analyzes the impact of using three different PBL depth definitions in this calculation. Two definitions are based on the scalar eddy diffusion coefficient and the third is based on the bulk Richardson number. Over land, the bulk Richardson number definition estimates shallower nocturnal PBLs than the other estimates while over water this definition generally produces deeper PBLs. The near surface wind velocity, temperature, and specific humidity responses to the change in turbulence are spatially and temporally heterogeneous, resulting in changes to tracer transport and concentrations. Near surface wind speed increases in the bulk Richardson number experiment cause Saharan dust increases on the order of 1 × 10−4 kg m−2 downwind over the Atlantic Ocean. Carbon monoxide (CO) surface concentrations are modified over Africa during boreal summer, producing differences on the order of 20 ppb, due to the model's treatment of emissions from biomass burning. While differences in carbon dioxide (CO2) are small in the time mean, instantaneous differences are on the order of 10 ppm and these are especially prevalent at high latitude during boreal winter. Understanding the sensitivity of trace gas and aerosol concentration estimates to PBL depth is important for studies seeking to calculate surface fluxes based on near-surface concentrations and to studies projecting future concentrations.

scholarly journals Stable Boundary Layer in Complex Terrain. Part I: Linking Fluxes and Intermittency to an Average Stability Index

2014 ◽  
Vol 53 (9) ◽  
pp. 2196-2215 ◽  
Author(s):  
Luiz E. Medeiros ◽  
David R. Fitzjarrald
Keyword(s):  
Richardson Number ◽  
Prediction Models ◽  
Weather Prediction ◽  
Stability Index ◽  
Bulk Richardson Number ◽  
Hudson Valley ◽  
Continuous Mixing ◽  
Heterogeneous Landscape ◽  
Turbulence Intermittency ◽  
Nearly Continuous

AbstractAverage heat and momentum fluxes observed by a network of surface stations during the Hudson Valley Ambient Meteorology Study (HVAMS) were found as functions of a spatially representative bulk Richardson number Ribr. Preferential sites were identified for the occurrence of strong turbulence under mesoscale stability conditions common to all stations. Locally sensed turbulence intermittency depends on the mesoscale flow stability. Nearly continuous turbulence with few long-lived intermittent events occurs when Ribr < Ricr, the critical gradient Richardson number. Less-continuous mixing associated with a larger number of events occurs when Ricr < Ribr < 5, with the weakest turbulence and fewer events observed for Ribr ≫ Ricr. It was found that the need to allow for extra mixing above the conventional critical bulk Richardson number in numerical weather prediction models is primarily a consequence of spatial averaging in a heterogeneous landscape and is secondarily the result of turbulence above Ricr at locations with “nonideal fetch.”

scholarly journals Near-Surface Temperature Inversion Growth Rate during the Onset of the Stable Boundary Layer

2017 ◽  
Vol 74 (10) ◽  
pp. 3433-3449 ◽  
Author(s):  
Ivo G. S. van Hooijdonk ◽  
Herman J. H. Clercx ◽  
Carsten Abraham ◽  
Amber M. Holdsworth ◽  
Adam H. Monahan ◽  
...  
Keyword(s):  
Boundary Layer ◽  
Growth Rate ◽  
Decay Rate ◽  
Stable Boundary Layer ◽  
Numerical Solutions ◽  
Weather Prediction ◽  
Temperature Inversion ◽  
Near Surface ◽  
Stable Boundary ◽  
History Effects

Abstract This study aims to find the typical growth rate of the temperature inversion during the onset of the stable boundary layer around sunset. The sunset transition is a very challenging period for numerical weather prediction, since neither accepted theories for the convective boundary layer nor those for the stable boundary layer appear to be applicable. To gain more insight in this period, a systematic investigation of the temperature inversion growth rate is conducted. A statistical procedure is used to analyze almost 16 years of observations from the Cabauw observational tower, supported by observations from two additional sites (Dome C and Karlsruhe). The results show that, on average, the growth rate of the temperature inversion (normalized by the maximum inversion during the night) weakly declines with increasing wind speed. The observed growth rate is quantitatively consistent among the sites, and it appears insensitive to various other parameters. The results were also insensitive to the afternoon decay rate of the net radiation except when this decay rate was very weak. These observations are compared to numerical solutions of three models with increasing complexity: a bulk model, an idealized single-column model (SCM), and an operational-level SCM. It appears only the latter could reproduce qualitative features of the observations using a first-order closure. Moreover, replacing this closure with a prognostic TKE scheme substantially improved the quantitative performance. This suggests that idealized models assuming instantaneous equilibrium flux-profile relations may not aid in understanding this period, since history effects may qualitatively affect the dynamics.

Observational Support for the Stability Dependence of the Bulk Richardson Number Across the Stable Boundary Layer

2013 ◽  
Vol 150 (3) ◽  
pp. 515-523 ◽  
Author(s):  
S. Basu ◽  
A. A. M. Holtslag ◽  
L. Caporaso ◽  
A. Riccio ◽  
G.-J. Steeneveld
Keyword(s):  
Boundary Layer ◽  
Stable Boundary Layer ◽  
Richardson Number ◽  
Bulk Richardson Number ◽  
Stable Boundary ◽  
The Stability

Boundary Layer and Surface Verification of the High-Resolution Rapid Refresh, Version 3

2020 ◽  
Vol 35 (6) ◽  
pp. 2255-2278
Author(s):  
Robert G. Fovell ◽  
Alex Gallagher
Keyword(s):  
Boundary Layer ◽  
High Resolution ◽  
Wind Speed ◽  
Prediction Models ◽  
Weather Prediction ◽  
Potential Temperature ◽  
Ground Level ◽  
Systematic Errors ◽  
Radiosonde Data ◽  
Near Surface

AbstractWhile numerical weather prediction models have made considerable progress regarding forecast skill, less attention has been paid to the planetary boundary layer. This study leverages High-Resolution Rapid Refresh (HRRR) forecasts on native levels, 1-s radiosonde data, and (primarily airport) surface observations across the conterminous United States. We construct temporally and spatially averaged composites of wind speed and potential temperature in the lowest 1 km for selected months to identify systematic errors in both forecasts and observations in this critical layer. We find near-surface temperature and wind speed predictions to be skillful, although wind biases were negatively correlated with observed speed and temperature biases revealed a robust relationship with station elevation. Above ≈250 m above ground level, below which radiosonde wind data were apparently contaminated by processing, biases were small for wind speed and potential temperature at the analysis time (which incorporates sonde data) but became substantial by the 24-h forecast. Wind biases were positive through the layer for both 0000 and 1200 UTC, and morning potential temperature profiles were marked by excessively steep lapse rates that persisted across seasons and (again) exaggerated at higher elevation sites. While the source or cause of these systematic errors are not fully understood, this analysis highlights areas for potential model improvement and the need for a continued and accessible archive of the data that make analyses like this possible.

scholarly journals Impact of planetary boundary layer turbulence on model climate and tracer transport

2015 ◽  
Vol 15 (13) ◽  
pp. 7269-7286 ◽  
Author(s):  
E. L. McGrath-Spangler ◽  
A. Molod ◽  
L. E. Ott ◽  
S. Pawson
Keyword(s):  
Boundary Layer ◽  
Planetary Boundary Layer ◽  
Richardson Number ◽  
Surface Wind ◽  
Specific Humidity ◽  
Boreal Summer ◽  
Tracer Transport ◽  
Bulk Richardson Number ◽  
Turbulent Length Scale ◽  
Near Surface

Abstract. Planetary boundary layer (PBL) processes are important for weather, climate, and tracer transport and concentration. One measure of the strength of these processes is the PBL depth. However, no single PBL depth definition exists and several studies have found that the estimated depth can vary substantially based on the definition used. In the Goddard Earth Observing System (GEOS-5) atmospheric general circulation model, the PBL depth is particularly important because it is used to calculate the turbulent length scale that is used in the estimation of turbulent mixing. This study analyzes the impact of using three different PBL depth definitions in this calculation. Two definitions are based on the scalar eddy diffusion coefficient and the third is based on the bulk Richardson number. Over land, the bulk Richardson number definition estimates shallower nocturnal PBLs than the other estimates while over water this definition generally produces deeper PBLs. The near-surface wind velocity, temperature, and specific humidity responses to the change in turbulence are spatially and temporally heterogeneous, resulting in changes to tracer transport and concentrations. Near-surface wind speed increases in the bulk Richardson number experiment cause Saharan dust increases on the order of 1 × 10−4 kg m−2 downwind over the Atlantic Ocean. Carbon monoxide (CO) surface concentrations are modified over Africa during boreal summer, producing differences on the order of 20 ppb, due to the model's treatment of emissions from biomass burning. While differences in carbon dioxide (CO2) are small in the time mean, instantaneous differences are on the order of 10 ppm and these are especially prevalent at high latitude during boreal winter. Understanding the sensitivity of trace gas and aerosol concentration estimates to PBL depth is important for studies seeking to calculate surface fluxes based on near-surface concentrations and for studies projecting future concentrations.

The Very Stable Boundary Layer on Nights with Weak Low-Level Jets

2007 ◽  
Vol 64 (9) ◽  
pp. 3068-3090 ◽  
Author(s):  
Robert M. Banta ◽  
Larry Mahrt ◽  
Dean Vickers ◽  
Jielun Sun ◽  
Ben B. Balsley ◽  
...  
Keyword(s):  
Boundary Layer ◽  
Time Series ◽  
Stable Boundary Layer ◽  
Weather Prediction ◽  
Density Currents ◽  
Lower Boundary Condition ◽  
Near Surface ◽  
Stable Boundary ◽  
Low Level ◽  
Low Level Jets

Abstract The light-wind, clear-sky, very stable boundary layer (vSBL) is characterized by large values of bulk Richardson number. The light winds produce weak shear, turbulence, and mixing, and resulting strong temperature gradients near the surface. Here five nights with weak-wind, very stable boundary layers during the Cooperative Atmosphere–Surface Exchange Study (CASES-99) are investigated. Although the winds were light and variable near the surface, Doppler lidar profiles of wind speed often indicated persistent profile shapes and magnitudes for periods of an hour or more, sometimes exhibiting jetlike maxima. The near-surface structure of the boundary layer (BL) on the five nights all showed characteristics typical of the vSBL. These characteristics included a shallow traditional BL only 10–30 m deep with weak intermittent turbulence within the strong surface-based radiation inversion. Above this shallow BL sat a layer of very weak turbulence and negligible turbulent mixing. The focus of this paper is on the effects of this quiescent layer just above the shallow BL, and the impacts of this quiescent layer on turbulent transport and numerical modeling. High-frequency time series of temperature T on a 60-m tower showed that 1) the amplitudes of the T fluctuations were dramatically suppressed at levels above 30 m in contrast to the relatively larger intermittent T fluctuations in the shallow BL below, and 2) the temperature at 40- to 60-m height was nearly constant for several hours, indicating that the very cold air near the surface was not being mixed upward to those levels. The presence of this quiescent layer indicates that the atmosphere above the shallow BL was isolated and detached both from the surface and from the shallow BL. Although some of the nights studied had modestly stronger winds and traveling disturbances (density currents, gravity waves, shear instabilities), these disturbances seemed to pass through the region without having much effect on either the SBL structure or on the atmosphere–surface decoupling. The decoupling suggests that under very stable conditions, the surface-layer lower boundary condition for numerical weather prediction models should act to decouple and isolate the surface from the atmosphere, for example, as a free-slip, thermally insulated layer. A multiday time series of ozone from an air quality campaign in Tennessee, which exhibited nocturnal behavior typical of polluted air, showed the disappearance of ozone on weak low-level jets (LLJ) nights. This behavior is consistent with the two-stratum structure of the vSBL, and with the nearly complete isolation of the surface and the shallow BL from the rest of the atmosphere above, in contrast to cases with stronger LLJs, where such coupling was stronger.

scholarly journals On the computation of planetary boundary layer height using the bulk Richardson number method

2014 ◽  
Vol 7 (3) ◽  
pp. 4045-4079 ◽  
Author(s):  
Y. Zhang ◽  
Z. Gao ◽  
D. Li ◽  
Y. Li ◽  
N. Zhang ◽  
...  
Keyword(s):  
Boundary Layer ◽  
Stable Boundary Layer ◽  
Richardson Number ◽  
Bulk Richardson Number ◽  
Stable Boundary ◽  
Boundary Layer Height ◽  
Unstable Boundary Layer ◽  
Richardson Method ◽  
Field Campaigns ◽  
Weakly Stable

Abstract. Experimental data from four intensive field campaigns are used to explore the variability of the critical bulk Richardson number, which is a key parameter for calculating the planetary boundary layer height (PBLH) in numerical weather and climate models with the bulk Richardson method. First, the PBLHs of three different thermally-stratified boundary layers (i.e., strongly stable boundary layer, weakly stable boundary layer, and unstable boundary layer) from the four field campaigns are determined using the turbulence method, the potential temperature gradient method, the low-level jet method, or the modified parcel method. Then for each type of boundary layer, an optimal critical Richardson numbers is obtained through linear fitting and statistical error minimization methods so that the bulk Richardson method with this optimal critical bulk Richardson number yields similar estimates of PBLHs as the methods mentioned above. We find that the optimal critical bulk Richardson number increases as the atmosphere becomes more unstable: 0.24 for strongly stable boundary layer, 0.31 for weakly stable boundary layer, and 0.39 for unstable boundary layer. Compared with previous schemes that use a single value of critical bulk Richardson number for calculating the PBLH, the new values of critical bulk Richardson number that proposed by this study yield more accurate estimate of PBLH.

Sign in / Sign up

Export Citation Format

Share Document