Evaluating the Effects of Raindrop Motion on the Accuracy of the Precipitation Inversion Algorithm by X-SAR

: Precipitation has a profound impact on both human life and the natural environment. X-band synthetic aperture radar (X-SAR) utilizes high-resolution microwave remote-sensing technology, providing opportunities for global precipitation measurements. The current precipitation inversion algorithms from X-SAR measurements assume that precipitation particles remain relatively stationary with the ground. However, the motion of raindrops could potentially reduce the accuracy of these algorithms. In this study, we first established a functional relationship between raindrop motion and SAR echoes based on the standard deviation of the raindrop Doppler velocity spectrum. Secondly, an exploratory algorithm was proposed to retrieve rainfall distribution under the raindrop motion error model (RMM) and quantitatively calculate the precipitation inversion error caused by raindrop motion. In comparison to conditions where the atmosphere is stationary, when the standard deviation of the Doppler velocity spectrum of raindrops is 1.1 m/s, the relative error of the retrieved surface rain rate increases from 2.1% to 35.8%. Numerical simulations show that SAR echoes are sensitive to changes in the standard deviation of the Doppler velocity spectrum, and the impact of raindrop motion on the accuracy of X-SAR precipitation measurements cannot be neglected.


Introduction
In the atmospheric, hydrological, oceanic, and environmental sciences, the study of rainfall is a very important component.Precipitation consists of a collection of liquid or solid water particles, known as hydrometeors, that begin within a cloud.The measurement of global precipitation is a major goal in environmental climate change research [1].International satellite missions, such as the tropical rainfall measuring mission (TRMM) [2] and the global precipitation mission (GPM) [3], enhance the advancements in satellite-based global precipitation measurements.
The X-band spaceborne synthetic aperture radar (X-SAR) satellite TerraSAR-X was successfully launched in 2007.It can acquire high-resolution images of any desired area, regardless of weather conditions, with a resolution of up to 1 m [4].With the rapid development of numerous satellite platforms equipped with high spatial resolution SAR technology, the application of SAR high-resolution ground data in various fields, such as water body detection, optical image enhancement, target detection, and terrain recognition, is becoming increasingly widespread [5][6][7][8].The traditional methods of estimating rainfall using rain gauges and land-based radars do not cover the entire world, particularly ocean areas [9].Satellite remote-sensing technology, however, can provide more comprehensive precipitation distribution information with its broad coverage and high-frequency measurements.The numerous inherent advantages of SAR provide the possibility for more convenient, nearly global, continuous, and high-precision precipitation inversion [10][11][12].Consequently, exploring how to utilize SAR for global precipitation monitoring holds Atmosphere 2024, 15, 265 2 of 20 significant implications for complementing the theory of precipitation measurements and meeting practical needs.
Research shows that the X band (10 GHz) near the Ku band (14 GHz), used by precipitation radar (PR), is more sensitive to rainfall effects than the L and C bands [13,14].Moore et al. (1997) observed SAR images and found that rain clouds were evident in both C-band and X-band images.For the first time, they used echo signals and image shadows to determine precipitation, showcasing the potential of spaceborne X-SAR for precipitation inversion [14].Pichugin et al. (1991) assumed a uniform distribution of precipitation in the vertical direction under a single-layer rainfall model [15].By introducing a particular function, they cleverly converted the non-linear expression for rainfall attenuation into a linear Volterra integral equation (VIE) of the second kind.This transformation made the task of rainfall retrieval straightforward.The VIE algorithm achieves high inversion accuracy but cannot be directly applied to models containing two or more meteorological particles.When this method is applied to vertically non-uniform models, the computational complexity of the analytical solution significantly increases.Marzano et al. (2008) proposed the model-oriented statistical (MOS) inversion algorithm [16].This method can retrieve precipitation distribution for a double-layer rainfall model with a completely separated snowfall layer and rainfall layer.However, it requires a large number of empirical parameters.MOS has high accuracy in rainstorms but performs poorly in light to moderate rain.Weinman et al. (2009) used precipitation data from Bangladesh and the Amazon region for inversion.The relative errors of the MOS and VIE algorithms are both within 20% [17], which makes them applicable for rainfall inversion in real-life scenarios.Zheng et al. (2020) proposed a method for detecting rainfall based on the correlation characteristics of the sea clutter and the difference in the correlation coefficients between sea wave images with and without rain.This method can detect rainfall from X-SAR images with uneven precipitation distribution [18].Zhao et al. (2021) used textural information to detect and segment rainfall on sea surface SAR images [19].
However, there are still several issues with spaceborne X-SAR precipitation measurement technology, one of which is the potential error caused by neglecting the movement of raindrops.Numerous studies show links between wind fields, raindrop motion, and SAR images.Draper and others conducted experiments that provided evidence of rainfall causing a significant attenuation in the C-band SAR backscattering cross section of the sea surface wind [20].Satake et al. analyzed the polarimetric data of rainfall obtained from an airborne X-SAR and explored the relationship between raindrop motion and SAR images [21].Yu found that heavy rain could reduce the accuracy of wind speed measurement by SAR up to 30%, especially under typhoon conditions, becoming more complex with increasing rainfall intensity [22].Studies focusing on rain effects on wind field retrieval have developed methods to recognize heavy rain-affected areas and reconstruct SAR-retrieved wind speed in those areas [23,24].However, many interesting rain-related phenomena revealed by SAR images are still not fully understood due to poor theoretical modeling of the rain-wind-wave interactions [25].
Our knowledge of the effects of wind-induced raindrop motion on SAR data limits the current precipitation retrieval methods.When SAR detects distributed targets such as precipitation, it assumes that the precipitation particles remain relatively stationary with the ground.In the actual rainfall environment, raindrop motion can impact the inversion model and introduce errors in the inversion of precipitation intensity.However, there is currently no publicly available quantitative analysis of the error caused by raindrop motion in SAR precipitation retrieval in the existing literature.The research on SAR precipitation measurements needs to include the content of raindrop motion error.To address the above problems, this paper establishes a raindrop motion error model and studies the impact of the standard deviation of the Doppler velocity spectrum of precipitation particles on SAR backscattering echoes.Subsequently, we propose an algorithm based on the raindrop motion error model (RMM).Combined with the MOS inversion method, the RMM algorithm can quantify the wind-induced raindrop motion's impact on precipitation retrieval accuracy.

SAR Rain Measurement Principles
The radar cross section (RCS) in m 2 is used to describe the backscattering characteristics of a target [26]: where the dimension of σ is the area, θ is the slant off-nadir angle, and φ is the azimuth angle, respectively.S s is the backscattering power density, and S i is the incident power density in the (θ, φ) direction.
The three-dimensional side view of precipitation detection by SAR is shown in Figure 1.The x-axis is the cross-track direction, the y-axis is the along-track direction, and the z-axis is the altitude direction.
Atmosphere 2024, 15, x FOR PEER REVIEW 3 of 22 raindrop motion error model (RMM).Combined with the MOS inversion method, the RMM algorithm can quantify the wind-induced raindrop motion's impact on precipitation retrieval accuracy.

SAR Rain Measurement Principles
The radar cross section (RCS) in m is used to describe the backscattering characteristics of a target [26]: where the dimension of  is the area,  is the slant off-nadir angle, and  is the azimuth angle, respectively. is the backscattering power density, and  is the incident power density in the ,  direction.The three-dimensional side view of precipitation detection by SAR is shown in Figure 1.The x-axis is the cross-track direction, the y-axis is the along-track direction, and the zaxis is the altitude direction.Assuming that the backscattering echoes of precipitation particles within a radar beam are incoherent, the energy of the total backscattering echoes is equal to the sum of the echo energy of each particle.When SAR detects precipitation in strip scan mode, the power of the backscattering echoes is as follows [26]: where  is the radar received power,  is the radar transmitted power,  is the radar antenna gain,  is the operating wavelength,  is the distance between the measured target and the radar,  is the total number of precipitation particles within the effective illuminated volume of the radar beam that simultaneously produce scattering echo energy and return to the radar antenna, and  is the backscattering cross section of an individual precipitation particle.
The radar volume reflectivity of precipitation is defined as the sum of the backscattering cross section of all precipitation particles per unit volume, denoted as  [27]: The volume resolution cell represents the minimum separable precipitation unit that SAR can distinguish during motion.It can be considered as the differential element of the effective illuminated volume.The volume of the cell is related to the three-dimensional resolutions of SAR [28]: Assuming that the backscattering echoes of precipitation particles within a radar beam are incoherent, the energy of the total backscattering echoes is equal to the sum of the echo energy of each particle.When SAR detects precipitation in strip scan mode, the power of the backscattering echoes is as follows [26]: where P r is the radar received power, P t is the radar transmitted power, G is the radar antenna gain, λ is the operating wavelength, r is the distance between the measured target and the radar, N is the total number of precipitation particles within the effective illuminated volume of the radar beam that simultaneously produce scattering echo energy and return to the radar antenna, and σ i is the backscattering cross section of an individual precipitation particle.The radar volume reflectivity of precipitation is defined as the sum of the backscattering cross section of all precipitation particles per unit volume, denoted as η [27]: The volume resolution cell represents the minimum separable precipitation unit that SAR can distinguish during motion.It can be considered as the differential element of the effective illuminated volume.The volume of the cell is related to the three-dimensional resolutions of SAR [28]: where r R is the slant-range resolution, r a is the azimuth resolution, and r c is the resolution corresponding to the beam width in the pitch plane.
In practical applications, SAR precipitation measurements need to consider the impact of precipitation particles on the attenuation of radar echoes.L is the one-way attenuation factor.Therefore, the backscattering echoes power for detecting precipitation with SAR within the volume resolution cell can be expressed as follows: Therefore, the sum of the radar backscattering cross sections of all precipitation particles within the volume resolution cell is as follows: The observation range of spaceborne SAR is typically extensive with significant variations in σ.Instead, we generally use the dimensionless normalized radar cross section (NRCS), denoted as σ SAR , which is obtained by comparing it with the area.The transformation relationship between σ and σ SAR will be discussed in detail in Section 4.1.
The absorption and scattering of radar waves by precipitation particles between the radar and the target cause signal attenuation.In radar meteorology, the precipitation particles' attenuation and scattering effects on electromagnetic waves are typically described by the power law relationship [29,30].The relationship between the rainfall attenuation coefficient k(x, z) in km −1 , radar reflectivity η(x, z) in km −1 , and precipitation R(x, z) in mm/h is summarized as follows: where |K 0 | 2 is a function of the refractive index of water when neglecting the influence of temperature; the value of |K 0 | 2 is approximately 0.93 for rain and 0.19 for snow [13].In this paper, the relevant empirical parameters a, b, i, and j adopt the typical values statistically calculated by Ulaby et al.For rain, a is 2.6 × 10 −3 , b is 1.11, i is 300, and j is 1.35; for snow, a is 5.6 × 10 −3 , b is 1.60, i is 182, and j is 1.60 [31].
From Equations ( 5), (7), and ( 8), the precipitation distribution information can be retrieved using SAR echo data, which is the fundamental principle of SAR precipitation measurements.

Raindrop Motion Error Model
When SAR detects distributed targets like precipitation, it assumes that the precipitation particles remain relatively stationary with the ground.In the actual precipitation environment, the echo signal received by the radar is influenced by various factors, such as raindrop motion, atmospheric turbulence, and the effect of surface horizontal winds.The signal is the comprehensive result of the backscattering of all scattering particles within the effective illuminated volume.Such a comprehensive result leads to the standard deviation of the Doppler velocity spectrum, indicating that the echo signal contains multiple frequency components.Raindrop motion affects the azimuthal resolution of SAR precipitation measurement that, in turn, affects the precipitation volume resolution cell, echo intensity, and other factors [23].
Spaceborne SAR is a two-dimensional, high-resolution microwave imaging radar.It achieves high resolution in range by transmitting large bandwidth signals and achieves high resolution in azimuth by synthetic aperture technology [32].Under the impact of raindrop motion, the equivalent azimuthal resolution is as follows [33]: where u is the platform speed, r is the distance between SAR and the measured target, and σ v is the standard deviation of the Doppler speed spectrum of raindrops.The Doppler speed spectrum width characterizes the extent to which various sizes of Doppler speeds within the effective illuminated volume deviate from their average value, which is caused by the different radial velocities of the scattering particles.When SAR measures distributed targets such as raindrops, σ v is variable under different meteorological conditions [34].It is approximately 0.6 m/s for snowfall in a stationary atmosphere and 10 m/s for rainfall in strong wind shear.In practical calculations, σ v takes 1 m/s as the typical value for rainfall in a stationary atmosphere [26].
Based on the analysis above, a positive linear relationship exists between the standard deviation of the raindrop Doppler velocity spectrum and the azimuthal resolution of SAR precipitation measurements.On this basis, we establish the raindrop motion error model, as shown in Figure 2. First, consider calculating the SAR backscattering cross section under the ideal rain cloud with a uniform distribution of precipitation in both the horizontal and vertical directions.We analyze the propagation of electromagnetic waves and backscattering situations to derive a calculation method for the SAR backscattering cross section in the raindrop motion model in the following.
Atmosphere 2024, 15, x FOR PEER REVIEW 5 of 22 high resolution in azimuth by synthetic aperture technology [32].Under the impact of raindrop motion, the equivalent azimuthal resolution is as follows [33]: where  is the platform speed,  is the distance between SAR and the measured target, and  is the standard deviation of the Doppler speed spectrum of raindrops.The Doppler speed spectrum width characterizes the extent to which various sizes of Doppler speeds within the effective illuminated volume deviate from their average value, which is caused by the different radial velocities of the scattering particles.When SAR measures distributed targets such as raindrops,  is variable under different meteorological conditions [34].It is approximately 0.6 m/s for snowfall in a stationary atmosphere and 10 m/s for rainfall in strong wind shear.In practical calculations,  takes 1 m/s as the typical value for rainfall in a stationary atmosphere [26].
Based on the analysis above, a positive linear relationship exists between the standard deviation of the raindrop Doppler velocity spectrum and the azimuthal resolution of SAR precipitation measurements.On this basis, we establish the raindrop motion error model, as shown in Figure 2. First, consider calculating the SAR backscattering cross section under the ideal rain cloud with a uniform distribution of precipitation in both the horizontal and vertical directions.We analyze the propagation of electromagnetic waves and backscattering situations to derive a calculation method for the SAR backscattering cross section in the raindrop motion model in the following.In Figure 2,  is the incident angle,  is the cloud top height, and  is the distance in the rain for attenuation for a differential element on the x-z plane with a maximum length of  in km. is the distance from the surface to a differential volume for rain scatter on the x-z plane with a maximum length of  in km.The beam below the cloud top is divided into resolution cells corresponding to the slant range resolution  .The light blue area in the center of Figure 2 shows one of the rain resolution cells.The thickness of the cell is the slant range resolution  .The width of the cell is the azimuth resolution  .It is noteworthy that, on the incident plane, the length of the rain resolution cell is  , perpendicular to the direction from the radar to the rain, while the length of the illuminated volume defined by the three-dimensional In Figure 2, θ is the incident angle, z t is the cloud top height, and R is the distance in the rain for attenuation for a differential element on the x-z plane with a maximum length of R max in km.L is the distance from the surface to a differential volume for rain scatter on the x-z plane with a maximum length of L max in km.
The beam below the cloud top is divided into resolution cells corresponding to the slant range resolution r R .The light blue area in the center of Figure 2 shows one of the rain resolution cells.The thickness of the cell is the slant range resolution r R .The width of the cell is the azimuth resolution r a .It is noteworthy that, on the incident plane, the length of the rain resolution cell is L max , perpendicular to the direction from the radar to the rain, while the length of the illuminated volume defined by the three-dimensional resolutions corresponds to the differential amount dl.From Figure 2, we can derive the following relationship: Then, we can obtain the following: The volume of the rain resolution cell from the bottom and the ground area coupled to the rain cell can be calculated as follows [14]: The total radar backscattering cross section is the sum of the two components: The ideal rain cloud assumes a uniform distribution of precipitation in both the horizontal and vertical directions.From Equations ( 7) and ( 8), we can infer that a certain amount of precipitation corresponds to fixed values of the attenuation coefficient k and the radar reflectivity η.Without taking into account the attenuation of electromagnetic waves along the propagation path, the surface backscattering cross section σ sr f and the precipitation backscattering cross section σ vol can be expressed as follows: where σ 0 is the average backscattering cross section for vertically polarized radiation incident on land at 30 • , which is approximately −7 dB [35].Now consider a more realistic situation where different parts of each rain resolution cell experience different attenuation paths for backscattering.The closer to the top of the rain cloud, the smaller the attenuation.Near the top of the rain cloud, little attenuation occurs.The attenuation factors K sr f and K vol are used to characterize the attenuation effect on the ground surface and precipitation backscattering echoes in the actual precipitation environment.Equations ( 17) and ( 18) can be rewritten as follows: Assume that the effective volume of the rain cell in the actual rainfall environment is V e f f .Calculating the integral in Equation (20) and substituting in the definition of V e f f , we obtain the following: In summary, the surface backscattering cross section and the volume backscattering cross section in the actual precipitation environment, considering the raindrop motion, can be expressed as follows: where r a is the equivalent azimuthal resolution, which varies with the standard deviation of the Doppler velocity spectrum of raindrops.σ v generally ranges from 0.6 to 8.0 m/s.Under the raindrop motion model established in this paper, we derive formulas for calculating the effective surface and precipitation backscattering under different standard deviations of the Doppler velocity spectrum.This allows for further analysis of the SAR radar scattering cross section changes caused by raindrop motion.

RMM Algorithm
The RMM algorithm merges the proposed raindrop motion error model with the MOS inversion algorithm, enabling a quantitative evaluation of the impact of raindrop motion on the accuracy of MOS.The overall RMM algorithm is schematically shown in Figure 3, and each step is illustrated in the following.
Atmosphere 2024, 15, x FOR PEER REVIEW 7 of 22 In summary, the surface backscattering cross section and the volume backscattering cross section in the actual precipitation environment, considering the raindrop motion, can be expressed as follows: where  is the equivalent azimuthal resolution, which varies with the standard deviation of the Doppler velocity spectrum of raindrops. generally ranges from 0.6 to 8.0 m/s.Under the raindrop motion model established in this paper, we derive formulas for calculating the effective surface and precipitation backscattering under different standard deviations of the Doppler velocity spectrum.This allows for further analysis of the SAR radar scattering cross section changes caused by raindrop motion.

RMM Algorithm
The RMM algorithm merges the proposed raindrop motion error model with the MOS inversion algorithm, enabling a quantitative evaluation of the impact of raindrop motion on the accuracy of MOS.The overall RMM algorithm is schematically shown in Figure 3, and each step is illustrated in the following.

RCS and NRCS
The RMM algorithm uses simulated SAR echo data with raindrop motion errors as equivalent SAR echo data.Under the raindrop motion error model, we can calculate the RCS with varying levels of raindrop motion error by controlling the standard deviation of the raindrop Doppler velocity spectrum.Since the MOS algorithm uses NRCS data as the input, we first need to normalize the RCS to NRCS to quantitatively analyze the raindrop motion's impact on the MOS accuracy.
The X-SAR precipitation measurement model used to calculate the NRCS is shown in Figure 4.This model neglects the earth's curvature, does not consider beam width, and simplifies spherical waves into plane waves for analysis.

RCS and NRCS
The RMM algorithm uses simulated SAR echo data with raindrop motion errors as equivalent SAR echo data.Under the raindrop motion error model, we can calculate the RCS with varying levels of raindrop motion error by controlling the standard deviation of the raindrop Doppler velocity spectrum.Since the MOS algorithm uses NRCS data as the input, we first need to normalize the RCS to NRCS to quantitatively analyze the raindrop motion's impact on the MOS accuracy.
The X-SAR precipitation measurement model used to calculate the NRCS is shown in Figure 4.This model neglects the earth's curvature, does not consider beam width, and simplifies spherical waves into plane waves for analysis.The NRCS is derived from the surface backscattering cross section from the ground and the volume backscattering cross section from precipitation [10], as follows: where  and  and  and  are, respectively, the upper and lower attenuation limits of backscattering echoes from the surface and the volume resolution cell. and  are the upper and lower limits of the resolution cell in the rainfall area. is the top height of the rain layer, and  is the top height of the snow layer.The rain cloud width  extends from the left corner  to the right corner  .
Based on the analysis in Section 3, we can similarly derive the radar backscattering cross section for the double-layer rain cloud shown in Figure 4: Thus, we can also obtain the conversion relationship between NRCS and RCS as follows: where  denotes the NRCS data; the surface backscattering cross section  and the volume backscattering cross section  are the two components of the RCS.Next, we use a simple example to illustrate how errors are introduced into the MOS in the RMM algorithm.Assuming  = 1 m/s and  = 2 m/s, the relationship can be obtained from Equation ( 9): The original MOS algorithm does not consider raindrop motion and only applies to a windless situation with  = 1 m/s.When the actual rainfall environment is windy with The NRCS is derived from the surface backscattering cross section from the ground and the volume backscattering cross section from precipitation [10], as follows: where z s2 and z s1 and z r4 and z r3 are, respectively, the upper and lower attenuation limits of backscattering echoes from the surface and the volume resolution cell.z r2 and z r1 are the upper and lower limits of the resolution cell in the rainfall area.z 0 is the top height of the rain layer, and z t is the top height of the snow layer.The rain cloud width w extends from the left corner x L to the right corner x R .
Based on the analysis in Section 3, we can similarly derive the radar backscattering cross section for the double-layer rain cloud shown in Figure 4: Thus, we can also obtain the conversion relationship between NRCS and RCS as follows: where σ SAR denotes the NRCS data; the surface backscattering cross section σ sr f and the volume backscattering cross section σ vol are the two components of the RCS.
Next, we use a simple example to illustrate how errors are introduced into the MOS in the RMM algorithm.Assuming σ v1 = 1 m/s and σ v2 = 2 m/s, the relationship can be obtained from Equation (9): The original MOS algorithm does not consider raindrop motion and only applies to a windless situation with σ v = 1 m/s.When the actual rainfall environment is windy with σ v = 2 m/s, the NRCS, used as the input data to the algorithm, contains a certain degree of motion error.We use σ mSAR to represent the NRCS data with the raindrop motion error as follows: Atmosphere 2024, 15, 265 9 of 20 A conclusion can be drawn about the conversion relationship.For the double-layer rain cloud model in Figure 4, the sizes of the two sets of NRCS data satisfy a proportional relationship when the standard deviation of the raindrop Doppler velocity spectrum σ v is taken as 1 m/s and 2 m/s, respectively.The difference between the two is 3 dB when converted into decibels.

Retrieval of Precipitation Distribution
The RMM mainly uses the MOS inversion method to obtain precipitation distribution.The spatial distribution of precipitation can be categorized into different types [36].To simplify the problem, it assumes that the two-dimensional precipitation distribution R(x, z) can be divided into horizontal and vertical distributions.The two distributions are considered uncorrelated and can be retrieved separately: The horizontal weighting factor H(x) varies linearly and symmetrically with the distance from the edges over a thickness d: There are 18 existing types of H(x) [37].The three most commonly used are rectangular (d = 0), trapezoidal (0 < d < w/2), and triangular (d = w/2) horizontal distributions.In this paper, the parameters with superscripts represent the retrieved values, while the parameters without superscripts represent the given values.By applying the concept of maximum likelihood classification to the σ SAR data analysis, we can obtain Ĥ(x) [38].The MOS inversion method performs statistical calculations on simulated NRCS data, yielding the mean vector m cSAR and the covariance matrix C cSAR .The typical parameter vector derived from the actual equivalent NRCS data acquired by SAR is denoted as x mSAR [39].Comparing simulation and observation, we calculate the likelihood distance as follows: where c H is the rain cloud shape class.Different likelihood distances are associated with distinct horizontal rainfall distributions.The desired horizontal distribution is the one with the minimum likelihood distance.
A widely accepted functional form of V(z) is the following: where p represents the frozen particle distribution coefficient.z 0 and z t can be obtained from temperature soundings.
The philosophy behind RMM is to statistically analyze NRCS data and primarily identify the horizontal distribution H(x).Then, after classifying the shape, we eventually invert the vertical distribution V(z) based on the parameterized model given in Equations ( 32) and (33).The details of the inversion method for the main rainfall parameters are shown in Appendix A.
From Equations ( 32) and ( 33), it can be observed that the accuracy of the retrieved surface rain rate V(0) is crucial for the inversion accuracy of the vertical precipitation distribution.During the inversion process of V(z) by RMM, V(0) is given by the following: where xmin is the x-axis position corresponding to the minimum σ mSAR .The values 1.13, 21.62, 2.58, and 23.3 are all empirical values, only applicable to rain clouds with θ = 30 • , z t = 13 km, and z 0 = 4.5 km.If the height of the rain cloud or the incidence angle of X-SAR changes, all the empirical values mentioned above need to be recalculated, which will require a large amount of weather data as assistance [40].

RCS under Raindrop Motion Error Model
In case 1, we set h = 514 km, λ = 3.1 cm, σ 0 = −7 dB, z t = 13 km, θ = 30 • , and u = 7300 m/s.The precipitation rate R varies from 0 to 100 mm/h, and the standard deviation of the Doppler velocity spectrum of raindrops σ v takes values of 0.6 m/s, 1.0 m/s, 2.0 m/s, and 8.0 m/s, respectively.Using Equations ( 16), (22), and ( 23), the simulated curves of the surface backscattering cross section σ sr f , the volume backscattering cross section σ vol , and the total radar backscattering cross section σ with varying R and σ v are shown in Figure 5.
Equations (32) and (33).The details of the inversion method for the main rainfall parameters are shown in Appendix A.
From Equations ( 32) and ( 33), it can be observed that the accuracy of the retrieved surface rain rate  0 is crucial for the inversion accuracy of the vertical precipitation distribution.During the inversion process of   by RMM,  0 is given by the following: where  is the x-axis position corresponding to the minimum  .The values 1.13, 21.62, 2.58, and 23.3 are all empirical values, only applicable to rain clouds with  = 30°,  = 13 km, and  0 = 4.5 km.If the height of the rain cloud or the incidence angle of X-SAR changes, all the empirical values mentioned above need to be recalculated, which will require a large amount of weather data as assistance [40].

RCS under Raindrop Motion Error Model
In case 1, we set ℎ = 514 km,  = 3.1 cm,  = −7 dB,  = 13 km,  = 30°, and  = 7300 m/s.The precipitation rate  varies from 0 to 100 mm/h, and the standard deviation of the Doppler velocity spectrum of raindrops  takes values of 0.6 m/s, 1.0 m/s, 2.0 m/s, and 8.0 m/s, respectively.Using Equations ( 16), (22), and ( 23), the simulated curves of the surface backscattering cross section  , the volume backscattering cross section  , and the total radar backscattering cross section  with varying  and  are shown in Figure 5.  Figure 5 shows that the RCS is sensitive to changes in the standard deviation of the raindrop Doppler velocity spectrum.As σ v increases under the same precipitation rate, both the area and volume within the effective beam illuminated range increase, leading to a rise in σ sr f and σ vol .Additionally, when σ v is 1 m/s and 2 m/s, respectively, the difference between the two sets of σ sr f and σ vol curves remains constant at 3 dB.The difference is determined by the coefficients of the two parts on the right side of Equation ( 26), confirming the accuracy of the theoretical analysis.Based on this result, we can conduct the experiments below to examine the impact of raindrop motion on the accuracy of MOS inversion results.

Retrieval Results of Precipitation Distribution
In this section, the RMM algorithm is used to simulate the inversion results of surface rainfall rate, horizontal distribution, and two-dimensional precipitation distribution under windy and windless conditions, respectively.By comparing these two sets of results, we can evaluate the impact of raindrop motion on inversion accuracy.

Simulation without Raindrop Motion
When σ v = 1 m/s, it means that the raindrop motion error model in Figure 3 is not considered, and the RMM algorithm is simplified to the original MOS algorithm.We first simulate and analyze the retrieval results of the original MOS algorithm with various combinations of rainfall parameters in cases 2-4.We set the relevant parameters of the rainfall area as follows: w = 6 km, z 0 = 4.5 km, z t = 13 km, θ = 30 • , and σ 0 = −7 dB.It is assumed that R(x, 0) consists of a horizontal distribution H(x) and a surface rain rate V(0): R(x, 0) = H(x)V(0) Case 2: Taking a rectangular rain cloud (d = 0 km), the surface rain rate is 30 mm/h.The NRCS data and the inversion results are shown in Figure 6.The surface rain rate retrieved by RMM of σ v = 1 m/s is 29.36 mm/h.
In this section, the RMM algorithm is used to simulate the inversion results of surface rainfall rate, horizontal distribution, and two-dimensional precipitation distribution under windy and windless conditions, respectively.By comparing these two sets of results, we can evaluate the impact of raindrop motion on inversion accuracy.

Simulation without Raindrop Motion
When  = 1 m/s, it means that the raindrop motion error model in Figure 3 is not considered, and the RMM algorithm is simplified to the original MOS algorithm.We first simulate and analyze the retrieval results of the original MOS algorithm with various combinations of rainfall parameters in cases 2-4.We set the relevant parameters of the rainfall area as follows:  = 6 km,  = 4.5 km,  = 13 km,  = 30°, and  = −7 dB.It is assumed that  , 0 consists of a horizontal distribution   and a surface rain rate  0 :   Case 5: Assume that the model surface rain rate varies from 5 to 100 mm/h, and the horizontal precipitation distribution takes rectangular ( = 0 km), trapezoidal ( = 1.5 km), and triangular ( = 3 km) shapes, respectively.Three sets of the retrieved surface rates, each consisting of 20 values, are shown in Figure 9.
Case 6: Assume that the model surface rain rate varies from 20 to 40 mm/h, and the horizontal precipitation distribution takes rectangular ( = 0 km), trapezoidal ( = 1.5 km), and triangular ( = 3 km) shapes, respectively.Three sets of the retrieved surface rates, each consisting of 21 values, are shown in Figure 10.    Figure 6e shows the absolute error between the retrieved two-dimensional precipitation distribution  ,  and the given two-dimensional precipitation distribution  ,  .More significant absolute errors occur at the edges of the rainfall area.This is mainly because of the error in the retrieved rainfall area width  .Furthermore, the MOS can be applied to rain clouds of various shapes, as shown in Figures 6-8.This method can also be applied to various rainfall rates, as shown in Figures 9 and 10.The above results prove that the retrieval results using the MOS are in good agreement with the given rain Figure 6e shows the absolute error between the retrieved two-dimensional precipitation distribution R(x, z) and the given two-dimensional precipitation distribution R(x, z).More significant absolute errors occur at the edges of the rainfall area.This is mainly because of the error in the retrieved rainfall area width ŵ.Furthermore, the MOS can be applied to rain clouds of various shapes, as shown in Figures 6-8.This method can also be applied to various rainfall rates, as shown in Figures 9 and 10.The above results prove that the retrieval results using the MOS are in good agreement with the given rain cloud model in windless conditions, where the impact of raindrop motion error is negligible.
In Figure 9, the MOS exhibits a better inversion effect during moderate to heavy rainfall.Consequently, our subsequent research will primarily focus on experiments where the surface rain rate ranges from 20 to 40 mm/h.

Simulation with Raindrop Motion
In Figure 5, the experimental results show the RCS with raindrop motion error that, according to Equation (26), can be converted to NRCS data.These data are then used for precipitation inversion.For a fixed surface rain rate, we use the RMM algorithm to evaluate the inversion error caused by raindrop motion in cases 7-9.In cases 10-12, we compare inversion results under a range of surface rain rates, consisting of 21 values, to better illustrate the motion error's impact on the retrieved surface rain rate.Case 13 further explores the variances in retrieved surface rain rates for different standard deviations of the raindrop Doppler velocity spectrum.Case 10: Taking a rectangular rain cloud, the surface rain rate varies from 10 to 50 mm/h, and  takes 1 m/s and 1.1 m/s, respectively.The surface rain rates retrieved by RMM are shown in Figure 14.
Case 11: Taking a trapezoidal rain cloud, the surface rain rate varies from 10 to 50 mm/h, and  takes 1 m/s and 1.1 m/s, respectively.The surface rain rates retrieved by RMM are shown in Figure 15.Several groups of the retrieved surface rainfall data are selected for an error analysis.The relative error (RE) of the surface rain rate is as follows: The root mean square error (RMS) of the surface rain rate is as follows: Several groups of the retrieved surface rainfall data are selected for an error analysis.The relative error (RE) of the surface rain rate is as follows: The root mean square error (RMS) of the surface rain rate is as follows: The attributes of the six selected cases are shown in Table 1, and Table 2 shows the relative error for these cases.Table 3 shows the root mean square error of the cases 10-12.Taking cases 2-4 as the control group, cases 7-9 can provide insight into the inversion errors caused by raindrop motion in the MOS algorithm.Taking the rectangular rain cloud as an example, it can be observed from Figures 6a and 11a that the NRCS curve corresponding to  = 1.1 m/s is approximately 0.4 dB higher than the original NRCS curve.Comparing Figures 6e and 11e, it can be observed that raindrop motion also introduces significant errors to the retrieved two-dimensional precipitation distribution  ,  , which are beyond the acceptable error range.For the inversion of surface rain rate, under a rectangular rain cloud with  = 1.1 m/s and  0 = 30 mm/h, the relative error of  0 increases from the original 2.13% to 35.80% in Table 2. Similar conclusions also apply to trapezoidal and triangular rain clouds.The  under the three different shapes of the rain clouds are very close to each other in Table 3.The  of  = 1.0 m/s is approximately 0.05, but the  of  = 1.1 m/s increases to approximately 0.32.The results above show that the inversion accuracy of MOS is significantly affected after introducing different raindrop motion errors.
In Figures 14-16, the curve of retrieved surface rain rates  0 noticeably decreases when  is greater than 1 m/s.Figure 17 provides a clearer view of how retrieved surface rain rates change with  under three different precipitation distributions.Taking  = 1 m/s as the dividing line, the two sides of the horizontal axis represent the introduction of varying levels of raindrop motion errors.The value of  0 increases when  is less than 1 m/s and decreases when  is larger than 1 m/s.Moreover, the degree of deviation from  0 in  0 is related to the extent of deviation from 1 m/s in  .

Conclusions
In Section 3, we establish the raindrop motion error model and analyze the RCS under different raindrop Doppler velocity spectrum standard deviations to identify the RCS's variation pattern.In this way, the relationship between the wind field, raindrop motion, and SAR echoes is established.When  takes 1 m/s and 2 m/s, respectively, the difference between the two groups of  and  remains constant with a total RCS difference of 3 dB.Next, the RMM algorithm is proposed, based on the raindrop motion error model.This algorithm allows for us to examine the impact of raindrop motion on the original MOS inversion model and accurately calculate the resulting error in precipitation retrieval.In the RMM algorithm, introducing the raindrop motion error increases the relative error of the retrieved surface rain rate  0 from the original 2.13% to 35.80% Taking cases 2-4 as the control group, cases 7-9 can provide insight into the inversion errors caused by raindrop motion in the MOS algorithm.Taking the rectangular rain cloud as an example, it can be observed from Figures 6a and 11a that the NRCS curve corresponding to σ v = 1.1 m/s is approximately 0.4 dB higher than the original NRCS curve.Comparing Figures 6e and 11e, it can be observed that raindrop motion also introduces significant errors to the retrieved two-dimensional precipitation distribution R(x, z), which are beyond the acceptable error range.For the inversion of surface rain rate, under a rectangular rain cloud with σ v = 1.1 m/s and V(0) = 30 mm/h, the relative error of V(0) increases from the original 2.13% to 35.80% in Table 2. Similar conclusions also apply to trapezoidal and triangular rain clouds.The RMS 21 under the three different shapes of the rain clouds are very close to each other in Table 3.The RMS 21 of σ v = 1.0 m/s is approximately 0.05, but the RMS 21 of σ v = 1.1 m/s increases to approximately 0.32.The results above show that the inversion accuracy of MOS is significantly affected after introducing different raindrop motion errors.
In Figures 14-16, the curve of retrieved surface rain rates V(0) noticeably decreases when σ v is greater than 1 m/s.Figure 17 provides a clearer view of how retrieved surface rain rates change with σ v under three different precipitation distributions.Taking σ v = 1 m/s as the dividing line, the two sides of the horizontal axis represent the introduction of varying levels of raindrop motion errors.The value of V(0) increases when σ v is less than 1 m/s and decreases when σ v is larger than 1 m/s.Moreover, the degree of deviation from V(0) in V(0) is related to the extent of deviation from 1 m/s in σ v .

Conclusions
In Section 3, we establish the raindrop motion error model and analyze the RCS under different raindrop Doppler velocity spectrum standard deviations to identify the RCS's variation pattern.In this way, the relationship between the wind field, raindrop motion, and SAR echoes is established.When σ v takes 1 m/s and 2 m/s, respectively, the difference between the two groups of σ sr f and σ vol remains constant with a total RCS difference of 3 dB.Next, the RMM algorithm is proposed, based on the raindrop motion error model.This algorithm allows for us to examine the impact of raindrop motion on the original MOS inversion model and accurately calculate the resulting error in precipitation retrieval.In the RMM algorithm, introducing the raindrop motion error increases the relative error of the retrieved surface rain rate V(0) from the original 2.13% to 35.80% under a rectangular rain cloud with σ v = 1.1 m/s and V(0) = 30 mm/h.The degree to which the retrieved surface rain rate differs from the given one is strongly correlated with the value of σ v .The inversion error significantly increases, indicating that the impact of raindrop motion on the accuracy of the MOS algorithm cannot be neglected.The current inversion algorithms, such as the MOS, require modifications to the relevant retrieval models when a wind field is present in the rainfall environment.
The original MOS algorithm does not account for wind fields influencing raindrops in actual rainfall environments.Therefore, in conjunction with the analysis and simulation in this paper, we suggest adding a motion compensation module to further improve the accuracy and applicability of the MOS algorithm.The problem becomes tractable when we obtain the standard deviation of the Doppler velocity spectrum in actual precipitation measurements.Adding a data preprocessing step to the MOS can rectify the impact of wind-induced raindrop motion on SAR rainfall retrieval accuracy.
This study provides a preliminary understanding of the theoretical framework regarding how raindrop motion impacts the MOS inversion error, using the standard deviation of the Doppler velocity spectrum of raindrops σ v .We have made some simplifications to the rain motion model and the double-layer rain cloud model.Naturally, the variation in σ v itself, influenced by environmental wind field velocities, is a highly complex process with many error factors.In the actual atmosphere, wind speed and direction vary with altitude [41].In turbulent air, particles of a certain diameter within the effective illuminated volume not only fall at speeds influenced by environmental wind and gravity but also move with the surrounding atmospheric turbulent [42].Further research could consider several aspects, such as the observation and modelling of the background wind field, the motion characteristics of raindrops within the wind field, and the acquisition of raindrop spectrum distribution data.Therefore, a more comprehensive study needs to focus on developing a refined correlation between wind field, wind speed, and the standard deviation of the Doppler velocity spectrum of raindrops and studying intricate rainfall distributions in more realistic situations.This paper's research findings can provide a certain reference for rainfall inversion, taking into account the error of raindrop motion.

Figure 2 .
Figure 2. Raindrop motion error model.(Below) The geometry for calculating attenuation of the returns from different parts of the rain resolution cell; (Upper right) The three-dimensional structure of the rain resolution cell.

Figure 2 .
Figure 2. Raindrop motion error model.(Below) The geometry for calculating attenuation of the returns from different parts of the rain resolution cell; (Upper right) The three-dimensional structure of the rain resolution cell.

Figure 3 .
Figure 3.The process of RMM.

Figure 3 .
Figure 3.The process of RMM.

Figure 4 .
Figure 4. Schematic view of the model used to compute the NRCS from a double-layer rain cloud.

Figure 4 .
Figure 4. Schematic view of the model used to compute the NRCS from a double-layer rain cloud.

Figure 5 .
Figure 5. Simulation results of RCS.(a) The surface backscattering cross section; (b) The volume backscattering cross section; (c) The total radar backscattering cross section.

Figure 5 .
Figure 5. Simulation results of RCS.(a) The surface backscattering cross section; (b) The volume backscattering cross section; (c) The total radar backscattering cross section.

Figure 9 .
Figure 9. Surface rain rate of 5 to 100 mm/h retrieved by RMM of  = 1 m/s.

Figure 9 .
Figure 9. Surface rain rate of 5 to 100 mm/h retrieved by RMM of σ v = 1 m/s.

Figure 9 .
Figure 9. Surface rain rate of 5 to 100 mm/h retrieved by RMM of  = 1 m/s.

Figure 10 .
Figure 10.Surface rain rate of 20 to 40 mm/h retrieved by RMM of  = 1 m/s.

Figure 10 .
Figure 10.Surface rain rate of 20 to 40 mm/h retrieved by RMM of σ v = 1 m/s.

Case 7 :
Taking a rectangular rain cloud, σ v is 1.1 m/s, and the surface rain rate is 30 mm/h.The inversion results of RMM are shown in Figure 11.The surface rain rate retrieved by RMM is 19.26 mm/h.Case 8: Taking a trapezoidal rain cloud, σ v is 1.1 m/s, and the surface rain rate is 30 mm/h.The inversion results of RMM are shown in Figure 12.The surface rain rate retrieved by RMM is 20.79 mm/h.Case 9: Taking a triangular rain cloud, σ v is 1.1 m/s, and the surface rain rate is 30 mm/h.The inversion results of RMM are shown in Figure 13.The surface rain rate retrieved by RMM is 20.80 mm/h.

Figure 13 . 22 Case 12 :
Figure 13.Triangular distribution diagrams retrieved by RMM of σ v = 1.1 m/s.(a) The normalized radar cross section; (b) Simulated precipitation distribution R(x, 0).Case 10: Taking a rectangular rain cloud, the surface rain rate varies from 10 to 50 mm/h, and σ v takes 1 m/s and 1.1 m/s, respectively.The surface rain rates retrieved by RMM are shown in Figure14.

Figure 14 .
Figure 14.Scatterplot of surface rain rates retrieved by RMM under rectangular horizontal distribution.

Figure 15 .
Figure 15.Scatterplot of surface rain rates retrieved by RMM under trapezoidal horizontal distribution.

Figure 16 .
Figure 16.Scatterplot of surface rain rates retrieved by RMM under triangular horizontal distribution.

Figure 14 .
Figure 14.Scatterplot of surface rain rates retrieved by RMM under rectangular horizontal distribution.Case 11: Taking a trapezoidal rain cloud, the surface rain rate varies from 10 to 50 mm/h, and σ v takes 1 m/s and 1.1 m/s, respectively.The surface rain rates retrieved by RMM are shown in Figure 15.

Figure 14 .
Figure 14.Scatterplot of surface rain rates retrieved by RMM under rectangular horizontal distribution.

Figure 15 .
Figure 15.Scatterplot of surface rain rates retrieved by RMM under trapezoidal horizontal distribution.

Figure 15 .
Figure 15.Scatterplot of surface rain rates retrieved by RMM under trapezoidal horizontal distribution.Case 12: Taking a triangular rain cloud, the surface rain rate varies from 10 to 50 mm/h, and σ v takes 1 m/s and 1.1 m/s, respectively.The surface rain rates retrieved by RMM are shown in Figure16.

Figure 15 .
Figure 15.Scatterplot of surface rain rates retrieved by RMM under trapezoidal horizontal distribution.

Figure 16 .
Figure 16.Scatterplot of surface rain rates retrieved by RMM under triangular horizontal distribution.

Figure 16 .
Figure 16.Scatterplot of surface rain rates retrieved by RMM under triangular horizontal distribution.

22 Figure 17 .
Figure 17.Surface rain rates of  ranging from 0.6 to 1.2 m/s retrieved by RMM.

Figure 17 .
Figure 17.Surface rain rates of σ v ranging from 0.6 to 1.2 m/s retrieved by RMM.

Table 1 .
The attributes of the 6 cases.

Table 2 .
Analysis of relative error of the cases.

Table 3 .
Root mean square error of the cases 10-12.Assume that the standard deviation of the raindrop Doppler velocity spectrum ranges from 0.6 to 1.2 m/s, the surface rain rate is 30 mm/h, and the horizontal distribution of precipitation