Assimilation Of Precipitation Data Research Articles

A Sequential Non-Gaussian Approach for Precipitation Data Assimilation

AbstractIn two recent studies, the authors presented a new data assimilation (DA) method for precipitation observations that does not require Gaussianity or linearity assumptions. The method, called localized ensemble mosaic assimilation (LEMA), initializes the new ensemble forecast by relaxing the background ensemble (prior) toward a single analysis composed of different column states taken from the ensemble members with the lowest error in the precipitation forecast. However, a limitation of the LEMA method is that relaxing the background ensemble toward that analysis severely reduces the spread of the ensemble, thus, limiting its usefulness for cycled DA applications. This study presents a new version of LEMA, called localized ensemble mosaic assimilation sequence (LEMAS), suitable for cycled DA operations. LEMAS constructs an ensemble of analysis mosaics using a small group of members closer to the observations instead of only the closest one. The new ensemble forecast is then initialized by recentering the prior ensemble around the mean of the analysis ensemble while scaling the original background perturbations to match the spread of the analysis mosaics. A series of ideal and real DA experiments are used to evaluate the potential of LEMAS for the assimilation of hourly accumulation observations. A comparison of LEMAS with the local ensemble transform Kalman filter (LETKF) using idealized experiments shows that LEMAS produces similar or slightly better forecast quality than the LETKF in temperature, water vapor, winds, and precipitation. Extending this comparison to real DA experiments assimilating Stage-IV precipitation observations shows that both methods produce precipitation forecasts of comparable quality.

A Heuristic Approach for Precipitation Data Assimilation: Effect of Forecast Errors and Assimilation of NCEP Stage IV Precipitation Analyses

Abstract Recently, Pérez Hortal et al. introduced a simple data assimilation (DA) technique named localized ensemble mosaic assimilation (LEMA) for the assimilation of radar-derived precipitation observations. The method constructs an analysis by assigning to each model grid point the information from the ensemble member that is locally closest to the precipitation observations. This study explores the effects of the forecasts errors in the performance of the method using a series of observing system simulation experiments (OSSEs) with different magnitudes of forecast errors employing a small ensemble of 20 members. The ideal experiments show that LEMA is able to produce forecasts with considerable and long-lived error reductions in the fields of precipitation, temperature, humidity, and wind. Nonetheless, the quality of the analysis deteriorates with increasing forecast errors beyond the spread of the ensemble. To overcome this limitation, we expand the spread of the ensemble used to construct the analysis mosaic by considering states at different times and states from forecasts initialized at different times (lagged forecasts). The ideal experiments show that the additional information in the expanded ensemble improves the performance of LEMA, producing larger and long-lived improvements in the state variables and in the precipitation forecast quality. Finally, the potential of LEMA is explored in real DA experiments using actual Stage IV precipitation observations. When LEMA uses only the background members, the quality of the precipitation forecast shows small or no improvements. However, the expanded ensemble improves the LEMA’s effectiveness, producing larger and more persistent improvements in precipitation forecasts.

Stably numerical solving inverse boundary value problem for data assimilation

In remote sensing data assimilation, inverse methods are often used to retrieve the boundary conditions from some of the observations in the interior scanning region. Needless to say, the inverse boundary value problem (IBVP) is ill-conditioned in general. Ultimately, the data assimilation must include mechanisms that enable them to overcome the numerical instability for solving IBVP. In this paper, we begin by studying the behavior of ill-conditioned IBVP, and find that the condition number varies with the number of the observations and distribution locations of the observations. Next, we define the number of equivalent independent equations as a novel measurement of ill-conditioning of a problem, which can measure the degree of bad conditions. Furthermore, the novel measure can answer how many additional observations are needed to stabilize the retrieving problem, and where additional observations are fixed up. Finally, we illustrate the proposed methodology by applying it to the study of precipitation data assimilation, with a particular emphasis on the analysis of the effect of the number of observations and their distribution locations. The new methodology appears to be particularly efficient in tackling the instability of the retrieving problem in data assimilation.

A Heuristic Approach for Precipitation Data Assimilation: Characterization Using OSSEs

Abstract We introduce a new technique for the assimilation of precipitation observations, the localized ensemble mosaic assimilation (LEMA). The method constructs an analysis by selecting, for each vertical column in the model, the ensemble member with precipitation at the ground that is locally closest to the observed values. The proximity between the modeled and observed precipitation is determined by the mean absolute difference of precipitation intensity, converted to reflectivity and measured over a spatiotemporal window centered at each grid point of the model. The underlying hypothesis of the approach is that the ensemble members that are locally closer to the observed precipitation are more probable to be closer to the “truth” in the state variables than the other members. The initial conditions for the new forecast are obtained by nudging the background states toward the mosaic of the closest ensemble members (analysis) over a 30 min time interval, reducing the impacts of the imbalances at the boundaries between the different selected members. The potential of the method is studied using observing system simulation experiments (OSSEs) employing a small ensemble of 20 members. The ensemble is produced by the WRF Model, run at a horizontal grid spacing of 20 km. The experiments lend support to the validity of the hypothesis and allow the determination of the optimal parameters for the approach. In the context of OSSE, this new data assimilation technique is able to produce forecasts with considerable and long-lived error reductions in the fields of precipitation, temperature, humidity, and wind.

Precipitation Data Assimilation System Based on a Neural Network and Case-Based Reasoning System

There are several methods to forecast precipitation, but none of them is accurate enough since predicting precipitation is very complicated and influenced by many factors. Data assimilation systems (DAS) aim to increase the prediction result by processing data from different sources in a general way, such as a weighted average, but have not been used for precipitation prediction until now. A DAS that makes use of mathematical tools is complex and hard to carry out. In our paper, machine learning techniques are introduced into a precipitation data assimilation system. After summarizing the theoretical construction of this method, we take some practical weather forecasting experiments and the results show that the new system is effective and promising.

Combined Assimilation of Satellite Precipitation and Soil Moisture: A Case Study Using TRMM and SMOS Data

This paper presents a framework that enables simultaneous assimilation of satellite precipitation and soil moisture observations into the coupled Weather Research and Forecasting (WRF) and Noah land surface model through variational approaches. The authors tested the framework by assimilating precipitation data from the Tropical Rainfall Measuring Mission (TRMM) and soil moisture data from the Soil Moisture Ocean Salinity (SMOS) satellite. The results show that assimilation of both TRMM and SMOS data can effectively improve the forecast skills of precipitation, top 10-cm soil moisture, and 2-m temperature and specific humidity. Within a 2-day time window, impacts of precipitation data assimilation on the forecasts remain relatively constant for forecast lead times greater than 6 h, while the influence of soil moisture data assimilation increases with lead time. The study also demonstrates that the forecast skill of precipitation, soil moisture, and near-surface temperature and humidity are further improved when both the TRMM and SMOS data are assimilated. In particular, the combined data assimilation reduces the prediction biases and root-mean-square errors, respectively, by 57% and 6% (for precipitation); 73% and 27% (for soil moisture); 17% and 9% (for 2-m temperature); and 33% and 11% (for 2-m specific humidity).

Assimilating the global satellite mapping of precipitation data with the Nonhydrostatic Icosahedral Atmospheric Model (NICAM)

AbstractThis study aims to propose two new approaches to improve precipitation forecasts from numerical weather prediction (NWP) models through effective data assimilation of satellite‐derived precipitation. The assimilation of precipitation data is known to be very difficult mainly because of highly non‐Gaussian statistics of precipitation variables. Following Lien et al., this study addresses the non‐Gaussianity issue by applying the Gaussian transformation (GT) based on the empirical cumulative distribution function (CDF) of precipitation. We propose a method that constructs the CDF with only recent 1 month samples, without using a long period of samples needed previously. We also propose a method to use the inverse GT, with which we can obtain realistic precipitation fields from biased NWP model outputs. We assimilate the Japan Aerospace eXploration Agency's Global Satellite Mapping of Precipitation (GSMaP) data into the Nonhydrostatic Icosahedral Atmospheric Model (NICAM) at 112 km horizontal resolution. Assimilating the GSMaP data results in improved weather forecasts compared to the control experiment assimilating only rawinsonde data. We find that horizontal observation thinning is necessary, probably due to the horizontal observation‐error correlations in the GSMaP data. We also obtained precipitation fields similar to GSMaP from the NICAM precipitation forecasts by using the inverse GT, leading to an improved precipitation forecast.

Precipitation data assimilation in WRFDA 4D-Var: implementation and application to convection-permitting forecasts over United States

Precipitation data assimilation has been developed in the Weather Research and Forecasting model data assimilation system (WRFDA) using four-dimensional variational (4D-Var) approach. Unlike other conventional data, precipitation is an integral quantity and it is not included as a control variable in WRFDA. A simplified Kessler scheme is used in tangent linear and adjoint model. Precipitation data are directly assimilated in WRFDA 4D-Var, and the assimilation of precipitation will have feedback to all the control variables via the constraint of the linearized physics package. Firstly, we present single observation experiments to exhibit how dynamic, thermodynamic and moisture fields are adjusted by assimilating rainfall information. Then, the National Centers for Environmental Prediction Stage IV precipitation data are assimilated for a heavy rainfall case on 9 June 2010 at a convection-permitting model setting (4-km). Finally, we conducted one-week experiments to further validate the robustness of the results for precipitation assimilation. Results show that precipitation assimilation has a positive impact on model fields, particularly on the low-level humidity. For the impact on precipitation forecasts, it indicates that precipitation assimilation reduces spin-up time efficiently, removes false alarms and produces model forecast precipitation closer to the observations through the changes in temperature, moisture and wind imposed in the analyses. The impact from precipitation assimilation persists up to three hours on average.

On the performance of the new NWP nowcasting system at the Danish Meteorological Institute during a heavy rain period

At the Danish Meteorological Institute, the NWP nowcasting system has been enhanced to include assimilation of 2D precipitation rates derived from weather radar observations. The assimilation is performed using a nudging-based technique. Here the rain rates are used to estimate the changes in the vertical profile of horizontal divergence needed to induce the observed rain rate. Verification of precipitation forecasts for a 17-day period in August 2010 based on the NWP nowcasting system is presented and compared to a reference without assimilation of precipitation data. In Denmark, this period was particularly rainy, with several heavy precipitation events. Three of these events are studied in detail. The verification is mainly based on scatter plots and fractions skill scores, which give scale-dependant indicators of the spatial skill of the forecasts. The study shows that the inclusion of precipitation observations has a positive impact on the spatial skill of the forecasts. This positive impact is the largest in the first hour, and then gradually decreases. On the average, the forecasts with assimilation of precipitation are skilful after 4 h on scales down to a few tens of kilometers. For the events studied, the assimilation improves the forecasted frequencies of heavy and light precipitation relative to the control, while there is some tendency to overpredict intermediate precipitation levels.

Dynamical Precipitation Downscaling for Hydrologic Applications Using WRF 4D-Var Data Assimilation: Implications for GPM Era

Abstract The objective of this study is to develop a framework for dynamically downscaling spaceborne precipitation products using the Weather Research and Forecasting (WRF) Model with four-dimensional variational data assimilation (4D-Var). Numerical experiments have been conducted to 1) understand the sensitivity of precipitation downscaling through point-scale precipitation data assimilation and 2) investigate the impact of seasonality and associated changes in precipitation-generating mechanisms on the quality of spatiotemporal downscaling of precipitation. The point-scale experiment suggests that assimilating precipitation can significantly affect the precipitation analysis, forecast, and downscaling. Because of occasional overestimation or underestimation of small-scale summertime precipitation extremes, the numerical experiments presented here demonstrate that the wintertime assimilation produces downscaled precipitation estimates that are in closer agreement with the reference National Centers for Environmental Prediction stage IV dataset than similar summertime experiments. This study concludes that the WRF 4D-Var system is able to effectively downscale a 6-h precipitation product with a spatial resolution of 20 km to hourly precipitation with a spatial resolution of less than 10 km in grid spacing—relevant to finescale hydrologic applications for the era of the Global Precipitation Measurement mission.

Assimilating Satellite Observations of Clouds and Precipitation into NWP Models

No Abstract available.

Dynamic Downscaling of the North American Monsoon with the NCEP–Scripps Regional Spectral Model from the NCEP CFS Global Model

Abstract The June–September (JJAS) 2000–07 NCEP coupled Climate Forecasting System (CFS) global hindcasts are downscaled over the North and South American continents with the NCEP–Scripps Regional Spectral Model (RSM) with anomaly nesting (AN) and without bias correction (control). A diagnosis of the North American monsoon (NAM) in CFS and RSM hindcasts is presented here. RSM reduces errors caused by coarse resolution but is unable to address larger-scale CFS errors even with bias correction. CFS has relatively weak Great Plains and Gulf of California low-level jets. Low-level jets are strengthened from downscaling, especially after AN bias correction. The RSM NAM hydroclimate shares similar flaws with CFS, with problematic diurnal and seasonal variability. Flaws in both diurnal and monthly variability are forced by erroneous convection-forced divergence outside the monsoon core region in eastern and southern Mexico. NCEP reanalysis shows significant seasonal variability errors, and AN shows little improvement in regional-scale flow errors. The results suggest that extreme caution must be taken when the correction is applied relative to reanalyses. Analysis also shows that North American Regional Reanalysis (NARR) NAM seasonal variability has benefited from precipitation data assimilation, but many questions remain concerning NARR’s representation of NAM.

Mesoscale Data Assimilation Experiment in the WWRP B08RDP

A mesoscale data assimilation experiment was performed for the Beijing 2008 Olympics Research and Development Project (B08RDP) under the World Weather Research Program (WWRP) conducted during the period around the Beijing 2008 Olympic Games.In this experiment, the Japan Meteorological Agency (JMA) hydrostatic mesoscale 4D-Var analysis system was modified and utilized to produce accurate initial fields over the China area. In addition to the conventional observations, precipitation data observed by rain-gauges were assimilated as well as data produced from a nowcasting system of the Australian Bureau of Meteorology.The analysis system with rainfall data outperformed other experiments, i.e., the mesoscale analysis without precipitation data and the JMA global analysis, with its quantitative precipitation forecast. This shows that the assimilation of precipitation data can have a positive impact on subsequent model forecasts and that it would be effective even for areas where successive rainfall observations are sparse.

Using Precipitation Observations in a Mesoscale Short-Range Ensemble Analysis and Forecasting System

Abstract A simple method to assimilate precipitation data from a synthesis of radar and gauge data is developed to operate alongside an ensemble Kalman filter that assimilates hourly surface observations. The mesoscale ensemble forecast system consists of 25 members with 30-km grid spacing and incorporates variability in both initial and boundary conditions and model physical process schemes. The precipitation assimilation method only incorporates information on when and where rainfall is observed. Model temperature and water vapor mixing ratio profiles at each grid point are modified if rainfall is observed but not predicted, or if rainfall is predicted but not observed. These modifications act to either increase or decrease, respectively, the likelihood that precipitation develops at that grid point. Two cases are examined in which this technique is applied to assimilate precipitation data every 15 min from 1200 to 1800 UTC, while hourly surface observations are also assimilated at the same time using the more sophisticated ensemble Kalman filter approach. Results show that the simple method for assimilating precipitation data helps the model develop precipitation where it is observed, resulting in the precipitation area being reproduced more accurately than in the run without precipitation-data assimilation, while not negatively influencing the positive results from the surface data assimilation. Improvement is also seen in the reliability of precipitation probabilities for a 1 mm h−1 threshold after the assimilation period, indicating that assimilating precipitation data may provide improved forecasts of the mesoscale environment for a few hours.

Adjoint Sensitivity of Surface Precipitation to Initial Conditions

Abstract The adjoint version of the Global Environmental Multiscale model including a comprehensive package of simplified and linearized physical processes (large-scale condensation, deep moist convection, vertical diffusion, and subgrid-scale orographic effects) is used to evaluate the sensitivity of surface precipitation to initial conditions for up to 24 h for two meteorological systems: a midlatitude front and a tropical cyclone. Such diagnostics are useful to improve the understanding on variational assimilation of precipitation data. In agreement with a similar study, the largest sensitivity is found with respect to the temperature field for both stratiform and convective precipitation. Close to the observation time and for stratiform precipitation, the sensitivity with respect to specific humidity is rather large, which corroborates conclusions from previous one-dimensional variational data assimilation experimentations. The sensitivity is then reduced significantly after the observation time. The sensitivities of surface precipitation to the wind components and to specific humidity are comparable and are at a maximum at the observation time. The sensitivity to the surface pressure is always much smaller than the sensitivity to the other variables. In general, sensitivities are largest at the observation time and then decrease. However, for the midlatitude perturbation, the sensitivity is enhanced after 12 h for stratiform precipitation and also for convective precipitation using a scheme based on the moisture convergence closure. This results from a dynamical coupling upstream of the area of interest through baroclinic instability as evidenced by vertically backward-tilted sensitivities. Such enhancement is not observed for the tropical case. The tangent-linear approximation remains acceptable for accumulated precipitation up to 24 h but is rather poor for instantaneous rain rates. The variational assimilation of accumulated precipitation should thus be favored.

“1D+4DVAR” Assimilation of NCEP Stage-IV Radar and Gauge Hourly Precipitation Data at ECMWF

Abstract The one- plus four-dimensional variational data assimilation (“1D+4DVAR”) method currently run in operations at ECMWF with rain-affected radiances from the Special Sensor Microwave Imager is used to study the potential impact of assimilating NCEP stage-IV analyses of hourly accumulated surface precipitation over the U.S. mainland. These data are a combination of rain gauge measurements and observations from the high-resolution Doppler Next-Generation Weather Radars. Several 1D+4DVAR experiments have been run over a month in spring 2005. First, the quality of the precipitation forecasts in the control experiment is assessed. Then, it is shown that the impact of the assimilation of the additional rain observations on global scores of dynamical fields and temperature is rather neutral, while precipitation scores are improved for forecast ranges up to 12 h. Additional 1D+4DVAR experiments in which all moisture-affected observations are removed over the United States demonstrate that the NCEP stage-IV precipitation data on their own can clearly be beneficial to the analyses and subsequent forecasts of the moisture field. This result suggests that the potential impact of precipitation observations is overshadowed by the influence of other high-quality humidity observations, in particular, radiosondes. It also confirms that the assimilation of precipitation observations has the ability to improve the quality of moisture analyses and forecasts in data-sparse regions. Finally, the limitations inherent in the current assimilation of precipitation data, their implications for the future, and possible ways of improvement are discussed.

Use of Scale Recursive Estimation for assimilation of precipitation data from TRMM (PR and TMI) and NEXRAD

The paper shows an application of Scale Recursive Estimation (SRE) used to assimilate rainfall rates estimated during a storm event from three remote sensing devices. These are the TMI radiometer and the PR radar, carried on board of the TRMM satellite and the KNQA Memphis Weather Surveillance radar, belonging to the NEXRAD network, each one providing rain rate estimates at a different spatial scale. The variability of rain rate process in scales is modeled as a multiplicative random cascade, including spatial intermittence. The observational noise in the estimates is modeled according to a multiplicative error. System estimation, including process and observational noise, is carried out using Maximum Likelihood Estimation implemented by a scale recursive Expectation Maximization (EM) algorithm. As a result, new rainfall rate estimates are obtained that feature decreased estimation error as compared to those coming from each device alone. The performance of the SRE-EM approach is compared with that of the latest methods proposed for data fusion of multisensor estimates. The proposed approach improves the current methods adopted for SRE and provides an alternative for data fusion in the field of precipitation.

A Numerical Study of the 1996 Saguenay Flood Cyclone: Effect of Assimilation of Precipitation Data on Quantitative Precipitation Forecasts

Abstract A one-dimensional variational (1DVAR) technique is applied to assimilate rain gauge precipitation data to extend the predictability of the Saguenay flood cyclone associated with a trough-merger event on 19–21 July 1996 in the Saguenay-Lac-Saint-Jean region of Quebec, Canada. Two 60-h simulations initialized at 0000 UTC 19 July were performed with the Canadian Mesoscale Compressible Community (MC2) model. The control (CTL) and NCEP simulations were initialized with the enhanced temperature and moisture profiles obtained from the 1DVAR scheme and the NCEP reanalysis data, respectively. Compared to observations, the CTL simulation reasonably reproduced the observed mass and wind fields and showed a marked improvement in the threat scores for heavy precipitation. The CTL run captured the observed spatial and temporal distribution of precipitation but overpredicted the area of precipitation. Sensitivity experiments showed that the threat (bias) scores are less (somewhat) sensitive to the specification of the observation error of the precipitation data. Of the four precipitation systems present at model initial time, the systems in the vicinity of the southern trough had the biggest impact on the threat score. Potential vorticity diagnostics of the CTL simulation suggested that the initial temperature and moisture field near the southern trough decreased the condensational heating relative to NCEP. This resulted in a stronger zonal wind component in the upper levels associated with the southern trough in CTL that retarded the eastward propagation of the northern trough, resulting in a correct placement of the surface precipitation and an improvement in the threat scores relative to NCEP.

Assimilation of Precipitation Data to the JMA Mesoscale Model with a Four-dimensional Variational Method and its Impact on Precipitation Forecasts

Two sets of forecast-analysis cycle experiments were performed with and without direct assimilation of precipitation amounts using JMA mesoscale four-dimensional variational data assimilation system (Meso4D-Var). With a devised cost function of precipitation observation, which is derived from the exponential distribution, Meso 4D-Var successfully assimilated precipitation data in various weather situations throughout a one-month experiment period. The result of experiments shows that the precipitation assimilation by Meso 4D-Var improves the model forecasts of both weak and moderate precipitation and it ameliorates a spin-up problem of the model precipitation.

A Nudging Scheme for the Assimilation of Precipitation Data into a Mesoscale Model

Abstract A nudging procedure for the assimilation of rainfall data into a mesoscale model [the Bologna Limited Area Model (BOLAM)] has been developed in order to improve short-range forecasting. The scheme modifies the model specific humidity profiles at every time step, according to the difference between observed and forecast precipitation. Different relaxation procedures are applied depending on the precipitation type (large-scale or convective rain), as estimated by the model itself. Optimizations of the nudging parameters and assessment of the scheme's performance was carried out in an idealized framework [(Observing System Simulation Experiment) OSSE-type strategy], implementing a lagged forecast scheme. Two events were selected for this purpose, both characterized by heavy precipitation in the Mediterranean Basin. The first was a severe orographic rainfall event, associated with the passage of a frontal system over the Alps during the Mesoscale Alpine Programme (MAP) field phase, in September 1999....