Role of land state in a high resolution mesoscale model for simulating the Uttarakhand heavy rainfall event over India

In 2013, Indian summer monsoon witnessed a very heavy rainfall event (>30 cm/day) over Uttarakhand in north India, claiming more than 5000 lives and property damage worth approximately 40 billion USD. This event was associated with the interaction of two synoptic systems, i.e., intensified subtropical westerly trough over north India and north-westward moving monsoon depression formed over the Bay of Bengal. The event had occurred over highly variable terrain and land surface characteristics. Although global models predicted the large scale event, they failed to predict realistic location, timing, amount, intensity and distribution of rainfall over the region. The goal of this study is to assess the impact of land state conditions in simulating this severe event using a high resolution mesoscale model. The land conditions such as multi-layer soil moisture and soil temperature fields were generated from High Resolution Land Data Assimilation (HRLDAS) modelling system. Two experiments were conducted namely, (1) CNTL (Control, without land data assimilation) and (2) LDAS, with land data assimilation (i.e., with HRLDAS-based soil moisture and temperature fields) using Weather Research and Forecasting (WRF) modelling system. Initial soil moisture correlation and root mean square error for LDAS is 0.73 and 0.05, whereas for CNTL it is 0.63 and 0.053 respectively, with a stronger heat low in LDAS. The differences in wind and moisture transport in LDAS favoured increased moisture transport from Arabian Sea through a convectively unstable region embedded within two low pressure centers over Arabian Sea and Bay of Bengal. The improvement in rainfall is significantly correlated to the persistent generation of potential vorticity (PV) in LDAS. Further, PV tendency analysis confirmed that the increased generation of PV is due to the enhanced horizontal PV advection component rather than the diabatic heating terms due to modified flow fields. These results suggest that, two different synoptic systems merged by the strong interaction of moving PV columns resulted in the strengthening and further amplification of the system over the region in LDAS. This study highlights the importance of better representation of the land surface fields for improved prediction of localized anomalous weather event over India.     

Comparison of AMSR-2 wind speed and sea surface temperature with moored buoy observations over the Northern Indian Ocean

The Northern Indian Ocean (NIO) is unique due to seasonal reversal of wind patterns, the formation of vortices and eddies which make satellite observations arduous. The veracity of sea surface wind (SSW) and sea surface temperature (SST) products of sun-synchronous AMSR-2 satellite are compared with high-temporal moored buoy observations over the NIO. The two year-long (2013–2014) comparisons reveal that the root-mean-square-error (RMSE) of AMSR-2 SST and SSW is \(<0.4{^{\circ }}\hbox {C}\) and \(<1.5\hbox { ms}^{-1}\), respectively, which are within the error range prescribed for the AMSR-2 satellite (\(\pm 0.8{^{\circ }}\hbox {C}\), \(\pm 1.5\hbox { ms}^{-1})\). The SST–wind relation is analyzed using data both from the buoy and satellite. As a result, the low-SST is associated with low-wind condition (positive slope) in the northern part of the Bay of Bengal (BoB), while low SST values are associated with high wind conditions (negative slope) over the southern BoB. Moreover, the AMSR-2 displayed larger slope for SST–wind relation and could be mainly due to overestimation of SST and underestimation of wind as compared to the buoy. The AMSR-2 SSW exhibited higher error during post-monsoon followed by monsoon season and could be attributed to the high wind conditions associated with intense oceanic vortices. The study suggests that the AMSR-2 products are reliable and can be used in tropical air–sea interactions, meso-scale features, and weather and climate studies.     

Impact of Doppler weather radar data on thunderstorm simulation during STORM pilot phase—2009

This study assesses the impact of Doppler weather radar (DWR) data (reflectivity and radial wind) assimilation on the simulation of severe thunderstorms (STS) events over the Indian monsoon region. Two different events that occurred during the Severe Thunderstorms Observations and Regional Modeling (STORM) pilot phase in 2009 were simulated. Numerical experiments—3DV (assimilation of DWR observations) and CNTL (without data assimilation)—were conducted using the three-dimensional variational data assimilation technique with the Advanced Research Weather Research and Forecasting model (WRF-ARW). The results show that consistent with prior studies the 3DV experiment, initialized by assimilation of DWR observations, performed better than the CNTL experiment over the Indian region. The enhanced performance was a result of improved representation and simulation of wind and moisture fields in the boundary layer at the initial time in the model. Assimilating DWR data caused higher moisture incursion and increased instability, which led to stronger convective activity in the simulations. Overall, the dynamic and thermodynamic features of the two thunderstorms were consistently better simulated after ingesting DWR data, as compared to the CNTL simulation. In the 3DV experiment, higher instability was observed in the analyses of thermodynamic indices and equivalent potential temperature (θ _e) fields. Maximum convergence during the mature stage was also noted, consistent with maximum vertical velocities in the assimilation experiment (3DV). In addition, simulated hydrometeor (water vapor mixing ratio, cloud water mixing ratio, and rain water mixing ratio) structures improved with the 3DV experiment, compared to that of CNTL. From the higher equitable threat scores, it is evident that the assimilation of DWR data enhanced the skill in rainfall prediction associated with the STS over the Indian monsoon region. These results add to the body of evidence now which provide consistent and notable improvements in the mesoscale model results over the Indian monsoon region after assimilating DWR fields.     

Improvement of Monsoon Depressions Forecast with Assimilation of Indian DWR Data Using WRF-3DVAR Analysis System

An attempt is made to evaluate the impact of Doppler Weather Radar (DWR) radial velocity and reflectivity in Weather Research and Forecasting (WRF)-3D variational data assimilation (3DVAR) system for prediction of Bay of Bengal (BoB) monsoon depressions (MDs). Few numerical experiments are carried out to examine the individual impact of the DWR radial velocity and the reflectivity as well as collectively along with Global Telecommunication System (GTS) observations over the Indian monsoon region. The averaged 12 and 24 h forecast errors for wind, temperature and moisture at different pressure levels are analyzed. This evidently explains that the assimilation of radial velocity and reflectivity collectively enhanced the performance of the WRF-3DVAR system over the Indian region. After identifying the optimal combination of DWR data, this study has also investigated the impact of assimilation of Indian DWR radial velocity and reflectivity data on simulation of the four different summer MDs that occurred over BoB. For this study, three numerical experiments (control no assimilation, with GTS and GTS along with DWR) are carried out to evaluate the impact of DWR data on simulation of MDs. The results of the study indicate that the assimilation of DWR data has a positive impact on the prediction of the location, propagation and development of rain bands associated with the MDs. The simulated meteorological parameters and tracks of the MDs are reasonably improved after assimilation of DWR observations as compared to the other experiments. The root mean square errors (RMSE) of wind fields at different pressure levels, equitable skill score and frequency bias are significantly improved in the assimilation experiments mainly in DWR assimilation experiment for all MD cases. The mean Vector Displacement Errors (VDEs) are significantly decreased due to the assimilation of DWR observations as compared to the CNTL and 3DV_GTS experiments. The study clearly suggests that the performance of the model simulation for the intense convective system which influences the large scale monsoonal flow is significantly improved after assimilation of the Indian DWR data from even one coastal locale within the MDs track.     

Customization of WRF-ARW model with physical parameterization schemes for the simulation of tropical cyclones over North Indian Ocean

The convection and planetary boundary layer (PBL) processes play significant role in the genesis and intensification of tropical cyclones (TCs). Several convection and PBL parameterization schemes incorporate these processes in the numerical weather prediction models. Therefore, a systematic intercomparison of performance of parameterization schemes is essential to customize a model. In this context, six combinations of physical parameterization schemes (2 PBL Schemes, YSU and MYJ, and 3 convection schemes, KF, BM, and GD) of WRF-ARW model are employed to obtain the optimum combination for the prediction of TCs over North Indian Ocean. Five cyclones are studied for sensitivity experiments and the out-coming combination is tested on real-time prediction of TCs during 2008. The tracks are also compared with those provided by the operational centers like NCEP, ECMWF, UKMO, NCMRWF, and IMD. It is found that the combination of YSU PBL scheme with KF convection scheme (YKF) provides a better prediction of intensity, track, and rainfall consistently. The average RMSE of intensity (13?hPa in CSLP and 11?m?s^?1 in 10-m wind), mean track, and landfall errors is found to be least with YKF combination. The equitable threat score (ETS) of YKF combination is more than 0.2 for the prediction of 24-h accumulated rainfall up to 125?mm. The vertical structural characteristics of cyclone inner core also recommend the YKF combination for Indian seas cyclones. In the real-time prediction of 2008 TCs, the 72-, 48-, and 24-h mean track errors are 172, 129, and 155?km and the mean landfall errors are 125, 73, and 66?km, respectively. Compared with the track of leading operational agencies, the WRF model is competing in 24?h (116?km error) and 72?h (166?km) but superior in 48-h (119?km) track forecast.     

Probability of occurrence of monthly and seasonal winter precipitation over Northwest India based on antecedent-monthly precipitation

Theoretical and Applied Climatology - Winter (December, January, and February (DJF)) precipitation over northwest India (NWI) is mainly associated with the eastward moving mid-latitude synoptic...     

Simulation of heavy rainfall events over Indian monsoon region using WRF-3DVAR data assimilation system

We present the results of the impact of the 3D variational data assimilation (3DVAR) system within the Weather Research and Forecasting (WRF) model to simulate three heavy rainfall events (25–28 June 2005, 29–31 July 2004, and 7–9 August 2002) over the Indian monsoon region. For each event, two numerical experiments were performed. In the first experiment, namely the control simulation (CNTL), the low-resolution global analyses are used as the initial and boundary conditions of the model. In the second experiment (3DV-ANA), the model integration was carried out by inserting additional observations in the model’s initial conditions using the 3DVAR scheme. The 3DVAR used surface weather stations, buoy, ship, radiosonde/rawinsonde, and satellite (oceanic surface wind, cloud motion wind, and cloud top temperature) observations obtained from the India Meteorological Department (IMD). After the successful inclusion of additional observational data using the 3DVAR data assimilation technique, the resulting reanalysis was able to successfully reproduce the structure of convective organization as well as prominent synoptic features associated with the mid-tropospheric cyclones (MTC). The location and intensity of the MTC were better simulated in the 3DV-ANA as compared to the CNTL. The results demonstrate that the improved initial conditions of the mesoscale model using 3DVAR enhanced the location and amount of rainfall over the Indian monsoon region. Model verification and statistical skill were assessed with the help of available upper-air sounding data. The objective verification further highlighted the efficiency of the data assimilation system. The improvements in the 3DVAR run are uniformly better as compared to the CNTL run for all the three cases. The mesoscale 3DVAR data assimilation system is not operational in the weather forecasting centers in India and a significant finding in this study is that the assimilation of Indian conventional and non-conventional observation datasets into numerical weather forecast models can help improve the simulation accuracy of meso-convective activities over the Indian monsoon region. Results from the control experiments also highlight that weather and regional climate model simulations with coarse analysis have high uncertainty in simulating heavy rain events over the Indian monsoon region and assimilation approaches, such as the 3DVAR can help reduce this uncertainty.     

A Study on Simulation of Heavy Rainfall Events Over Indian Region with ARW-3DVAR Modeling System

An attempt is made to evaluate the impact of the three dimensional variational (3DVAR) data assimilation within the Weather Research Forecasting (WRF) modeling system to simulate two heavy rainfall events which occured on 26–27 July 2005 and 27–30 July 2006. During the 26–27 July 2005 event, the unprecedented localized intense rainfall 90–100 cm was recorded over the northeast parts of Mumbai city; however, southern parts received only 10 cm. Model simulation with the data assimilation experiment is reasonably well predicted for the rainfall intensity (800 mm) in 24 h and with accurate location over Mumbai agreeing with observation. Divergence, vorticity, vertical velocity and moisture parameters are evaluated during the various stages of the event. It is noticed that maximum convergence and vorticity during the mature stage; at the same time the vertical velocity also follows a similar trend during the period in the assimilation experiment. Vorticity budget terms over the location of heavy rainfall revealed that the contribution of the positive tilting term produced positive vorticity which triggered the convection and negative contribution to vorticity from the tilting term to precede the dissipation of the system. Model simulations from the second rain event, the off-shore trough at sea level along the west coast of India, is well represented after assimilation of observations during day-1 and day-2 as compared to the control simulations; the orientation of the off-shore trough is well matched with that of the observed. The intensity and spatial distribution of the rainfall has considerably improved in the assimilation simulation. The statistical skill scores also revealed that the precipitation forecast during the period has appreciably improved due to assimilation of observations. The results of this study indicate a positive impact of the 3DVAR assimilation on the simulation of heavy rainfall events.     

Simulation of Bay of Bengal Tropical Cyclones with WRF Model: Impact of Initial and Boundary Conditions

An attempt is made to delineate the relative performances and credentials of GFS, FNL, and NCMRWF global analyses/forecast products as initial and boundary conditions (IBCs) to the WRF-ARW model in the simulation of four Bay of Bengal tropical cyclones (TCs). The results suggest that FNL could simulate horizontal advection of vorticity maxima at 850 hPa; hence, the tracks are more realistic with least errors as compared to GFS and NCMRWF. The mean landfall errors for 24-, 48-, and 72-hour forecasts are 73, 41, and 72 km, respectively. The TC intensity is well captured by NCMRWF IBCs, as it could predict 850 hPa vorticity maxima. The 24-hour accumulated rainfall is well simulated with FNL, and equitable threat score is more than 0.2 up to 100 mm with minimum bias.