Abstract

Model performance and sensitivity to model physics options are studied with the Weather Research and Forecasting model (version 3.1.1) over Delhi region in India for surface and upper air meteorological parameters in summer and winter seasons. A case study with the model has been performed with different configurations, and the best physics options suited for this region have been, determined. Comparison between estimated and observed data was carried out through standard statistical measures. Generally, the combination of Pleim-Xiu land surface model, Pleim surface layer scheme, and Asymmetric Convective Model has been found to produce better estimates of temperature and relative humidity for Delhi region. Wind speed and direction estimations were observed best for MM5 similarity surface layer along with Yonsei University boundary layer scheme. Nested domains with higher resolutions were not helpful in improving the simulation results as per the current availability of the data. Overall, the present case study shows that the model has performed reasonably well over the subtropical region of Delhi.

1. Introduction

Numerical weather prediction models have different sets of physical parameterization configuration options. As the model becomes more sophisticated, greater number of physical processes can be incorporated into it, and there is a range of physical schemes which can be used to simulate them. Thus, one of the essential steps in numerical weather simulation is to choose the best set of physics options for the region and time period under consideration. Numerous sensitivity studies for one or more physics options of WRF model have been undertaken in different parts of the world such as Spain [1], Japan and Korea [2], Alaska [3], South America [4], western United States [5], southern United States [6], West Africa [7], and others. In India, some sensitivity studies for WRF have been undertaken which mainly focus on extreme events like thunderstorm [8], tropical cyclone [9], and heavy precipitation [1012]. The present work is a case study in which an attempt is made to apply several physics options of the Weather Research and Forecasting model (WRF v 3.1.1) to a subtropical region, namely, Delhi in India, for examining model sensitivity and evaluating the model’s performance as suited to this region.

2. Model and Configurations

2.1. WRF Modeling System

The Weather Research and Forecasting model is developed for mesoscale modeling and is a supported “community model”, that is, a free and shared resource with distributed development and centralized support. Its development is led by NCAR, NOAA/ESRL, and NOAA/NCEP/EMC with partnerships and collaborations with universities and other government agencies in the US and overseas. The Advanced Research WRF (ARW) dynamical core has an equation set which is fully compressible, Eulerian and nonhydrostatic with a run-time hydrostatic option [13]. It is conservative for scalar variables. The model uses terrain-following, hydrostatic-pressure vertical coordinate with the top of the model being a constant pressure surface. The horizontal grid is the Arakawa-C grid. The time integration scheme in the model uses the third-order Runge-Kutta scheme, and the spatial discretization employs 2nd- to 6th-order schemes. The model supports both idealized and real-data applications with various lateral boundary condition options. The model also supports one-way, two-way, and moving nest options [14].

2.2. Model Configurations

There are various options for physical parameterisations of key boundary layer phenomena in this model for (i) microphysics, (ii) cumulus parameterization (CP), (iii) surface layer (SL), (iv) land-surface model (LSM), and (v) planetary boundary layer (PBL). The present study captures all of the above physical parameterisations in various options representing all of these five categories. The model experiments based on these are detailed in Table 1. The following subsections describe the physical parameterisations in each of these categories [15] and their consideration for model simulations in this study.

2.2.1. Microphysics

Microphysics includes explicitly resolved water vapor, cloud, and precipitation processes. Out of the several microphysics schemes of the model, the options suitable for high resolution-domain are chosen as the innermost domain is of a resolution of 2 km. These include Purdue Lin scheme [16] and WRF Single-Moment 6-class scheme (WSM6) [17]. However, a relatively simpler Eta microphysics option [18] has also been chosen for comparison.

2.2.2. Cumulus Parameterization

Cumulus parameterization schemes are responsible for the subgrid-scale effects of convective and/or shallow clouds but are valid only for coarse-grid sizes such as greater than 10 km. Grell 3d ensemble cumulus scheme is most suitable for the present analysis as it is applicable for higher-resolution domains allowing for subsidence in neighboring columns. Apart from this, Kain-Fritsch scheme [19] is also considered for model experiments since this is the most widely exercised option for Cumulus parameterization.

2.2.3. Surface Layer and Planetary Boundary Layer

The surface layer schemes calculate friction velocities and exchange coefficients that enable the calculation of surface heat and moisture fluxes by the land-surface models and surface stress in the planetary boundary layer (PBL) scheme. The PBL schemes determine the flux profiles within the well-mixed boundary layer and the stable layer, and thus provide atmospheric tendencies of temperature, moisture (including clouds), and horizontal momentum in the entire atmospheric column. In some cases, a surface layer option is tied to a particular boundary layer option. In the present study, all surface layer schemes applicable in WRF, that is, MM5 similarity [20], Eta similarity [21], Pleim-Xiu surface layer [22], and QNSE surface layer schemes have been used. The default boundary layer scheme was Mellor-Yamada-Janjic scheme [21] while the other options that were tried were Yonsei University scheme [23] (tied with MM5 similarity for surface layer), asymmetric convective model (ACM2) scheme (tied with Pleim-Xiu surface layer) [22], and quasinormal scale elimination (QNSE) PBL. Other PBL schemes such as BouLaC PBL and LES PBL are applicable to their relevant special cases and thus have not been tested in the present study.

2.2.4. Land Surface Model (LSM)

Land surface models (LSMs) play an important role in building soil temperature and moisture profiles and canopy properties. Consequently, the land surface schemes work as a fuel for atmospheric processes or in other words provide the forcing from the ground below to the lower part of the boundary layer which subsequently gets transported to the rest of the atmosphere. LSMs use atmospheric information from the surface layer scheme, radiative forcing from the radiation scheme, and precipitation forcing from the microphysics and convective schemes, together with internal information on the land’s state variables and land-surface properties, to provide heat and moisture fluxes over land points and sea-ice points.

These fluxes provide a lower boundary condition for the vertical transport to the PBL. Out of all the LSMs applicable in WRF, Noah Land Surface Model (default), RUC Land Surface Model, and Pleim-Xiu Land Surface Model have been applied in the simulations. The simple 5-layer thermal diffusion scheme has not been tested as it includes only soil temperature. Also, since the area under analysis is a subtropical region, fractional sea-ice scheme (which requires snow cover data) has not been considered.

In the present study, a total of nine simulations for each season were carried out. A default set of schemes was fixed for one simulation which was named as Case  1. On the basis of the more widely used options, the default case was fixed as Purdue Lin Scheme [16] for Microphysics, Kain-Fritsch scheme [19] for Cumulus parameterization, and Eta Similarity theory [21] for Surface Layer, Noah Land Surface Model [24, 25], and Mellor-Yamada-Janjic Planetary Boundary Layer. Thereafter, different options were adopted for a given physics category keeping options in other physics categories as constant. Table 1 lists out the various options that were used in the simulations. In all, apart from the configuration in default case simulation (simulation no. 1), two simulations for microphysics variations (simulation no. 6, 8), three simulations for surface layer schemes (simulation no. 2, 4, 7), two simulations for LSMs variations (simulation no. 3, 7), three simulations for surface layer schemes (simulation no. 2, 4, 7), three simulations for PBL schemes (Simulation no. 4, 5, 7), and one simulation for cumulus parameterization (simulation no. 9) have been tested. Rapid radiative transfer model (RRTM) longwave radiation scheme [26] and Dudhia [27] shortwave radiation scheme have been applied to all simulations.

3. Model Domain and Simulation Period

Figure 1 shows the domains used for modeling in the present study. The parent domain (D1) covers the Indian subcontinent with a spatial resolution of 18 km. The domain is a peninsular region flanked by Himalayan mountains in the north. The first nested domain (D2) constitutes northern India region with resolution of 6 km. The second nested domain (D3), which is the area under analysis, covers the city of National Capital Region of Delhi and surroundings with a domain resolution of 2 km centered at 28.52°N and 77.12°E. D3 lies in the subtropical climate zone (Keppen classification: Cwa). Geologically, this region is bounded by the Indo-Gangetic alluvial plains in the North and East, by Thar Desert in the West and by old Aravalli hill ranges in the South. There is a ridge trending along NNE-SSW direction which constitutes a small area of Delhi’s terrain which is otherwise generally flat [28]. Seasonally, the year can be divided into four main periods. Summer is experienced in the months of March–June followed by monsoon months of July, August, and September. Postmonsoon months are October and November while the period of December–February constitutes the winter season. The maximum temperature ranges from 41 to 45°C in peak summer season and the minimum temperature in winter season is in the range of 3–6°C in coldest period of Dec.-Jan. For the present study, summer and winter conditions have been simulated for about two weeks in each season. The period from 29 April 2009 00 : 00 UTC to 14 May 2009 00 : 00 UTC was chosen for summer simulations while for the winter season the simulation period was from 25 December 2009 00 : 00 UTC to 9 January 2010 00 : 00 UTC. The chosen simulation periods are representative of peak periods of their respective seasons with no recorded case of precipitation. The first day of each simulation was considered as spin-up period, and the remaining 14 days have been analysed. The vertical resolution of the model is 33 levels with pressure ranging from 976 mbar to 54 mbar. Up to a height of 2.2 km, there are 10 levels spaced at approximately 0.24, 0.32 0.42, 0.55, 0.7, 0.91, 1.16, 1.49, 1.87, and 2.26 km.

4. Data and Analysis

USGS 24 classification category land-use data was used for interpolating topography and land use with spatial resolution of 30′′. NCEP final analysis data (FNL) of resolution 1° was used as input for initial and boundary conditions to the model. The sources of these datasets are available for download from WRF model website (http://www.mmm.ucar.edu/wrf/users/downloads.html). The output has been generated from the model at every 30 minutes. The data used for validation of surface variables has been derived from http://www.wunderground.com/. This website collects weather conditions directly from weather stations located in various countries which are owned by government agencies and international airports [29]. In Delhi, there are two aeronautical weather stations of Indian Meteorological Department which are part of WMO World Weather Information Service. First is the aeronautical weather station located at Indira Gandhi International Airport in Delhi (28.57°N, 77.12°E) which will be referred to as “station 1” hereafter in this text. Within Delhi region, this is the only location for which data is available in public domain at an interval of 30-60 minutes. The other weather station is at Safdarjung Airport (28.58°N 77.2°E) and surface data for this site is available at 3 hour interval. This station will be referred to as “station 2” hereafter in this text. The observations are taken using standards techniques and practices as laid out by WMO [30]. Statistical measures such as correlation score ( ), root mean squared error (RMSE), fractional bias (FB), and number of estimated values within a factor of 2 (FAC2) [31, 32] have been used for evaluation of model performance. Values of these measures were estimated for all key meteorological variables at the surface, namely, temperature, wind speed, and relative humidity for each simulation.

Station 2 is also the source of radiosonde measurements which has been used for analysis of upper air variables such as temperature and wind speed. Upper air data is available for two soundings per day. The observed upper air data has been accessed from Integrated Global Radiosonde Archive of National Climatic Data Centre [33].

5. Results and Discussions

5.1. Surface Variables
5.1.1. Temperature and Relative Humidity

Statistical evaluation for both stations 1 and 2 is presented in Tables 2, 3, 4, and 5 where the best values obtained for each case are written in bold. It can be noted that simulation 7 works satisfactorily for temperature in terms of statistical evaluation reporting least RMSEs and fractional bias and highest correlations with observed data for temperature. Figure 2 displays diurnal variation of observed and estimated surface temperature by simulation 7 averaged over each hour for the entire study period comprising 14 days each of both summer and winter seasons for both the observation points. Top panel shows hourly variation at station 1 and bottom panel displays 3-hourly variation for station 2. By and large, there is a general tendency towards overprediction of daytime temperature averages while underprediction for night-time temperatures during summer period. However, a consistent overprediction is shown for the winter period; however night-time prediction is far closer to observations than the daytime predictions. Different studies in past have observed both underprediction [6, 34] and overestimation [7, 35] for temperature, but the result of present study is closest to the one noted by Hanna et al. [36] which was carried out in Oklahoma and Kansas (US) that observed underestimation during night-time and overestimation in daytime. In India, Das et al. [37] observed underprediction for temperature during monsoon season.

Simulation 7 is an experimental case with Pleim-Xiu (PX) land surface model, Pleim surface layer scheme, and Asymmetric Convective Model for PBL. The PX LSM coupled with the ACM2 PBL and Pleim surface layer scheme has been found to be well suited for extended simulations (i.e., more than a week such as in present study) where the indirect soil moisture and temperature nudging scheme leads to more accurate near-surface meteorology [22].

Overall, relative humidity is underpredicted by all simulations in the model as indicated by positive fractional bias in Tables 25. As with temperature, for relative humidity also, some studies have observed overestimation such as in Mediterranean regions [1] while some others have observed underprediction such as in dry western Africa [7].

For relative humidity also, simulation 7 is the best-performing simulation for both summer and winter seasons with consistently least RMSE and FB and highest and FAC2. The PX LSM simulates the evolution of soil moisture and temperature in two layers (0-1 and 1–100 cm); canopy moisture consists of algorithms for nudging of soil moisture and soil temperature. The algorithms of the coupled ACM2 PBL allow for smooth transition from eddy diffusion in stable conditions to the combined local and nonlocal transport in unstable conditions. This makes this combination well suited for consistent transport of moisture fluxes within the PBL.

From the above discussions, it may be concluded that simulation 7 performs best for both temperature and relative humidity.

5.1.2. Wind Speed and Direction

Wind direction can be better represented through wind rose diagram as all directions together with wind speed can be clearly displayed. Wind roses have been constructed for summer and winter simulation period in Figures 3 and 4 for station 1. Since hourly data is not available for station 2, it has not been included for wind rose analysis. However, statistical analysis of wind speed has been performed for estimations at both stations 1 and 2 in Tables 25. A general tendency by the model to overpredict wind speed is observed by all simulations which has also been observed in many earlier studies [1, 2, 4, 36]. As noted in earlier studies, wind speed is significantly affected by local fluctuations especially in highly unstable conditions such as in summer [6], and thus wind sensitivities tend to have more “noise” [38]. Nonetheless, for present study, simulation 4 clearly shows the best performance amongst all the 9 cases for both summer and winter conditions at both stations. Major difference between default case and simulation 4 lies in the two PBL schemes used, that is, the local and nonlocal approach, respectively. The YSU scheme (used in simulation 4) takes into account momentum and mass transfer from large-scale eddies, whereas the MYJ scheme (in default case simulation) determines diffusion and turbulence locally [7, 39].

Wind direction, being an angular quantity, has not been analysed statistically because statistical measures are based on arithmetic difference between observed and estimated quantities which is not applicable for vector parameters like wind direction. From wind roses, however, it can be noted that simulation 4 is the closest to observed pattern in terms of both direction and speed for summer and winter season (Figures 4 and 5).

5.2. Impact of Domain Resolution

WRF is suitable for use across scales ranging from meters to thousands of kilometers (http://www.mmm.ucar.edu/wrf/users/). Simulations over higher-resolution domains are computationally expensive. Thus, comparison of simulations at different domain resolutions was undertaken in the present study. Figures 5 and 6 display time series of the difference between observed values and best-performing simulations for each variable, that is, temperature (simulation 7), wind speed (simulation 4), and relative humidity (simulation 7) at different domain resolutions (d01: 18 km, d02: 6 km, d03: 2 km) for summer season at station 1 and 2, respectively. From the time series plots, it is clear that a resolution of 2 km does not significantly improve the simulation results when compared with 18 and 6 km resolution cases. At 18 km resolution, the linear distance between observation station point and model grid point was about 932 m while at 2 km resolution, this distance was 117 m. The difference between observed and simulated values at these stations is quite similar even when the output is obtained at 18 km resolution or 6 km resolution. Similar results were obtained for winter period also, though the plots are not shown here. Higher resolutions are needed for certain studies like local ambient temperature or air quality studies, and improved model inputs should be provided for such cases. Improved model simulations at higher resolutions may perhaps require better model inputs at finer resolutions such as meteorological parameters, land-use data, and assimilation with more local data sources.

5.3. Upper Air Variables
5.3.1. Temperature and Wind Speed Profile in Boundary Layer

Figure 7 depicts height versus temperature in the boundary layer for two days of summer (simulation 7) at 0530 and 1730 hours (local time) of May 1 and May 8, 2009. Estimated profile is quite close to the observed profile for all the considered cases.

As explained in Section 5.2, wind speed is greatly influenced by local factors and thus not suitable for validation with point observation. However, to have a general idea of performance of the model for vertical profile of wind speed, profile was obtained from simulation 4 for two summer days (Figure 8). Deviations of simulated profile from the observed profile increase with increase in height. This tendency has also been observed in some earlier studies [2, 36].

6. Conclusions

As discussed in the previous sections, there is no single combination of physics options that performs best for all the meteorological parameters. However, the present study suggests the suitable options for different meteorological variables under winter and summer conditions in a subtropical environment. Temperature and relative humidity are best estimated by coupled Pleim-Xiu ACM scheme while wind speed and direction are best simulated by MM5 surface layer LSM and YSU PBL schemes both for summer and winter seasons. Therefore, the model configuration should be chosen from the viewpoint of the objective of the study undertaken. The best-performing options for temperature and relative humidity are likely to work better for fog-related studies while best performing options for wind velocity can perform better for dispersion studies. The main conclusions of this study are as follows.(i)A case study for different physics options applicable in WRF (version 3.1.1) model for surface and upper air meteorological parameters for summer and winter seasons over subtropical region of Delhi has been performed.(ii)In general, the selection of surface layer, land surface model, and planetary boundary layer scheme was found to have more impact on output in comparison to microphysics and cumulus parameterization for both seasons.(iii)Identification of model options for a given parameter and season is helpful for a particular model application. MM5-YSU combination works well for wind speed, and Pleim-ACM combination gives better temperature and moisture estimations.(iv)Nested domains (18, 6 and 2 km) with higher resolutions were not helpful in improving the simulation results as per the current availability of the data. (v)Overall, the present case study shows that the model performance is satisfactory over the subtropical region of Delhi.

Acknowledgments

The present study is a part of research project “Atmospheric Chemical Modelling in relation to Photochemical Smog over Delhi” funded by the Department of Science and Technology (DST), Government of India (DST No: SR/S4/AS : 249/05). The second author acknowledges Junior Research Fellowship support from DST.