Search Results
You are looking at 1 - 10 of 28 items for
- Author or Editor: Pedro Jimenez x
- Refine by Access: All Content x
Abstract
The ability of the Weather Research and Forecasting (WRF) model to reproduce the surface wind direction over complex terrain is examined. A simulation spanning a winter season at a high horizontal resolution of 2 km is compared with wind direction records from a surface observational network located in the northeastern Iberian Peninsula. A previous evaluation has shown the ability of WRF to reproduce the wind speed over the region once the effects of the subgrid-scale topography are parameterized. Hence, the current investigation complements the previous findings, providing information about the model's ability to reproduce the direction of the surface flow. The differences between the observations and the model are quantified in terms of scores explicitly designed to handle the circular nature of the wind direction. Results show that the differences depend on the wind speed. The larger the wind speed is, the smaller are the wind direction differences. Areas with more complex terrain show larger systematic differences between model and observations; in these areas, a statistical correction is shown to help. The importance of the grid point selected for the comparison with observations is also analyzed. A careful selection is relevant to reducing comparative problems over complex terrain.
Abstract
The ability of the Weather Research and Forecasting (WRF) model to reproduce the surface wind direction over complex terrain is examined. A simulation spanning a winter season at a high horizontal resolution of 2 km is compared with wind direction records from a surface observational network located in the northeastern Iberian Peninsula. A previous evaluation has shown the ability of WRF to reproduce the wind speed over the region once the effects of the subgrid-scale topography are parameterized. Hence, the current investigation complements the previous findings, providing information about the model's ability to reproduce the direction of the surface flow. The differences between the observations and the model are quantified in terms of scores explicitly designed to handle the circular nature of the wind direction. Results show that the differences depend on the wind speed. The larger the wind speed is, the smaller are the wind direction differences. Areas with more complex terrain show larger systematic differences between model and observations; in these areas, a statistical correction is shown to help. The importance of the grid point selected for the comparison with observations is also analyzed. A careful selection is relevant to reducing comparative problems over complex terrain.
Abstract
The wind stress formulation in an atmospheric model over shallow waters is investigated using year-long observations of the wind profile within the first 100 m of the atmosphere and mesoscale simulations. The model experiments use a range of planetary boundary layer parameterizations to quantify the uncertainty related to the turbulent closure assumptions and thus to isolate the dominant influence of the surface roughness formulation. Results indicate that a positive wind speed bias exists when common open-ocean formulations for roughness are adopted for a region with a water depth of 30 m. Imposition of a wind stress formulation that is consistent with previous shallow-water estimates is necessary to reconcile model wind speeds with observations, providing modeling evidence that supports the increase of surface drag over shallow waters. The possibility of including water depth in the parameterization of roughness length is examined.
Abstract
The wind stress formulation in an atmospheric model over shallow waters is investigated using year-long observations of the wind profile within the first 100 m of the atmosphere and mesoscale simulations. The model experiments use a range of planetary boundary layer parameterizations to quantify the uncertainty related to the turbulent closure assumptions and thus to isolate the dominant influence of the surface roughness formulation. Results indicate that a positive wind speed bias exists when common open-ocean formulations for roughness are adopted for a region with a water depth of 30 m. Imposition of a wind stress formulation that is consistent with previous shallow-water estimates is necessary to reconcile model wind speeds with observations, providing modeling evidence that supports the increase of surface drag over shallow waters. The possibility of including water depth in the parameterization of roughness length is examined.
Abstract
The Weather Research and Forecasting (WRF) model presents a high surface wind speed bias over plains and valleys that constitutes a limitation for the increasing use of the model for several applications. This study attempts to correct for this bias by parameterizing the effects that the unresolved topographic features exert over the momentum flux. The proposed parameterization is based on the concept of a momentum sink term and makes use of the standard deviation of the subgrid-scale orography as well as the Laplacian of the topographic field. Both the drag generated by the unresolved terrain and the possibility of an increase in the speed of the flow over the mountains and hills, where it is herein shown that WRF presents a low wind speed bias, are considered in the scheme. The surface wind simulation over a complex-terrain region that is located in the northeast of the Iberian Peninsula is improved with the inclusion of the new parameterization. In particular, the underestimation of the wind speed spatial variability resulting from the mentioned biases is corrected. The importance of selecting appropriate grid points to compare with observations is also examined. The wind speed from the nearest grid point is not always the most appropriate one for this comparison, nearby ones being more representative. The new scheme not only improves the climatological winds but also the intradiurnal variations at the mountains, over which the default WRF shows limitations in reproducing the observed wind behavior. Some advantages of the proposed formulation for wind-resource evaluation are also discussed.
Abstract
The Weather Research and Forecasting (WRF) model presents a high surface wind speed bias over plains and valleys that constitutes a limitation for the increasing use of the model for several applications. This study attempts to correct for this bias by parameterizing the effects that the unresolved topographic features exert over the momentum flux. The proposed parameterization is based on the concept of a momentum sink term and makes use of the standard deviation of the subgrid-scale orography as well as the Laplacian of the topographic field. Both the drag generated by the unresolved terrain and the possibility of an increase in the speed of the flow over the mountains and hills, where it is herein shown that WRF presents a low wind speed bias, are considered in the scheme. The surface wind simulation over a complex-terrain region that is located in the northeast of the Iberian Peninsula is improved with the inclusion of the new parameterization. In particular, the underestimation of the wind speed spatial variability resulting from the mentioned biases is corrected. The importance of selecting appropriate grid points to compare with observations is also examined. The wind speed from the nearest grid point is not always the most appropriate one for this comparison, nearby ones being more representative. The new scheme not only improves the climatological winds but also the intradiurnal variations at the mountains, over which the default WRF shows limitations in reproducing the observed wind behavior. Some advantages of the proposed formulation for wind-resource evaluation are also discussed.
Abstract
As wind farms grow in number and size worldwide, it is important that their potential impacts on the environment are studied and understood. The Fitch parameterization implemented in the Weather Research and Forecasting (WRF) Model since version 3.3 is a widely used tool today to study such impacts. We identified two important issues related to the way the added turbulent kinetic energy (TKE) generated by a wind farm is treated in the WRF Model with the Fitch parameterization. The first issue is a simple “bug” in the WRF code, and the second issue is the excessive value of a coefficient, called C TKE, that relates TKE to the turbine electromechanical losses. These two issues directly affect the way that a wind farm wake evolves, and they impact properties like near-surface temperature and wind speed at the wind farm as well as behind it in the wake. We provide a bug fix and a revised value of C TKE that is one-quarter of the original value. This 0.25 correction factor is empirical; future studies should examine its dependence on parameters such as atmospheric stability, grid resolution, and wind farm layout. We present the results obtained with the Fitch parameterization in the WRF Model for a single turbine with and without the bug fix and the corrected C TKE and compare them with high-fidelity large-eddy simulations. These two issues have not been discovered before because they interact with one another in such a way that their combined effect is a somewhat realistic vertical TKE profile at the wind farm and a realistic wind speed deficit in the wake. All WRF simulations that used the Fitch wind farm parameterization are affected, and their conclusions may need to be revisited.
Abstract
As wind farms grow in number and size worldwide, it is important that their potential impacts on the environment are studied and understood. The Fitch parameterization implemented in the Weather Research and Forecasting (WRF) Model since version 3.3 is a widely used tool today to study such impacts. We identified two important issues related to the way the added turbulent kinetic energy (TKE) generated by a wind farm is treated in the WRF Model with the Fitch parameterization. The first issue is a simple “bug” in the WRF code, and the second issue is the excessive value of a coefficient, called C TKE, that relates TKE to the turbine electromechanical losses. These two issues directly affect the way that a wind farm wake evolves, and they impact properties like near-surface temperature and wind speed at the wind farm as well as behind it in the wake. We provide a bug fix and a revised value of C TKE that is one-quarter of the original value. This 0.25 correction factor is empirical; future studies should examine its dependence on parameters such as atmospheric stability, grid resolution, and wind farm layout. We present the results obtained with the Fitch parameterization in the WRF Model for a single turbine with and without the bug fix and the corrected C TKE and compare them with high-fidelity large-eddy simulations. These two issues have not been discovered before because they interact with one another in such a way that their combined effect is a somewhat realistic vertical TKE profile at the wind farm and a realistic wind speed deficit in the wake. All WRF simulations that used the Fitch wind farm parameterization are affected, and their conclusions may need to be revisited.
Abstract
A new one-dimensional 1.5-order planetary boundary layer (PBL) scheme, based on the K–ε turbulence closure applied to the Reynolds-averaged Navier–Stokes (RANS) equations, is developed and implemented within the Weather Research and Forecasting (WRF) Model. The new scheme includes an analytic solution of the coupled equations for turbulent kinetic energy and dissipation rate. Different versions of the PBL scheme are proposed, with increasing levels of complexity, including a model for the calculation of the Prandtl number, a correction to the dissipation rate equation, and a prognostic equation for the temperature variance. Five different idealized cases are tested: four of them explore convective conditions, and they differ in initial thermal stratification and terrain complexity, while one simulates the very stable boundary layer case known as GABLS. For each case study, an ensemble of different large-eddy simulations (LES) is taken as reference for the comparison with the novel PBL schemes and other state-of-the-art 1- and 1.5-order turbulence closures. Results show that the new PBL K–ε scheme brings improvements in all the cases tested in this study. Specifically, the more significant are obtained with the turbulence closure including a prognostic equation for the temperature variance. Moreover, the largest benefits are obtained for the idealized cases simulating a typical thermal circulation within a two-dimensional valley. This suggests that the use of prognostic equations for dissipation rate and temperature variance, which take into account their transport and history, is particularly important with the increasing complexity of PBL dynamics.
Abstract
A new one-dimensional 1.5-order planetary boundary layer (PBL) scheme, based on the K–ε turbulence closure applied to the Reynolds-averaged Navier–Stokes (RANS) equations, is developed and implemented within the Weather Research and Forecasting (WRF) Model. The new scheme includes an analytic solution of the coupled equations for turbulent kinetic energy and dissipation rate. Different versions of the PBL scheme are proposed, with increasing levels of complexity, including a model for the calculation of the Prandtl number, a correction to the dissipation rate equation, and a prognostic equation for the temperature variance. Five different idealized cases are tested: four of them explore convective conditions, and they differ in initial thermal stratification and terrain complexity, while one simulates the very stable boundary layer case known as GABLS. For each case study, an ensemble of different large-eddy simulations (LES) is taken as reference for the comparison with the novel PBL schemes and other state-of-the-art 1- and 1.5-order turbulence closures. Results show that the new PBL K–ε scheme brings improvements in all the cases tested in this study. Specifically, the more significant are obtained with the turbulence closure including a prognostic equation for the temperature variance. Moreover, the largest benefits are obtained for the idealized cases simulating a typical thermal circulation within a two-dimensional valley. This suggests that the use of prognostic equations for dissipation rate and temperature variance, which take into account their transport and history, is particularly important with the increasing complexity of PBL dynamics.
Abstract
Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δx) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δx ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation.
Significance Statement
The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.
Abstract
Generating accurate weather forecasts of planetary boundary layer (PBL) properties is challenging in many geographical regions, oftentimes due to complex topography or horizontal variability in, for example, land characteristics. While recent advances in high-performance computing platforms have led to an increase in the spatial resolution of numerical weather prediction (NWP) models, the horizontal gridcell spacing (Δx) of many regional-scale NWP models currently fall within or are beginning to approach the gray zone (i.e., Δx ≈ 100–1000 m). At these gridcell spacings, three-dimensional (3D) effects are important, as the most energetic turbulent eddies are neither fully parameterized (as in traditional mesoscale simulations) nor fully resolved [as in traditional large-eddy simulations (LES)]. In light of this modeling challenge, we have implemented a 3D PBL parameterization for high-resolution mesoscale simulations using the Weather Research and Forecasting Model. The PBL scheme, which is based on the algebraic model developed by Mellor and Yamada, accounts for the 3D effects of turbulence by calculating explicitly the momentum, heat, and moisture flux divergences in addition to the turbulent kinetic energy. In this study, we present results from idealized simulations in the gray zone that illustrate the benefit of using a fully consistent turbulence closure framework under convective conditions. While the 3D PBL scheme reproduces the evolution of convective features more appropriately than the traditional 1D PBL scheme, we highlight the need to improve the turbulent length scale formulation.
Significance Statement
The spatial resolution of weather models continues to increase at a rapid rate in accordance with the enhancement of computing power. As a result, smaller-scale atmospheric features become more explicitly resolved. However, most numerical models still ignore the impact of horizontal weather variations on boundary layer flows, which becomes more important at these smaller spatial scales. To address this issue, we have implemented a new modeling approach, using fundamental principles, which accounts for horizontal variability. Our results show that including three-dimensional effects of turbulence is necessary to achieve realistic boundary layer characteristics. This novel technique may be useful for many applications including complex terrain flows, pollutant dispersion, and surface–atmosphere interaction studies.
Abstract
Meteorological data of good quality are important for understanding both global and regional climates. In this respect, great efforts have been made to evaluate temperature- and precipitation-related records. This study summarizes the evaluations made to date of the quality of wind speed and direction records acquired at 41 automated weather stations in the northeast of the Iberian Peninsula. Observations were acquired from 1992 to 2005 at a temporal resolution of 10 and 30 min. A quality assurance system was imposed to screen the records for 1) manipulation errors associated with storage and management of the data, 2) consistency limits to ensure that observations are within their natural limits of variation, and 3) temporal consistency to assess abnormally low/high variations in the individual time series. In addition, the most important biases of the dataset are analyzed and corrected wherever possible. A total of 1.8% wind speed and 3.7% wind direction records was assumed invalid, pointing to specific problems in wind measurement. The study not only tries to contribute to the science with the creation of a wind dataset of improved quality, but it also reports on potential errors that could be present in other wind datasets.
Abstract
Meteorological data of good quality are important for understanding both global and regional climates. In this respect, great efforts have been made to evaluate temperature- and precipitation-related records. This study summarizes the evaluations made to date of the quality of wind speed and direction records acquired at 41 automated weather stations in the northeast of the Iberian Peninsula. Observations were acquired from 1992 to 2005 at a temporal resolution of 10 and 30 min. A quality assurance system was imposed to screen the records for 1) manipulation errors associated with storage and management of the data, 2) consistency limits to ensure that observations are within their natural limits of variation, and 3) temporal consistency to assess abnormally low/high variations in the individual time series. In addition, the most important biases of the dataset are analyzed and corrected wherever possible. A total of 1.8% wind speed and 3.7% wind direction records was assumed invalid, pointing to specific problems in wind measurement. The study not only tries to contribute to the science with the creation of a wind dataset of improved quality, but it also reports on potential errors that could be present in other wind datasets.
Abstract
This study presents a systematic analysis of convective parameterizations performance with interactive radiation, microphysics, and surface on an idealized day with shallow convection. To this end, we analyze a suite of mesoscale numerical experiments (i.e., with parameterized turbulence). In the first set, two different convection schemes represent shallow convection at a 9-km resolution. These experiments are then compared with model results omitting convective parameterizations at 9- and 3-km horizontal resolution (gray zone). Relevant in our approach is to compare the results against two simulations by different large-eddy simulation (LES) models. Results show that the mesoscale experiments, including the 3-km resolution, are unable to adequately represent the timing, intensity, height, and extension of the shallow cumulus field. The main differences with LES experiments are the following: a too late onset, too high cloud base, and a too early transport of moisture too high, overestimating the second cloud layer. Related to this, both convective parameterizations produce warm and dry biases of up to 2 K and 2 g kg−1, respectively, in the cloud layer. This misrepresentation of the cloud dynamics leads to overestimated shortwave radiation variability, both spacewise and timewise. Domain-averaged shortwave radiation at the surface, however, compares satisfactorily with LES. The shortwave direct and diffuse partition is misrepresented by the convective parameterizations with an underestimation (overestimation) of diffuse (direct) radiation both locally and, by a relative 40% (10%), of the domain average.
Abstract
This study presents a systematic analysis of convective parameterizations performance with interactive radiation, microphysics, and surface on an idealized day with shallow convection. To this end, we analyze a suite of mesoscale numerical experiments (i.e., with parameterized turbulence). In the first set, two different convection schemes represent shallow convection at a 9-km resolution. These experiments are then compared with model results omitting convective parameterizations at 9- and 3-km horizontal resolution (gray zone). Relevant in our approach is to compare the results against two simulations by different large-eddy simulation (LES) models. Results show that the mesoscale experiments, including the 3-km resolution, are unable to adequately represent the timing, intensity, height, and extension of the shallow cumulus field. The main differences with LES experiments are the following: a too late onset, too high cloud base, and a too early transport of moisture too high, overestimating the second cloud layer. Related to this, both convective parameterizations produce warm and dry biases of up to 2 K and 2 g kg−1, respectively, in the cloud layer. This misrepresentation of the cloud dynamics leads to overestimated shortwave radiation variability, both spacewise and timewise. Domain-averaged shortwave radiation at the surface, however, compares satisfactorily with LES. The shortwave direct and diffuse partition is misrepresented by the convective parameterizations with an underestimation (overestimation) of diffuse (direct) radiation both locally and, by a relative 40% (10%), of the domain average.
