Browse
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.
Abstract
This paper investigates the surface-layer processes associated with the morning transition from nighttime downslope winds to daytime upslope winds over a semi-isolated massif. It provides an insight into the characteristics of the transition and its connection with the processes controlling the erosion of the temperature inversion at the foot of the slope. First, a criterion for the identification of days prone to the development of purely thermally driven slope winds is proposed and adopted to select five representative case studies. Then, the mechanisms leading to different patterns of erosion of the nocturnal temperature inversion at the foot of the slope are analyzed. Three main patterns of erosion are identified: the first is connected to the growth of the convective boundary layer at the surface, the second is connected to the descent of the inversion top, and the third is a combination of the previous two. The first pattern is linked to the initiation of the morning transition through surface heating, and the second pattern is connected to the top-down dilution mechanism and so to mixing with the above air. The discriminating factor in the determination of the erosion pattern is identified in the partitioning of turbulent sensible heat flux at the surface.
Significance Statement
The purpose of this study is to improve our understanding of the thermally driven slope circulations with a focus on the unsteady processes associated with the morning transition and the erosion patterns of the nocturnal temperature inversion, so far in the literature less investigated and understood than the evening transition. Understanding this diurnal process will advance our abilities to model it and to improve the accuracy of weather forecasting in complex terrain.
Abstract
This paper investigates the surface-layer processes associated with the morning transition from nighttime downslope winds to daytime upslope winds over a semi-isolated massif. It provides an insight into the characteristics of the transition and its connection with the processes controlling the erosion of the temperature inversion at the foot of the slope. First, a criterion for the identification of days prone to the development of purely thermally driven slope winds is proposed and adopted to select five representative case studies. Then, the mechanisms leading to different patterns of erosion of the nocturnal temperature inversion at the foot of the slope are analyzed. Three main patterns of erosion are identified: the first is connected to the growth of the convective boundary layer at the surface, the second is connected to the descent of the inversion top, and the third is a combination of the previous two. The first pattern is linked to the initiation of the morning transition through surface heating, and the second pattern is connected to the top-down dilution mechanism and so to mixing with the above air. The discriminating factor in the determination of the erosion pattern is identified in the partitioning of turbulent sensible heat flux at the surface.
Significance Statement
The purpose of this study is to improve our understanding of the thermally driven slope circulations with a focus on the unsteady processes associated with the morning transition and the erosion patterns of the nocturnal temperature inversion, so far in the literature less investigated and understood than the evening transition. Understanding this diurnal process will advance our abilities to model it and to improve the accuracy of weather forecasting in complex terrain.
Abstract
In the realm of boundary layer flows in complex terrain, low-level jets (LLJs) have received considerable attention, although little literature is available for double-nosed LLJs that remain not well understood. To this end, we use the Mountain Terrain Atmospheric Modeling and Observations (MATERHORN) dataset to demonstrate that double-nosed LLJs developing within the planetary boundary layer (PBL) are common during stable nocturnal conditions and present two possible mechanisms responsible for their formation. It is observed that the onset of a double-nosed LLJ is associated with a temporary shape modification of an already-established LLJ. The characteristics of these double-nosed LLJs are described using a refined version of identification criteria proposed in the literature, and their formation is classified in terms of two driving mechanisms. The wind-driven mechanism encompasses cases where the two noses are associated with different air masses flowing one on top of the other. The wave-driven mechanism involves the vertical momentum transport by an inertial–gravity wave to generate the second nose. The wave-driven mechanism is corroborated by the analysis of nocturnal double-nosed LLJs, where inertial–gravity waves are generated close to the ground by a sudden flow perturbation.
Abstract
In the realm of boundary layer flows in complex terrain, low-level jets (LLJs) have received considerable attention, although little literature is available for double-nosed LLJs that remain not well understood. To this end, we use the Mountain Terrain Atmospheric Modeling and Observations (MATERHORN) dataset to demonstrate that double-nosed LLJs developing within the planetary boundary layer (PBL) are common during stable nocturnal conditions and present two possible mechanisms responsible for their formation. It is observed that the onset of a double-nosed LLJ is associated with a temporary shape modification of an already-established LLJ. The characteristics of these double-nosed LLJs are described using a refined version of identification criteria proposed in the literature, and their formation is classified in terms of two driving mechanisms. The wind-driven mechanism encompasses cases where the two noses are associated with different air masses flowing one on top of the other. The wave-driven mechanism involves the vertical momentum transport by an inertial–gravity wave to generate the second nose. The wave-driven mechanism is corroborated by the analysis of nocturnal double-nosed LLJs, where inertial–gravity waves are generated close to the ground by a sudden flow perturbation.
Abstract
This study examines the sensitivity of numerical simulations of near-surface atmospheric conditions to the initial surface albedo and snow depth during an observed ice fog event in the Heber Valley of northern Utah. Numerical simulation results from the mesoscale community Weather Research and Forecasting (WRF) Model are compared with observations from the Mountain Terrain Atmospheric Modeling and Observations (MATERHORN) Program fog field program. It is found that near-surface cooling during the nighttime is significantly underestimated by the WRF Model, resulting in the failure of the model to reproduce the observed fog episode. Meanwhile, the model also overestimates the temperature during the daytime. Nevertheless, these errors could be reduced by increasing the initial surface albedo and snow depth, which act to cool the near-surface atmosphere by increasing the reflection of downward shortwave radiation and decreasing the heating effects from the soil layer. Overall results indicate the important effects of snow representation on the simulation of near-surface atmospheric conditions and highlight the need for snow measurements in the cold season for improved model physics parameterizations.
Abstract
This study examines the sensitivity of numerical simulations of near-surface atmospheric conditions to the initial surface albedo and snow depth during an observed ice fog event in the Heber Valley of northern Utah. Numerical simulation results from the mesoscale community Weather Research and Forecasting (WRF) Model are compared with observations from the Mountain Terrain Atmospheric Modeling and Observations (MATERHORN) Program fog field program. It is found that near-surface cooling during the nighttime is significantly underestimated by the WRF Model, resulting in the failure of the model to reproduce the observed fog episode. Meanwhile, the model also overestimates the temperature during the daytime. Nevertheless, these errors could be reduced by increasing the initial surface albedo and snow depth, which act to cool the near-surface atmosphere by increasing the reflection of downward shortwave radiation and decreasing the heating effects from the soil layer. Overall results indicate the important effects of snow representation on the simulation of near-surface atmospheric conditions and highlight the need for snow measurements in the cold season for improved model physics parameterizations.
Abstract
Topographic effects on radiation, including both topographic shading and slope effects, are included in the Weather Research and Forecasting (WRF) Model, and here they are made compatible with the immersed boundary method (IBM). IBM is an alternative method for representing complex terrain that reduces numerical errors over sloped terrain, thus extending the range of slopes that can be represented in WRF simulations. The implementation of topographic effects on radiation is validated by comparing land surface fluxes, as well as temperature and velocity fields, between idealized WRF simulations both with and without IBM. Following validation, the topographic shading implementation is tested in a semirealistic simulation of flow over Granite Mountain, Utah, where topographic shading is known to affect downslope flow development in the evening. The horizontal grid spacing is 50 m and the vertical grid spacing is approximately 8–27 m near the surface. Such a case would fail to run in WRF with its native terrain-following coordinates because of large local slope values reaching up to 55°. Good agreement is found between modeled surface energy budget components and observations from the Mountain Terrain Atmospheric Modeling and Observations (MATERHORN) program at a location on the east slope of Granite Mountain. In addition, the model captures large spatiotemporal inhomogeneities in the surface sensible heat flux that are important for the development of thermally driven flows over complex terrain.
Abstract
Topographic effects on radiation, including both topographic shading and slope effects, are included in the Weather Research and Forecasting (WRF) Model, and here they are made compatible with the immersed boundary method (IBM). IBM is an alternative method for representing complex terrain that reduces numerical errors over sloped terrain, thus extending the range of slopes that can be represented in WRF simulations. The implementation of topographic effects on radiation is validated by comparing land surface fluxes, as well as temperature and velocity fields, between idealized WRF simulations both with and without IBM. Following validation, the topographic shading implementation is tested in a semirealistic simulation of flow over Granite Mountain, Utah, where topographic shading is known to affect downslope flow development in the evening. The horizontal grid spacing is 50 m and the vertical grid spacing is approximately 8–27 m near the surface. Such a case would fail to run in WRF with its native terrain-following coordinates because of large local slope values reaching up to 55°. Good agreement is found between modeled surface energy budget components and observations from the Mountain Terrain Atmospheric Modeling and Observations (MATERHORN) program at a location on the east slope of Granite Mountain. In addition, the model captures large spatiotemporal inhomogeneities in the surface sensible heat flux that are important for the development of thermally driven flows over complex terrain.
Abstract
In numerical weather prediction and in reanalysis, robust approaches for observation bias correction are necessary to approach optimal data assimilation. The success of bias correction can be limited by model errors. Here, simultaneous estimation of observation and model biases, and the model state for an analysis, is explored with ensemble data assimilation and a simple model. The approach is based on parameter estimation using an augmented state in an ensemble adjustment Kalman filter. The observation biases are modeled with a linear term added to the forward operator. A bias is introduced in the forcing term of the model, leading to a model with complex errors that can be used in imperfect-model assimilation experiments.
Under a range of model forcing biases and observation biases, accurate observation bias estimation and correction are possible when the model forcing bias is simultaneously estimated and corrected. In the presence of both model error and observation biases, estimating one and ignoring the other harms the assimilation more than not estimating any errors at all, because the biases are not correctly attributed. Neglecting a large model forcing bias while estimating observation biases results in filter divergence; the observation bias parameter absorbs the model forcing bias, and recursively and incorrectly increases the increments. Neglecting observation bias results in suboptimal assimilation, but the model forcing bias parameter estimate remains stable because the model dynamics ensure covariance between the parameter and the model state.
Abstract
In numerical weather prediction and in reanalysis, robust approaches for observation bias correction are necessary to approach optimal data assimilation. The success of bias correction can be limited by model errors. Here, simultaneous estimation of observation and model biases, and the model state for an analysis, is explored with ensemble data assimilation and a simple model. The approach is based on parameter estimation using an augmented state in an ensemble adjustment Kalman filter. The observation biases are modeled with a linear term added to the forward operator. A bias is introduced in the forcing term of the model, leading to a model with complex errors that can be used in imperfect-model assimilation experiments.
Under a range of model forcing biases and observation biases, accurate observation bias estimation and correction are possible when the model forcing bias is simultaneously estimated and corrected. In the presence of both model error and observation biases, estimating one and ignoring the other harms the assimilation more than not estimating any errors at all, because the biases are not correctly attributed. Neglecting a large model forcing bias while estimating observation biases results in filter divergence; the observation bias parameter absorbs the model forcing bias, and recursively and incorrectly increases the increments. Neglecting observation bias results in suboptimal assimilation, but the model forcing bias parameter estimate remains stable because the model dynamics ensure covariance between the parameter and the model state.
Abstract
Large temperature fluctuations (LTFs), defined as a drop of the near-surface temperature of at least 3°C in less than 30 min followed by a recovery of at least half of the initial drop, were frequently observed during the Mountain Terrain Atmospheric Modeling and Observations (MATERHORN) program. Temperature time series at over 100 surface stations were examined in an automated fashion to identify and characterize LTFs. LTFs occur almost exclusively at night and at locations elevated 50–100 m above the basin floors, such as the east slope of the isolated Granite Mountain (GM). Temperature drops associated with LTFs were as large as 13°C and were typically greatest at heights of 4–10 m AGL. Observations and numerical simulations suggest that LTFs are the result of complex flow interactions of stably stratified flow with a mountain barrier and a leeside cold-air pool (CAP). An orographic wake forms over GM when stably stratified southwesterly nocturnal flow impinges on GM and is blocked at low levels. Warm crest-level air descends in the lee of the barrier, and the generation of baroclinic vorticity leads to periodic development of a vertically oriented vortex. Changes in the strength or location of the wake and vortex cause a displacement of the horizontal temperature gradient along the slope associated with the CAP edge, resulting in LTFs. This mechanism explains the low frequency of LTFs on the west slope of GM as well as the preference for LTFs to occur at higher elevations later at night, as the CAP depth increases.
Abstract
Large temperature fluctuations (LTFs), defined as a drop of the near-surface temperature of at least 3°C in less than 30 min followed by a recovery of at least half of the initial drop, were frequently observed during the Mountain Terrain Atmospheric Modeling and Observations (MATERHORN) program. Temperature time series at over 100 surface stations were examined in an automated fashion to identify and characterize LTFs. LTFs occur almost exclusively at night and at locations elevated 50–100 m above the basin floors, such as the east slope of the isolated Granite Mountain (GM). Temperature drops associated with LTFs were as large as 13°C and were typically greatest at heights of 4–10 m AGL. Observations and numerical simulations suggest that LTFs are the result of complex flow interactions of stably stratified flow with a mountain barrier and a leeside cold-air pool (CAP). An orographic wake forms over GM when stably stratified southwesterly nocturnal flow impinges on GM and is blocked at low levels. Warm crest-level air descends in the lee of the barrier, and the generation of baroclinic vorticity leads to periodic development of a vertically oriented vortex. Changes in the strength or location of the wake and vortex cause a displacement of the horizontal temperature gradient along the slope associated with the CAP edge, resulting in LTFs. This mechanism explains the low frequency of LTFs on the west slope of GM as well as the preference for LTFs to occur at higher elevations later at night, as the CAP depth increases.
Abstract
Weather Research and Forecasting (WRF) Model simulations of the autumn 2012 and spring 2013 Mountain Terrain Atmospheric Modeling and Observations Program (MATERHORN) field campaigns are validated against observations of components of the surface energy balance (SEB) collected over contrasting desert-shrub and playa land surfaces of the Great Salt Lake Desert in northwestern Utah. Over the desert shrub, a large underprediction of sensible heat flux and an overprediction of ground heat flux occurred during the autumn campaign when the model-analyzed soil moisture was considerably higher than the measured soil moisture. Simulations that incorporate in situ measurements of soil moisture into the land surface analyses and use a modified parameterization for soil thermal conductivity greatly reduce these errors over the desert shrub but exacerbate the overprediction of latent heat flux over the playa. The Noah land surface model coupled to WRF does not capture the many unusual playa land surface processes, and simulations that incorporate satellite-derived albedo and reduce the saturation vapor pressure over the playa only marginally improve the forecasts of the SEB components. Nevertheless, the forecast of the 2-m temperature difference between the playa and desert shrub improves, which increases the strength of the daytime off-playa breeze. The stronger off-playa breeze, however, does not substantially reduce the mean absolute errors in overall 10-m wind speed and direction. This work highlights some deficiencies of the Noah land surface model over two common arid land surfaces and demonstrates the importance of accurate land surface analyses over a dryland region.
Abstract
Weather Research and Forecasting (WRF) Model simulations of the autumn 2012 and spring 2013 Mountain Terrain Atmospheric Modeling and Observations Program (MATERHORN) field campaigns are validated against observations of components of the surface energy balance (SEB) collected over contrasting desert-shrub and playa land surfaces of the Great Salt Lake Desert in northwestern Utah. Over the desert shrub, a large underprediction of sensible heat flux and an overprediction of ground heat flux occurred during the autumn campaign when the model-analyzed soil moisture was considerably higher than the measured soil moisture. Simulations that incorporate in situ measurements of soil moisture into the land surface analyses and use a modified parameterization for soil thermal conductivity greatly reduce these errors over the desert shrub but exacerbate the overprediction of latent heat flux over the playa. The Noah land surface model coupled to WRF does not capture the many unusual playa land surface processes, and simulations that incorporate satellite-derived albedo and reduce the saturation vapor pressure over the playa only marginally improve the forecasts of the SEB components. Nevertheless, the forecast of the 2-m temperature difference between the playa and desert shrub improves, which increases the strength of the daytime off-playa breeze. The stronger off-playa breeze, however, does not substantially reduce the mean absolute errors in overall 10-m wind speed and direction. This work highlights some deficiencies of the Noah land surface model over two common arid land surfaces and demonstrates the importance of accurate land surface analyses over a dryland region.
Abstract
Operational Weather Research and Forecasting (WRF) Model forecasts run over Dugway Proving Ground (DPG) in northwest Utah, produced by the U.S. Army Test and Evaluation Command Four-Dimensional Weather System (4DWX), underpredict the amplitude of the diurnal temperature cycle during September and October. Mean afternoon [2000 UTC (1300 LST)] and early morning [1100 UTC (0400 LST)] 2-m temperature bias errors evaluated against 195 surface stations using 6- and 12-h forecasts are –1.37° and 1.66°C, respectively. Bias errors relative to soundings and 4DWX-DPG analyses illustrate that the afternoon cold bias extends from the surface to above the top of the planetary boundary layer, whereas the early morning warm bias develops in the lowest model levels and is confined to valleys and basins. These biases are largest during mostly clear conditions and are caused primarily by a regional overestimation of near-surface soil moisture in operational land surface analyses, which do not currently assimilate in situ soil moisture observations. Bias correction of these soil moisture analyses using data from 42 North American Soil Moisture Database stations throughout the Intermountain West reduces both the afternoon and early morning bias errors and improves forecasts of upper-level temperature and stability. These results illustrate that the assimilation of in situ and remotely sensed soil moisture observations, including those from the recently launched NASA Soil Moisture Active Passive (SMAP) mission, have the potential to greatly improve land surface analyses and near-surface temperature forecasts over arid regions.
Abstract
Operational Weather Research and Forecasting (WRF) Model forecasts run over Dugway Proving Ground (DPG) in northwest Utah, produced by the U.S. Army Test and Evaluation Command Four-Dimensional Weather System (4DWX), underpredict the amplitude of the diurnal temperature cycle during September and October. Mean afternoon [2000 UTC (1300 LST)] and early morning [1100 UTC (0400 LST)] 2-m temperature bias errors evaluated against 195 surface stations using 6- and 12-h forecasts are –1.37° and 1.66°C, respectively. Bias errors relative to soundings and 4DWX-DPG analyses illustrate that the afternoon cold bias extends from the surface to above the top of the planetary boundary layer, whereas the early morning warm bias develops in the lowest model levels and is confined to valleys and basins. These biases are largest during mostly clear conditions and are caused primarily by a regional overestimation of near-surface soil moisture in operational land surface analyses, which do not currently assimilate in situ soil moisture observations. Bias correction of these soil moisture analyses using data from 42 North American Soil Moisture Database stations throughout the Intermountain West reduces both the afternoon and early morning bias errors and improves forecasts of upper-level temperature and stability. These results illustrate that the assimilation of in situ and remotely sensed soil moisture observations, including those from the recently launched NASA Soil Moisture Active Passive (SMAP) mission, have the potential to greatly improve land surface analyses and near-surface temperature forecasts over arid regions.
Abstract
Expansion in the availability of relocatable near-surface atmospheric observing sensors introduces the question of where placement maximizes gain in forecast accuracy. As one possible method of addressing observation placement, the performance of ensemble sensitivity analysis (ESA) is examined for high-resolution (Δx = 4 km) predictions in complex terrain and during weak flow. ESA can be inaccurate when the underlying assumptions of linear dynamics (and Gaussian statistics) are violated, or when the sensitivity cannot be robustly sampled. A case study of a fog event at Salt Lake City International Airport (KSLC) in Utah provides a useful basis for examining these issues, with the additional influence of complex terrain. A realistic upper-air observing network is used in perfect-model ensemble data assimilation experiments, providing the statistics for ESA. Results show that water vapor mixing ratios over KSLC are sensitive to potential temperature on the first model layer tens of kilometers away, 6 h prior to verification and prior to the onset of fog. Potential temperatures indicate inversion strength in the Salt Lake basin; the ESA predicts southerly flow and strengthened inversions will increase water vapor over KSLC. Linearity tests show that the nonlinear response is about twice the expected response. Experiments with smaller ensembles show that qualitatively similar conclusions about the sensitivity pattern can be reached with ensembles as small as 48 members, but smaller ensembles do not produce accurate sensitivity estimates. Taken together, the results motivate a closer look at the fundamental characteristics of ESA when dynamics (and therefore correlations) are weak.
Abstract
Expansion in the availability of relocatable near-surface atmospheric observing sensors introduces the question of where placement maximizes gain in forecast accuracy. As one possible method of addressing observation placement, the performance of ensemble sensitivity analysis (ESA) is examined for high-resolution (Δx = 4 km) predictions in complex terrain and during weak flow. ESA can be inaccurate when the underlying assumptions of linear dynamics (and Gaussian statistics) are violated, or when the sensitivity cannot be robustly sampled. A case study of a fog event at Salt Lake City International Airport (KSLC) in Utah provides a useful basis for examining these issues, with the additional influence of complex terrain. A realistic upper-air observing network is used in perfect-model ensemble data assimilation experiments, providing the statistics for ESA. Results show that water vapor mixing ratios over KSLC are sensitive to potential temperature on the first model layer tens of kilometers away, 6 h prior to verification and prior to the onset of fog. Potential temperatures indicate inversion strength in the Salt Lake basin; the ESA predicts southerly flow and strengthened inversions will increase water vapor over KSLC. Linearity tests show that the nonlinear response is about twice the expected response. Experiments with smaller ensembles show that qualitatively similar conclusions about the sensitivity pattern can be reached with ensembles as small as 48 members, but smaller ensembles do not produce accurate sensitivity estimates. Taken together, the results motivate a closer look at the fundamental characteristics of ESA when dynamics (and therefore correlations) are weak.
