A new discrete multiplicative random cascade model for downscaling intermittent rainfall fields

Spatial downscaling of rainfall fields is a challenging mathematical problem for which many different types of methods have been proposed. One popular solution consists of redistributing rainfall amounts over smaller and smaller scales by means of a discrete multiplicative random cascade (DMRCs). This works well for slowly varying homogeneous rainfall fields but often fails in the presence of intermittency (i.e., large amounts of zero rainfall values). The most common workaround in this case is to use two separate cascade models, namely one for the occurrence and another for the intensity. In this paper, a new and simpler approach based on the notion of equal-volume areas (EVAs) is proposed. Unlike classical cascades where rainfall amounts are redistributed over grid cells of equal size, the EVA cascade splits grid cells into areas of different sizes, with each of them containing exactly half of the original amount of water. The relative areas of the subgrid cells are determined by drawing random values from a logit-normal cascade generator model with scale and intensity-dependent standard deviation (SD). The process ends when the amount of water in each subgrid cell is smaller than a fixed-bucket capacity, at which point the output of the cascade can be resampled over a regular Cartesian mesh. The present paper describes the implementation of the EVA cascade model and gives some first results for 100 selected events in the Netherlands. Performance is assessed by comparing the outputs of the EVA model to bilinear interpolation and to a classical DMRC model based on fixed grid cell sizes. Results show that, on average, the EVA cascade outperforms the classical method, producing fields with more realistic distributions, small-scale extremes and spatial structures. Improvements are mostly credited to the higher robustness of the EVA model in the presence of intermittency and to the lower variance of its generator. However, both approaches have their advantages and weaknesses. For example, while the classical cascade tends to overestimate small-scale variability and extremes, the EVA model tends to produce fields that are slightly too smooth and block shaped compared to the observations. The complementary nature of the two approaches, and the fact that they produce errors of opposite signs, opens up new possibilities for quality control and bias corrections of downscaled fields.

A new discrete multiplicative random cascade model for downscaling intermittent rainfall fields

Spatial downscaling of rainfall fields is a challenging mathematical problem for which many different types of methods have been proposed. One popular solution consists in redistributing rainfall amounts over smaller and smaller scales by means of a discrete multiplicative random cascade (DMRC). This works well for slowly varying, homogeneous rainfall fields but often fails in the presence of intermittency (i.e., large amounts of zero rainfall values). The most common workaround in this case is to use two separate cascade models, one for the occurrence and another for the intensity. In this paper, a new and simpler approach based on the notion of equal-volume areas (EVAs) is proposed. Unlike classical cascades where rainfall amounts are redistributed over grid cells of equal size, the EVA cascade splits grid cells into areas of different sizes, each of them containing exactly half of the original amount of water. The relative areas of the sub-grid cells are determined by drawing random values from a logit-normal cascade generator model with scale and intensity dependent standard deviation. The process ends when the amount of water in each sub-grid cell is smaller than a fixed bucket capacity, at which point the output of the cascade can be re-sampled over a regular Cartesian mesh. The present paper describes the implementation of the EVA cascade model and gives some first results for 100 selected events in the Netherlands. Performance is assessed by comparing the outputs of the EVA model to bilinear interpolation and to a classical DMRC model based on fixed grid cell sizes. Results show that on average, the EVA cascade outperforms the classical method, producing fields with more realistic distributions, small-scale extremes and spatial structures. Improvements are mostly credited to the higher robustness of the EVA model to the presence of intermittency and to the lower variance of its generator. However, improvements are not systematic and both approaches have their advantages and weaknesses. For example, while the classical cascade tends to overestimate small-scale extremes and variability, the EVA model tends to produce fields that are slightly too smooth and blocky compared with observations.

Multiplicative cascade models for fine spatial downscaling of rainfall: parameterization with rain gauge data

Capturing the spatial distribution of high-intensity rainfall over short-time intervals is critical for accurately assessing the efficacy of urban stormwater drainage systems. In a stochastic simulation framework, one method of generating realistic rainfall fields is by multiplicative random cascade (MRC) models. Estimation of MRC model parameters has typically relied on radar imagery or, less frequently, rainfall fields interpolated from dense rain gauge networks. However, such data are not always available. Furthermore, the literature is lacking estimation procedures for spatially incomplete datasets. Therefore, we proposed a simple method of calibrating an MRC model when only data from a moderately dense network of rain gauges is available, rather than from the full rainfall field. The number of gauges needs only be sufficient to adequately estimate the variance in the ratio of the rain rate at the rain gauges to the areal average rain rate across the entire spatial domain. In our example for Warsaw, Poland, we used 25 gauges over an area of approximately 1600 km2. MRC models calibrated using the proposed method were used to downscale 15-min rainfall rates from a 20 by 20 km area to the scale of the rain gauge capture area. Frequency distributions of observed and simulated 15-min rainfall at the gauge scale were very similar. Moreover, the spatial covariance structure of rainfall rates, as characterized by the semivariogram, was reproduced after allowing the probability density function of the random cascade generator to vary with spatial scale.

Multiplicative cascade models for fine spatial downscaling of rainfall: parameterization with rain gauge data

Capturing the spatial distribution of high-intensity rainfall over short-time intervals is critical for accurately assessing the efficacy of urban stormwater drainage systems. In a stochastic simulation framework, one method of generating realistic rainfall fields is by multiplicative random cascade (MRC) models. Estimation of MRC model parameters has typically relied on radar imagery or, less frequently, rainfall fields interpolated from dense rain gauge networks. However, such data are not always available. Furthermore, the literature is lacking estimation procedures for spatially incomplete datasets. Therefore, we proposed a simple method of calibrating an MRC model when only data from a moderately dense network of rain gauges are available, rather than from the full rainfall field. The number of gauges need only be sufficient to adequately estimate the variance in the ratio of the rain rate at the rain gauges to the areal average rain rate across the entire spatial domain. In our example for Warsaw, Poland, we used 25 gauges over an area of approximately 1600 km2. MRC models calibrated using the proposed method were used to downscale 15-min rainfall rates from a 20 by 20 km area to the scale of the rain gauge capture area. Frequency distributions of observed and simulated 15-min rainfall at the gauge scale were very similar. Moreover, the spatial covariance structure of rainfall rates, as characterized by the semivariogram, was reproduced after allowing the probability density of the random cascade generator to vary with spatial scale.

THE MAXIMUM OF MULTIFRACTAL CASCADES: EXACT DISTRIBUTION AND APPROXIMATIONS

We study the distribution of the maximum M of multifractal measures using discrete cascade representations. For such discrete cascades, the exact distribution of M can be found numerically. We evaluate the sensitivity of the distribution of M to simplifying approximations, including independence of the measure among the cascade tiles and replacement of the dressing factor by a random variable with the same distribution type as the cascade generator. We also examine how the distribution of M varies with the dimensionality of the support and the multiplicity of the cascade. Of these factors, dependence of the measure among different cascade tiles has the highest effect on the distribution of M. This effect comes mainly from long-range dependence. We use these findings to propose a simple approximation to the distribution of M and give charts to implement the approximation for beta-lognormal cascades.