Improvement of X-Band Polarization Radar Melting Layer Recognition by the Bayesian Method and ITS Impact on Hydrometeor Classification

Using melting layer(ML)and non-melting layer(NML)data observed with the X-band dual linear polarization Doppler weather radar(X-POL)in Shunyi,Beijing,the reflectivity(ZH),differential reflectivity(ZDR),and correlation coefficient(CC)in the ML and NML are obtained in several stable precipitation processes.The prior probability density distributions(PDDs)of the ZH,ZDR and CC are calculated first,and then the probabilities of ZH,ZDR and CC at each radar gate are determined(PBB in the ML and PNB in the NML)by the Bayesian method.When PBB>PNB the gate belongs to the ML,and when PBB<PNB the gate belongs to the NML.The ML identification results with the Bayesian method are contrasUsing melting layer(ML)and non-melting layer(NML)data observed with the X-band dual linear polarization Doppler weather radar(X-POL)in Shunyi,Beijing,the reflectivity(ZH),differential reflectivity(ZDR),and correlation coefficient(CC)in the ML and NML are obtained in several stable precipitation processes.The prior probability density distributions(PDDs)of the ZH,ZDR and CC are calculated first,and then the probabilities of ZH,ZDR and CC at each radar gate are determined(PBB in the ML and PNB in the NML)by the Bayesian method.When PBB>PNB the gate belongs to the ML,and when PBB<PNB the gate belongs to the NML.The ML identification results with the Bayesian method are contrasted under the conditions of the independent PDDs and joint PDDs of the ZH,ZDR and CC.The results suggest that MLs can be identified effectively,although there are slight differences between the two methods.Because the values of the polarization parameters are similar in light rain and dry snow,it is difficult for the polarization radar to distinguish them.After using the Bayesian method to identify the ML,light rain and dry snow can be effectively separated with the X-POL observed data.ted under the conditions of the independent PDDs and joint PDDs of the ZH,ZDR and CC.The results suggest that MLs can be identified effectively,although there are slight differences between the two methods.Because the values of the polarization parameters are similar in light rain and dry snow,it is difficult for the polarization radar to distinguish them.After using the Bayesian method to identify the ML,light rain and dry snow can be effectively separated with the X-POL observed data.

Modified Method for Estimating Organic Carbon Density in Discontinuous Karst Soil Using Ground-Penetrating Radar and Geostatistics

The conventional method which assumes the soil distribution is continuous was unsuitable for estimating soil organic carbon density(SOCD) in Karst areas because of its discontinuous soil distribution. The accurate estimation of SOCD in Karst areas is essential for carbon sequestration assessment in China. In this study, a modified method,which considers the vertical proportion of soil area in the profile when calculating the SOCD, was developed to estimate the SOCD in a typical Karst peak-cluster depression area in southwest China. In the modified method, ground-penetrating radar(GPR) technology was used to detect the distribution and thickness of soil. The accuracy of the method was confirmed through comparison with the data obtained using a validation method, in which the soil thickness was measured by excavation. In comparison with the conventional method and average-soil-depth method,the SOCD estimated using the GPR method showed the minimum relative error with respect to that obtained using the validation method. At a regional scale, the average SOCDs at depths of 0-20 cm and 0-100 cm, which were interpolated by ordinary kriging,were 1.49(ranging from 0.03-5.65) and 2.26(0.09-11.60) kgm-2based on GPR method in our study area(covering 393.6 hm2), respectively. Therefore, the modified method can be applied on the accurate estimation of SOCD in discontinuous soil areas such as Karst regions.

Finite difference time domain method forward simulation of complex geoelectricity ground penetrating radar model

The ground penetrating radar(GPR) forward simulation all aims at the singular and regular models, such as sandwich model, round cavity, square cavity, and so on, which are comparably simple. But as to the forward of curl interface underground or “v” figure complex model, it is difficult to realize. So it is important to forward the complex geoelectricity model. This paper takes two Maxwell’s vorticity equations as departure point, makes use of the principles of Yee’s space grid model theory and the basic principle finite difference time domain method, and deduces a GPR forward system of equation of two dimensional spaces. The Mur super absorbed boundary condition is adopted to solve the super strong reflection on the interceptive boundary when there is the forward simulation. And a self-made program is used to process forward simulation to two typical geoelectricity model.

Modified OMP method for multi-target parameter estimation in frequency-agile distributed MIMO radar

Introducing frequency agility into a distributed multipleinput multiple-output(MIMO)radar can significantly enhance its anti-jamming ability.However,it would cause the sidelobe pedestal problem in multi-target parameter estimation.Sparse recovery is an effective way to address this problem,but it cannot be directly utilized for multi-target parameter estimation in frequency-agile distributed MIMO radars due to spatial diversity.In this paper,we propose an algorithm for multi-target parameter estimation according to the signal model of frequency-agile distributed MIMO radars,by modifying the orthogonal matching pursuit(OMP)algorithm.The effectiveness of the proposed method is then verified by simulation results.

Three-Dimensional Wind Field Retrieved from Dual-Doppler Radar Based on a Variational Method:Refinement of Vertical Velocity Estimates

In this paper,a scheme of dual-Doppler radar wind analysis based on a three-dimensional variational method is proposed and performed in two steps.First,the horizontal wind field is simultaneously recovered through minimizing a cost function defined as a radial observation term with the standard conjugate gradient method,avoiding a weighting parameter specification step.Compared with conventional dual-Doppler wind synthesis approaches,this variational method minimizes errors caused by interpolation from radar observation to analysis grid in the iterative solution process,which is one of the main sources of errors.Then,through the accelerated Liebmann method,the vertical velocity is further reestimated as an extra step by solving the Poisson equation with impermeable conditions imposed at the ground and near the tropopause.The Poisson equation defined by the second derivative of the vertical velocity is derived from the mass continuity equation.Compared with the method proposed by O’Brien,this method is less sensitive to the uncertainty of the boundary conditions and has better stability and reliability.Furthermore,the method proposed in this paper is applied to Doppler radar observation of a squall line process.It is shown that the retrieved vertical wind profile agrees well with the vertical profile obtained with the velocity–azimuth display(VAD)method,and the retrieved radial velocity as well as the analyzed positive and negative velocity centers and horizontal wind shear of the squall line are in accord with radar observations.There is a good correspondence between the divergence field of the derived wind field and the vertical velocity.And,the horizontal and vertical circulations within and around the squall line,as well as strong updrafts,the associated downdrafts,and associated rear inflow of the bow echo,are analyzed well.It is worth mentioning that the variational method in this paper can be applied to simultaneously synthesize the three-dimensional wind field from multiple-Doppler radar observations.

Practical calibration method for instantaneous polarization radar systems utilizing a metal sphere

A practical calibration method is proposed for instantaneous polarization radar systems.The method only needs one measurement by using a metal sphere.The distortions of system and the actual polarization scattering matrix（PSM）of target can be obtained.First,an instantaneous polarization radar system is presented.The system can obtain PSM by a single pulse echo.The dual-polarization antenna can transmit and receive two orthogonal polarization waves.The multilayer micro-strip patch antenna is adopted for this kind of radar system.Second,based on the multi-port network theory,the operation and system errors of instantaneous polarization radar system are analyzed.By making assumption on the cross-talk factors of antenna,distortion matrices of R and Tare derived.Finally,the calibration method based on instantaneous polarization measurement is introduced.Simulation results show the performance of this calibration method.The values of calibrated PSM are in agreement with the actual ones after calibration.

An element-free Galerkin method for ground penetrating radar numerical simulation

An element-free Galerkin method(EFGM) is used to solve the two-dimensional(2D) ground penetrating radar(GPR)modelling problems, due to its simple pre-processing, the absence of elements and high accuracy. Different from element-based numerical methods, this approach makes nodes free from the elemental restraint and avoids the explicit mesh discretization. First, we derived the boundary value problem for the 2D GPR simulation problems. Second, a penalty function approach and a boundary condition truncated method were used to enforce the essential and the absorbing boundary conditions, respectively. A three-layered GPR model was used to verify our element-free approach. The numerical solutions show that our solutions have an excellent agreement with solutions of a finite element method(FEM). Then, we used the EFGM to simulate one more complex model to show its capability and limitations. Simulation results show that one obvious advantage of EFGM is the absence of element mesh, which makes the method very flexible. Due to the use of MLS fitting, a key feature of EFM, is that both the dependent variable and its gradient are continuous and have high precision.

RADAR CROSS-SECTION COMPUTATIONS OF ARBITRARILY COMPLICATED TARGETS BY APPLYING THE PANEL METHOD

The present paper deals with the method for the radar cross-section (RCS)computations of arbitrarily complicated targets based on the work by D. Klement et al.(1988).This method is convenient in use, fast in operation and precise in calculating RCS of a complicatedtarget. With this method, the RCS of classic scatterers, for example, a cone and a cylinder, arecomputed with the result of good agreement with experimental data. Furthermore, the RCS’of an aircraft model at various attitudes are calculated with the result of good agreement withexperimental data also.

Efficient Calculation of the Radar Cross Section at the Interior Resonance by the Inverse Power Method

It is well known that the incorrect results will be given using either the electric or magnetic field integral equation to calculate the radar cross section (RCS) of a closed body at the interior resonance. In this paper, an effective iterative technique is used to correct the calculated surface current density from the electric field integral equation. The radar cross section is computed for an infinite conducting circular cylinder at the interior resonance, and the obtained results are in good agreement with the analytical results. The backscattering cross section of an infinite triangular cylinder in the vicinity of a resonant frequency is also calculated. It is shown that the presence method is efficient and accurate.

Solving multi-object radar cross section based on wide-angle parabolic equation method

Based on a Pade approximation, a wide-angle parabolic equation method is introduced for computing the multiobject radar cross section （RCS） for the first time. The method is a paraxial version of the scalar wave equation, which solves the field by marching them along the paraxial direction. Numerical results show that a single wide-angle parabofic equation run can compute multi-object RCS efficiently for angles up to 45 ° . The method provides anew and efficient numerical method for computation electromagnetics.

Scaled radar cross section measurement method for lossy targets via dynamically matching reflection coefficients in THz band

In the terahertz band,the dispersive characteristic of dielectric material is one of the major problems in the scaled radar cross section(RCS)measurement,which is inconsistent with the electrodynamics similitude deducted according to the Maxwell’s equations.Based on the high-frequency estimation method of physical optics(PO),a scaled RCS measurement method for lossy objects is proposed through dynamically matching the reflection coefficients according to the distribution of the object facets.Simulations of the model of SLICY are conducted,and the inversed RCS of the lossy prototype is obtained using the proposed method.Comparing the inversed RCS with the calculated results,the validity of the proposed method is demonstrated.The proposed method provides an effective solution to the scaled RCS measurement for lossy objects in the THz band.

A Mo LC+Mo M-based G^0 distribution parameter estimation method with application to synthetic aperture radar target detection

The accuracy of background clutter model is a key factor which determines the performance of a constant false alarm rate(CFAR) target detection method. G0 distribution is one of the optimal statistic models in the synthetic aperture radar(SAR) image background clutter modeling and can accurately model various complex background clutters in the SAR images. But the application of the distribution is greatly limited by its disadvantages that the parameter estimation is complex and the local detection threshold is difficult to be obtained. In order to solve the above-mentioned problems, an synthetic aperture radar CFAR target detection method using the logarithmic cumulant(Mo LC) + method of moment(Mo M)-based G0 distribution clutter model is proposed. In the method, G0 distribution is used for modeling the background clutters, a new Mo LC+Mo M-based parameter estimation method coupled with a fast iterative algorithm is used for estimating the parameters of G0 distribution and an exquisite dichotomy method is used for obtaining the local detection threshold of CFAR detection, which greatly improves the computational efficiency, detection performance and environmental adaptability of CFAR detection. Experimental results show that the proposed SAR CFAR target detection method has good target detection performance in various complex background clutter environments.

The Adjustment Analysis Method of the Active Surface Antenna Based on Convolutional Neural Network

Active surface technique is one of the key technologies to ensure the reflector accuracy of the millimeter/submillimeter wave large reflector antenna.The antenna is complex,large-scale,and high-precision equipment,and its active surfaces are affected by various factors that are difficult to comprehensively deal with.In this paper,based on the advantage of the deep learning method that can be improved through data learning,we propose the active adjustment value analysis method of large reflector antenna based on deep learning.This method constructs a neural network model for antenna active adjustment analysis in view of the fact that a large reflector antenna consists of multiple panels spliced together.Based on the constraint that a single actuator has to support multiple panels(usually 4),an autonomously learned neural network emphasis layer module is designed to enhance the adaptability of the active adjustment neural network model.The classical 8-meter antenna is used as a case study,the actuators have a mean adjustment error of 0.00252 mm,and the corresponding antenna surface error is0.00523 mm.This active adjustment result shows the effectiveness of the method in this paper.

A New Algorithm of Rain Type Classification for GPM Dual-Frequency Precipitation Radar in Summer Tibetan Plateau

In this study, a new rain type classification algorithm for the Dual-Frequency Precipitation Radar(DPR) suitable over the Tibetan Plateau(TP) was proposed by analyzing Global Precipitation Measurement(GPM) DPR Level-2 data in summer from 2014 to 2020. It was found that the DPR rain type classification algorithm(simply called DPR algorithm) has mis-identification problems in two aspects in summer TP. In the new algorithm of rain type classification in summer TP,four rain types are classified by using new thresholds, such as the maximum reflectivity factor, the difference between the maximum reflectivity factor and the background maximum reflectivity factor, and the echo top height. In the threshold of the maximum reflectivity factors, 30 d BZ and 18 d BZ are both thresholds to separate strong convective precipitation, weak convective precipitation and weak precipitation. The results illustrate obvious differences of radar reflectivity factor and vertical velocity among the three rain types in summer TP, such as the reflectivity factor of most strong convective precipitation distributes from 15 d BZ to near 35 d BZ from 4 km to 13 km, and increases almost linearly with the decrease in height. For most weak convective precipitation, the reflectivity factor distributes from 15 d BZ to 28 d BZ with the height from 4 km to 9 km. For weak precipitation, the reflectivity factor mainly distributes in range of 15–25 d BZ with height within 4–10 km. It is also shows that weak precipitation is the dominant rain type in summer TP, accounting for 40%–80%,followed by weak convective precipitation(25%–40%), and strong convective precipitation has the least proportion(less than 30%).

Sensitivity Analysis of Electromagnetic Scattering from Dielectric Targets with Polynomial Chaos Expansion and Method of Moments

In this paper,an adaptive polynomial chaos expansion method(PCE)based on the method of moments(MoM)is proposed to construct surrogate models for electromagnetic scattering and further sensitivity analysis.The MoM is applied to accurately solve the electric field integral equation(EFIE)of electromagnetic scattering from homogeneous dielectric targets.Within the bistatic radar cross section(RCS)as the research object,the adaptive PCE algorithm is devoted to selecting the appropriate order to construct the multivariate surrogate model.The corresponding sensitivity results are given by the further derivative operation,which is compared with those of the finite difference method(FDM).Several examples are provided to demonstrate the effectiveness of the proposed algorithm for sensitivity analysis of electromagnetic scattering from homogeneous dielectric targets.

An empirical method for joint inversion of wave and wind parameters based on SAR and wave spectrometer data

Synthetic aperture radar(SAR)and wave spectrometers,crucial in microwave remote sensing,play an essential role in monitoring sea surface wind and wave conditions.However,they face inherent limitations in observing sea surface phenomena.SAR systems,for instance,are hindered by an azimuth cut-off phenomenon in sea surface wind field observation.Wave spectrometers,while unaffected by the azimuth cutoff phenomenon,struggle with low azimuth resolution,impacting the capture of detailed wave and wind field data.This study utilizes SAR and surface wave investigation and monitoring(SWIM)data to initially extract key feature parameters,which are then prioritized using the extreme gradient boosting(XGBoost)algorithm.The research further addresses feature collinearity through a combined analysis of feature importance and correlation,leading to the development of an inversion model for wave and wind parameters based on XGBoost.A comparative analysis of this model with ERA5 reanalysis and buoy data for of significant wave height,mean wave period,wind direction,and wind speed reveals root mean square errors of 0.212 m,0.525 s,27.446°,and 1.092 m/s,compared to 0.314 m,0.888 s,27.698°,and 1.315 m/s from buoy data,respectively.These results demonstrate the model’s effective retrieval of wave and wind parameters.Finally,the model,incorporating altimeter and scatterometer data,is evaluated against SAR/SWIM single and dual payload inversion methods across different wind speeds.This comparison highlights the model’s superior inversion accuracy over other methods.

Direction finding of bistatic MIMO radar in strong impulse noise

For bistatic multiple-input multiple-output(MIMO)radar,this paper presents a robust and direction finding method in strong impulse noise environment.By means of a new lower order covariance,the method is effective in suppressing impulse noise and achieving superior direction finding performance using the maximum likelihood(ML)estimation method.A quantum equilibrium optimizer algorithm(QEOA)is devised to resolve the corresponding objective function for efficient and accurate direc-tion finding.The results of simulation reveal the capability of the presented method in success rate and root mean square error over existing direction-finding methods in different application situations,e.g.,locating coherent signal sources with very few snapshots in strong impulse noise.Other than that,the Cramér-Rao bound(CRB)under impulse noise environment has been drawn to test the capability of the presented method.

Method of Monitoring Three-dimensional Mining Surface Deformation Based on InSAR

In order to solve the problems of small monitoring range,long time and high cost of existing sedimentation observation methods,based on two-view sentinel No.1 radar images of Guqiao mining area in Huainan City from November 4,2017 to November 28,2017,surface change information was obtained in combination with D-InSAR,and the three-dimensional surface deformation was monitored by two-pass method and single line of sight D-InSAR method.The results show that during the research period of 24 d,the maximum deformation of the mining area reached 71 mm,and the southern subsidence was the most obvious,which was in line with the mining subsidence law.The maximum displacement from the north to the south was about 250 mm,while the maximum displacement from the east to the west was about 80 mm,and the maximum subsidence in the center was 110 mm.It is concluded that D-InSAR technique has a good effect on the inversion of the mining subsidence,and this method is suitable for three-dimensional surface monitoring in areas with similar geological conditions.The monitoring results have certain reference value.

Application of geophysical methods in fine detection of urban concealed karst:A case study of Wuhan City,China

The construction of modern livable cities faces challenges in karst areas,including ground collapse and engineering problems.Wuhan,with a population of 13.74×10^(6) and approximately 1161 km^(2)of soluble rocks in the urban area of 8569.15 km^(2),predominantly consists of concealed karst areas where occasional ground collapse events occur,posing significant threats to underground engineering projects.To address these challenges,a comprehensive geological survey was conducted in Wuhan,focusing on major karstrelated issues.Geophysical methods offer advantages over drilling in detecting concealed karst areas due to their efficiency,non-destructiveness,and flexibility.This paper reviewed the karst geological characteristics in Wuhan and the geophysical exploration methods for karst,selected eight effective geophysical methods for field experimentation,evaluated their suitability,and proposed method combinations for different karst scenarios.The results show that different geophysical methods have varying applicability for karst detection in Wuhan,and combining multiple methods enhances detection effectiveness.The specific recommendations for method combinations provided in this study serve as a valuable reference for karst detection in Wuhan.

PHASE AUTOFOCUSING ALGORITHM FOR COMPRESSED INVERSE SYNTHETIC APERTURE RADAR IMAGING

As same as the conventional inverse synthetic aperture radar（ISAR）, the compressed ISAR also requires the echo signal based motion compensation, which consists of the range alignment and the phase autofoeusing. A phase autofocusing algorithm for compressed ISAR imaging is presented. In the algorithm, phase autofocusing for the sparse ISAR echoes is accomplished using the eigenvector method. Experimental results validate the effectiveness of the algorithm.