Browse
Abstract
Modification of grasslands into irrigated and non-irrigated agriculture in the Great Plains results in significant impacts on weather and climate. However, there has been lack of observational data-based studies solely focused on impacts of irrigation on the PBL and convective conditions. The Great Plains Irrigation Experiment (GRAINEX) during the 2018 growing season collected data over irrigated and non-irrigated land uses over Nebraska to understand these impacts. Specifically, the objective was to determine whether the impacts of irrigation are sustained throughout the growing season.
The data analyzed include latent and sensible heat flux, air temperature, dew point temperature, equivalent temperature (moist enthalpy), PBL height, lifting condensation level (LCL), level of free convection (LFC), and PBL mixing ratio. Results show increased partitioning of energy into latent heat compared to sensible heat over irrigated areas while average maximum air was decreased and dewpoint temperature was increased from the early to peak growing season. Radiosonde data suggest reduced planetary boundary layer (PBL) heights at all launch sites from the early to peak growing season. However, reduction of PBL height was much greater over irrigated areas compared to non-irrigated croplands. Compared to the early growing period, LCL and LFC heights were also lower during the peak growing period over irrigated areas. Results note, for the first time, that the impacts of irrigation on PBL evolution and convective environment can be sustained throughout the growing season and regardless of background atmospheric conditions. These are important findings and applicable to other irrigated areas in the world.
Abstract
Modification of grasslands into irrigated and non-irrigated agriculture in the Great Plains results in significant impacts on weather and climate. However, there has been lack of observational data-based studies solely focused on impacts of irrigation on the PBL and convective conditions. The Great Plains Irrigation Experiment (GRAINEX) during the 2018 growing season collected data over irrigated and non-irrigated land uses over Nebraska to understand these impacts. Specifically, the objective was to determine whether the impacts of irrigation are sustained throughout the growing season.
The data analyzed include latent and sensible heat flux, air temperature, dew point temperature, equivalent temperature (moist enthalpy), PBL height, lifting condensation level (LCL), level of free convection (LFC), and PBL mixing ratio. Results show increased partitioning of energy into latent heat compared to sensible heat over irrigated areas while average maximum air was decreased and dewpoint temperature was increased from the early to peak growing season. Radiosonde data suggest reduced planetary boundary layer (PBL) heights at all launch sites from the early to peak growing season. However, reduction of PBL height was much greater over irrigated areas compared to non-irrigated croplands. Compared to the early growing period, LCL and LFC heights were also lower during the peak growing period over irrigated areas. Results note, for the first time, that the impacts of irrigation on PBL evolution and convective environment can be sustained throughout the growing season and regardless of background atmospheric conditions. These are important findings and applicable to other irrigated areas in the world.
Abstract
This study investigates the connection between the arrival of dry stratospheric air with the Soberanes Fire (2016). The Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT) and Goddard Earth Observing System Forward Processing model (GEOS-FP) are used for back-trajectories and offshore deep stratospheric intrusion (SI) in conjunction with the ignition and outbreak of the fire. The back-trajectory analysis indicates most air reaching the vertical column was critically dry, exhibiting relative humidity values below 10%. As the fire ignited, dry air arrived from due west at heights of 1-3 km about 24 hours prior. During the overnight fire growth, dry air arrived from the northwest to north-northwest at heights of 3.5-5.5 km 48-72 hours prior. The synoptic and the GEOS-FP analysis demonstrate offshore mid-to-low stratospheric intrusion. On July 21, 2016, an enclosed upper-level low approached the California/Oregon border along the northwesterly subtropical jet stream hours before the fire outbreak. The GEOS-FP results of potential vorticity, specific humidity, and ozone along the back-trajectories to the west and northwest of the fire suggest a stratospheric intrusion event into the mid-to-low troposphere at the back-trajectory start points, and vertical velocity indicates sinking motion. The specific humidity analyzed at the arrival time shows the transport of the abnormally dry air to the Soberanes Fire. Results suggest a connection between dry stratospheric air transported to the Soberanes Fire at ignition and overnight accelerated growth, supported by a dark bank in satellite water vapor imagery. The prediction of low-level transport of dry stratospheric air to the coastal communities could help predict the occurrence of wildfire outbreaks, or periods of accelerated fire growth.
Abstract
This study investigates the connection between the arrival of dry stratospheric air with the Soberanes Fire (2016). The Hybrid Single-Particle Lagrangian Integrated Trajectory model (HYSPLIT) and Goddard Earth Observing System Forward Processing model (GEOS-FP) are used for back-trajectories and offshore deep stratospheric intrusion (SI) in conjunction with the ignition and outbreak of the fire. The back-trajectory analysis indicates most air reaching the vertical column was critically dry, exhibiting relative humidity values below 10%. As the fire ignited, dry air arrived from due west at heights of 1-3 km about 24 hours prior. During the overnight fire growth, dry air arrived from the northwest to north-northwest at heights of 3.5-5.5 km 48-72 hours prior. The synoptic and the GEOS-FP analysis demonstrate offshore mid-to-low stratospheric intrusion. On July 21, 2016, an enclosed upper-level low approached the California/Oregon border along the northwesterly subtropical jet stream hours before the fire outbreak. The GEOS-FP results of potential vorticity, specific humidity, and ozone along the back-trajectories to the west and northwest of the fire suggest a stratospheric intrusion event into the mid-to-low troposphere at the back-trajectory start points, and vertical velocity indicates sinking motion. The specific humidity analyzed at the arrival time shows the transport of the abnormally dry air to the Soberanes Fire. Results suggest a connection between dry stratospheric air transported to the Soberanes Fire at ignition and overnight accelerated growth, supported by a dark bank in satellite water vapor imagery. The prediction of low-level transport of dry stratospheric air to the coastal communities could help predict the occurrence of wildfire outbreaks, or periods of accelerated fire growth.
Abstract
Measurements from the South Dakota School of Mines and Technology T-28 hail-penetrating aircraft are analyzed using recently developed data processing techniques with the goals of identifying where the large hail is found relative to vertical motion and improving the detection of hail microphysical properties from radar. Hail particle size distributions (PSD) and environmental conditions (temperature, relative humidity, liquid water content, air vertical velocity) were digitally collected by the T28 between 1995 and 2003 and synthesized by Detwiler et al. The PSD were forward modeled by Cecchini et al. to simulate the radar reflectivity of the PSD at multiple radar wavelengths. The T-28 penetrated temperatures primarily between 0° and −10°C. The largest hailstones were sampled near the updraft/downdraft interface. Liquid water contents were highest in the updraft cores, whereas total (liquid + frozen) water contents were highest near the updraft/downdraft interface. The fitted properties of the PSD (intercept and slope) are directly related to each other but do not show any dependence on the region of the hailstorm where sampled. The PSD measurements and the radar reflectivity calculations at multiple radar wavelengths facilitated the development of relationships between the PSD bulk properties—hail kinetic energy and kinetic energy flux—and the radar reflectivity. Rather than using the oft-assumed sphericity and solid ice physical properties, actual measurements of hail properties are used in the analysis. Results from the maximum estimated size of hail (MESH) and vertical integrated liquid water (VIL) algorithms are evaluated based on this analysis.
Significance Statement
Hailstorms in the United States have caused over $10 billion in damage for each of the last 14 years, according to insurance industry estimates (). Algorithms have been developed to identify the presence and size of hail from radar. Numerical simulations of hailstorms have improved significantly since the 1970s, and further improvements will provide better resolution and more accurate estimates of the sizes of hailstones falling to the ground. Measurements of the properties of hailstones—their mass and terminal velocities—have been improved in recent years but in general are not incorporated in the algorithms developed for radar estimates of hail sizes or for the properties of hail used in the model simulations. This study synthesizes in situ aircraft data, computed radar backscatter cross sections, together with recent estimates of the physical characteristics of hailstones to improve the representation of hail in numerical models and quantitative assessment hail properties in storms using weather radar.
Abstract
Measurements from the South Dakota School of Mines and Technology T-28 hail-penetrating aircraft are analyzed using recently developed data processing techniques with the goals of identifying where the large hail is found relative to vertical motion and improving the detection of hail microphysical properties from radar. Hail particle size distributions (PSD) and environmental conditions (temperature, relative humidity, liquid water content, air vertical velocity) were digitally collected by the T28 between 1995 and 2003 and synthesized by Detwiler et al. The PSD were forward modeled by Cecchini et al. to simulate the radar reflectivity of the PSD at multiple radar wavelengths. The T-28 penetrated temperatures primarily between 0° and −10°C. The largest hailstones were sampled near the updraft/downdraft interface. Liquid water contents were highest in the updraft cores, whereas total (liquid + frozen) water contents were highest near the updraft/downdraft interface. The fitted properties of the PSD (intercept and slope) are directly related to each other but do not show any dependence on the region of the hailstorm where sampled. The PSD measurements and the radar reflectivity calculations at multiple radar wavelengths facilitated the development of relationships between the PSD bulk properties—hail kinetic energy and kinetic energy flux—and the radar reflectivity. Rather than using the oft-assumed sphericity and solid ice physical properties, actual measurements of hail properties are used in the analysis. Results from the maximum estimated size of hail (MESH) and vertical integrated liquid water (VIL) algorithms are evaluated based on this analysis.
Significance Statement
Hailstorms in the United States have caused over $10 billion in damage for each of the last 14 years, according to insurance industry estimates (). Algorithms have been developed to identify the presence and size of hail from radar. Numerical simulations of hailstorms have improved significantly since the 1970s, and further improvements will provide better resolution and more accurate estimates of the sizes of hailstones falling to the ground. Measurements of the properties of hailstones—their mass and terminal velocities—have been improved in recent years but in general are not incorporated in the algorithms developed for radar estimates of hail sizes or for the properties of hail used in the model simulations. This study synthesizes in situ aircraft data, computed radar backscatter cross sections, together with recent estimates of the physical characteristics of hailstones to improve the representation of hail in numerical models and quantitative assessment hail properties in storms using weather radar.
Abstract
We use a simple risk model for U.S. hurricane wind and surge economic damage to investigate the impact of projected changes in the frequencies of hurricanes of different intensities due to climate change. For average annual damage, we find that changes in the frequency of category-4 storms dominate. For distributions of annual damage, we find that changes in the frequency of category-4 storms again dominate for all except the shortest return periods. Sensitivity tests show that accounting for landfall, uncertainties, and correlations leads to increases in damage estimates. When we propagate the distributions of uncertain frequency changes to give a best estimate of the changes in damage, the changes are moderate. When we pick individual scenarios from within the distributions of frequency changes, we find a significant probability of much larger changes in damage. The inputs on which our study depends are highly uncertain, and our methods are approximate, leading to high levels of uncertainty in our results. Also, the damage changes we consider are only part of the total possible change in hurricane damage due to climate change. Total damage change estimates would also need to include changes due to other factors, including possible changes in genesis, tracks, size, forward speed, sea level, rainfall, and exposure. Nevertheless, we believe that our results give important new insights into U.S. hurricane risk under climate change.
Significance Statement
We investigate how changes in the frequencies of hurricanes of different intensities as a result of climate change may contribute to changes in U.S. economic damage due to wind and surge. We find that economic damage will likely increase as a result of projected increases in the frequency of landfalling hurricanes. Analysis of our results shows that increases in the frequency of category-4 storms are the main driver of the changes. Our best estimate results, based on a multimodel ensemble, give modest increases in damage, but within the ensemble there are individual scenarios that give much larger increases in damage. The large range of individual damage estimates is a motivation for continuing efforts to reduce the uncertainty around hurricane projections under climate change.
Abstract
We use a simple risk model for U.S. hurricane wind and surge economic damage to investigate the impact of projected changes in the frequencies of hurricanes of different intensities due to climate change. For average annual damage, we find that changes in the frequency of category-4 storms dominate. For distributions of annual damage, we find that changes in the frequency of category-4 storms again dominate for all except the shortest return periods. Sensitivity tests show that accounting for landfall, uncertainties, and correlations leads to increases in damage estimates. When we propagate the distributions of uncertain frequency changes to give a best estimate of the changes in damage, the changes are moderate. When we pick individual scenarios from within the distributions of frequency changes, we find a significant probability of much larger changes in damage. The inputs on which our study depends are highly uncertain, and our methods are approximate, leading to high levels of uncertainty in our results. Also, the damage changes we consider are only part of the total possible change in hurricane damage due to climate change. Total damage change estimates would also need to include changes due to other factors, including possible changes in genesis, tracks, size, forward speed, sea level, rainfall, and exposure. Nevertheless, we believe that our results give important new insights into U.S. hurricane risk under climate change.
Significance Statement
We investigate how changes in the frequencies of hurricanes of different intensities as a result of climate change may contribute to changes in U.S. economic damage due to wind and surge. We find that economic damage will likely increase as a result of projected increases in the frequency of landfalling hurricanes. Analysis of our results shows that increases in the frequency of category-4 storms are the main driver of the changes. Our best estimate results, based on a multimodel ensemble, give modest increases in damage, but within the ensemble there are individual scenarios that give much larger increases in damage. The large range of individual damage estimates is a motivation for continuing efforts to reduce the uncertainty around hurricane projections under climate change.
Abstract
We assess site-specific surface short-wave radiation forecasts from two high resolution configurations of the South African Weather Service numerical weather prediction model, at 4 km and 1.5 km. The models exhibit good skill overall in forecasting surface short-wave radiation, with zero median error for all radiation components. This information is relevant to support a growing Renewable Energy sector in South Africa, particularly for photovoltaics. Further model performance analysis has shown an imbalance between cloud and solar radiation forecasting errors. In addition, cloud over-prediction does not necessarily equate to under-estimating solar radiation. Overcast cloud regimes are predicted too often with an associated positive mean radiation bias, whereas the relative abundance of partly cloudy regimes is under-predicted by the models with mixed radiation biases. Challenges highlighted by the misrepresentation of partly cloudy regimes in solar radiation error attribution may be used to inform improvements to the numerical core, namely the cloud and radiation schemes.
Abstract
We assess site-specific surface short-wave radiation forecasts from two high resolution configurations of the South African Weather Service numerical weather prediction model, at 4 km and 1.5 km. The models exhibit good skill overall in forecasting surface short-wave radiation, with zero median error for all radiation components. This information is relevant to support a growing Renewable Energy sector in South Africa, particularly for photovoltaics. Further model performance analysis has shown an imbalance between cloud and solar radiation forecasting errors. In addition, cloud over-prediction does not necessarily equate to under-estimating solar radiation. Overcast cloud regimes are predicted too often with an associated positive mean radiation bias, whereas the relative abundance of partly cloudy regimes is under-predicted by the models with mixed radiation biases. Challenges highlighted by the misrepresentation of partly cloudy regimes in solar radiation error attribution may be used to inform improvements to the numerical core, namely the cloud and radiation schemes.
Abstract
Wet bulb globe temperature (WBGT) is used to assess environmental heat stress and accounts for the influences of air temperature, humidity, wind speed, and radiation on heat stress. Measurements of WBGT are highly sensitive to slight changes in environmental conditions and can vary several degrees Celsius across small distances (10s to 100s of meters). Compared to observations with an ISO-compliant WBGT meter, this work assesses the accuracy of WBGT measurements made with a popular handheld meter (the Kestrel 5400 Heat Stress Tracker) and WBGT estimates. Measurements were made during the summers of 2019-2021 in a variety of suburban and urban environments in North Carolina, including three high school campuses. WBGT can be estimated from standard weather station variables, and many of these stations report cloud cover in lieu of solar radiation. Therefore, this work also evaluates the accuracy of clear-sky radiation estimates and adjustments to those estimates based on cloud cover. WBGT estimated with the method from Liljegren et al. (2008) from a weather station were on average 0.2°C warmer than observed WBGT, while the Kestrel 5400 WBGT was 0.7°C warmer. Large variations in WBGT were observed across surfaces and shade conditions, with differences of 0.9°C (0.3–1.4°C) between a tennis court and a neighboring grass field. The method for estimating clear-sky radiation in Ryan & Stolzenbach (1972) was most accurate and the clear-sky radiation modified by percentage cloud cover was found to be within 75 w/m2 of observations on average.
Abstract
Wet bulb globe temperature (WBGT) is used to assess environmental heat stress and accounts for the influences of air temperature, humidity, wind speed, and radiation on heat stress. Measurements of WBGT are highly sensitive to slight changes in environmental conditions and can vary several degrees Celsius across small distances (10s to 100s of meters). Compared to observations with an ISO-compliant WBGT meter, this work assesses the accuracy of WBGT measurements made with a popular handheld meter (the Kestrel 5400 Heat Stress Tracker) and WBGT estimates. Measurements were made during the summers of 2019-2021 in a variety of suburban and urban environments in North Carolina, including three high school campuses. WBGT can be estimated from standard weather station variables, and many of these stations report cloud cover in lieu of solar radiation. Therefore, this work also evaluates the accuracy of clear-sky radiation estimates and adjustments to those estimates based on cloud cover. WBGT estimated with the method from Liljegren et al. (2008) from a weather station were on average 0.2°C warmer than observed WBGT, while the Kestrel 5400 WBGT was 0.7°C warmer. Large variations in WBGT were observed across surfaces and shade conditions, with differences of 0.9°C (0.3–1.4°C) between a tennis court and a neighboring grass field. The method for estimating clear-sky radiation in Ryan & Stolzenbach (1972) was most accurate and the clear-sky radiation modified by percentage cloud cover was found to be within 75 w/m2 of observations on average.
Abstract
This study focuses on the application of two standard inflow turbulence generation methods for growing convective boundary layer (CBL) simulations: the recycle–rescale (R-R) and the digital filter–based (DF) methods, which are used in computational fluid dynamics. The primary objective of this study is to expand the applicability of the R-R method to simulations of thermally driven CBLs. This method is called the extended R-R method. However, in previous studies, the DF method has been extended to generate potential temperature perturbations. This study investigated whether the extended DF method can be applied to simulations of growing thermally driven CBLs. In this study, idealized simulations of growing thermally driven CBLs using the extended R-R and DF methods were performed. The results showed that both extended methods could capture the characteristics of thermally driven CBLs. The extended R-R method reproduced turbulence in thermally driven CBLs better than the extended DF method in the spectrum and histogram of vertical wind speed. However, the height of the thermally driven CBL was underestimated in about 100 m compared with the extended DF method. Sensitivity experiments were conducted on the parameters used in the extended DF and R-R methods. The results showed that underestimation of the length scale in the extended DF method causes a shortage of large-scale turbulence components. The other point suggested by the results of the sensitivity experiments is that the length of the driver region in the extended R-R method should be sufficient to reproduce the spanwise movement of the roll vortices.
Significance Statement
Inflow turbulence generation methods for large-eddy simulation (LES) models are crucial for the better downscaling of meteorological mesoscale models (RANS models) to microscale models (LES models). Various CFD methods have been developed, but few have been applied to simulations of thermally driven convective boundary layers (CBLs). To address this problem, we focused on a method that recycles turbulence [the recycle–rescale (R-R) method] and another method that synthetically generates turbulence [the digital filter–based (DF) method]. This study extends the R-R method to manage turbulence in thermally driven CBLs. In addition, this study investigated the applicability of the DF method to thermally driven CBL simulations. Both extended methods are effective for downscaling experiments and capture the characteristics of thermally driven CBLs.
Abstract
This study focuses on the application of two standard inflow turbulence generation methods for growing convective boundary layer (CBL) simulations: the recycle–rescale (R-R) and the digital filter–based (DF) methods, which are used in computational fluid dynamics. The primary objective of this study is to expand the applicability of the R-R method to simulations of thermally driven CBLs. This method is called the extended R-R method. However, in previous studies, the DF method has been extended to generate potential temperature perturbations. This study investigated whether the extended DF method can be applied to simulations of growing thermally driven CBLs. In this study, idealized simulations of growing thermally driven CBLs using the extended R-R and DF methods were performed. The results showed that both extended methods could capture the characteristics of thermally driven CBLs. The extended R-R method reproduced turbulence in thermally driven CBLs better than the extended DF method in the spectrum and histogram of vertical wind speed. However, the height of the thermally driven CBL was underestimated in about 100 m compared with the extended DF method. Sensitivity experiments were conducted on the parameters used in the extended DF and R-R methods. The results showed that underestimation of the length scale in the extended DF method causes a shortage of large-scale turbulence components. The other point suggested by the results of the sensitivity experiments is that the length of the driver region in the extended R-R method should be sufficient to reproduce the spanwise movement of the roll vortices.
Significance Statement
Inflow turbulence generation methods for large-eddy simulation (LES) models are crucial for the better downscaling of meteorological mesoscale models (RANS models) to microscale models (LES models). Various CFD methods have been developed, but few have been applied to simulations of thermally driven convective boundary layers (CBLs). To address this problem, we focused on a method that recycles turbulence [the recycle–rescale (R-R) method] and another method that synthetically generates turbulence [the digital filter–based (DF) method]. This study extends the R-R method to manage turbulence in thermally driven CBLs. In addition, this study investigated the applicability of the DF method to thermally driven CBL simulations. Both extended methods are effective for downscaling experiments and capture the characteristics of thermally driven CBLs.
Abstract
We investigated 9 indices of Spatio-temporal extreme precipitation events over Nicaragua during 2001-2016, from GPCC, CHIRPS, and IMERG, and their correlation with teleconnection patterns. The main objectives were to evaluate the variability of extreme precipitation events, to know the performance of IMERG and CHIRPS in the characterization of these extreme events, using GPCC and 4 rain gauges as references, and finally to determine the teleconnection patterns that have the highest correlation with these indices. The spatial coverage of the area with the highest number of consecutive days with daily precipitation less than 1 mm corresponds to the Pacific region, with annual mean values of up to 120 continuous days. Some extreme precipitation event indices (RR, RX1day, and RX5day) show a decreasing trend, suggesting that the study area has been experiencing a reduction of extreme precipitation indices in terms of intensities and duration throughout the study period. In addition, it was observed that CHIRPS shows a better fit when dealing with precipitation events that do not exceed certain thresholds and IMERG improves when describing intense precipitation event patterns. We found that EOFPAC, NIÑA 3.4, PACWARM, and SOI have a greater influence on extreme precipitation events, these results suggest that they are being controlled by ENSO episodes, providing a better understanding of the climate configuration, as a prediction and forecasting potential, useful for agriculture, land use and risk management.
Abstract
We investigated 9 indices of Spatio-temporal extreme precipitation events over Nicaragua during 2001-2016, from GPCC, CHIRPS, and IMERG, and their correlation with teleconnection patterns. The main objectives were to evaluate the variability of extreme precipitation events, to know the performance of IMERG and CHIRPS in the characterization of these extreme events, using GPCC and 4 rain gauges as references, and finally to determine the teleconnection patterns that have the highest correlation with these indices. The spatial coverage of the area with the highest number of consecutive days with daily precipitation less than 1 mm corresponds to the Pacific region, with annual mean values of up to 120 continuous days. Some extreme precipitation event indices (RR, RX1day, and RX5day) show a decreasing trend, suggesting that the study area has been experiencing a reduction of extreme precipitation indices in terms of intensities and duration throughout the study period. In addition, it was observed that CHIRPS shows a better fit when dealing with precipitation events that do not exceed certain thresholds and IMERG improves when describing intense precipitation event patterns. We found that EOFPAC, NIÑA 3.4, PACWARM, and SOI have a greater influence on extreme precipitation events, these results suggest that they are being controlled by ENSO episodes, providing a better understanding of the climate configuration, as a prediction and forecasting potential, useful for agriculture, land use and risk management.
Abstract
We investigated the ability of three planetary boundary layer (PBL) schemes in the Weather Research and Forecasting (WRF) model to simulate boundary layer turbulence in the “grey zone” (i.e. resolutions from 100 m to 1 km). The three schemes chosen are the well-established MYNN PBL scheme and the two newest PBL schemes added to WRF: the SMS-3DTKE scheme and the EEPS scheme. The SMS-3DTKE scheme is designed to be scale-aware and avoid the double-counting of turbulent kinetic energy (TKE) in simulations within the grey zone. We evaluated their performance using aircraft measurements obtained during three research flights immediately downwind of Manhattan, New York City. The MYNN PBL scheme simulates TKE best, despite not being scale-aware and slightly underestimating TKE from observations, while the SMS-3DTKE scheme appears to be overly scale-aware for the three flights examined, in particular when combined with the MM5 surface layer scheme. The EEPS scheme significantly underestimates TKE, mostly in the elevated layers of the boundary layer. Additionally, we examined the impact of flow over tall buildings on observed TKE and found that only the windiest day showed a significant increase in TKE directly downwind of Manhattan. This impact was, however, not reproduced by any of the model configurations, regardless of the land use data selected, although the better resolved NLCD land use led to a slight improvement of the spatial distribution of TKE, implying that more explicit representation of the impact of tall buildings may be needed to fully capture their impact on boundary layer turbulence.
Abstract
We investigated the ability of three planetary boundary layer (PBL) schemes in the Weather Research and Forecasting (WRF) model to simulate boundary layer turbulence in the “grey zone” (i.e. resolutions from 100 m to 1 km). The three schemes chosen are the well-established MYNN PBL scheme and the two newest PBL schemes added to WRF: the SMS-3DTKE scheme and the EEPS scheme. The SMS-3DTKE scheme is designed to be scale-aware and avoid the double-counting of turbulent kinetic energy (TKE) in simulations within the grey zone. We evaluated their performance using aircraft measurements obtained during three research flights immediately downwind of Manhattan, New York City. The MYNN PBL scheme simulates TKE best, despite not being scale-aware and slightly underestimating TKE from observations, while the SMS-3DTKE scheme appears to be overly scale-aware for the three flights examined, in particular when combined with the MM5 surface layer scheme. The EEPS scheme significantly underestimates TKE, mostly in the elevated layers of the boundary layer. Additionally, we examined the impact of flow over tall buildings on observed TKE and found that only the windiest day showed a significant increase in TKE directly downwind of Manhattan. This impact was, however, not reproduced by any of the model configurations, regardless of the land use data selected, although the better resolved NLCD land use led to a slight improvement of the spatial distribution of TKE, implying that more explicit representation of the impact of tall buildings may be needed to fully capture their impact on boundary layer turbulence.
Abstract
This study investigates the impacts of grid spacing and station network on surface analyses and forecasts including temperature, humidity, and winds in Beijing Winter Olympic complex terrain. The high-resolution analyses are generated by a rapid-refresh integrated system that includes a topographic downscaling procedure. Results show that surface analyses are more accurate with a higher targeted grid spacing. In particular, the average analysis errors of surface temperature, humidity, and winds are all significantly reduced when the grid size is increased. This improvement is mainly attributed to a more realistic simulation of the topographic effects in the integrated system because the topographic downscaling at higher grid spacing can add more details in a complex mountain region. From 1 km to 100 m, 1–12-h forecasts of temperature and humidity are also largely improved, while the wind only shows a slight improvement for 1–6-h forecasts. The influence of station network on the surface analyses is further examined. Results show that the spatial distributions of temperature and humidity at a 100-m space scale are more realistic and accurate when adding an intensive automatic weather station network, as more observational information can be absorbed. The adding of a station network can also reduce forecast errors, which can last for about 6 h. However, although surface winds display better analysis skill when more stations are added, the wind at the mountaintop region sometimes encounters a marginally worse effect for both analysis and forecast. The results are helpful to improve the analysis and forecast products in complex terrain and have some implications for downscaling from a coarse grid size to a finer grid.
Abstract
This study investigates the impacts of grid spacing and station network on surface analyses and forecasts including temperature, humidity, and winds in Beijing Winter Olympic complex terrain. The high-resolution analyses are generated by a rapid-refresh integrated system that includes a topographic downscaling procedure. Results show that surface analyses are more accurate with a higher targeted grid spacing. In particular, the average analysis errors of surface temperature, humidity, and winds are all significantly reduced when the grid size is increased. This improvement is mainly attributed to a more realistic simulation of the topographic effects in the integrated system because the topographic downscaling at higher grid spacing can add more details in a complex mountain region. From 1 km to 100 m, 1–12-h forecasts of temperature and humidity are also largely improved, while the wind only shows a slight improvement for 1–6-h forecasts. The influence of station network on the surface analyses is further examined. Results show that the spatial distributions of temperature and humidity at a 100-m space scale are more realistic and accurate when adding an intensive automatic weather station network, as more observational information can be absorbed. The adding of a station network can also reduce forecast errors, which can last for about 6 h. However, although surface winds display better analysis skill when more stations are added, the wind at the mountaintop region sometimes encounters a marginally worse effect for both analysis and forecast. The results are helpful to improve the analysis and forecast products in complex terrain and have some implications for downscaling from a coarse grid size to a finer grid.