Radar rainfall uncertainty modelling influenced by wind

Radar‐based estimates of rainfall are affected by many sources of uncertainties, which would propagate through the hydrological model when radar rainfall estimates are used as input or initial conditions. An elegant solution to quantify these uncertainties is to model the empirical relationship between radar measurements and rain gauge observations (as the ‘ground reference’). However, most current studies only use a fixed and uniform model to represent the uncertainty of radar rainfall, without consideration of its variation under different synoptic regimes. Wind is such a typical weather factor, as it not only induces error in rain gauge measurements but also causes the raindrops observed by weather radar to drift when they reach the ground. For this reason, as a first attempt, this study introduces the wind field into the uncertainty model and designs the radar rainfall uncertainty model under different wind conditions. We separate the original dataset into three subsamples according to wind speed, which are named as WDI (0–2 m/s), WDII (2–4 m/s) and WDIII (>4 m/s). The multivariate distributed ensemble generator is introduced and established for each subsample. Thirty typical events (10 at each wind range) are selected to explore the behaviours of uncertainty under different wind ranges. In each time step, 500 ensemble members are generated, and the values of 5th to 95th percentile values are used to produce the uncertainty bands. Two basic features of uncertainty bands, namely dispersion and ensemble bias, increase significantly with the growth of wind speed, demonstrating that wind speed plays a considerable role in influencing the behaviour of the uncertainty band. On the basis of these pieces of evidence, we conclude that the radar rainfall uncertainty model established under different wind conditions should be more realistic in representing the radar rainfall uncertainty. This study is only a start in incorporating synoptic regimes into rainfall uncertainty analysis, and a great deal of more effort is still needed to build a realistic and comprehensive uncertainty model for radar rainfall data. Copyright © 2014 John Wiley & Sons, Ltd.     

Modelling radar-rainfall estimation uncertainties using elliptical and Archimedean copulas with different marginal distributions

Abstract

Given that radar-based rainfall has been broadly applied in hydrological studies, quantitative modelling of its uncertainty is critically important, as the error of input rainfall is the main source of error in hydrological modelling. Using an ensemble of rainfall estimates is an elegant solution to characterize the uncertainty of radar-based rainfall and its spatial and temporal variability. This paper has fully formulated an ensemble generator for radar precipitation estimation based on the copula method. Each ensemble member is a probable realization that represents the unknown true rainfall field based on the distribution of radar rainfall (RR) error and its spatial error structure. An uncertainty model consisting of a deterministic component and a random error factor is presented based on the distribution of gauge rainfall conditioned on the radar rainfall (GR|RR). Two kinds of copulas (elliptical and Archimedean copulas) are introduced to generate random errors, which are imposed by the deterministic component. The elliptical copulas (e.g. Gaussian and t-copula) generate the random errors based on the multivariate distribution, typically of decomposition of the error correlation matrix using the LU decomposition algorithm. The Archimedean copulas (e.g. Clayton and Gumbel) utilize the conditional dependence between different radar pixels to obtain random errors. Based on those, a case application is carried out in the Brue catchment located in southwest England. The results show that the simulated uncertainty bands of rainfall encompass most of the reference raingauge measurements with good agreement between the simulated and observed spatial dependences. This indicates that the proposed scheme is a statistically reliable method in ensemble radar rainfall generation and is a useful tool for describing radar rainfall uncertainty.

Editor D. Koutsoyiannis; Associate editor S. Grimaldi     

Probabilistic radar rainfall nowcasts using empirical and theoretical uncertainty models

Weather radar has a potential to provide accurate short‐term (0–3 h) forecasts of rainfall (i.e. radar nowcasts), which are of great importance in warnings and risk management for hydro‐meteorological events. However, radar nowcasts are affected by large uncertainties, which are not only linked to limitations in the forecast method but also because of errors in the radar rainfall measurement. The probabilistic quantitative precipitation nowcasting approach attempts to quantify these uncertainties by delivering the forecasts in a probabilistic form. This study implements two forms of probabilistic quantitative precipitation nowcasting for a hilly area in the south of Manchester, namely, the theoretically based scheme [ensemble rainfall forecasts (ERF)‐TN] and the empirically based scheme (ERF‐EM), and explores which one exhibits higher predictive skill. The ERF‐TN scheme generates ensemble forecasts of rainfall in which each ensemble member is determined by the stochastic realisation of a theoretical noise component. The so‐called ERF‐EM scheme proposed and applied for the first time in this study, aims to use an empirically based error model to measure and quantify the combined effect of all the error sources in the radar rainfall forecasts. The essence of the error model is formulated into an empirical relation between the radar rainfall forecasts and the corresponding ‘ground truth’ represented by the rainfall field from rain gauges measurements. The ensemble members generated by the two schemes have been compared with the rain gauge rainfall. The hit rate and the false alarm rate statistics have been computed and combined into relative operating characteristic curves. The comparison of the performance scores for the two schemes shows that the ERF‐EM achieves better performance than the ERF‐TN at 1‐h lead time. The predictive skills of both schemes are almost identical when the lead time increases to 2 h. In addition, the relation between uncertainty in the radar rainfall forecasts and lead time is also investigated by computing the dispersion of the generated ensemble members. Copyright © 2013 John Wiley & Sons, Ltd.     

Analyses of the warm season rainfall climatology of the northeastern US using regional climate model simulations and radar rainfall fields

We examine the warm season (April-September) rainfall climatology of the northeastern US through analyses of high-resolution radar rainfall fields from the Hydro-NEXRAD system and regional climate model simulations using the weather research and forecasting (WRF) model. Analyses center on the 5-year period from 2003 to 2007 and the study area includes the New York-New Jersey metropolitan region covered by radar rainfall fields from the Fort Dix, NJ WSR-88D. The objective of this study is to develop and test tools for examining rainfall climatology, with a special focus on heavy rainfall. An additional emphasis is on rainfall climatology in regions of complex terrain, like the northeastern US, which is characterized by land-water boundaries, large heterogeneity in land use and cover, and mountainous terrain in the western portion of the region. We develop a 5-year record of warm season radar rainfall fields for the study region using the Hydro-NEXRAD system. We perform regional downscaling simulations for the 5-year study period using the WRF model. Radar rainfall fields are used to characterize the interannual, seasonal and diurnal variation of rainfall over the study region and to examine spatial heterogeneity of rainfall. Regional climate model simulations are characterized by a wet bias in the rainfall fields, with the largest bias in the high-elevation regions of the model domain. We show that model simulations capture broad features of the interannual, seasonal, and diurnal variation of rainfall. Model simulations do not capture spatial gradients in radar rainfall fields around the New York metropolitan region and land-water boundaries to the east. The model climatology of convective available potential energy (CAPE) is used to interpret the regional distribution of warm season rainfall and the seasonal and diurnal variability of rainfall. We use hydrologic and meteorological observations from July 2007 to examine the interactions of land surface processes and rainfall from a regional perspective.     

Ensemble rainfall forecasting with numerical weather prediction and radar‐based nowcasting models

The overall objective of this study is to improve the forecasting accuracy of the precipitation in the Singapore region by means of both rainfall forecasting and nowcasting. Numerical Weather Predication (NWP) and radar‐based rainfall nowcasting are two important sources for quantitative precipitation forecast. In this paper, an attempt to combine rainfall prediction from a high‐resolution mesoscale weather model and a radar‐based rainfall model was performed. Two rainfall forecasting methods were selected and examined: (i) the weather research and forecasting model (WRF); and (ii) a translation model (TM). The WRF model, at a high spatial resolution, was run over the domain of interest using the Global Forecast System data as initializing fields. Some heavy rainfall events were selected from data record and used to test the forecast capability of WRF and TM. Results obtained from TM and WRF were then combined together to form an ensemble rainfall forecasting model, by assigning weights of 0.7 and 0.3 weights to TM and WRF, respectively. This paper presented results from WRF and TM, and the resulting ensemble rainfall forecasting; comparisons with station data were conducted as well. It was shown that results from WRF are very useful as advisory of anticipated heavy rainfall events, whereas those from TM, which used information of rain cells already appearing on the radar screen, were more accurate for rainfall nowcasting as expected. The ensemble rainfall forecasting compares reasonably well with the station observation data. Copyright © 2012 John Wiley & Sons, Ltd.     

A comparison of gauge and radar precipitation data for simulating an extreme hydrological event in the Severn Uplands,UK

This paper provides a comparison of gauge and radar precipitation data sources during an extreme hydrological event. November–December 2006 was selected as a time period of intense rainfall and large river flows for the Severn Uplands, an upland catchment in the United Kingdom. A comparison between gauge and radar precipitation time‐series records for the event indicated discrepancies between data sources, particularly in areas of higher elevation. The HEC‐HMS rainfall‐runoff model was selected to assess the accuracy of the precipitation to simulate river flows for the extreme event. Gauge, radar and gauge‐corrected radar rainfall were used as model inputs. Universal cokriging was used to geostatistically interpolate gauge data with radar and elevation data as covariates. This interpolated layer was used to calculate the mean‐field bias and correct the radar composites. Results indicated that gauge‐ and gauge‐corrected radar‐driven models replicated flows adequately for the extreme event. Gauge‐corrected flow predictions produced an increase in flow prediction accuracy when compared with the raw radar, yet predictions were comparative in accuracy to those using the interpolated gauge network. Subsequent investigations suggested this was due to an adequate spatial and temporal resolution of the precipitation gauge network within the Severn Uplands. Results suggested that the six rain gauges could adequately represent precipitation variability of the Severn Uplands to predict flows at an approximately equal accuracy to that obtained by radar. Temporally, radar produced an increase in flow prediction accuracy in mountainous reaches once the gauge time step was in excessive of an hourly interval. Copyright © 2010 John Wiley & Sons, Ltd.     

Copula‐based uncertainty modelling: application to multisensor precipitation estimates

The multisensor precipitation estimates (MPE) data, available in hourly temporal and 4 km × 4 km spatial resolution, are produced by the National Weather Service and mosaicked as a national product known as Stage IV. The MPE products have a significant advantage over rain gauge measurements due to their ability to capture spatial variability of rainfall. However, the advantages are limited by complications related to the indirect nature of remotely sensed precipitation estimates. Previous studies confirm that efforts are required to determine the accuracy of MPE and their associated uncertainties for future use in hydrological and climate studies. So far, various approaches and extensive research have been undertaken to develop an uncertainty model. In this paper, an ensemble generator is presented for MPE products that can be used to evaluate the uncertainty of rainfall estimates. Two different elliptical copula families, namely, Gaussian and t‐copula are used for simulations. The results indicate that using t‐copula may have significant advantages over the well‐known Gaussian copula particularly with respect to extremes. Overall, the model in which t‐copula was used for simulation successfully generated rainfall ensembles with similar characteristics to those of the ground reference measurements. Copyright © 2010 John Wiley & Sons, Ltd.     

Experimental Seasonal Forecast of Monsoon 2005 Using T170L42 AGCM on PARAM Padma

As a part of the Experimental Extended Range Monsoon Prediction Experiment, ensemble mode seasonal runs for the monsoon season of 2005 were made using the National Centre for Environmental Prediction (NCEP), T170L42 AGCM. The seasonal runs were made using six initial atmospheric conditions based on the NCEP operational analysis and with forecast monthly sea-surface temperature (SST) of the NCEP Coupled forecast system (CFS). These simulations were carried out on the PARAM Padma supercomputer of Centre for Development of Advanced Computing (C-DAC), India. The model climatology was prepared by integrating the model for ten years using climatological SST as the lower boundary. The climatology of the model compares well with the observed, in terms of the spatial distribution of rainfall over the Indian land mass. The model-simulated rainfall compares well with the Tropical Rainfall Measuring Mission (TRMM) estimates for the 2005 monsoon season. Compared to the model climatology (7.81 mm/day), the model had simulated a normal rainfall (7.75 mm/day) for the year 2005 which is in agreement with the observations (99% of long-term mean). However, the model could not capture the observed increase in September rainfall from that of a low value in August 2005. The circulation patterns simulated by the model are also comparable to the observed patterns. The ensemble mean onset is found to be nearer to the observed onset date within one pentad.     

Interannual climate variability in South America: impacts on seasonal precipitation, extreme events, and possible effects of climate change

Interannual variability is an important modulator of synoptic and intraseasonal variability in South America. This paper seeks to characterize the main modes of interannual variability of seasonal precipitation and some associated mechanisms. The impact of this variability on the frequency of extreme rainfall events and the possible effect of anthropogenic climate change on this variability are reviewed. The interannual oscillations of the annual total precipitation are mainly due to the variability in austral autumn and summer. While autumn is the dominant rainy season in the northern part of the continent, where the variability is highest (especially in the northeastern part), summer is the rainy season over most of the continent, thanks to a summer monsoon regime. In the monsoon season, the strongest variability occurs near the South Atlantic Convergence Zone (SACZ), which is one of the most important features of the South American monsoon system. In all seasons but summer, the most important source of variability is ENSO (El Ni?o Southern Oscillation), although ENSO shows a great contribution also in summer. The ENSO impact on the frequency of extreme precipitation events is also important in all seasons, being generally even more significant than the influence on seasonal rainfall totals. Climate change associated with increasing emission of greenhouse gases shows potential to impact seasonal amounts of precipitation in South America, but there is still great uncertainty associated with the projected changes, since there is not much agreement among the models’ outputs for most regions in the continent, with the exception of southeastern South America and southern Andes. Climate change can also impact the natural variability modes of seasonal precipitation associated with ENSO.     

Conditional simulation of remotely sensed rainfall data using a non-Gaussian v-transformed copula

Quantification of rainfall and its spatial and temporal variability is extremely important for reliable hydrological and meteorological modeling. While rain gauge measurements do not provide reasonable areal representation of rainfall, remotely sensed precipitation estimates offer much higher spatial resolution. However, uncertainties associated with remotely sensed rainfall estimates are not well quantified. This issue is important considering the fact that uncertainties in input rainfall are the main sources of error in hydrologic processes. Using an ensemble of rainfall estimates that resembles multiple realizations of possible true rainfall, one can assess uncertainties associated with remotely sensed rainfall data. In this paper, ensembles are generated by imposing rainfall error fields over remotely sensed rainfall estimates. A non-Gaussian copula-based model is introduced for simulation of rainfall error fields. The v-transformed copula is employed to describe the dependence structure of rainfall error estimates without the influence of the marginal distribution. Simulations using this model can be performed unconditionally or conditioned on ground reference measurements such that rain gauge data are honored at their locations. The presented model is implemented for simulation of rainfall ensembles across the Little Washita watershed, Oklahoma. The results indicate that the model generates rainfall fields with similar spatio-temporal characteristics and stochastic properties to those of observed rainfall data.     

Integration of gauge and radar rainfall to enable best simulation of hydrological parameters

This study is about use of spatially distributed rain in physically based hydrological models. In recent years, spatially distributed radar rainfall data have become available. The distributed radar rain is used to precisely model hydrologic processes and it is more realistic than the past practice of distribution methods like Thiessen polygons. Radar provides a highly accurate spatial distribution of rainfall and greatly improves the basin average rainfall estimates. However, quantification of the exact amount of rainfall from radar observation is relatively difficult. Thus, the fundamental idea of this study is to apply hourly gauge and radar rainfall data in a distributed hydrological model to simulate hydrological parameters. Hence the comparison is made between the outcomes of the WetSpa model from radar rainfall distribution and gauge rainfall distributed by the Thiessen polygon technique. The comparative plots of the hydrograph and the results of hydrological components such as evapotranspiration, surface runoff, soil moisture, recharge and interflow, reflect the spatially distributed radar input performing well for model outflow simulation.

EDITOR D. Koutsoyiannis

ASSOCIATE EDITOR F. Pappenberger     

Towards flash-flood prediction in the dry Dead Sea region utilizing radar rainfall information

Flash-flood warning models can save lives and protect various kinds of infrastructure. In dry climate regions, rainfall is highly variable and can be of high-intensity. Since rain gauge networks in such areas are sparse, rainfall information derived from weather radar systems can provide useful input for flash-flood models. This paper presents a flash-flood warning model which utilizes radar rainfall data and applies it to two catchments that drain into the dry Dead Sea region. Radar-based quantitative precipitation estimates (QPEs) were derived using a rain gauge adjustment approach, either on a daily basis (allowing the adjustment factor to change over time, assuming available real-time gauge data) or using a constant factor value (derived from rain gauge data) over the entire period of the analysis. The QPEs served as input for a continuous hydrological model that represents the main hydrological processes in the region, namely infiltration, flow routing and transmission losses. The infiltration function is applied in a distributed mode while the routing and transmission loss functions are applied in a lumped mode. Model parameters were found by calibration based on the 5 years of data for one of the catchments. Validation was performed for a subsequent 5-year period for the same catchment and then for an entire 10-year record for the second catchment. The probability of detection and false alarm rates for the validation cases were reasonable. Probabilistic flash-flood prediction is presented applying Monte Carlo simulations with an uncertainty range for the QPEs and model parameters. With low probability thresholds, one can maintain more than 70% detection with no more than 30% false alarms. The study demonstrates that a flash-flood warning model is feasible for catchments in the area studied.     

Impact of Doppler Radar Wind in Simulating the Intensity and Propagation of Rainbands Associated with Mesoscale Convective Complexes Using MM5-3DVAR System

Pre-monsoon rainfall around Kolkata (northeastern part of India) is mostly of convective origin as 80% of the seasonal rainfall is produced by Mesoscale Convective Systems (MCS). Accurate prediction of the intensity and structure of these convective cloud clusters becomes challenging, mostly because the convective clouds within these clusters are short lived and the inaccuracy in the models initial state to represent the mesoscale details of the true atmospheric state. Besides the role in observing the internal structure of the precipitating systems, Doppler Weather Radar (DWR) provides an important data source for mesoscale and microscale weather analysis and forecasting. An attempt has been made to initialize the storm-scale numerical model using retrieved wind fields from single Doppler radar. In the present study, Doppler wind velocities from the Kolkata Doppler weather radar are assimilated into a mesoscale model, MM5 model using the three-dimensional variational data assimilation (3DVAR) system for the prediction of intense convective events that occurred during 0600 UTC on 5 May and 0000 UTC on 7 May, 2005. In order to evaluate the impact of the DWR wind data in simulating these severe storms, three experiments were carried out. The results show that assimilation of Doppler radar wind data has a positive impact on the prediction of intensity, organization and propagation of rain bands associated with these mesoscale convective systems. The assimilation system has to be modified further to incorporate the radar reflectivity data so that simulation of the microphysical and thermodynamic structure of these convective storms can be improved.     

Adjustment of wind-drift effect for real-time systematic error correction in radar rainfall data

An effective bias correction procedure using gauge measurement is a significant step for radar data processing to reduce the systematic error in hydrological applications. In these bias correction methods, the spatial matching of precipitation patterns between radar and gauge networks is an important premise. However, the wind-drift effect on radar measurement induces an inconsistent spatial relationship between radar and gauge measurements as the raindrops observed by radar do not fall vertically to the ground. Consequently, a rain gauge does not correspond to the radar pixel based on the projected location of the radar beam. In this study, we introduce an adjustment method to incorporate the wind-drift effect into a bias correlation scheme. We first simulate the trajectory of raindrops in the air using downscaled three-dimensional wind data from the weather research and forecasting model (WRF) and calculate the final location of raindrops on the ground. The displacement of rainfall is then estimated and a radar–gauge spatial relationship is reconstructed. Based on this, the local real-time biases of the bin-average radar data were estimated for 12 selected events. Then, the reference mean local gauge rainfall, mean local bias, and adjusted radar rainfall calculated with and without consideration of the wind-drift effect are compared for different events and locations. There are considerable differences for three estimators, indicating that wind drift has a considerable impact on the real-time radar bias correction. Based on these facts, we suggest bias correction schemes based on the spatial correlation between radar and gauge measurements should consider the adjustment of the wind-drift effect and the proposed adjustment method is a promising solution to achieve this.     

Spatial reconstruction of rainfall fields from rain gauge and radar data

Rainfall is a phenomenon difficult to model and predict, for the strong spatial and temporal heterogeneity and the presence of many zero values. We deal with hourly rainfall data provided by rain gauges, sparsely distributed on the ground, and radar data available on a fine grid of pixels. Radar data overcome the problem of sparseness of the rain gauge network, but are not reliable for the assessment of rain amounts. In this work we investigate how to calibrate radar measurements via rain gauge data and make spatial predictions for hourly rainfall, by means of Monte Carlo Markov Chain algorithms in a Bayesian hierarchical framework. We use zero-inflated distributions for taking zero-measurements into account. Several models are compared both in terms of data fitting and predictive performances on a set of validation sites. Finally, rainfall fields are reconstructed and standard error estimates at each prediction site are shown via easy-to-read spatial maps.     

Hydrological model sensitivity to parameter and radar rainfall estimation uncertainty

Radar estimates of rainfall are being increasingly applied to flood forecasting applications. Errors are inherent both in the process of estimating rainfall from radar and in the modelling of the rainfall–runoff transformation. The study aims at building a framework for the assessment of uncertainty that is consistent with the limitations of the model and data available and that allows a direct quantitative comparison between model predictions obtained by using radar and raingauge rainfall inputs. The study uses radar data from a mountainous region in northern Italy where complex topography amplifies radar errors due to radar beam occlusion and variability of precipitation with height. These errors, together with other error sources, are adjusted by applying a radar rainfall estimation algorithm. Radar rainfall estimates, adjusted and not, are used as an input to TOPMODEL for flood simulation over the Posina catchment (116 km²). Hydrological model parameter uncertainty is explicitly accounted for by use of the GLUE (Generalized Likelihood Uncertainty Estimation). Statistics are proposed to evaluate both the wideness of the uncertainty limits and the percentage of observations which fall within the uncertainty bounds. Results show the critical importance of proper adjustment of radar estimates and the use of radar estimates as close to ground as possible. Uncertainties affecting runoff predictions from adjusted radar data are close to those obtained by using a dense raingauge network, at least for the lowest radar observations available. Copyright © 2004 John Wiley & Sons, Ltd.     

Time–frequency characterization of sub-divisional scale seasonal rainfall in India using the Hilbert–Huang transform

Time–frequency characterization is useful in understanding the nonlinear and non-stationary signals of the hydro-climatic time series. The traditional Fourier transform, and wavelet transform approaches have certain limitations in analyzing non-linear and non-stationary hydro-climatic series. This paper presents an effective approach based on the Hilbert–Huang transform to investigate time–frequency characteristics, and the changing patterns of sub-divisional rainfall series in India, and explored the possible association of monsoon seasonal rainfall with different global climate oscillations. The proposed approach integrates the complete ensemble empirical mode decomposition with adaptive noise algorithm and normalized Hilbert transform method for analyzing the spectral characteristics of two principal seasonal rainfall series over four meteorological subdivisions namely Assam-Meghalaya, Kerala, Orissa and Telangana subdivisions in India. The Hilbert spectral analysis revealed the dynamic nature of dominant time scales for two principal seasonal rainfall time series. From the trend analysis of instantaneous amplitudes of multiscale components called intrinsic mode functions (IMFs), it is found that both intra and inter decadal modes are responsible for the changes in seasonal rainfall series of different subdivisions and significant changes are noticed in the amplitudes of inter decadal modes of two seasonal rainfalls in the four subdivisions since 1970s. Further, the study investigated the links between monsoon rainfall with the global climate oscillations such as Quasi Bienniel Oscillation (QBO), El Nino Southern Oscillation (ENSO), Sunspot Number (SN), Atlantic Multidecadal Oscillation (AMO) etc. The study noticed that the multiscale components of rainfall series IMF1, IMF2, IMF3, IMF4 and IMF5 have similar periodic structure of QBO, ENSO, SN, tidal forcing and AMO respectively. As per the seasonal rainfall patterns is concerned, the results of the study indicated that for Assam-Meghalaya subdivision, there is a likelihood of extreme rare events at ~0.2 cycles per year, and both monsoon and pre-monsoon rainfall series have decreasing trends; for Kerala subdivision, extreme events can be expected during monsoon season with shorter periodicity (~2.5 years), and monsoon rainfall has statistically significant decreasing trend and post-monsoon rainfall has a statistically significant increasing trend; and for Orissa subdivision, there are chances of extremes rainfall events in monsoon season and a relatively stable rainfall pattern during post-monsoon period, but both monsoon and post-monsoon rainfall series showed an overall decreasing trend; for Telangana subdivision, there is a likelihood of extreme events during monsoon season with a periodicity of ~4 years, but both monsoon and post-monsoon rainfall series showed increasing trends. The results of correlation analysis of IMF components of monsoon rainfall and five climate indices indicated that the association is expressed well only for low frequency modes with similar evolution of trend components.     

Temporal precipitation estimation from nearby radar reflectivity using dynamic factor analysis in the mountainous watershed – a case during Typhoon Morakot

The radar reflectivity (Z)–rain intensity (R) relationship fluctuates in both temporal and spatial scales. The dynamic factor analysis (DFA) and min/max autocorrelation factor analysis (MAFA) was specifically designed for considering various space–time integrations of gauge rainfall and radar reflectivity. We detect representative radar reflectivity observed around rainfall stations that were most responsible for rainfall intensity and identify the crucial patterns of the radar reflectivity in the Kaoping River watershed during Typhoon Morakot. Result shows that the MAFA and DFA can reduce the uncertainty of the dynamic Z‐R relationship effectively. The MAFA separates an entire area into two subareas (southern and northern areas) according to the relationships between the radar reflectivity and min/max autocorrelation factor (MAF) axes. For both areas, the different extents of temporal rainfall correlated with the radar reflectivity were determined using DFA. Especially in the northern area, the radar reflectivity was significantly related to the rainfall intensity for most stations without mountain blockage. Mountain blockages associated with the presence of terrain and wind direction were inferred the major factors that affected the relationship between radar reflectivity and rainfall intensity in the mountainous watershed. Further study can consider the terrain effect and meteorological information, such as wind speed and direction in the DFA model, with the dominant radar reflectivity to estimate the temporal rainfall patterns. Copyright © 2012 John Wiley & Sons, Ltd.     

The isotopic signature of monsoon conditions,cloud modes,and rainfall type

This work provides a comprehensive physically based framework for the interpretation of the north Australian rainfall stable isotope record (δ¹⁸O and δ²H). Until now, interpretations mainly relied on statistical relationships between rainfall amount and isotopic values on monthly timescales. Here, we use multiseason daily rainfall stable isotope and high resolution (10 min) ground‐based C‐band polarimetric radar data and show that the five weather types (monsoon regimes) that constitute the Australian wet season each have a characteristic isotope ratio. The data suggest that this is not only due to changes in regional rainfall amount during these regimes but, more importantly, is due to different rain and cloud types that are associated with the large scale circulation regimes. Negative (positive) isotope anomalies occurred when stratiform rainfall fractions were large (small) and the horizontal extent of raining areas were largest (smallest). Intense, yet isolated, convective conditions were associated with enriched isotope values whereas more depleted isotope values were observed when convection was widespread but less intense. This means that isotopic proxy records may record the frequency of which these typical wet season regimes occur. Positive anomalies in paleoclimatic records are most likely associated with periods where continental convection dominates and convection is sea‐breeze forced. Negative anomalies may be interpreted as periods when the monsoon trough is active, convection is of the oceanic type, less electric, and stratiform areas are wide spread. This connection between variability of rainfall isotope anomalies and the intrinsic properties of convection and its large‐scale environment has important implications for all fields of research that use rainfall stable isotopes.     

A non-parametric automatic blending methodology to estimate rainfall fields from rain gauge and radar data

Quantitative estimation of rainfall fields has been a crucial objective from early studies of the hydrological applications of weather radar. Previous studies have suggested that flow estimations are improved when radar and rain gauge data are combined to estimate input rainfall fields. This paper reports new research carried out in this field. Classical approaches for the selection and fitting of a theoretical correlogram (or semivariogram) model (needed to apply geostatistical estimators) are avoided in this study. Instead, a non-parametric technique based on FFT is used to obtain two-dimensional positive-definite correlograms directly from radar observations, dealing with both the natural anisotropy and the temporal variation of the spatial structure of the rainfall in the estimated fields. Because these correlation maps can be automatically obtained at each time step of a given rainfall event, this technique might easily be used in operational (real-time) applications. This paper describes the development of the non-parametric estimator exploiting the advantages of FFT for the automatic computation of correlograms and provides examples of its application on a case study using six rainfall events. This methodology is applied to three different alternatives to incorporate the radar information (as a secondary variable), and a comparison of performances is provided. In particular, their ability to reproduce in estimated rainfall fields (i) the rain gauge observations (in a cross-validation analysis) and (ii) the spatial patterns of radar fields are analyzed. Results seem to indicate that the methodology of kriging with external drift [KED], in combination with the technique of automatically computing 2-D spatial correlograms, provides merged rainfall fields with good agreement with rain gauges and with the most accurate approach to the spatial tendencies observed in the radar rainfall fields, when compared with other alternatives analyzed.