Abstract

Various configurations of the intermittent data assimilation procedure for Level-II Weather Surveillance Radar-1988 Doppler radar data are examined for the analysis and prediction of a tornadic thunderstorm that occurred on 8 May 2003 near Oklahoma City, Oklahoma. Several tornadoes were produced by this thunderstorm, causing extensive damages in the south Oklahoma City area. Within the rapidly cycled assimilation system, the Advanced Regional Prediction System three-dimensional variational data assimilation (ARPS 3DVAR) is employed to analyze conventional and radar radial velocity data, while the ARPS complex cloud analysis procedure is used to analyze cloud and hydrometeor fields and adjust in-cloud temperature and moisture fields based on reflectivity observations and the preliminary analysis of the atmosphere. Forecasts for up to 2.5 h are made from the assimilated initial conditions. Two one-way nested grids at 9- and 3-km grid spacings are employed although the assimilation configuration experiments are conducted for the 3-km grid only while keeping the 9-km grid configuration the same. Data from the Oklahoma City radar are used. Different combinations of the assimilation frequency, in-cloud temperature adjustment schemes, and the length and coverage of the assimilation window are tested, and the results are discussed with respect to the length and evolution stage of the thunderstorm life cycle. It is found that even though the general assimilation method remains the same, the assimilation settings can significantly impact the results of assimilation and the subsequent forecast. For this case, a 1-h-long assimilation window covering the entire initial stage of the storm together with a 10-min spinup period before storm initiation works best. Assimilation frequency and in-cloud temperature adjustment scheme should be set carefully to add suitable amounts of potential energy during assimilation. High assimilation frequency does not necessarily lead to a better result because of the significant adjustment during the initial forecast period. When a short assimilation window is used, covering the later part of the initial stage of storm and using a high assimilation frequency and a temperature adjustment scheme based on latent heat release can quickly build up the storm and produce a reasonable analysis and forecast. The results also show that when the data from a single Doppler radar are assimilated with properly chosen assimilation configurations, the model is able to predict the evolution of the 8 May 2003 Oklahoma City tornadic thunderstorm well for up to 2.5 h. The implications of the choices of assimilation settings for real-time applications are discussed.

Abstract

When assessed using the difference between urban and rural air temperatures, the urban heat island (UHI) is most prominent during the nighttime. Typically, nocturnal UHI intensity is maintained throughout the night. The UHI intensity over Dallas–Fort Worth (DFW), Texas, however, experienced frequent “collapses” (sudden decreases) around midnight during August 2011, while the region was experiencing an intense heat wave. Observational and modeling studies were conducted to understand this unique phenomenon. Sea-breeze passage was found to be ultimately responsible for the collapses of the nocturnal UHI. Sea-breeze circulation developed along the coast of the Gulf of Mexico during the daytime. During the nighttime, the sea-breeze circulation was advected inland (as far as ~400 km) by the low-level jet-enhanced southerly flow, maintaining the characteristics of sea-breeze fronts, including the enhanced wind shear and vertical mixing. Ahead of the front, surface radiative cooling enhanced the near-surface temperature inversion in rural areas through the night with calm winds. During the frontal passage (around midnight at DFW), the enhanced vertical mixing at the leading edge of the fronts brought warmer air to the surface, leading to rural surface warming events. In contrast, urban effects led to a nearly neutral urban boundary layer. The enhanced mechanical mixing associated with sea-breeze fronts, therefore, did not increase urban surface temperature. The different responses to the sea-breeze frontal passages between rural (warming) and urban areas (no warming) led to the collapse of the UHI. The inland penetration of sea-breeze fronts at such large distances from the coast and their effects on UHI have not been documented in the literature.

Abstract

In this two-part paper, the impact of level-II Weather Surveillance Radar-1988 Doppler (WSR-88D) reflectivity and radial velocity data on the prediction of a cluster of tornadic thunderstorms in the Advanced Regional Prediction System (ARPS) model are studied. Radar reflectivity data are used primarily in a cloud analysis procedure that retrieves the amount of hydrometeors and adjusts in-cloud temperature, moisture, and cloud fields, while radial velocity data are analyzed through a three-dimensional variational (3DVAR) scheme that contains a mass divergence constraint in the cost function. In Part I, the impact of the cloud analysis and modifications to the scheme are examined while Part II focuses on the impact of radial velocity and the mass divergence constraint.

The case studied is that of the 28 March 2000 Fort Worth, Texas, tornado outbreaks. The same case was studied by Xue et al. using the ARPS Data Analysis System (ADAS) and an earlier version of the cloud analysis procedure with WSR-88D level-III data. Since then, several modifications to the cloud analysis procedure, including those to the in-cloud temperature adjustment and the analysis of precipitation species, have been made. They are described in detail with examples.

The assimilation and predictions use a 3-km grid nested inside a 9-km one. The level-II reflectivity data are assimilated, through the cloud analysis, at 10-min intervals in a 1-h period that ends a little over 1 h preceding the first tornado outbreak. Experiments with different settings within the cloud analysis procedure are examined. It is found that the experiment using the improved cloud analysis procedure with reflectivity data can capture the important characteristics of the main tornadic thunderstorm more accurately than the experiment using the early version of cloud analysis. The contributions of different modifications to the above improvements are investigated.

Abstract

In this two-part paper, the impact of level-II Weather Surveillance Radar-1988 Doppler (WSR-88D) radar reflectivity and radial velocity data on the prediction of a cluster of tornadic thunderstorms in the Advanced Regional Prediction System (ARPS) model is studied. Radar reflectivity data are used primarily in a cloud analysis procedure that retrieves the amount of hydrometeors and adjusts in-cloud temperature, moisture, and cloud fields, while radial velocity data are analyzed through a three-dimensional variational (3DVAR) data assimilation scheme that contains a 3D mass divergence constraint in the cost function. In Part I, the impact of the cloud analysis and modifications to the scheme are discussed. In this part, the impact of radial velocity data and the mass divergence constraint in the 3DVAR cost function are studied.

The case studied is that of the 28 March 2000 Fort Worth tornadoes. The addition of the radial velocity improves the forecasts beyond that experienced with the cloud analysis alone. The prediction is able to forecast the morphology of individual storm cells on the 3-km grid up to 2 h; the rotating supercell characteristics of the storm that spawned two tornadoes are well captured; timing errors in the forecast are less than 15 min and location errors are less than 10 km at the time of the tornadoes.

When forecasts were made with radial velocity assimilation but not reflectivity, they failed to predict nearly all storm cells. Using the current 3DVAR and cloud analysis procedure with 10-min intermittent assimilation cycles, reflectivity data are found to have a greater positive impact than radial velocity. The use of radial velocity does improve the storm forecast when combined with reflectivity assimilation, by, for example, improving the forecasting of the strong low-level vorticity centers associated with the tornadoes. Positive effects of including a mass divergence constraint in the 3DVAR cost function are also documented.

Abstract

Since the 1950s, a countergradient flux term has been added to some K-profile-based first-order PBL schemes, allowing them to simulate the slightly statically stable upper part of the convective boundary layer (CBL) observed in a limited number of aircraft soundings. There is, however, substantial uncertainty in inferring detailed CBL structure, particularly the level of neutral stability (z _n ), from such a limited number of soundings. In this study, composite profiles of potential temperature are derived from multiyear early afternoon radiosonde data over Beijing, China. The CBLs become slightly stable above z _n ~ 0.31–0.33z _i , where z _i is the CBL depth. These composite profiles are used to evaluate two K-profile PBL schemes, the Yonsei University (YSU) and Shin–Hong (SH) schemes, and to optimize the latter through parameter calibration. In one-dimensional simulations using the WRF Model, YSU simulates a stable CBL above z _n ~ 0.24z _i , while default SH simulates a thick superadiabatic lower CBL with z _n ~ 0.45z _i . Experiments with the analytic solution of a K-profile PBL model show that adjusting the countergradient flux profile leads to significant changes in the thermal structure of CBL, informing the calibration of SH. The SH scheme replaces the countergradient heat flux term in its predecessor YSU scheme with a three-layer nonlocal heating profile, with f _nl specifying the peak value and $z_{*\mathrm{SL}}$ specifying the height of this peak value. Increasing f _nl to 1.1 lowers z _n , but to too low a value, while simultaneously increasing $z_{*\mathrm{SL}}$ to 0.4 leads to a more appropriate z _n ~ 0.36z _i . The calibrated SH scheme performs better than YSU and default SH for real CBLs.

Abstract

A 50-m-grid-spacing Advanced Regional Prediction System (ARPS) simulation of the 8 May 2003 Oklahoma City tornadic supercell is examined. A 40-min forecast run on the 50-m grid produces two F3-intensity tornadoes that track within 10 km of the location of the observed long-track F4-intensity tornado.

The development of both simulated tornadoes is analyzed to determine the processes responsible for tornadogenesis. Trajectory-based analyses of vorticity components and their time evolution reveal that tilting of low-level frictionally generated horizontal vorticity plays a dominant role in the development of vertical vorticity near the ground. This result represents the first time that such a mechanism has been shown to be important for generating near-surface vertical vorticity leading to tornadogenesis.

A sensitivity simulation run with surface drag turned off was found to be considerably different from the simulation with drag included. A tornado still developed in the no-drag simulation, but it was much shorter lived and took a substantially different track than the observed tornadoes as well as the simulated tornadoes in the drag simulation. Tilting of baroclinic vorticity in an outflow surge may have played a role in tornadogenesis in the no-drag simulation.

Abstract

The 8 May 2003 Oklahoma City, Oklahoma, tornadic supercell is predicted with the Advanced Regional Prediction System (ARPS) model using four nested grids with 9-km, 1-km, 100-m, and 50-m grid spacings. The Oklahoma City Weather Surveillance Radar-1988 Doppler (WSR-88D) radial velocity and reflectivity data are assimilated through the ARPS three-dimensional variational data assimilation (3DVAR) and cloud analysis on the 1-km grid to generate a set of initial conditions that includes a well-analyzed supercell and associated low-level mesocyclone. Additional 1-km experiments show that the use of radial velocity and the proper use of a divergence constraint in the 3DVAR play an important role in the establishment of the low-level mesocyclone during the assimilation and forecast. Assimilating reflectivity data alone failed to predict the mesocyclone intensification. The 100-m grid starts from the interpolated 1-km control initial conditions, while the further nested 50-m grid starts from the 20-min forecast on the 100-m grid. The forecasts on both grids cover the entire period of the observed tornado outbreak and successfully capture the development of tornadic vortices. The intensity of a tornado on the 50-m grid reaches the high end of category 3 on the Fujita scale (F3), while the corresponding simulated tornado on the 100-m grid reaches F2 intensity. The timing of the tornadogenesis on both grids agrees with the observations very well, although the predicted tornado was slightly weaker and somewhat shorter lived. The predicted tornado track parallels the observed damage track although it is displaced northward by about 8 km. The predicted tornado vortices have realistic structures similar to those documented in previous theoretical, idealized modeling and some observational studies. The prediction of an observed tornado in a supercell with a similar degree of realism has not been achieved before.

Abstract

This study evaluated the Weather Research and Forecasting (WRF) Model sensitivity to different planetary boundary layer (PBL) schemes (the YSU and MYJ schemes) and urban schemes including the bulk scheme (BULK), single-layer urban canopy model (UCM), multilayer building environment parameterization (BEP) model, and multilayer building energy model (BEM). Daily reinitialization simulations were conducted over Dallas–Fort Worth during a dry summer month (July 2011) and a wet summer month (July 2015) with weaker (stronger) daytime (nocturnal) UHI in 2011 than 2015. All urban schemes overestimated the urban daytime 2-m temperature in both summers, but BEP and BEM still reproduced the daytime urban cool island in the dry summer. All urban schemes reproduced the nocturnal urban heat island, with BEP producing the weakest one due to its unrealistic urban cooling. BULK and UCM overestimated the urban canopy wind speed, while BEP and BEM underestimated it. The urban schemes showed prominent impact on daytime PBL profiles. UCM + MYJ showed a superior performance than other configurations. The relatively large (small) aspect ratio between building height and road width in UCM (BEM) was responsible for the overprediction (underprediction) of urban canopy temperature. The relatively low (high) building height in UCM (BEM) was responsible for the overprediction (underprediction) of urban canopy wind speed. Improving urban schemes and providing realistic urban parameters were critical for improving urban canopy simulation.

Abstract

The magnitude of water vapor content within the near-storm inflow can either support or deter the storm’s upscale growth and maintenance. However, the heterogeneity of the moisture field near storms remains poorly understood because the operational observation network lacks detail. This observational study illustrates that near-storm inflow water vapor environments are both significantly heterogeneous and different than the far-inflow storm environment. This study also depicts the importance of temporal variation of water vapor mixing ratio (WVMR) to instability during the peak tornadic seasons in the U.S. Southeast and Great Plains regions during the Verification of the Origins of Rotation in Tornadoes Experiment Southeast 2018 (VSE18) campaign and the Targeted Observation by Radar and UAS of Supercells (TORUS) campaign, respectively. VSE18 results suggest that the surface processes control WVMR variation significantly in lower levels, with the highest WVMR mainly located near the surface in inflows in the southeast region. In contrast, TORUS results show more vertically homogeneous WVMR profiles and rather uniform water vapor distribution variation occurring in deep, moist stratified inflows in the Great Plains region. Temporal water vapor variations within 5-min periods could lead to over 1000 J kg⁻¹ CAPE changes in both VSE18 and TORUS, which represent significant potential buoyancy perturbations for storms to intensify or decay. These temporal water vapor and instability evolutions of moving storms remain difficult to capture via radiosondes and fixed in situ or profiling instrumentation, yet may exert a strong impact on storm evolution. This study suggests that improving observations of the variability of near-storm inflow moisture can accurately refine a potential severe weather threat.

Abstract

Many studies have investigated urban heat island (UHI) intensity for cities around the world, which is normally quantified as the temperature difference between urban location(s) and rural location(s). A few open questions still remain regarding the UHI, such as the spatial distribution of UHI intensity, temporal (including diurnal and seasonal) variation of UHI intensity, and the UHI formation mechanism. A dense network of atmospheric monitoring sites, known as the Oklahoma City (OKC) Micronet (OKCNET), was deployed in 2008 across the OKC metropolitan area. This study analyzes data from OKCNET in 2009 and 2010 to investigate OKC UHI at a subcity spatial scale for the first time. The UHI intensity exhibited large spatial variations over OKC. During both daytime and nighttime, the strongest UHI intensity is mostly confined around the central business district where land surface roughness is the highest in the OKC metropolitan area. These results do not support the roughness warming theory to explain the air temperature UHI in OKC. The UHI intensity of OKC increased prominently around the early evening transition (EET) and stayed at a fairly constant level throughout the night. The physical processes during the EET play a critical role in determining the nocturnal UHI intensity. The near-surface rural temperature inversion strength was a good indicator for nocturnal UHI intensity. As a consequence of the relatively weak near-surface rural inversion, the strongest nocturnal UHI in OKC was less likely to occur in summer. Other meteorological factors (e.g., wind speed and cloud) can affect the stability/depth of the nighttime boundary layer and can thus modulate nocturnal UHI intensity.