scholarly journals Sensitivity of an Idealized Tropical Cyclone to the Configuration of the Global Forecast System–Eddy Diffusivity Mass Flux Planetary Boundary Layer Scheme

Atmosphere ◽  
2021 ◽  
Vol 12 (2) ◽  
pp. 284
Author(s):  
Evan A. Kalina ◽  
Mrinal K. Biswas ◽  
Jun A. Zhang ◽  
Kathryn M. Newman

The intensity and structure of simulated tropical cyclones (TCs) are known to be sensitive to the planetary boundary layer (PBL) parameterization in numerical weather prediction models. In this paper, we use an idealized version of the Hurricane Weather Research and Forecast system (HWRF) with constant sea-surface temperature (SST) to examine how the configuration of the PBL scheme used in the operational HWRF affects TC intensity change (including rapid intensification) and structure. The configuration changes explored in this study include disabling non-local vertical mixing, changing the coefficients in the stability functions for momentum and heat, and directly modifying the Prandtl number (Pr), which controls the ratio of momentum to heat and moisture exchange in the PBL. Relative to the control simulation, disabling non-local mixing produced a ~15% larger storm that intensified more gradually, while changing the coefficient values used in the stability functions had little effect. Varying Pr within the PBL had the greatest impact, with the largest Pr (~1.6 versus ~0.8) associated with more rapid intensification (~38 versus 29 m s−1 per day) but a 5–10 m s−1 weaker intensity after the initial period of strengthening. This seemingly paradoxical result is likely due to a decrease in the radius of maximum wind (~15 versus 20 km), but smaller enthalpy fluxes, in simulated storms with larger Pr. These results underscore the importance of measuring the vertical eddy diffusivities of momentum, heat, and moisture under high-wind, open-ocean conditions to reduce uncertainty in Pr in the TC PBL.

2017 ◽  
Vol 30 (17) ◽  
pp. 6661-6682 ◽  
Author(s):  
Shira Raveh-Rubin

Dry-air intrusions (DIs) are dry, deeply descending airstreams from the upper troposphere toward the planetary boundary layer (PBL). The significance of DIs spans a variety of aspects, including the interaction with convection, extratropical cyclones and fronts, the PBL, and extreme surface weather. Here, a Lagrangian definition for DI trajectories is used and applied to ECMWF interim reanalysis (ERA-Interim) data. Based on the criterion of a minimum descent of 400 hPa during 48 h, a first global Lagrangian climatology of DI trajectories is compiled for the years 1979–2014, allowing quantitative understanding of the occurrence and variability of DIs, as well as the dynamical and thermodynamical interactions that determine their impact. DIs occur mainly in winter. While traveling equatorward from 40°–50° latitude, DIs typically reach the lower troposphere (with maximum frequencies of ~10% in winter) in the storm-track regions, as well as over the Mediterranean Sea, Arabian Sea, and eastern North Pacific, off the western coast of South America, South Africa, and Australia, and across the Antarctic coast. The DI descent is nearly adiabatic, with a mean potential temperature decrease of 3 K in two days. Relative humidity drops strongly during the first descent day and increases in the second day, because of mixing into the moist PBL. Significant destabilization of the lower levels occurs beneath DIs, accompanied by increased 10-m wind gusts, intense surface heat and moisture fluxes, and elevated PBL heights. Interestingly, only 1.2% of all DIs are found to originate from the stratosphere.


Atmosphere ◽  
2020 ◽  
Vol 11 (12) ◽  
pp. 1317
Author(s):  
Tito Maldonado ◽  
Jorge A. Amador ◽  
Erick R. Rivera ◽  
Hugo G. Hidalgo ◽  
Eric J. Alfaro

Hurricane Otto (2016) was characterised by remarkable meteorological features of relevance for the scientific community and society. Scientifically, among the most important attributes of Otto is that it underwent a rapid intensification (RI) process. For society, this cyclone severely impacted Costa Rica and Nicaragua, leaving enormous economic losses and many fatalities. In this study, a set of three numerical simulations are performed to examine the skill of model estimations in reproducing RI and trajectory of Hurricane Otto by comparing the results of a global model to a regional model using three different planetary boundary layer parameterizations (PBL). The objective is to set the basis for future studies that analyse the physical reasons why a particular simulation (associated with a certain model setup) performs better than others in terms of reproducing RI and trajectory. We use the regional model Weather Research and Forecasting—Advanced Research WRF (WRF-ARW) with boundary and initial conditions provided by the Global Forecast System (GFS) analysis (horizontal resolution of 0.5 degrees). The PBL used are the Medium Range Forecast, the Mellor-Yamada-Janjic (MYJ), and the Yonsei University (YSU) parameterizations. The regional model is run in three static domains with horizontal grid spacing of 27, 9 and 3 km, the latter covering the spacial extent of Otto during the simulation period. WRF-ARW results improve the GFS forecast, in almost every aspect evaluated in this study, particularly, the simulated trajectories in WRF-ARW show a better representation of the cyclone path and movement compared to GFS. Even though the MYJ experiment was the only one that exhibited an abrupt 24-h change in the storm’s surface wind, close to the 25-knot threshold, the YSU scheme presented the fastest intensification, closest to reality.


2021 ◽  
Author(s):  
Andrey Debolskiy ◽  
Evgeny Mortikov ◽  
Andrey Glazunov ◽  
Christof Lüpkes

<p>According to the Monin-Obukhov similarity theory (MOST), in the stratified surface layer of the atmosphere, the mean vertical velocity and scalars gradients are related to the turbulent fluxes of these quantities and to the distance z from the surface in a universal manner. The stability parameter ζ=z/L, where L is the Obukhov turbulent length scale, is the only dimensionless parameter that determines the flux-gradient relationships. This imposes a dependency of the dimensionless velocity and buoyancy gradients on ζ in form of universal nondimensional stability functions for  the surface layer. Over the decades a number of them were proposed and derived mostly from extensive field campaigns of measurements in the ABL. The stability functions differ from each other by both open coefficients and functional dependence on  ζ.  They have a limited range of applicability, which is often extended by incorporating the assumption about their asymptotic behavior.</p><p>           A generalization of MOST by considering the dependence of the dimensionless gradients on the local stability parameter z/Λ  in the framework of first order closures allows the extension of  the universal stability functions from the surface layer to most of the ABL. However, because of applicability constraints, differences in the asymptotic behavior and in other implied assumptions, it is not immediately obvious, which set of stability functions will perform best. In this study we analyze a set of stability functions which are implemented in a uniform manner into a one-dimensional first-order closure.  The latter applies a turbulent mixing length with generalized local MOST scaling which fits to a surface schemes employing corresponding functions for consistency. We use two numerical experiment setups accompanied with LES data for validation which correspond to the weakly stable GABLES1 case and to LES simulations of the very stable ABL based on measurements at the Antarctic station DOME-C (van der Linden et al. 2019). We also focus on the sensitivity of the 1D model results to coarser grids with respect to both the used  surface flux schemes and  the ABL turbulence closures since their are meant to be used in climate models because of numerical efficiency.</p><p>Authors want to aknowledge partial funding by Russian Foundation for Basic Research (RFBR project N 20-05-00776), sensitivity analysis and closure development were performed with support  of Russian Science Foundation (RSF No 20-17-00190). Steven van der Linden for providing LES data of DOME-C based experiments.</p><p>References:</p><p>van der Linden S.J. et al. Large-Eddy Simulations of the Steady Wintertime Antarctic Boundary Layer // Boundary Layer Meteorology 173.2 (2019): 165-192.</p>


2021 ◽  
Author(s):  
Vladimir Gryanik ◽  
Christof Luepkes ◽  
Andrey Grachev ◽  
Dmitry Sidorenko

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


2017 ◽  
Vol 145 (4) ◽  
pp. 1413-1426 ◽  
Author(s):  
Jun A. Zhang ◽  
Robert F. Rogers ◽  
Vijay Tallapragada

Abstract This study evaluates the impact of the modification of the vertical eddy diffusivity (Km) in the boundary layer parameterization of the Hurricane Weather Research and Forecasting (HWRF) Model on forecasts of tropical cyclone (TC) rapid intensification (RI). Composites of HWRF forecasts of Hurricanes Earl (2010) and Karl (2010) were compared for two versions of the planetary boundary layer (PBL) scheme in HWRF. The results show that using a smaller value of Km, in better agreement with observations, improves RI forecasts. The composite-mean, inner-core structures for the two sets of runs at the time of RI onset are compared with observational, theoretical, and modeling studies of RI to determine why the runs with reduced Km are more likely to undergo RI. It is found that the forecasts with reduced Km at the RI onset have a shallower boundary layer with stronger inflow, more unstable near-surface air outside the eyewall, stronger and deeper updrafts in regions farther inward from the radius of maximum wind (RMW), and stronger boundary layer convergence closer to the storm center, although the mean storm intensity (as measured by the 10-m winds) is similar for the two groups. Finally, it is found that the departure of the maximum tangential wind from the gradient wind at the eyewall, and the inward advection of angular momentum outside the eyewall, is much larger in the forecasts with reduced Km. This study emphasizes the important role of the boundary layer structure and dynamics in TC intensity change, supporting recent studies emphasizing boundary layer spinup mechanism, and recommends further improvement to the HWRF PBL physics.


2010 ◽  
Vol 10 (6) ◽  
pp. 1129-1149 ◽  
Author(s):  
M. Milelli ◽  
M. Turco ◽  
E. Oberto

Abstract. The forecast in areas of very complex topography, as for instance the Alpine region, is still a challenge even for the new generation of numerical weather prediction models which aim at reaching the km-scale. The problem is enhanced by a general lack of standard observations, which is even more evident over the southern side of the Alps. For this reason, it would be useful to increase the performance of the mathematical models by locally assimilating non-conventional data. Since in ARPA Piemonte there is the availability of a great number of non-GTS stations, it has been decided to assimilate the 2 m temperature, coming from this dataset, in the very-high resolution version of the COSMO model, which has a horizontal resolution of about 3 km, more similar to the average resolution of the thermometers. Four different weather situations have been considered, ranging from spring to winter, from cloudy to clear sky. The aim of the work is to investigate the effects of the assimilation of non-GTS data in order to create an operational very high-resolution analysis, but also to test the option of running in the future a very short-range forecast starting from these analyses (RUC or Rapid Update Cycle). The results, in terms of Root Mean Square Error, Mean Error and diurnal cycle of some surface variables such as 2 m temperature, 2 m relative humidity and 10 m wind intensity show a positive impact during the assimilation cycle which tends to dissipate a few hours after the end of it. Moreover, the 2 m temperature assimilation has a slightly positive or neutral impact on the vertical profiles of temperature, eventhough some calibration is needed for the precipitation field which is too much perturbed during the assimilation cycle, while it is unaffected in the forecast period. So the stability of the planetary boundary layer, on the one hand, has not been particularly improved by the new-data assimilation, but, on the other hand, it has not been destroyed. It has to be pointed out that a correct description of the planetary boundary layer, even only the lowest part of it, could be helpful to the forecasters and, in general, to the users, in order to deal with meteorological hazards such as snow (in particular snow/rain limit definition), or fog (description of temperature inversions).


2008 ◽  
Vol 47 (3) ◽  
pp. 835-852 ◽  
Author(s):  
Jérôme Patoux ◽  
Ralph C. Foster ◽  
Robert A. Brown

Abstract Oceanic surface pressure fields are derived from the NASA Quick Scatterometer (QuikSCAT) surface wind vector measurements using a two-layer similarity planetary boundary layer model in the midlatitudes and a mixed layer planetary boundary layer model in the tropics. These swath-based surface pressure fields are evaluated using the following three methods: 1) a comparison of bulk pressure gradients with buoy pressure measurements in the North Pacific and North Atlantic Oceans, 2) a least squares difference comparison with the European Centre for Medium-Range Weather Forecasts (ECMWF) surface pressure analyses, and 3) a parallel spectral analysis of the QuikSCAT and ECMWF surface pressure fields. The correlation coefficient squared between scatterometer-derived pressure fields and buoys is found to be R2 = 0.936. The average root-mean-square difference between the scatterometer-derived and the ECMWF pressure fields ranges from 1 to 3 hPa, depending on the latitude and season, and decreases after the assimilation of QuikSCAT winds in the ECMWF numerical weather prediction model. The spectral components of the scatterometer-derived pressure fields are larger than those of ECMWF surface analyses at all scales in the midlatitudes and only at shorter wavelengths in the tropics.


2016 ◽  
Author(s):  
Katherine McCaffrey ◽  
Laura Bianco ◽  
Paul Johnston ◽  
James M. Wilczak

Abstract. Observations of turbulence in the planetary boundary layer are critical for developing and evaluating boundary layer parameterizations in mesoscale numerical weather prediction models. These observations, however, are expensive, and rarely profile the entire boundary layer. Using optimized configurations for 449 MHz and 915 MHz wind profiling radars during the eXperimental Planetary boundary layer Instrumentation Assessment, improvements have been made to the historical methods of measuring vertical velocity variance through the time series of vertical velocity, as well as the Doppler spectral width. Using six heights of sonic anemometers mounted on a 300-m tower, correlations of up to R2 = 0.74 are seen in measurements of the large-scale variances from the radar time series, and R2 = 0.79 in measurements of small-scale variance from radar spectral widths. The total variance, measured as the sum of the small- and large-scales agrees well with sonic anemometers, with R2 = 0.79. Correlation is higher in daytime, convective boundary layers than nighttime, stable conditions when turbulence levels are smaller. With the good agreement with the in situ measurements, highly-resolved profiles up to 2 km can be accurately observed from the 449 MHz radar, and 1 km from the 915 MHz radar. This optimized configuration will provide unique observations for the verification and improvement to boundary layer parameterizations in mesoscale models.


2016 ◽  
Vol 73 (5) ◽  
pp. 2021-2038 ◽  
Author(s):  
Xuelong Chen ◽  
Bojan Škerlak ◽  
Mathias W. Rotach ◽  
Juan A. Añel ◽  
Zhonbgo Su ◽  
...  

Abstract The planetary boundary layer (PBL) over the Tibetan Plateau (with a mean elevation about 4 km above sea level) reaches an unmatched height of 9515 m above sea level. The proximity of this height to the tropopause facilitates an exchange between the stratosphere and the boundary layer. However, the underlying mechanisms responsible for this unique PBL have remained uncertain. Here, the authors explore these mechanisms and their relative importance using measurements of the PBL, the associated surface fluxes, and single-column and regional numerical simulations, as well as global reanalysis data. Results indicate that the dry conditions of both ground soil and atmosphere in late winter cannot explain the special PBL alone. Rather, the results from a single-column model demonstrate the key influence of the stability of the free atmosphere upon the growth of extremely deep PBLs over the Tibetan Plateau. Simulations with the numerical weather prediction model Consortium for Small-Scale Modelling (COSMO) exhibit good correspondence with the observed mean PBL structure and realistic turbulent kinetic energy distributions throughout the PBL. Using ERA-Interim, the authors furthermore find that weak atmospheric stability and the resultant deep PBLs are associated with higher upper-level potential vorticity (PV) values, which in turn correspond to a more southerly jet position and higher wind speeds. Upper-level PV structures and jet position thus influence the PBL development over the Tibetan Plateau.


Sign in / Sign up

Export Citation Format

Share Document