## 1. Introduction

It is well known that nonuniform beamfilling (NUBF) is a major error source in quantitative rainfall remote measurement with a spaceborne rain radar if the footprint size of the radar is comparable to or larger than a convective cell size. The effect of NUBF depends also upon the algorithm for estimating rain rate. It has been shown that the use of the path-integrated attenuation (PIA) derived from the surface return method [surface reference technique (SRT)] suffers more from the NUBF effect than the use of radar reflectivity factor with a *Z–R* relation for rain-rate estimation (Nakamura 1991; Amayenc et al. 1993). However, SRT is recognized as a very important method to achieve a quantitative retrieval of heavy rain (about 20 mm h^{−1} or higher) because it provides a reference to stabilize the rain attenuation correction that often becomes unstable if we know only a measured and attenuated *Z*-factor profile (Meneghini and Kozu 1990). It is therefore crucial to develop a method to judge the existence of significant NUBF (to issue a “warning” flag), and, if possible, to develop a method to correct the SRT-derived path attenuation for the NUBF effect. Since the attenuation coefficient (in dB per unit distance) has an approximately linear relationship with rain rate, obtaining accurate PIA is also useful to estimate path-integrated rain rate.

We have studied a method for correcting the SRT-measured PIA for the NUBF effect (Kozu and Iguchi 1996, 1997) and have shown that the method has a potential to reduce bias errors in PIA estimation based upon case studies using shipborne and ground-based radar data from the Tropical Ocean Global Atmosphere (TOGA) Coupled Ocean–Atmosphere Response Experiment (COARE) (Short et al. 1995). This method is based upon a characteristic of spatial variability of rain intensity, which has been used for the SSM/I rain-rate retrieval algorithm by Kummerow and Giglio (1994). They showed a fairly good correlation between the normalized standard deviations (NSDs) of rain at 1-km resolution and the NSD at 12.5-km resolution, the latter being obtained from the 85-GHz channel of SSM/I, and made an NUBF correction using the estimated 1-km resolution NSD. In the case of radar measurement, however, the spatial scale of interest is approximately equal to the correlation length of individual rain cells, which is much smaller than those of SSM/I. Therefore, we may find characteristics different from the result of Kummerow and Giglio (1994). Graves (1993) also pointed out that statistical properties of rain field tend to be unstable as the size of instantaneous field of view (IFOV) becomes small, making a statistical NUBF correction difficult in his statistical modeling of NUBF for microwave radiometer measurements of oceanic rainfall.

This paper, which is the result of a study following the work of Kozu and Iguchi (1996, 1997), describes comprehensive statistical analyses of shipborne radar data obtained between November 1992 and February 1993 over the tropical Pacific. In the following sections, we describe 1) basic formulation of the problem, 2) the NUBF correction method, 3) statistical characteristics of rainfall, and 4) a simulation of the NUBF correction.

## 2. Rain and observation models

### a. PDF and its parameters

We assume a vertically uniform, horizontally nonuniform rain model in an IFOV as shown in Fig. 1. This model is valid when we consider the SRT because it deals with only the path-integrated values. In this situation, we assume that we can measure a PIA with the SRT, a storm height, and a *Z*-factor profile that includes propagation loss due to rain.

*R,*

*R*

*R, σ*

_{R}(i.e., we assume “uniform” statistics within the IFOV), then the PDF,

*f*(

*R*), is expressed as

*μ*

_{R}and

*ξ*

^{2}

_{R}

*R,*which are related to

*R*

*σ*

_{R}as follows:

*A*

_{SRT}, in decibels, is given by

*A*

*R*

*LaR*

^{b}

*L*is the storm depth, and

*a*and

*b*are the coefficient and exponent of the

*k–R*relation (attenuation coefficient versus rain-rate relation). We define the uniform PIA,

*A*

_{u}, as the PIA that corresponds to

*R*

*k–R*relation that can be obtained when rain is uniform within the IFOV:

*A*

_{u}

*La*

*R*

^{b}

*R*

*σ*

_{R}as parameters, we can relate

*A*

_{SRT}to

*A*

_{u}. The relationship between those quantities assuming the lognormal PDF and a rain depth

*L*of 5 km is shown in Fig. 2.

*k–R*and

*Z–R*relations are generally expressed as power laws,

*k*

*aR*

^{b}

*Z*

*αR*

^{β}

*k*and

*Z*are also lognormal and the means and variances of ln

*k*(

*μ*

_{k},

*ξ*

^{2}

_{k}

*Z*(

*μ*

_{Z},

*ξ*

^{2}

_{Z}

*k*(

*k*

*σ*

_{k}) and those of

*Z*(

*Z*

*σ*

_{Z}) are expressed as

*k*and

*Z*are different from the NSD of rain rate but can easily be converted. Please note that the parameter

*σ*

_{R}in Fig. 2 is the NSD of rain rate (not of attenuation coefficient). In the calculation of

*A*

_{u}and

*A*

_{SRT}, the difference in PDF parameters is incorporated.

### b. Relation between point and area-averaged quantities

Before describing the NUBF correction method, we consider the relationship between a “point” quantity *Q* and an “area-averaged” counterpart 〈*Q*〉 that is a basis of the current method. A similar problem has been discussed in the study of methods to predict PIA from a point measurement of rain rate with a rain gauge (Morita and Higuti 1978). They derived expressions for the mean and variance of PIA using a “normalized spatial autocorrelation function” (hereafter, ACF) of rain rate expressed by the form of exp(−*ζ**r**r* is the distance between two points. In this paper, we derive an expression of *Q* (vertically averaged or integrated quantity with high spatial resolution) and 〈*Q*〉 (area-averaged value of *Q* over the IFOV). Note that the difference between the notation *Q**Q*〉 is that the former is the expectation and the latter is still a statistical random variable.

*ρ*

*r*

*ζr*

^{η}

*r*is the distance between two points on a two-dimensional rain field. As we see later in this paper, this model can well fit the “averaged” 1D ACF that is generated from the observed 2D ACF by making an “iso-range” circular averaging. Of course, there may be a wide range of directional properties in the 2D ACF, depending on storm type and orientation. Nevertheless, the approximation used here would be useful to obtain approximate properties of the relation between NSDs in connection to the ACF of rain field. Let

*Q*(

*x, y*) be the point value at (

*x, y*) and 〈

*Q*(

*x*

_{0},

*y*

_{0})〉 be the spatially averaged

*Q*over a circular area having the diameter

*D*centered at (

*x*

_{0},

*y*

_{0}). In addition, let

*Q*

*σ*

^{′}

_{Q}

*Q*. We add (′) to

*σ*

_{Q}

*Q*, which may be different from the variability of

*Q*within an IFOV. Then we have the following relations:

*ρ*(

*r*) is independent of the location (

*x, y*), {

*Q*〉

^{2}

*Q*

^{2}} can be written as

*P*(

*λ, z*) is the incomplete gamma function, and the double integral over

*θ*

_{1}and

*r*

_{2}is represented by a function

*B*(

*D, ζ, η*). Note that the above 2D integration is made over the circular region having a diameter

*D.*A point

*p*(

*r*

_{2},

*θ*

_{2}) is first specified in the circle. The variable

*θ*

_{1}is a polar angle about

*p*(

*r*

_{2},

*θ*

_{2}), and

*r*

_{1}is defined as a distance from

*p*(

*r*

_{2},

*θ*

_{2}) to the edge of the circle in the direction of

*θ*

_{1}. Thus the upper limit of the integration over

*r*

_{1}(

*R*

_{m}) is a function of

*θ*

_{1}and

*r*

_{2}. Thus, the relation between

*σ*

^{′2}

_{Q}

*σ*

^{2}

_{〈Q〉}

*σ*

^{′2}

_{Q}

*σ*

^{2}

_{〈Q〉}

*B*(

*D, ζ, η*). It should be noted that the PDF form has nothing to do with Eq. (16), but the PDF does affect the correction factor to estimate

*A*

_{u}(see Fig. 2).

The difference between *σ*_{Q}*σ*^{′}_{Q}

## 3. Description of correction method

As we discussed earlier, our final goal is to estimate *A*_{u} from radar-measurable quantities. The first step is to estimate parameters of the rain-rate PDF within an IFOV of interest, specifically *R**σ*_{R}. In this section we describe formulations using measurable quantities but do not include measurement errors to simplify the problem. The mean and NSD in this section represent the sample mean and sample NSD, respectively, which is different from the formulations in section 2, where we considered hypothetical random rain field with parent PDFs. To avoid a confusion, we denote a sample NSD “*σ̂,**σ*” as was used in section 2.

Since we cannot directly measure 〈*R*〉, we cannot apply the direct relations connecting PDF parameters [Eqs. (2), (7)–(12)]; however, a quantity 〈*Q*〉 we can measure such as *A*_{SRT} would serve a proxy of 〈*R*〉 to estimate the NSD of *R.* To obtain an estimate of *σ*〈*Q*〉, we use data at nine IFOVs (eight surrounding IFOVs and the IFOV at which we want to estimate *σ*_{R}), as shown in Fig. 3, which illustrates the geometry of spaceborne radar IFOVs related to this problem. Specifically, we estimate *σ*_{R} from SRT-derived PIA, *A*_{SRT}. Although *A*_{SRT} is not proportional to 〈*k*〉, we will use this as a radar-measurable quantity. In the following, we formulate a discrete-model relationship for actual observation conditions, which can also be used for the simulation we discuss later.

*R*〉

_{i}(

*i*= 1–9) be the IFOV-averaged path-averaged rain rate, which consists of

*n*rain-rate elements (

*R*

_{j,i},

*j*= 1,

*n*) within the IFOV:

*w*

_{j}is the weighting function representing the two-way antenna pattern. The PIA we should obtain when the rain is uniform in the IFOV,

*A*

_{u,i}(dB), is

*A*

_{u,i}

*La*

*R*〉

^{b}

_{i}

*La*

*R*

_{i}

*b*depends upon the frequency and drop size distribution, but it is fairly close to unity for the frequency range of interest (10–35 GHz) (Olsen et al. 1978). The last approximate equality in Eq. (5′) [see also Eq. (5)] represents this fact with a slight modification of the coefficient

*a*to

*a*′ to get the best fit. Therefore, the NSD of PIA is generally close to that of path-averaged rain rate, and the

*A*

_{SRT}–

*A*

_{u}relationship derived for a specific rain depth can be used for different rain depth cases. In other words, the variation of rain depth is not a problem if the rain rate averaged over a given rain depth follows the same type of PDF form as that assumed to generate the

*A*

_{SRT}–

*A*

_{u}relationship. Note that the “given” rain depth should be equal to or higher than the maximum rain depth of interest.

*R*

_{j,i}, the PIA element

*A*

_{j,i}is expressed as

*A*

_{j,i}

*Lk*

_{j}

*La*

*R*

^{b}

_{j,i}

*La*

^{′}

*R*

_{j,i}

*A*

_{SRT}) is given by

*i*th IFOV (

*σ̂*

_{R,i}

*σ*

_{〈Q〉}

*Q*〉:

*Q*〉

*Q*〉

_{1}–〈

*Q*〉

_{9}.

As the radar-measurable quantity, we will use the PIA measured with the SRT (*A*_{SRT}). In the procedure for correcting *A*_{SRT}, the coarse-resolution NSD, *σ̂*_{ASRT}*A*_{SRT}, where *σ̂*_{ASRT}*σ̂*_{〈Q〉}*Q*〉 = *A*_{SRT}. Next, *σ*_{R,5} is estimated from *σ̂*_{ASRT}*σ̂*_{ASRT}*σ̂*_{R}*A*_{u,5} is estimated from *A*_{SRT,5} and *σ*_{R,5} using the relation between *A*_{SRT} and *A*_{u} for a given *σ*_{R} assuming a lognormal PDF within the IFOV. In this case, the bias error in the sample NSD (*σ̂*_{R}*σ̂*_{R}–*σ̂*_{ASRT}

## 4. Statistical characteristics of rainfall over tropical Pacific

### a. Data source

To study the rainfall characteristics statistically, we use 42 constant-altitude plan position indicator (CAPPI) maps at 2-km altitude obtained from the Massachusetts Institute of Technology (MIT) radar on board the R/V *Vickers* during the TOGA COARE campaign from November 1992 to February 1993. Those rain scenes are listed in Table 1 with a “subjective” rain-type classification. Though subjective, the classification may be useful to get approximate information about the rain-type dependence of the statistical characteristics. Examples of CAPPI data for types 1–4 are shown in Fig. 4. Please note that the classification is based on the characteristics of the area having a radius of 40 km centered at the radar location. The parameters to process the CAPPI *Z*-factor maps are listed in Table 2.

### b. Spatial correlation properties of rain rate

As we discussed in section 3, the relation of variabilities between “point” rainfall quantity (*Q*) and area-averaged (〈*Q*〉) quantity can be estimated if we know the parameters of the ACF of *Q*, and the averaging area (represented by *D*). This is approximately true when we relate *R* and 〈*Q*〉, where *Q* is not necessarily equal to *R.* In other words, characterizing the ACF of *R* would be useful to estimate the *σ*^{′}_{R}*σ*_{〈Q〉}*ζ* and *η.* On average, *ζ* ≈ 0.135 and *η* ≈ 1.13 for convective types (types 1, 2, and 4) and *ζ* ≈ 0.064 and *η* ≈ 1.10 for the stratiform type (type 3). Using Eq. (16) with *D* = 4 km, the relation between *σ*_{〈R〉} and *σ*^{′}_{R}

*σ*

^{′}

_{R}

*σ*

_{R}. A theoretically rigorous treatment of this problem is beyond the scope of this paper; in the following, we consider this problem intuitively. Random rain-rate fluctuations, in which wavelengths are less than 4 km, clearly contribute to the variability of rain rate in a 4-km IFOV. The contribution of a fluctuation with spatial wavelength (

*λ*) to

*σ*

_{R}decreases when

*λ*becomes longer. A simple simulation evaluating the reduction factor versus

*λ*has indicated that the contribution is reduced to about 70% when

*λ*= 8 km. Thus, we assume the following approximate relation:

*σ*

_{R}

*σ*

^{′2}

_{R}

*σ*

^{2}

_{〈R〉(D=8km)}

^{1/2}

*σ*

^{2}

_{〈R〉(D=8km)}

*σ*

^{2}

_{〈R〉}

*c*

_{〈R〉}relating

*σ*

_{〈R〉}to

*σ*

_{R},

*σ*

_{R}

*c*

_{〈R〉}

*σ*

_{〈R〉}

*c*

_{〈R〉}is about 0.72. For

*ζ*≈ 0.135 and

*η*≈ 1.13 (mean values for convective rains),

*c*

_{〈R〉}is about 0.66. Equation (23) can be generalized as follows to estimate

*σ*

_{R}from any 〈

*Q*〉:

*σ*

_{R}

*c*

_{〈Q〉}

*σ*

_{〈Q〉}

*σ̂*

_{R}

*c*

_{〈Q〉}

*σ̂*

_{〈Q〉}

### c. Relation between finescale and coarse-scale rainfall variabilities

The relation Eq. (23a) can also be obtained directly from actual radar data by a regression analysis, although we have to notice that *c*_{〈Q〉}*σ̂*_{R}*σ̂*_{〈Q〉}*A*_{SRT} and 〈*R*〉; the former is used to obtain a statistical relation to estimate *σ*_{R} from the downlooking radar measurement, and the latter to compare with the theoretical calculation discussed above. The difference between 〈*R*〉 and *A*_{SRT} will be evaluated in section 4e. Figure 8 shows histograms of the coefficient *c*_{ASRT}, which corresponds to *c*_{〈Q〉}*σ̂*_{R}*σ̂*_{ASRT}*r*_{ASRT} for the convective (types 1, 2, and 4) and stratiform (type 3) cases. For the convective cases, histograms are given for two cases: (case 1) the minimum rain rate of 0.1 mm h^{−1} and (case 2) the 5 mm h^{−1} rain rate to investigate the dependence of the regression result on rain rate since the NUBF correction is more important for higher rain rates. Note that the correlation coefficient is about 0.5 and slightly higher at the stratiform case and that in the convective case the coefficient *c*_{ASRT} is generally greater than it is in the stratiform case. As for the rain-rate dependence, case 2 gives slightly lower and more stable *c*_{ASRT} values than case 1.

Let us compare the result of the regression analysis between *σ̂*_{R}*σ̂*_{〈R〉}*c*_{〈R〉} for 32 convective scenes. The mean value of 0.77 (case 2) compares fairly well with the calculation result (about 0.72) based on the ACF analysis shown in Fig. 7 with *ζ* = 0.135 and *η* = 1.13. On the other hand, the result of case 1 (mean = 0.96) is significantly higher than the calculation result. Although the sample NSDs should underestimate the NSD of the parent PDF by the factor of *N* − 1)/*N*)]*N* is the number of data, the effect of this bias error, which causes a slight positive bias, would be rather small (about 2%) considering that both *σ̂*_{R}*σ̂*_{〈R〉}*σ̂*_{R}*σ̂*_{〈R〉}*c*_{〈R〉} obtained from the *σ̂*_{R}*σ̂*_{〈R〉}*r*_{c} derived from ACF analyses is shown in Fig. 10, where *r*_{c} is defined as the distance at which the ACF is reduced to *e*^{−1}. It is found that there is a clear correlation between them, suggesting the possibility to relate the the ACF properties to the *σ*_{R}–*σ*_{ 〈R〉} relationship.

### d. Lognormality of rain-rate PDF

To validate the assumption of lognormal distribution of rain rate within an FOV, the lognormality of rain rate for the spatial domain of several kilometers squares is examined by a *χ*^{2} test. For this test, 49 rain-rate samples over a domain of 7 km by 7 km square are used, in which the *χ*^{2} value is calculated only when at least 45 samples out of 49 are recognized as “rainy.” An example of the results is shown in Fig. 11, where *χ*^{2} values are expressed as a cumulative distribution. In the distribution, *χ*^{2} greater than about 10 is recognized as nonlognormal with the level of significance of 0.05 since we used an eight-category histogram to calculate the *χ*^{2} value that should follow the *χ*^{2} distribution for five degrees of freedom. Note that about 50% of the samples do not follow a lognormal distribution. We should note that a no-rain IFOV is not used to calculate a *χ*^{2} value since there is an ambiguity in the treatment of “zero” rain rate. Figure 12 shows a summary of lognormality check of 42 rain scenes; the histogram of percentage of rain areas judged nonlognormal with 5% of significance from the *χ*^{2} test. It is found that about 50% of the rain area may not be expressed as a lognormally distributed random field.

Including the zero rain rate in the *χ*^{2} test should increases the *χ*^{2} value, that is, nonlognormality. This suggests that the nonlognormality may be significant at storm edges and small isolated convective cells. In such cases, significant no-rain areas may exist within an IFOV, and we need to assume a delta-function PDF for the no-rain areas and another PDF for rainy areas. Such a mixed model was used in the NUBF study for microwave radiometer measurements (Graves 1993), but in the case of the spaceborne radar measurement, a model including a nonrainy region in an IFOV would be more difficult to apply because of the saturation of *A*_{SRT} due to the dominance of echo power from the nonrainy region (Kozu and Iguchi 1996). It is a next-step study subject to find alternative PDF that can better describe the rain-rate variability, to develop a method to detect a mixed rainy/nonrainy IFOV, and to apply a special procedure to such a mixed IFOV.

In spite of the above problems, we will assume the lognormal PDF in the simulation (described in the next section) to relate *A*_{SRT} to *A*_{u} because of its simplicity.

### e. Radar measurable quantity other than A_{SRT}

*A*

_{SRT}as the measurable quantity 〈

*Q*〉. In cases where either rain attenuation is too small or normalized radar cross section (

*σ*

^{0}) of the surface varies considerably, we may not be able to use

*A*

_{SRT}. There may be other quantities that can be used to estimate PIA, however. One possibility is the PIA obtained from the sum of attenuations at each range gate estimated from the profile of

*Z*factor including attenuation (

*Z*

_{m}) with a

*k–Z*relation (

*A*

_{SU}or

*A*

_{ST}). The quantity

*Z*

_{m}(

*r*) is given by

*r*

_{0}is the range at the storm top, and

*γ*and

*δ*are the coefficient and the exponent of the

*k–Z*relation. PIA estimates (

*A*

_{SU}and

*A*

_{ST}) are expressed as follows:

*A*

_{SU}and

*A*

_{ST}is that the former uses

*Z*

_{m}factors at high altitudes only, so that it may better be used for relatively intense rain cases than the latter, although it does not use rain information at rain bottom. Since

*A*

_{SU}and

*A*

_{ST}depend on the vertical structure of rain rate, the values obtained from the vertically uniform rain rate model assumed in this paper would not be the same as those from natural rain. Nevertheless, investigating the properties of these quantities under the vertically uniform rain would be useful for obtaining an approximate idea about the feasibility of using the

*Z*

_{m}profile alone. To test if we can use any quantity as 〈

*Q*〉, depending on the rain and surface characteristics, correlations between

*c*

_{〈Q〉}

*R*〉 and those obtained with the use of

*A*

_{SRT},

*A*

_{SU}, or

*A*

_{ST}are plotted in Fig. 13, suggesting the validity of using

*A*

_{SRT},

*A*

_{SU}, and

*A*

_{ST}in place of 〈

*R*〉.

*ε*is a parameter to be adjusted by using

*A*

_{SRT}or by assuming the

*Z*-factor profile is vertically constant just above the surface. Although we do not test this type of PIA estimate, we can expect that

*ξ*is similar to the quantities we have tested here since

*ξ*is the same as

*A*

_{SRT}when SRT is reliable, and it should be generally a better estimate of PIA than

*A*

_{SU}or

*A*

_{ST}.

## 5. Simulation of the correction method

### a. Simulation procedure

For a test of the NUBF correction method, simulations are performed using the same radar dataset as those used in the above analyses. The procedure of the simulation is depicted in Fig. 14. For the preparation of the simulation, we generate a lookup table (actually a set of regression lines of the result shown in Fig. 2, see Table 3) to obtain *A*_{u} from *A*_{SRT} and *σ*_{R}.

We first convert the 1-km spatial resolution *Z*-factor map to the corresponding rain rate and PIA. Next, we calculate the coarse-resolution rain rate and PIA of 4-km resolution by using Eqs. (17), (5′), and (19). At the same time, we calculate *σ̂*_{R}*σ̂*_{ASRT}*σ̂*_{R}*σ̂*_{ASRT}*c*_{ASRT} for convective rain, case 2 (=0.723, see Fig. 8) is used for all types of rainfall.

With the above preparation, the simulation starts with an examination of the coarse-resolution radar measurement (*A*_{SRT}) of each IFOV. For a specific IFOV, we obtain *A*_{SRT} and *σ̂*_{ASRT}*σ̂*_{ASRT}*σ̂*_{R}*A*_{SRT} and the estimated *σ̂*_{R}*σ*_{R}, are then input to the lookup table to obtain an estimate of *A*_{u} (*Ã*_{u}).

Since we use the same dataset for the simulation as that used to derive the *σ̂*_{R}*σ̂*_{ASRT}*c*_{ASRT} to estimate *σ̂*_{R}*σ̂*_{ASRT}*A*_{u} from *A*_{SRT} and the estimated *σ*_{R} is a highly nonlinear, error-sensitive process. In addition, the validity of assuming the lognormal PDF can be tested with this simulation.

### b. Simulation result

Following the procedure described above, simulations of the NUBF correction are performed; the results are shown in Figs. 15–18. In the figures, three scattergrams are compared for each rain type; one is the correlation between *A*_{u} and *A*_{SRT} (without correction), the second is that between *A*_{u} and *Ã*_{u} (estimated *A*_{u} with the NUBF correction) using the estimated *σ*_{R}, and the third is the same as the second except that the “true” *σ*_{R} (strictly speaking, true *σ̂*_{R}*σ*_{R} significantly improves the tendency of overcorrection, suggesting the errors in estimating *σ*_{R}. On the other hand, in type 2, improvement by using the true *σ*_{R} is not significant, suggesting an improper PDF model causing an erroneous *A*_{SRT}–*A*_{u} relation. In the cases of types 3 and 4, the result of NUBF correction is better than types 1 and 2; in particular, the NUBF effect is very small in type 3, leading to stable NUBF corrections. From the above results it is clear that the “instantaneous” or “IFOV-based” NUBF correction can have considerable errors. Nevertheless, it can be concluded that the correction scheme tested here does work in general to reduce the significant negative bias errors due to the NUBF. Note that Durden et al. (1998) also obtained an encouraging result in testing an NUBF correction scheme using the data acquired by a KU-band airborne rain radar flown in the TOGA COARE field campaign.

The effectiveness of the NUBF correction can be seen more clearly from statistical comparisons. Figure 19 shows the ratio of cumulative probabilities of PIA for types 2 and 4. Each cumulative probability represents the percentage of occurrence of PIA that exceeds the given PIA value. In the figure, *A*_{SRT}/*A*_{u} and *Ã*_{u}/*A*_{u} are plotted. Note that the NUBF correction can significantly improve the estimation accuracy of the cumulative distribution of PIA.

## 6. Conclusions

We have studied a method to make a correction of the path-integrated attenuation derived from spaceborne radar measurement for nonuniform beamfilling. The key to this method is the estimation of the finescale rain-rate variability within an IFOV (normalized standard deviation; *σ*_{R}) from the local statistics of a radar-measurable quantity (〈*Q*〉), such as PIA derived from the surface reference technique. The estimated *σ*_{R} is then used to obtain a correction factor to estimate a rainfall quantity such as rain rate and PIA, which should be obtained when rain rate is uniform in the IFOV. In this paper, we have focused our attention to use SRT-derived PIA (*A*_{SRT}) to estimate the uniform PIA (*A*_{u}).

Statistical analyses have been made on spatial variabilities of these radar-measurable quantities using a shipborne radar dataset obtained from the TOGA COARE field campaign. The analyses include the spatial autocorrelation function, the finescale and coarse-scale rainfall variabilities, and lognormality of a rain-rate PDF. There are reasonably good correlations between the coarse-scale variability of *A*_{SRT} and the finescale variability in an IFOV, and the regression coefficient (slope) of these two quantities is fairly stable. Based on the statistical analyses, the method is tested with a simulation using the same dataset. The test result indicates that this method is effective in reducing bias errors in the estimation of rain-rate statistics. It is also effective in making the NUBF correction for individual IFOV basis; in this case, however, one must account for the variability of local rainfall statistical characteristics, which may cause significant errors in estimating *A*_{u}. In the actual implementation of this method, therefore, some limitation on the magnitude of the correction factor (to retrieve *A*_{u} from *A*_{SRT}) or the magnitude of the estimated finescale variability (*σ*_{R}) may be needed to avoid a significant overcorrection for the NUBF.

## Acknowledgments

We wish to thank Drs. C. Kummerow and D. A. Short of NASA/Goddard Space Flight Center (GSFC), Prof. K. Shimizu of Keio University, and Dr. N. Kashiwagi of the Institute of Statistical Mathematics for their valuable suggestions. Acknowledgment is also given to Mr. R. Okada, Remote Sensing Technology Center of Japan, who contributed to many of the statistical data analyses. The MIT radar data were kindly provided from the TOGA COARE project and TRMM Office, NASA/GSFC. This study was partly supported by National Space Development Agency of Japan (NASDA) under a joint research program between CRL and NASDA for the study of TRMM Precipitation Radar algorithms.

## REFERENCES

Amayenc, P., M. Marzoug, and J. Testud, 1993: Analysis of cross-beam resolution effects in rainfall rate profile retrieval from a spaceborne radar.

*IEEE Trans.,***GE-31,**417–425.Durden, S. L. Z. S. Haddad, A. Kitiyakara, and F. K. Li, 1998: Effects of nonuniform beam filling on rainfall retrieval for the TRMM precipitation radar.

*J. Atmos. Oceanic Technol.,***15,**635–646.Graves, C. E., 1993: A model for the beam-filling effect associated with the microwave retrieval of rain.

*J. Atmos. Oceanic Technol.,***10,**5–14.Iguchi, T., T. Kozu, R. Meneghini, and K. Okamoto, 1997: Rain profiling algorithm for the TRMM precipitation radar.

*Proc. International Geoscience and Remote Sensing Symp. (IGARSS’97),*Singapore, 1636–1638.Kozu, T., and T. Iguchi, 1996: A preliminary study of non-uniform beam filling correction for spaceborne rain radar measurement.

*IEICE Trans. Comm.,***E79-B,**763–769.——, and ——, 1997: Correction: A preliminary study of nonuniform beam filling correction for spaceborne rain radar measurement.

*IEICE Trans. Comm.,***E80-B,**989.Kummerow, C., and L. Giglio, 1994: A passive microwave technique for estimating rainfall and vertical structure information from space. Part I: Algorithm description.

*J. Appl. Meteor.,***33,**3–18.Meneghini, R., and T. Kozu, 1990:

*Spaceborne Weather Radar.*Artech House, 199 pp.Morita, K., and I. Higuti, 1978: Statistical studies on rain attenuation and site diversity effect on Earth to satellite links in microwave and millimeter wavebands.

*Trans. IEICE, Japan,***E61,**425–432.Nakamura, K., 1991: Biases of rain retrieval algorithms for spaceborne radar caused by nonuniformity of rain.

*J. Atmos. Oceanic Technol.,***8,**363–373.Olsen, R. L., D. V. Rogers, and D. B. Hodge, 1978: The

*aR**b*relationship in the calculation of rain attenuation.*IEEE Trans.,***AP-26,**318–329.Short, D. A., P. A. Kucera, B. S. Ferrier, and O. W. Thiele, 1995: COARE IOP rainfall from shipborne radars: 1. Rain mapping algorithms. Preprints,

*27th Conf. on Radar Meteorology,*Vail, CO, Amer. Meteor. Soc., 678–690.

Relationship between “uniform” PIA (*A*_{u}) and SRT-derived PIA (*A*_{SRT}) assuming the lognormal PDF and a rain depth *L* of 5 km as the normalized standard deviation of rain rate (*σ*_{R}) as a parameter.

Citation: Journal of Atmospheric and Oceanic Technology 16, 11; 10.1175/1520-0426(1999)016<1722:NBCFSR>2.0.CO;2

Relationship between “uniform” PIA (*A*_{u}) and SRT-derived PIA (*A*_{SRT}) assuming the lognormal PDF and a rain depth *L* of 5 km as the normalized standard deviation of rain rate (*σ*_{R}) as a parameter.

Citation: Journal of Atmospheric and Oceanic Technology 16, 11; 10.1175/1520-0426(1999)016<1722:NBCFSR>2.0.CO;2

Relationship between “uniform” PIA (*A*_{u}) and SRT-derived PIA (*A*_{SRT}) assuming the lognormal PDF and a rain depth *L* of 5 km as the normalized standard deviation of rain rate (*σ*_{R}) as a parameter.

Citation: Journal of Atmospheric and Oceanic Technology 16, 11; 10.1175/1520-0426(1999)016<1722:NBCFSR>2.0.CO;2

Concept of down-looking rain radar IFOV used in the NUBF correction scheme.

Concept of down-looking rain radar IFOV used in the NUBF correction scheme.

Concept of down-looking rain radar IFOV used in the NUBF correction scheme.

Example of MIT radar 2-km CAPPI images (radius = 150 km). Note that only the 80 km by 80 km square area centered at the radar site is used for the NUBF study.

Example of MIT radar 2-km CAPPI images (radius = 150 km). Note that only the 80 km by 80 km square area centered at the radar site is used for the NUBF study.

Example of MIT radar 2-km CAPPI images (radius = 150 km). Note that only the 80 km by 80 km square area centered at the radar site is used for the NUBF study.

Example of the averaged 1D ACF of rain rate.

Example of the averaged 1D ACF of rain rate.

Example of the averaged 1D ACF of rain rate.

Histograms of parameters of the averaged 1D ACF *ρ*(*r*) = exp(−*ζr*^{η}) for rain type 1. The parameters are obtained for the left and right sides of each rain scene.

Histograms of parameters of the averaged 1D ACF *ρ*(*r*) = exp(−*ζr*^{η}) for rain type 1. The parameters are obtained for the left and right sides of each rain scene.

Histograms of parameters of the averaged 1D ACF *ρ*(*r*) = exp(−*ζr*^{η}) for rain type 1. The parameters are obtained for the left and right sides of each rain scene.

Dependence of the coefficient *c*_{〈R〉} in *σ*_{R}–*σ*_{〈R〉} relation on the ACF parameters *ζ* and *η.*

Dependence of the coefficient *c*_{〈R〉} in *σ*_{R}–*σ*_{〈R〉} relation on the ACF parameters *ζ* and *η.*

Dependence of the coefficient *c*_{〈R〉} in *σ*_{R}–*σ*_{〈R〉} relation on the ACF parameters *ζ* and *η.*

(a) Histograms of correlation coefficient *r*_{ASRT} and the coefficient *c*_{ASRT} in *σ̂*_{R}*σ̂*_{ASRT}^{−1} (case 1) and 5 mm h^{−1} (case 2). (b) The same as (a) except for type 3.

(a) Histograms of correlation coefficient *r*_{ASRT} and the coefficient *c*_{ASRT} in *σ̂*_{R}*σ̂*_{ASRT}^{−1} (case 1) and 5 mm h^{−1} (case 2). (b) The same as (a) except for type 3.

(a) Histograms of correlation coefficient *r*_{ASRT} and the coefficient *c*_{ASRT} in *σ̂*_{R}*σ̂*_{ASRT}^{−1} (case 1) and 5 mm h^{−1} (case 2). (b) The same as (a) except for type 3.

Histogram of *c*_{〈R〉} for 32 convective scenes. Upper and lower panels show case 1 and case 2, respectively.

Histogram of *c*_{〈R〉} for 32 convective scenes. Upper and lower panels show case 1 and case 2, respectively.

Histogram of *c*_{〈R〉} for 32 convective scenes. Upper and lower panels show case 1 and case 2, respectively.

Scattergram between *C*_{〈R〉} obtained from the *σ̂*_{R}*σ̂*_{〈R〉}*r*_{c} derived from the ACF analysis.

Scattergram between *C*_{〈R〉} obtained from the *σ̂*_{R}*σ̂*_{〈R〉}*r*_{c} derived from the ACF analysis.

Scattergram between *C*_{〈R〉} obtained from the *σ̂*_{R}*σ̂*_{〈R〉}*r*_{c} derived from the ACF analysis.

An example of the *χ*^{2} test results, where *χ*^{2} values are expressed as a cumulative distribution.

An example of the *χ*^{2} test results, where *χ*^{2} values are expressed as a cumulative distribution.

An example of the *χ*^{2} test results, where *χ*^{2} values are expressed as a cumulative distribution.

A summary of lognormality check of 42 rain scenes; the histogram of percentage of rain areas judged nonlognormal with 5% of significance from the *χ*^{2} test.

A summary of lognormality check of 42 rain scenes; the histogram of percentage of rain areas judged nonlognormal with 5% of significance from the *χ*^{2} test.

A summary of lognormality check of 42 rain scenes; the histogram of percentage of rain areas judged nonlognormal with 5% of significance from the *χ*^{2} test.

Correlations between *c*_{〈Q〉}*R*〉 and those obtained with the use of *A*_{SRT}, *A*_{SU}, or *A*_{ST}.

Correlations between *c*_{〈Q〉}*R*〉 and those obtained with the use of *A*_{SRT}, *A*_{SU}, or *A*_{ST}.

Correlations between *c*_{〈Q〉}*R*〉 and those obtained with the use of *A*_{SRT}, *A*_{SU}, or *A*_{ST}.

Procedure of the simulation.

Procedure of the simulation.

Procedure of the simulation.

Scattergrams (a) between *A*_{u} and *A*_{SRT} (without correction), (b) between *A*_{u} and *Ã*_{u} (with NUBF correction) using the estimated *σ*_{R}, and (c) between *A*_{u} and *Ã*_{u} (with NUBF correction) using the“true” *σ*_{R} for type 1 rainfall. Solid and dotted lines represent the linear regression result and the 1:1 correspondence line, respectively.

Scattergrams (a) between *A*_{u} and *A*_{SRT} (without correction), (b) between *A*_{u} and *Ã*_{u} (with NUBF correction) using the estimated *σ*_{R}, and (c) between *A*_{u} and *Ã*_{u} (with NUBF correction) using the“true” *σ*_{R} for type 1 rainfall. Solid and dotted lines represent the linear regression result and the 1:1 correspondence line, respectively.

Scattergrams (a) between *A*_{u} and *A*_{SRT} (without correction), (b) between *A*_{u} and *Ã*_{u} (with NUBF correction) using the estimated *σ*_{R}, and (c) between *A*_{u} and *Ã*_{u} (with NUBF correction) using the“true” *σ*_{R} for type 1 rainfall. Solid and dotted lines represent the linear regression result and the 1:1 correspondence line, respectively.

The same as Fig. 15 except for type 2.

The same as Fig. 15 except for type 2.

The same as Fig. 15 except for type 2.

The same as Fig. 15 except for type 3.

The same as Fig. 15 except for type 3.

The same as Fig. 15 except for type 3.

The same as Fig. 15 except for type 4.

The same as Fig. 15 except for type 4.

The same as Fig. 15 except for type 4.

Ratio of cumulative probabilities of PIA for type 2 and type 4.

Ratio of cumulative probabilities of PIA for type 2 and type 4.

Ratio of cumulative probabilities of PIA for type 2 and type 4.

Summary of rain scenes used for the current study. For example, 921112__0301 represents 0301 UTC 12 Nov 1992.

Parameters for the analysis of MIT radar data.

Coefficients *a*_{0}, *a*_{1}, and *a*_{2} in log_{10}*A _{u}* =

*a*

_{0}+

*a*

_{1}(log

_{10}

*A*

_{SRT}) +

*a*

_{2}(log

_{10}

*A*

_{SRT})

^{2}.