Introduction
Arid areas are areas with dry climates, covering about 30% of the land area. Common features of arid areas are low precipitation, large temperature differences, evaporation may be much higher than precipitation, and water scarcity is the main factor limiting plant growth [Reference Guo, Qu, Wang, Xiao and Pang1,Reference Guan, Cheng, Guo and Ji2]. The frequent occurrence and long-term persistence of drought will not only bring huge losses to socio-economic, especially agricultural production, but also cause many adverse ecological and environmental impacts such as water shortage, increased desertification, and frequent sandstorms. At the same time, due to the abundance of heat, with reasonable irrigation and fertilization has the potential to become a high-yield area. For this reason, the study of arid areas, which focuses on various parameters such as water content, vegetation, sand and dust and temperature, can be studied and analyzed by radar remote sensing echoes [Reference Khanal, Klopfenstein, Kushal, Ramarao, Fulton, Douridas and Shearer3–Reference Sarala and Jacob6] and scattering echoes [Reference Li, Zhang, Jiang and Wei7–Reference Rowell and Stein9] to obtain relevant information about arid areas. By extracting and analyzing useful data, it can contribute to many aspects such as environmental detection, breeding and grazing, target detection, and precision strikes.
It can be noted from literatures that international scholars have done a lot of research on the characteristics of environmental echoes under different environments, but relatively little research has been done on arid regions. Tian [Reference Tian10] calculated the snow electromagnetic scattering coefficients using the small perturbation method, and generalized the effects of moisture and roughness on snow scattering. The study of [Reference Tian10] is targeted and its findings are not applicable to other environmental facets. Comite [Reference Comite and Pierdicca11] proposed a first-order small slope approximation (SSA) algorithm, established a reliable electromagnetic scattering model, studied the monostatic and bistatic scattering characteristics of cultivated land, and proposed a new understanding about anisotropic soil scattering, which is important for the inversion of soil water content. His study is very informative, but the effects of roughness and electromagnetic frequency are less studied. Johnson [Reference Johnson and Ouellette12] established an approximate and accurate surface scattering model, focusing on the locations where the minima occur in the scattering of the two polarization modes HH and VV, and found that the locations vary with the polarization mode and the roughness. Han [Reference Han13] focused on the soil types and compositions of saline, dry and wetlands, established a dielectric constant model, and proposed an inversion algorithm based on electromagnetic scattering. And the correctness of the inversion model was verified by case calculations. Han's research results can be used for soil moisture detection in saline soils, but its application is limited to some extent due to the specificity of its application.
In this paper, we focus on electromagnetic scattering in arid regions, and the experimental parameters are based on measured data from the southeastern region of Ejina Banner, Inner Mongolia [Reference Xie14]. Moreover, an accurate soil dielectric constant model is established and the effects of various factors, such as electromagnetic frequency and soil moisture, on the soil dielectric constant are summarized. In terms of research method, the traditional SSA method is modified in this paper to speed up the calculation process and reduce the error, and the correctness of the modified SSA is verified by simulation. Finally, it is found that the scattered echoes in arid areas have a large relationship with electromagnetic wave band, soil water content and ground roughness, and this connection is linear and regular.
Four component soil dielectric model
Soil dielectric constants are affected by various factors, and an accurate mathematical model to describe the dielectric properties of soils is necessary if highly accurate soil-related information is to be obtained. In this subsection, a four-component dielectric model for soils in arid regions was developed, and the effects of water content and electromagnetic band on the dielectric constant of soils were investigated. According to the four-component soil dielectric model, the empirical equations are as follows [Reference Wang, Guo and Wang15]
In the above formula, W p stands for soil moisture compression point, m t stands for critical body moisture. S(%) represents the sand content of the soil and C(%) represents the clay content. The dry soil density is represented by the symbol ρ b:
where ρ is the pore accumulation rate of the soil, ρ s represents the density of the rock in the soil and generally takes the value 2.65 g/cm3.
Soil moisture is represented by the symbol m v and is often used to calculate the relative dielectric constant of the soil.
In the above equation, the relative permittivity of ice, rock and air are represented by $\varepsilon _i$, $\varepsilon _r$and $\varepsilon _a$, respectively, whose values can be found in the literature [Reference Xie14]. The relative permittivity of pure water is
In the above Debye equation, $\varepsilon _{w0}$ denotes DC permittivity, τ w denotes relaxation time, and $T( ^\circ {\rm C})$ denotes soil temperature. They are related as
The effects of soil moisture and electromagnetic wave frequency on the relative permittivity were plotted in Fig. 1.
Rough surface modeling and scattering calculation
Two-dimensional rough surface modeling
In the study of this paper, in order to improve the accuracy of the model, the two-dimensional rough surface was generated with full consideration of the actual conditions of the soil ground and integrating the Monte Carlo method [Reference Caflisch16] and the Gaussian spectral function [Reference Bourlier, Bergine and Saillard17].
The surface height function of a two-dimensional rough surface is
In the formula, L x, L y represents the length of the rough surface in two dimensions, and b m,n is calculated as
N(0, 1) is the random sampling points on the normal distribution, N x and N y are the number of sampling points on the rough surface in two dimensions, W(k xm, k yn) represents
where h is the root-mean-square height of the Gaussian rough surface, l x and l y are the correlation lengths of the rough surface in two dimensions,k xm and k yn are the discrete points of the spatial frequencies in two dimensions, respectively.
The schematic diagram of the two-dimensional rough surfaces with different roughness generated by the method in this paper is shown in Fig. 2.
Conical incident wave
A finite truncation of the infinitely extended rough surface is required in the numerical calculation. In order to eliminate the effects such as reflection and edge bypassing caused by the abrupt truncation of the rough surface edge, a conical incident wave is used instead of plane wave incidence.
The conical wave requires attention to the beam width, which determines the irradiation width of the incident wave. When it is too large, the size of the rough surface needs to be increased accordingly to ensure the irradiated area; when it is too small, the rough surface and the target cannot be fully irradiated and a reasonable compound scattering field cannot be obtained. Reasonable calculation results and faster efficiency can be obtained when g and rough surface size L satisfy the following relationship
For TE incident waves, the incident electromagnetic field is
where E(k x, k y) is the incident spectrum, and the horizontal polarization direction ${\boldsymbol h}_i$ and vertical polarization direction ${\boldsymbol v}_i$ are denoted as
η 0 is the free space wave impedance, k x, k y, k z are the spatial spectral domains in the x, y, z directions, respectively, $k_\rho = \sqrt {k_x^2 + k_y^2 }$.
For TM incident waves, the incident electromagnetic field is
The angles of the incident waves are set to $\theta _i = 60^\circ$,$\varphi _i = 0^\circ$, and the two-dimensional normalized conical incident wave amplitudes are shown in Fig. 3. From the figure, it can be seen that the magnetic field is strongest at the center of the rough surface, and the farther away from the center, the weaker the magnetic field is, which slowly goes to zero at the edge. Therefore, the errors caused by the truncation of the rough surface can be avoided in the numerical simulation.
Modified small slope approximation
Based on the fact that the soil rough surface is anisotropic, the classical SSA method with appropriate modifications is used in this paper.
The core of the SSA method to solve the electromagnetic scattering problem is to solve for the scattering amplitude T. In the scattering field calculation, the solution equation for T is as follows
where the point on the scattering cross section S is labeled by the two-dimensional vector x, $P_{\boldsymbol k}( {\boldsymbol r})$ satisfying the boundary conditions and the Helmholtz equation.
Its complete solution can be obtained by integrating:
where r is the surface element vector and $G_k( {\boldsymbol r}-{\boldsymbol x})$ represents the Green's function.
According to the above equation, the scattering amplitude value T of the soil rough surface can be found, which in turn leads to σ suf:
Combine (23)(24) (25)(26)(27):
I(α) is the spectral integral of the rough surface and ${\boldsymbol k}-{\boldsymbol q} = ( Q_h, \;Q_z)$.Q h and Q z can be solved by the following equations.
It can be found that the soil rough surface scattering coefficient depends on the following integral.
where
In the above equation, J 0(y) is the first class zero-order Bessel function and B(y, α) is the calculated coefficient.
Combining the above equations, I(α) can be expanded into the following equation:
The summation based on the traditional SSA requires manual truncation, and when α is large, the formula converges quickly and the error generated by manual truncation is small; when α is small, the formula converges slowly, which will lead to the error generated by truncation beyond the engineering range, which is uncontrollable in the actual calculation, so it is unrealistic to apply equation (35) for the solution. Therefore, modifications based on the traditional SSA method can be considered to speed up the calculation process and reduce the error.
Based on the above idea, we can find that Γ( − z) has a remainder −( − 1)n/n! at a non-negative integer n. Obviously, by applying the residue theorem, the series terms in I(α) can be turned into:
Substituting (36) back into (35):
where
In the above equation, B(N) decays exponentially, and the decay accelerates as u increases.
Thus, the modified SSA method can be used to calculate the scattering coefficient of soil rough surface to ensure the accuracy and calculation speed at the same time, which meets the engineering requirements.
Validation of algorithm
The main work of this subsection is to verify the correctness of the SSA algorithm used in this paper. The validation algorithm uses the MOM algorithm, which is convincingly accurate and precise in its calculation. In the validation, both SSA and MOM methods were used to calculate the backscattering coefficients for the same rough surface with the following settings of the rough surface parameters with reference to the arid region: dimension L x × L y = 10λ × 10λ, correlation length l x × l y = 0.8λ × 0.8λ, root mean square height h = 0.1λ, dielectric constant (13.61, − 0.03) and incident wave frequency f = 300MHz.
The calculation results of the two methods are plotted in Fig. 4, and the two curves match well, indicating that the algorithm is effective in calculating the rough surface of Gaussian media.
Calculation results and analysis
The soil data in this paper were obtained from the measured data in the southeastern part of Ejina Banner, Inner Mongolia, which has typical characteristics of an arid region with little rainfall, high evaporation, sufficient sunshine, large temperature difference and sandy wind. The soil of this arid land is mainly composed of sand, chalk and clay, accounting for 66.97%, 20.78 and 12.25%, respectively, with a density of 1.173 g/cm3 and an average water content of 10.35%.
In the simulation, the effects of different soil roughness, different soil moisture content and different incident wavelengths on the backscattering characteristics in arid regions are mainly studied and discussed. Where, the size of the rough surface is L x × L y = 10λ × 10λ, and the correlation length is l x × l y = 0.8λ × 0.8λ.
Influence of different roughness
This subsection focuses on the effects of different soil roughness on the scattering characteristics in arid regions, and the backscattering coefficients were calculated for root-mean-square heights of 0.05, 0.10 and 0.15 m, respectively, and plotted in Fig. 5. The incident wave frequency is 300 MHz, the soil moisture content is 10.35%, and the HH and VV polarization results are as follows.
Observing Fig. 6, in general, the scattering intensity decreases linearly with the increase of radar incidence angle, and the scattering intensity is maximum when the incident wave is perpendicular to the ground. From the calculation results of different roughness, as the soil roughness becomes larger, the backward scattering increases, and they are positively correlated. According to this property, the roughness of the region can be obtained by comparison. And the roughness of arid areas is often affected by desertification - the more severe the desertification, the more sand and gravel fill the gullies, the flatter the land surface, the lower the roughness, so to some extent it can reflect the degree of desertification in the area. In another way, the intensity of the backscattering coefficient in a certain area can be continuously detected to study whether the desertification in the area has been aggravated or weakened, and to provide a reference for the subsequent decision of regional management.
Effect of soil moisture content
In arid areas where precipitation is scarce and evaporation is high, water scarcity is a major limiting factor for plant growth. However, arid areas have enough heat and they have the potential to become highly productive if they receive proper irrigation and fertilization. Therefore, this subsection calculates the scattering coefficients of soils with different moisture levels and discusses their scattering characteristics.
According to the actual measurement data in the southeastern part of Ejina Banner, Inner Mongolia, the soil moisture content varies greatly in different periods, with the highest value of 38.07% and the lowest value of 0.23%. In order to fully investigate the electromagnetic scattering characteristics in arid regions, the scattering results were calculated for two extreme water content conditions, as shown in Fig. 7.
It was found that for the HH polarization results, the two curves could be clearly distinguished in the (0°, 85°) range, and their overlap was very small. While in the VV polarization mode, the scattering coefficients of soils with different moisture contents become closer and closer with increasing angles, and finally almost overlap. In general, the soil moisture content has a strong influence on the scattering coefficients, and the differences are greatest when the incident wave is perpendicular to a rough surface, which can be easily distinguished with the naked eye. It can be inferred from the simulated results of soil water content under two extreme conditions in the region that soil water content is an important factor affecting the electromagnetic scattering characteristics in arid regions. The variation of soil water content will lead to fluctuations of electromagnetic scattering intensity on the rough surface. The results of this paper can be combined with remote sensing techniques to invert soil water content based on radar scattering coefficients, which is important for soil water content studies such as irrigation control, ecological studies, vorticity covariance, and slope stabilization climate science.
The frequency effect
Different wavelengths of electromagnetic waves have different penetrating abilities and different sensitivities to rough surface parameters. Therefore, it is important to study the electromagnetic scattering characteristics of arid areas under different electromagnetic spectrum environments and find the electromagnetic “optimal solution” for environmental protection and ecological monitoring.
In this subsection of the study, the soil moisture content was 10.35% and the root mean square height of the rough surface was 0.05 m. The calculated results were plotted in Fig. 8.
Through the curves we can find that, intuitively, the backscattering coefficient of the arid rough surface decreases continuously with the increase of the incident wave electromagnetic frequency. However, there are some differences between the results of HH polarization and VV polarization: ①For the VV polarization results, in the range of (0°, 90°), the three curves always do not intersect and are distinguished obviously; ②For the HH polarization results, in the angle of (0°, 48°), the three curves can be distinguished by the naked eye; in the range of (49°, 90°), the scattering coefficients of the three frequencies are extremely close to each other, and the magnitude cannot be determined by visual observation.
Therefore, when studying the effect of electromagnetic frequency on the scattering coefficient, better results can be obtained by choosing the incident wave with VV polarization.
Calculation
This paper focuses on the electromagnetic scattering characteristics in arid regions, and the parameters in the study are all derived from the measured data in the southeast of Ejinaqi, Inner Mongolia [Reference Xie14]. Soil water content, ground roughness and incident wave frequency were found to be important influencing factors in the study, and the scattering coefficients in arid lands subsequently showed regular variations. The reasons for the appearance of this regular characteristic were analyzed, and improved methods for monitoring soil moisture and desertification control were proposed, which have important significance in agricultural and ecological management. Besides, the results of this paper can be effectively combined with remote sensing technology to invert soil water content and study the ecological situation of vegetation, which is important for soil water content studies, such as irrigation control, ecological studies, vorticity covariance, slope stability climatology, etc. In the next research, it is planned to apply remote sensing technology to promote multi-faceted research on geomorphology, climate, and vegetation in arid regions.
Conflict of interest
There are no conflicts of interest to disclose for all the authors.
Fengyun Jiang was born in Jiangxi, China. She received the bachelor and master degree from Hohai University and Nanchang University, in 1999 and 2007, respectively. She is currently a lecturer of Yichun University, Yichun, Jiangxi, China. Her research interests include electromagnetic computing.