Impact of period and timescale of FDDA analysis nudging on the numerical simulation of tropical cyclones in the Bay of Bengal

In this study, the impact of four-dimensional data assimilation (FDDA) analysis nudging is examined on the prediction of tropical cyclones (TC) in the Bay of Bengal to determine the optimum period and timescale of nudging. Six TCs (SIDR: November 13–16, 2007; NARGIS: April 29–May 02, 2008; NISHA: November 25–28, 2008; AILA: May 23–26, 2009; LAILA: May 18–21, 2010; JAL: November 04–07, 2010) were simulated with a doubly nested Weather Research and Forecasting (WRF) model with a horizontal resolution of 9 km in the inner domain. In the control run for each cyclone, the National Centers for Environmental Prediction (NCEP) Global Forecast System (GFS) analysis and forecasts at 0.5° resolution are used for initial and boundary conditions. In the FDDA experiments available surface, upper air observations obtained from NCEP Atmospheric Data Project (ADP) data sets were used for assimilation after merging with the first guess through objective analysis procedure. Analysis nudging experiments with different nudging periods (6, 12, 18, and 24 h) indicated a period of 18 or 24 h of nudging during the pre-forecast stage provides maximum impact on simulations in terms of minimum track and intensity forecasts. To determine the optimum timescale of nudging, two cyclone cases (NARGIS: April 28–May 02, 2008; NISHA: November 25–28, 2008) were simulated varying the inverse timescales as 1.0e?4 to 5.0e?4 s^?1 in steps of 1.0e?4 s^?1. A positive impact of assimilation is found on the simulated characteristics with a nudging coefficient of either 3.0e?4 or 4.0e?4 s^?1 which corresponds to a timescale of about 1 h for nudging dynamic (u,v) and thermodynamical (t,q) fields.     

Real-time prediction of a severe cyclone ��Jal�� over Bay of Bengal using a high-resolution mesoscale model WRF (ARW)

Real-time predictions for the JAL severe cyclone formed in November 2010 over Bay of Bengal using a high-resolution Weather Research and Forecasting (WRF ARW) mesoscale model are presented. The predictions are evaluated with different initial conditions and assimilation of observations. The model is configured with two-way interactive nested domains and with fine resolution of 9?km for the region covering the Bay of Bengal. Simulations are performed with NCEP GFS 0.5° analysis and forecasts for initial/boundary conditions. To examine the impact of initial conditions on the forecasts, eleven real-time numerical experiments are conducted with model integration starting at 00, 06, 12, 18 UTC 4 Nov, 5?Nov and 00, 06, 12 UTC 6 Nov and all ending at 00 UTC 8 Nov. Results indicated that experiments starting prior to 18 UTC 04 Nov produced faster moving cyclones with higher intensity relative to the IMD estimates. The experiments with initial time at 18 UTC 04 Nov, 00 UTC 05 Nov and with integration length of 78?h and 72?h produced best prediction comparable with IMD estimates of the cyclone track and intensity parameters. To study the impact of observational assimilation on the model predictions FDDA, grid nudging is performed separately using (1) land-based automated weather stations (FDDAAWS), (2) MODIS temperature and humidity profiles (FDDAMODIS), and (3) ASCAT and OCEANSAT wind vectors (FDDAASCAT). These experiments reduced the pre-deepening period of the storm by 12?h and produced an early intensification. While the assimilation of AWS data has shown meagre impact on intensity, the assimilation of scatterometer winds produced an intermittent drop in intensity in the peak stage. The experiments FDDAMODIS and FDDAQSCAT produced minimum error in track and intensity estimates for a 90-h prediction of the storm.     

An HWRF-based ensemble assessment of the land surface feedback on the post-landfall intensification of Tropical Storm Fay (2008)

While tropical cyclones (TCs) usually decay after landfall, Tropical Storm Fay (2008) initially developed a storm central eye over South Florida by anomalous intensification overland. Unique to the Florida peninsula are Lake Okeechobee and the Everglades, which may have provided a surface feedback as the TC tracked near these features around the time of peak intensity. Analysis is done with the use of an ensemble model-based approach with the Developmental Testbed Center (DTC) version of the Hurricane WRF (HWRF) model using an outer domain and a storm-centered moving nest with 27- and 9-km grid spacing, respectively. Choice of land surface parameterization and small-scale surface features may influence TC structure, dictate the rate of TC decay, and even the anomalous intensification after landfall in model experiments. Results indicate that the HWRF model track and intensity forecasts are sensitive to three features in the model framework: land surface parameterization, initial boundary conditions, and the choice of planetary boundary layer (PBL) scheme. Land surface parameterizations such as the Geophysical Fluid Dynamics Laboratory (GFDL) Slab and Noah land surface models (LSMs) dominate the changes in storm track, while initial conditions and PBL schemes cause the largest changes in the TC intensity overland. Land surface heterogeneity in Florida from removing surface features in model simulations shows a small role in the forecast intensity change with no substantial alterations to TC track.     

Flood forecasts based on multi-model ensemble precipitation forecasting using a coupled atmospheric-hydrological modeling system

The recent improvement of numerical weather prediction (NWP) models has a strong potential for extending the lead time of precipitation and subsequent flooding. However, uncertainties inherent in precipitation outputs from NWP models are propagated into hydrological forecasts and can also be magnified by the scaling process, contributing considerable uncertainties to flood forecasts. In order to address uncertainties in flood forecasting based on single-model precipitation forecasting, a coupled atmospheric-hydrological modeling system based on multi-model ensemble precipitation forecasting is implemented in a configuration for two episodes of intense precipitation affecting the Wangjiaba sub-region in Huaihe River Basin, China. The present study aimed at comparing high-resolution limited-area meteorological model Canadian regional mesoscale compressible community model (MC2) with the multiple linear regression integrated forecast (MLRF), covering short and medium range. The former is a single-model approach; while the latter one is based on NWP models [(MC2, global environmental multiscale model (GEM), T213L31 global spectral model (T213)] integrating by a multiple linear regression method. Both MC2 and MLRF are coupled with Chinese National Flood Forecasting System (NFFS), MC2-NFFS and MLRF-NFFS, to simulate the discharge of the Wangjiaba sub-basin. The evaluation of the flood forecasts is performed both from a meteorological perspective and in terms of discharge prediction. The encouraging results obtained in this study demonstrate that the coupled system based on multi-model ensemble precipitation forecasting has a promising potential of increasing discharge accuracy and modeling stability in terms of precipitation amount and timing, along with reducing uncertainties in flood forecasts and models. Moreover, the precipitation distribution of MC2 is more problematic in finer temporal and spatial scales, even for the high resolution simulation, which requests further research on storm-scale data assimilation, sub-grid-scale parameterization of clouds and other small-scale atmospheric dynamics.     

Assessing the influence of assimilating radar-observed radial winds on the simulation of a tropical cyclone

Observations by Doppler weather radar are crucial for nowcasting and short-time forecasting of severe weather events as they bring in refined information of the atmosphere. However, due to the inevitable noises and non-meteorological signals, they cannot be assimilated straightforwardly into a numerical model. In the present study, assimilation of the radial component of wind velocity observed by two Doppler radars is performed in the numerical simulation of Supertyphoon Rammasun (2014) just before its landfall. After several quality-control steps, the radar-observed radial velocities are de-aliased, noise-reduced and assimilated into the model to improve initial conditions for the high-resolution simulation. Results show that only when using global background error covariance matrix can the observational increment be properly assimilated into the model, correcting large-scale background steering flow and yielding a simulated track close to the observed one. However, little improvement is found in simulating the TC core-scale structures by the assimilation of radar velocity as compared to the radar-observed flow, primarily due to the insufficient spatial resolution of the model that may lead to the incorrect representation of the TC core structure and the rejection of some core-region observations during the data assimilation procedure. Moreover, assimilation-induced asymmetries consume a certain portion of mean kinetic energy, preventing the simulated Rammasun from axisymmetrization and thus intensification as compared with the non-assimilated experiment.     

The effect of a surface data assimilation technique and the traditional four-dimensional data assimilation on the simulation of a monsoon depression over India using a mesoscale model

The objective of this study is to investigate the impact of a surface data assimilation (SDA) technique, together with the traditional four-dimensional data assimilation (FDDA), on the simulation of a monsoon depression that formed over India during the field phase of the 1999 Bay of Bengal Monsoon Experiment (BOBMEX). The SDA uses the analyzed surface data to continuously assimilate the surface layer temperature as well as the water vapor mixing ratio in the mesoscale model. The depression for the greater part of this study was offshore and since successful application of the SDA would require surface information, a method of estimating surface temperature and surface humidity using NOAA-TOVS satellites was used. Three sets of numerical experiments were performed using a coupled mesoscale model. The first set, called CONTROL, uses the NCEP (National Center for Environmental Prediction) reanalysis for the initial and lateral boundary conditions in the MM5 simulation. The second and the third sets implemented the SDA of temperature and moisture together with the traditional FDDA scheme available in the MM5 model. The second set of MM5 simulation implemented the SDA scheme only over the land areas, and the third set extended the SDA technique over land as well as sea. Both the second and third sets of the MM5 simulation used the NOAA-TOVS and QuikSCAT satellite and conventional upper air and surface meteorological data to provide an improved analysis. The results of the three sets of MM5 simulations are compared with one another and with the analysis and the BOBMEX 1999 buoy, ship, and radiosonde observations. The predicted sea level pressure of both the model runs with assimilation resembles the analysis closely and also captures the large-scale structure of the monsoon depression well. The central sea level pressures of the depression for both the model runs with assimilation were 2–4 hPa lower than the CONTROL. The results of both the model runs with assimilation indicate a larger spatial area as well as increased rainfall amounts over the coastal regions after landfall compared with the CONTROL. The impact of FDDA and SDA, the latter over land, resulted in reduced errors of the following: 1.45 K in temperature, 0.39 m s⁻¹ in wind speed, and 14° in wind direction compared with the BOBMEX buoy observation, and 1.43 m s⁻¹ in wind speed, 43° in wind direction, and 0.75% in relative humidity compared with the CONTROL. The impact of SDA over land and sea compared with SDA over land only showed a further marginal reduction of errors: 0.23 K in air temperature (BOBMEX buoy) and 1.33 m s⁻¹ in wind speed simulations.     

Development of extended WRF variational data assimilation system (WRFDA) for WRF non-hydrostatic mesoscale model

The paper intends to present the development of the extended weather research forecasting data assimilation (WRFDA) system in the framework of the non-hydrostatic mesoscale model core of weather research forecasting system (WRF-NMM), as an imperative aspect of numerical modeling studies. Though originally the WRFDA provides improved initial conditions for advanced research WRF, we have successfully developed a unified WRFDA utility that can be used by the WRF-NMM core, as well. After critical evaluation, it has been strategized to develop a code to merge WRFDA framework and WRF-NMM output. In this paper, we have provided a few selected implementations and initial results through single observation test, and background error statistics like eigenvalues, eigenvector and length scale among others, which showcase the successful development of extended WRFDA code for WRF-NMM model. Furthermore, the extended WRFDA system is applied for the forecast of three severe cyclonic storms: Nargis (27 April–3 May 2008), Aila (23–26 May 2009) and Jal (4–8 November 2010) formed over the Bay of Bengal. Model results are compared and contrasted within the analysis fields and later on with high-resolution model forecasts. The mean initial position error is reduced by 33% with WRFDA as compared to GFS analysis. The vector displacement errors in track forecast are reduced by 33, 31, 30 and 20% to 24, 48, 72 and 96 hr forecasts respectively, in data assimilation experiments as compared to control run. The model diagnostics indicates successful implementation of WRFDA within the WRF-NMM system.     

Impact of 3D Var GSI-ENKF hybrid data assimilation system

The hybrid two-way coupled 3DEnsVar assimilation system was tested with the NCMRWF global data assimilation forecasting system. At present, this system consists of T574L64 deterministic model and the grid-point statistical interpolation analysis scheme. In this experiment, the analysis system is modified with a two-way coupling with an 80 member Ensemble Kalman Filter of T254L64 resolution and runs are carried out in parallel to the operational system for the Indian summer monsoon season (June–September) for the year 2015 to study its impact. Both the assimilation systems are based on NCEP GFS system. It is found that hybrid assimilation marginally improved the quality of the forecasts of all variables over the deterministic 3D Var system, in terms of statistical skill scores and also in terms of circulation features. The impact of the hybrid system in prediction of extreme rainfall and cyclone track is discussed.     

Optimizing the routing of ship’s tropical cyclone avoidance based on the numerical forecasts

In this paper, an algorithm to design a shortest-time route for a ship to avoid a tropical cyclone (TC) is proposed. The proposed algorithm takes into account the influence of the changing winds and sea waves on ship’s speed and the corresponding risk using the forecasts from a numerical weather prediction model. Experimental results show that the new algorithm is able to save more time comparing with the traditional sector diagram typhoon avoidance method. In the application of the new algorithm to the navigation practice, the distance between adjacent alternative waypoints should be adjusted to meet the navigational needs, and the route should be updated simultaneously with TC forecasts from a numerical weather prediction model.     

Operational tropical cyclone intensity prediction—an empirical technique

One very specific operational requirement of the Tropical Cyclone (TC) Programme of the Regional Specialized Meteorological Centre, New Delhi is to provide 12-hourly forecasts valid up to 48 h (preferably 72 h) on the intensity of cyclones over the southern Indian Seas. In this paper, a simple empirical model for predicting the intensity of TCs occurring in the Bay of Bengal is proposed. The model parameter has been determined from a database assembled on 30 recent cyclones, and the model itself is based on the assumption that a TC intensifies exponentially. A method for correcting the forecast during subsequent observation hours (6- or 12-h intervals) is also presented. The results show that the forecast skill for forecasts of up to 48 h is reasonably good. The absolute mean errors are less than 12 knots for 48-h forecasts, with the forecast skill decreasing with time. With the incorporation of a correction procedure based on the latest observations, some improvement in the forecast skill can be obtained. The model is expected to be useful to operational forecasters.     

Application of the representer method for parameter estimation in numerical reservoir models

A data assimilation method was applied to estimate poorly known parameters (permeabilities) in a numerical reservoir model. Most variational methods for data assimilation are based on the assumption that the model is perfect except for the poorly known parameters. The representer method allows also for model errors, i.e. for uncertainties in the state variables (pressures and saturations). The method is based on minimizing a cost functional, assuming all the errors and parameters to be multivariate Gaussian random variables with given mean and covariances. The uncertain parameters and variables are expanded into a finite sum of basis functions called representers, and the gradients of the cost functional are obtained with an adjoint method. This approach gives an optimal parametrization in the sense that the final result is equal to the solution of the full inverse problem. The method was tested on a simple one-dimensional model to simulate two-phase (oil-water) flow through a heterogeneous reservoir. The results show that the method is able to provide an acceptable estimate of the permeability field. We used pressure measurements from a small number of observation wells in between the injection and production wells, but the representer method could be used equally well to assimilate data from other sources. The method appears to be a promising data assimilation tool for applications in reservoir engineering.     

High resolution numerical weather prediction over the Indian subcontinent

In this study, the Florida State University Global Spectral Model (FSUGSM), in association with a high-resolution nested regional spectral model (FSUNRSM), is used for short-range weather forecasts over the Indian domain. Three-day forecasts for each day of August 1998 were performed using different versions of the FSUGSM and FSUNRSM and were compared with the observed fields (analysis) obtained from the European Center for Medium Range Weather Forecasts (ECMWF). The impact of physical initialization (a procedure that assimilates observed rain rates into the model atmosphere through a set of reverse algorithms) on rainfall forecasts was examined in detail. A very high nowcasting skill for precipitation is obtained through the use of high-resolution physical initialization applied at the regional model level. Higher skills in wind and precipitation forecasts over the Indian summer monsoon region are achieved using this version of the regional model with physical initialization. A relatively new concept, called the ‘multimodel/multianalysis superensemble’ is described in this paper and is applied for the wind and precipitation forecasts over the Indian subcontinent. Large improvement in forecast skills of wind at 850 hPa level over the Indian subcontinent is shown possible through the use of the multimodel superensemble. The multianalysis superensemble approach that uses the latest satellite data from the Tropical Rainfall Measuring Mission (TRMM) and the Defense Meteorological Satellite Program (DMSP) has shown significant improvement in the skills of precipitation forecasts over the Indian monsoon region.     

Experimental superensemble forecasts of tropical cyclones over the Bay of Bengal

This study entails the implementation of an experimental real time forecast capability for tropical cyclones over the Bay of Bengal basin of North Indian Ocean. This work is being built on the experience gained from a number of recent studies using the concept of superensemble developed at the Florida State University (FSU). Real time hurricane forecasts are one of the major components of superensemble modeling at FSU. The superensemble approach of training followed by real time forecasts produces the best forecasts for tracks and intensity (up to 5 days) of Atlantic hurricanes and Pacific typhoons. Improvements in track forecasts of about 25–35% compared to current operational forecast models has been noted over the Atlantic Ocean basin. The intensity forecasts for hurricanes are only marginally better than the best models. In this paper, we address tropical cyclone forecasts over the Bay of Bengal for the years 1996–2000. The main result from this study is that the position and intensity errors for tropical cyclone forecasts over the Bay of Bengal from the multimodel superensemble are generally less than those of all of the participating models during 1- to 3-day forecasts. Some of the major tropical cyclones, such as the November 1996 Andhra Pradesh cyclone and October 1999 Orissa super cyclone were well handled by this superensemble approach. A conclusion from this study is that the proposed approach may be a viable way to construct improved forecasts of Bay of Bengal tropical cyclone positions and intensity.     

Observation targeting with a second-order adjoint method for increased predictability

The efficiency of current adjoint-based observations targeting strategies in variational data assimilation is closely determined by the underlying assumption of a linear propagation of initial condition errors into the model forecasts. A novel targeting strategy is proposed in the context of four-dimensional variational data assimilation (4D-Var) to account for nonlinear error growth as the forecast lead time increases. A quadratic error growth model is shown to maintain the accuracy in tracking the nonlinear evolution of initial condition perturbations, as compared to the first-order approximation. A second-order adjoint model is used to provide the derivative information that is necessary in the higher-order Taylor series approximation. The observation targeting approach relies on the dominant eigenvectors of the Hessian matrix associated with a specific forecast error aspect as an indicator of the directions of largest quadratic error growth. A comparative qualitative analysis between observation targeting based on first- and second-order adjoint information is presented in idealized 4D-Var experiments with a two-dimensional global shallow-water model. The results indicate that accounting for the quadratic error growth in the targeting strategy is of particular benefit as the forecast lead time increases.     

结合糙率校正的河网水情数据同化

为解决河网水动力模型重要参数糙率与水力状态量水位、流量的同步校正问题,以糙率和水力状态量作为河网非线性动态系统变量,采用扩展卡尔曼滤波,构建结合糙率动态校正的河网水情数据同化模型.通过算例计算,系统分析了水位动态噪声水平、糙率动态噪声水平、糙率初始值及测站个数对模型校正的影响.结果表明:模型能够有效用于水位状态量的实时校正;靠近测站的糙率校正值趋于真值,远离测站的糙率校正值趋于初始值;通过调整糙率动态噪声水平,可以有效控制糙率的修正量,防止糙率修正过大而引起计算失效问题.     

Forecasting of cyclone Viyaru and Phailin by NWP-based cyclone prediction system (CPS) of IMD – an evaluation

An objective NWP-based cyclone prediction system (CPS) was implemented for the operational cyclone forecasting work over the Indian seas. The method comprises of five forecast components, namely (a) Cyclone Genesis Potential Parameter (GPP), (b) Multi-Model Ensemble (MME) technique for cyclone track prediction, (c) cyclone intensity prediction, (d) rapid intensification, and (e) predicting decaying intensity after the landfall. GPP is derived based on dynamical and thermodynamical parameters from the model output of IMD operational Global Forecast System. The MME technique for the cyclone track prediction is based on multiple linear regression technique. The predictor selected for the MME are forecast latitude and longitude positions of cyclone at 12-hr intervals up to 120 hours forecasts from five NWP models namely, IMD-GFS, IMD-WRF, NCEP-GFS, UKMO, and JMA. A statistical cyclone intensity prediction (SCIP) model for predicting 12 hourly cyclone intensity (up to 72 hours) is developed applying multiple linear regression technique. Various dynamical and thermodynamical parameters as predictors are derived from the model outputs of IMD operational Global Forecast System and these parameters are also used for the prediction of rapid intensification. For forecast of inland wind after the landfall of a cyclone, an empirical technique is developed. This paper briefly describes the forecast system CPS and evaluates the performance skill for two recent cyclones Viyaru (non-intensifying) and Phailin (rapid intensifying), converse in nature in terms of track and intensity formed over Bay of Bengal in 2013. The evaluation of performance shows that the GPP analysis at early stages of development of a low pressure system indicated the potential of the system for further intensification. The 12-hourly track forecast by MME, intensity forecast by SCIP model and rapid intensification forecasts are found to be consistent and very useful to the operational forecasters. The error statistics of the decay model shows that the model was able to predict the decaying intensity after landfall with reasonable accuracy. The performance statistics demonstrates the potential of the system for improving operational cyclone forecast service over the Indian seas.     

Integrated geothermal modeling on different scales in the Northeast German basin

Geothermal modeling is an important part of large-scale basin studies. Based on a new 3D structural model of the Northeast German basin, the present day regional geothermal field is modeled. Range and regional trend of the modeled temperature values are in agreement with the published data. Due to the high spatial resolution, the calculated temperature distribution provides additional information with respect to areas where no measured data is available. The results are used as input and boundary parameters for small-scale models of geothermal energy production. In general, in many regions not enough data is available to define all necessary physical or chemical parameters for modeling. In this context, data obtained from the large-scale model help to constrain unknown parameters. Subsequently, the small-scale model is used to simulate various production schemes focusing on enhanced predictions with respect to the possible lifetime of such installations. The simulation results also show the need for elaborated models if reliable predictions of the temperature evolution are required.     

Impact of hybrid GSI analysis using ETR ensembles

Performance of a hybrid assimilation system combining 3D Var based NGFS (NCMRWF Global Forecast System) with ETR (Ensemble Transform with Rescaling) based Global Ensemble Forecast (GEFS) of resolution T-190L28 is investigated. The experiment is conducted for a period of one week in June 2013 and forecast skills over different spatial domains are compared with respect to mean analysis state. Rainfall forecast is verified over Indian region against combined observations of IMD and NCMRWF. Hybrid assimilation produced marginal improvements in overall forecast skill in comparison with 3D Var. Hybrid experiment made significant improvement in wind forecasts in all the regions on verification against mean analysis. The verification of forecasts with radiosonde observations also show improvement in wind forecasts with the hybrid assimilation. On verification against observations, hybrid experiment shows more improvement in temperature and wind forecasts at upper levels. Both hybrid and operational 3D Var failed in prediction of extreme rainfall event over Uttarakhand on 17 June, 2013.     

The Indian summer monsoon as revealed by the NCMRWF system

In this study, we present the mean seasonal features of the Indian summer monsoon circulation in the National Centre for Medium Range Weather Forecasting (NCMRWF) global data assimilation and forecast system. The large-scale budgets of heat and moisture are examined in the analyzed and model atmosphere. The daily operational analyses and forecasts (day 1 through day 5) produced for the summer seasons comprising June, July and August of 1995 and 1993 have been considered for the purpose. The principal aim of the study is two-fold. Primarily, to comprehend the influence of the systematic errors over the Indian summer monsoon, secondarily, to analyze the performance of the model in capturing the interseasonal variability. The heat and moisture balances show reduction in the influx of heat and moisture in the model forecasts compared to the analyzed atmosphere over the monsoon domain. Consequently, the diabatic heating also indicates reducing trend with increase in the forecast period. In effect, the strength of Indian summer monsoon, which essentially depends on these parameters, weakens considerably in the model forecasts. Despite producing feeble monsoon circulation, the model captures interseasonal variability realistically. Although, 1995 and 1993 are fairly normal monsoon seasons, the former received more rainfall compared to the latter in certain pockets of the monsoon domain. This is clearly indicated by the analyzed and model atmosphere in terms of energetics.     

Case studies of seasonal rainfall forecasts for Hong Kong and its vicinity using a regional climate model

Seasonal climate forecasts are one of the most promising tools for providing early warnings for natural hazards such as floods and droughts. Using two case studies, this paper documents the skill of a regional climate model in the seasonal forecasting of below normal rainfall in southern China during the rainy seasons of July–August–September 2003 and April–May–June 2004. The regional model is based on the Regional Spectral Model of the National Centers for Environmental Prediction of the United States. It is the first time that the model has been applied to a region dominated by the East Asian Monsoon. The article shows that the regional climate model, when being forced by reasonably good forecasts from a global model, can generate useful seasonal rainfall forecasts for the region, where it is dominated by the East Asia monsoon. The spatial details of the dry conditions obtained from the regional climate model forecast are also found to be comparable with the observed distribution.