Simulation of mesoscale structure of thunderstorm using ARPS Model

In this paper mesoscale structures of two thunderstorm events over Delhi have been simulated using ARPS (Advanced Regional Prediction System) model. Numerical experiments were carried out using radiosonde data of Delhi and applying a potential temperature perturbation for triggering convective activity. The simulation exercise demonstrates strong updrafts and downdrafts associated with the thunderstorm cells, indicating the presence of very strong localized convection. The development and evolution of thunderstorm and propagation of associated precipitation zone are clearly brought out in this simulation study.

Precipitation simulation of synoptic scale systems over western Himalayan region using Advanced Regional Prediction System (ARPS) model

Operational weather prediction over western Himalayan region is a challenging job due to scarcity of data and complex topography that interacts with approaching weather system. Accurate prediction of complex weather phenomena requires dense data network which is difficult to establish in mountain due to complex terrain and hostile weather conditions over Himalaya. The alternate method to overcome this problem is by ingesting three-dimensional meteorological variables from global model’s analysis and forecast values as initial and lateral boundary conditions in meso-scale numerical models. Simultaneously, data assimilation is a potential tool in which non-conventional [satellite, radar and Automatic Weather Station (AWS)] and conventional (surface and upper air observations) data are ingested in the numerical models to generate high resolution and accurate initial fields for the initialization of the mesoscale model. In the present study, Advanced Regional Prediction System (ARPS) model has been used for the prediction of synoptic weather system known as Western Disturbance (WD) that affects the weather of western and central Himalaya during winter period (November – April).The ARPS model has been selected for this study because the model has its own objective analysis and quality control system. It has the capacity to ingest the satellite, Doppler weather radar data and other types of observations. Its assimilation system can also be used to overcome the problem of data scarcity in Himalayan region. In this study, initial and lateral boundary fields are taken from the T-80 spectral global model operationally used at National Centre for Medium Range Prediction (NCMRWF), Noida (UP), India. The global model’s analysis was taken as the initial condition and 24 hour’s interval forecasts as lateral boundary conditions. The model has been used for the simulation of few WDs for 96 hours (Four days). The comparison of ARPS simulation with T-80 forecast shows that the ARPS model could produce better results in respect of the circulation of WDs and hence it can be utilized for the operational weather prediction over the Indian region.

Explicit Prediction of Hail in a Long-Lasting Multicellular Convective System in Eastern China Using Multimoment Microphysics Schemes

During the afternoon of 28 April 2015, a multicellular convective system swept southward through much of Jiangsu Province, China, over about 7 h, producing egg-sized hailstones on the ground. The hailstorm event is simulated using the Advanced Regional Prediction System (ARPS) at 1-km grid spacing. Different configurations of the Milbrandt–Yau microphysics scheme are used, predicting one, two, and three moments of the hydrometeor particle size distributions (PSDs). Simulated reflectivity and maximum estimated size of hail (MESH) derived from the simulations are verified against reflectivity observed by operational S-band Doppler radars and radar-derived MESH, respectively. Comparisons suggest that the general evolution of the hailstorm is better predicted by the three-moment scheme, and neighborhood-based MESH evaluation further confirms the advantage of the three-moment scheme in hail size prediction. Surface accumulated hail mass, number, and hail distribution characteristics within simulated storms are examined across sensitivity experiments. Results suggest that multimoment schemes produce more realistic hail distribution characteristics, with the three-moment scheme performing the best. Size sorting is found to play a significant role in determining hail distribution within the storms. Detailed microphysical budget analyses are conducted for each experiment, and results indicate that the differences in hail growth processes among the experiments can be mainly ascribed to the different treatments of the shape parameter within different microphysics schemes. Both the differences in size sorting and hail growth processes contribute to the simulated hail distribution differences within storms and at the surface.

The Impact of a Stochastically Perturbing Microphysics Scheme on an Idealized Supercell Storm

The concept of stochastic parameterization provides an opportunity to represent spatiotemporal errors caused by microphysics schemes that play important roles in supercell simulations. In this study, two stochastic methods, the stochastically perturbed temperature tendency from microphysics (SPTTM) method and the stochastically perturbed intercept parameters of microphysics (SPIPM) method, are implemented within the Lin scheme, which is based on the Advanced Regional Prediction System (ARPS) model, and are tested using an idealized supercell case. The SPTTM and SPIPM methods perturb the temperature tendency and the intercept parameters (IPs), respectively. Both methods use recursive filters to generate horizontally smooth perturbations and adopt the barotropic structure for the perturbation r, which is multiplied by tendencies or parameters from this parameterization. A double-moment microphysics scheme is used for the truth run. Compared to the multiparameter method, which uses randomly perturbed prescribed parameters, stochastic methods often produce larger ensemble spreads and better forecast the intensity of updraft helicity (UH). The SPTTM method better predicts the intensity by intensifying the midlevel heating with its positive perturbation r, whereas it performs worse in the presence of negative perturbation. In contrast, the SPIPM method can increase the intensity of UH by either positive or negative perturbation, which increases the likelihood for members to predict strong UH.

A study of the influence of forest gaps on fire–atmosphere interactions

Much uncertainty exists regarding the possible role that gaps in forest canopies play in modulating fire–atmosphere interactions in otherwise horizontally homogeneous forests. This study examines the influence of gaps in forest canopies on atmospheric perturbations induced by a low-intensity fire using the ARPS-CANOPY model, a version of the Advanced Regional Prediction System (ARPS) model with a canopy parameterization. A series of numerical experiments are conducted with a stationary low-intensity fire, represented in the model as a line of enhanced surface sensible heat flux. Experiments are conducted with and without forest gaps, and with gaps in different positions relative to the fire line. For each of the four cases considered, an additional simulation is performed without the fire to facilitate comparison of the fire-perturbed atmosphere and the background state. Analyses of both mean and instantaneous wind velocity, turbulent kinetic energy, air temperature, and turbulent mixing of heat are presented in order to examine the fire-perturbed atmosphere on multiple timescales. Results of the analyses indicate that the impact of the fire on the atmosphere is greatest in the case with the gap centered on the fire and weakest in the case with the gap upstream of the fire. It is shown that gaps in forest canopies have the potential to play a role in the vertical as well as horizontal transport of heat away from the fire. Results also suggest that, in order to understand how the fire will alter wind and turbulence in a heterogeneous forest, one needs to first understand how the forest heterogeneity itself influences the wind and turbulence fields without the fire.

A study of the influence of forest gaps on fire-atmosphere interactions

Much uncertainty exists regarding the possible role that gaps in forest canopies play in modulating fire-atmosphere interactions in otherwise horizontally homogeneous forests. This study examines the impact of forest gaps on fire-atmosphere interactions using the ARPS-CANOPY model, a version of the Advanced Regional Prediction System (ARPS) model with a canopy parameterization. A series of numerical experiments are conducted with a stationary low-intensity fire, represented in the model as a line of enhanced surface sensible heat flux. Experiments are conducted with and without forest gaps, and with gaps in different positions relative to the fireline. For each of the four cases considered, an additional simulation is performed without the fire to facilitate comparison of the fire-perturbed atmosphere and the background state. Analyses of both mean and instantaneous wind velocity, turbulent kinetic energy, air temperature, and turbulent mixing of heat are presented in order to examine the fire-perturbed atmosphere on multiple time scales. Results of the analyses indicate that the impact of the fire on the atmosphere is greatest in the case with the gap centered on the fire, and weakest in the case with the gap upstream of the fire. It is shown that gaps in forest canopies have the potential to play a substantial role in the vertical as well as horizontal transport of heat away from the fire. Results also suggest that in order to understand how the fire will alter wind and turbulence in a heterogeneous forest, one needs to first understand how the forest heterogeneity itself influences the wind and turbulence fields without the fire.

Mean and Turbulent Flow Downstream of a Low-Intensity Fire: Influence of Canopy and Background Atmospheric Conditions

This study examines the sensitivity of mean and turbulent flow in the planetary boundary layer and roughness sublayer to a low-intensity fire and evaluates whether the sensitivity is dependent on canopy and background atmospheric properties. The ARPS-CANOPY model, a modified version of the Advanced Regional Prediction System (ARPS) model with a canopy parameterization, is utilized for this purpose. A series of numerical experiments are conducted to evaluate whether the ability of the fire to alter downstream wind, temperature, turbulent kinetic energy (TKE), and vertical heat flux differs between forested and open areas, sparse and dense forests, weak and strong background flow, and neutral and convective background stability. Analysis of all experiments shows that, in general, mean and turbulent flow both prior to and during a low-intensity fire is damped in the presence of a canopy. Greater sensitivity to the fire is found in cases with strong ambient wind speed than in cases with quiescent or weak wind speed. Furthermore, sensitivity of downstream atmospheric conditions to the fire is shown to be strongest with a neutrally stratified background. An analysis of the TKE budget reveals that both buoyancy and wind shear contribute to TKE production during the period of time in which the fire conditions are applied to the model. On the basis of the results of the ARPS simulations, caution is advised when applying ARPS-simulation results to predictions of smoke transport and dispersion: smoke-model users should consider whether canopy impacts on the atmosphere are accounted for and whether neutral stratification is assumed.

Some Observing System Simulation Experiments with a Hybrid 3DEnVAR System for Storm-Scale Radar Data Assimilation

A hybrid three-dimensional ensemble–variational data assimilation (3DEnVAR) algorithm is developed based on the 3D variational data assimilation (3DVAR) and ensemble Kalman filter (EnKF) programs with the Advanced Regional Prediction System (ARPS). The method uses the extended control variable approach to combine the static and ensemble-derived flow-dependent forecast error covariances. The method is applied to the assimilation of simulated data from two radars for a supercell storm. Some sensitivity experiments are performed to answer questions about how flow-dependent covariance estimated from the forecast ensemble can be best used in the hybrid 3DEnVAR scheme. When the ensemble size is relatively small (with 5 or 10 ensemble members), it is found that experiments with a weaker weighting value for the ensemble covariance leads to better analysis results. Even when severe sampling errors exist, introducing ensemble-estimated covariances into the variational method still benefits the analysis. For reasonably large ensemble sizes (50–100 members), a stronger relative weighting (>0.8) for the ensemble covariance leads to better analyses from the hybrid 3DEnVAR. In addition, the sensitivity experiments also indicate that the best results are obtained when the number of the augmented control variables is a function of three spatial dimensions and ensemble members, and is the same for all analysis variables.

Assimilation of T-TREC-Retrieved Winds from Single-Doppler Radar with an Ensemble Kalman Filter for the Forecast of Typhoon Jangmi (2008)

A tropical cyclone (TC) circulation Tracking Radar Echo by Correlation technique (T-TREC) developed recently is applied to derive horizontal winds from single Doppler radar reflectivity Z data (combined with radial velocity Vr data when available). The typically much longer maximum range of Z observations compared to Vr data allows for much larger spatial coverage of the T-TREC-retrieved winds (VTREC) when a TC first enters the maximum range of a coastal radar. Retrieved using data from more than one scan volume, the T-TREC winds also contain valuable cross-beam wind information. The VTREC or Vr data at 30-min intervals are assimilated into the Advanced Regional Prediction System (ARPS) model at 3-km grid spacing using an ensemble Kalman filter, over a 2-h window shortly after Typhoon Jangmi (2008) entered the Vr coverage area of an operational weather radar of Taiwan. The assimilation of VTREC data produces analyses of the typhoon structure and intensity that more closely match observations than analyses produced using Vr data or the reference Global Forecast System (GFS) analysis. Subsequent 28-h forecasts of intensity, track, structure, and precipitation are also improved by assimilating VTREC data. Further sensitivity experiments show that assimilation of VTREC data can build up a reasonably strong vortex in 1 h, while a longer assimilation period is required to spin up the vortex when assimilating Vr. Although the difference between assimilating VTREC and Vr is smaller when the assimilation window is longer, the improvement from assimilating VTREC is still evident. Assimilating Z data in addition to Vr or VTREC results in little further improvement.

Numerical Prediction of the 8 May 2003 Oklahoma City Tornadic Supercell and Embedded Tornado Using ARPS with the Assimilation of WSR-88D Data

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.