Search

Objectives: In 2015, the Mw 7.8 earthquake in Nepal caused landslides and seriously impacted Nepal and its surrounding areas. Methods: In this paper, the Sindhupalchowk area, which is close to the epicenter of the earthquake, is selected as the study area. First, based on the L-band ALOS-2 and C-band Sentinel-1A synthetic aperture radar data, the detection and identification of post-earthquake landslides are carried out by using stacking-InSAR(interferometric synthetic aperture radar) technology. Second, the two-dimensional time series deformation characteristics of a typical landslide are obtained by combining the ascending and descending orbit data and multi-dimensional small baseline subset-InSAR technology. Results: Combined with optical images, 14 hidden points of landslide are delineated. The main deformation of the typical landslide occurs in the horizontal east-west direction, and the maximum deformation rate is −69 mm/a. Conclusions: By analyzing the trend and periodic term deformation signals in time series of the typical landslide, it is found that the earthquake has a significant acceleration effect on the landslide movement, and the increase of rainfall makes the landslide displacement show a periodic change from August to November every year. This paper can provide a reference for post-earthquake landslide monitoring re⁃ search. [ABSTRACT FROM AUTHOR]

Multi-epoch interferometric synthetic aperture radar (InSAR) is a highly effective technique for monitoring deformation in urban areas. However, interpreting InSAR deformation can be challenging due to various factors, including inherent geometric imaging distortion, the intricate structure and deformation properties of targets in urban scenes, and the multiple scattering of microwave signals between objects in urban scenes. This paper discusses the challenges involved in interpreting time-series InSAR deformation: (1) Precisely identifying the location of deformation signals and linking them to their corresponding objects, i. e., determining where the deformation signal occurs, (2) understanding the mechanisms and factors that cause the detected deformation signals, i.e., determining what the deformation signal represents, (3) establishing the connection among the detected deformation signals, the deformation events, and the scattering mechanisms. We suggest a parametric framework to improve the accurate interpretation of InSAR deformation. This framework includes several factors, including kinematic characteristics (deformation rate, cumulative deformation, deformation gradient, and deformation model), geometric parameters (position, size, structure, orientation, and roughness), semantic information (land cover type, terrain morphology, texture, and auxiliary information on natural and anthropogenic disturbance) and physical properties (scattering mechanism, penetrability, extensibility, conductivity, and thermal conductivity). Our approach aims to enhance the representation of coherent points for a better understanding of InSAR deformation. This paper offers a comprehensive overview of the advancements achieved in extracting parameters of InSAR coherent points and interpreting deformation based on geometric parameters, semantic information, and physical properties. High-precision 3D positioning is crucial for InSAR fine monitoring in urban areas. It helps determine the source of deformation signals and facilitates the analysis of deformation mechanisms. Semantic information, such as 3D models, high-resolution optical images, laser point cloud data, and land use data, can aid in interpreting InSAR deformation. By combining InSAR deformation data with a deep learning approach, there is an opportunity to interpret deformations effectively. In urban environments, the scattering mechanism of ground objects is complex. Multiple scattering signals can provide effective observations of deformation and information about the target's size. However, combining the scattering mechanism of synthetic aperture radar signals to carry out parameter inversion and deformation mechanism interpretation of urban target terrain remains a challenge. The framework, which considers the geometric parameters, semantic information, and physical attributes of InSAR coherent points, will be crucial for deformation interpretation and mechanism cognition. This framework will enable fine deformation monitoring, intelligent recognition, and application in future urban scenes. [ABSTRACT FROM AUTHOR]

Objectives: The key phenological information of oilseed rapeseed (Brassica napus L.) plays an important role in field management, viewing time prediction and yield estimation of the oilseed rape. It is also an important part of precision agriculture. Polarimetric synthetic aperture radar technology shows great potential in phenological phase identification with its all-weather monitoring capability and its sensitivity to the crop structural information. Methods: First, we identified the 5 phenological phases of the oilseed rape on the test area with 5 time series full-polarization Radarsat-2 data, which covers the whole growth period of the oilseed rape. 6 typical Stokes parameters are extracted and applied in the identification of oilseed rape phenological phases, the extracted Stokes parameters includ averaged intensity(g0), normalized average intensity(g0m), averaged degree of polarization(ρm), perimeter degree of zero orientation route(Pdor),inclination degree of zero aperture route(Idap), and arc asymmetry degree of zero aperture route(Aadap). Then, The phenological phases of oilseed rape is identified by the decision tree (DT) algorithm based on the comparative analysis of the dynamic response of the 6 special Stokes parameters to rape growth stages. Results and Conclusions: Among the extracted Stokes parameters applied in this study, except ρm and Aadap, other parameters show great sensitivity to the change of the oilseed rape phenological phases. The DT algorithm also perform well in the classification of the oilseed rape phenological phases. The classification results agree well with the field measured samples, and the overall classification accuracy is 87.4%, while the highest classifi⁃ cation accuracy of each phenological phase is 94.3%. [ABSTRACT FROM AUTHOR]

Objectives: Differential interferometric synthetic aperture radar (D-InSAR) has its own unique monitoring methods in acquiring surface deformation. In view of the difficulties in the application of D-InSAR technology in three-dimensional (3D) deformation monitoring of coal mine surface, this paper proposes a dynamic monitoring method based on fusion of single sight D-InSAR and Boltzmann-Knothe (BK) model. Methods: First, taking into account the characteristics of D-InSAR that it is difficult to obtain large gradient deformations, the Boltzmann function model with good boundary fit was selected, and BK model was constructed by the fusion of the Knothe time function. Second, due to the decoherence of the time-space baseline and the limitation of D-InSAR principle, it can often only obtain the short-term, one-dimensional deformation of the surface in line of sight (LOS). Therefore, based on the projection relationship between the 3D deformation and the LOS deformation, combined with the geological mining conditions of the target area, the fitness function with non-cumulative deformation was constructed. Finally, a dynamic monitoring method of surface 3D deformation of coal mine based on fusion of single sight D-InSAR and BK model was proposed by introducing the fireworks algorithm (FWA). Results and Conclusions: Simulation experiments show that the proposed method can reliably and accurately inverted all dynamic prediction parameters of mining subsidence, and the relative error of the solved parameters was between 0.11% and 7.51%. The subsidence and horizontal movement monitored by this proposed method were consistent with the real values, and the monitoring effect of the PIM-Knothe method was distorted. The results of the anti-error experiment show that under the influence of observation error or model parameter error, the overall effect of the solved parameters by the method in this paper was better, and the ability to resist errors was shown to a certain extent. The proposed method was applied to the 13121 working face of Gubei Mine in Huainan, and the surface 3D deformation of the coal mine from 2019-10-25 to 2019-12-12 was obtained. The conversion of short-term D-InSAR observation to 3D deformation dynamic monitoring based on the proposed method was realized, and the monitoring results were more consistent with the surface movement and deformation measured by the observation station. The reliability and scientificity of the dynamic prediction method for solving the surface 3D deformation of coal seam based on D-InSAR technology has been verified. [ABSTRACT FROM AUTHOR]

Objectives: Interferogram filtering is the key to the subsequent processing steps of interferometric synthetic aperture radar, such as phase unwrapping and geocoding. However, the existing filtering methods can't retain features in dense fringes and accurately estimate phase in low-coherence regions. The convolutional neural networks (CNN) is introduced to learn noise features and solve this problem. Methods: First, we selected a certain number of interferograms as samples, and divided them into training set, test set, and validation set. Second,we preprocessed the training set samples, and cut the preprocessed training set interferograms into small fixed-size block and randomly extract it as a model training sample, used the above steps to train the autoencoder filter model, after a certain number of iterations, the model was fitted. Results: Experiments were carried out on spaceborne imaging radar-C-band synthetic aperture radar data and Sentinel-1A data. Our proposed method was compared with Goldstein filter, mean filter, Lee filter, Frost filter, and improved denoising convolutional neural network (DnCNN). Goldstein filter can remove most of the noise while maintaining fringes edge, has good denoising ability. Mean filter performs well in high-coherence areas, but performs poorly in low-coherence areas, and can't filter noise well. Lee filter maintains it well image resolution, but the denoising effect is weak, and it can be seen from the filtering results that there is still a lot of noise. Frost filtering is weak in low-coherence areas and there is a lot of noise, but the fringe edges are well maintained in high-coherence areas. The improved DnCNN filter can significantly eliminate the noise, but it can't distinguish the fringe edge and the noise well. Our proposed method can suppress the noise very well, and it can maintain the fringe edge well in the low-coherence area and the high-coherence area. Conclusions: This proposed method can greatly improve the phase quality of the interferogram, suppress the noise to a greater extent, and restore more image details and maintain the edge continuity of the interference fringe. [ABSTRACT FROM AUTHOR]

Objectives: Multi-temporal interferometric synthetic aperture radar (InSAR) technology is widely used in deformation monitoring over wide areas due to its characteristics of large-scale, high spatial-temporal resolution and the capability to monitor deformation signal with millimeter precision even submillimeter. For the time-consuming problem in acquiring image and data preprocessing over research area, this paper uses looking inside the continents from space with synthetic aperture radar (LiCSAR) products and takes a new time series analysis LiCSBAS method, which greatly improve the computational efficiency. Methods: The quality check and the closure error of the phase loop are used to detect the unwrapping error in the two-dimensional phase unwrapping, and the interferogram with low coherence and the large error pixels are excluded. To minimize atmospheric artifacts, utilizing the generic atmospheric correction online service for In‑ SAR (GACOS) products and based on this, a method for joint estimation of deformation and atmospheric phase is proposed. Taking the Chengdu Plain and western mountainous area as the study area, totally 90 [ABSTRACT FROM AUTHOR]

Objectives: Machine learning classifier‑based water body extraction methods for polarimetric synthetic aperture radar (PolSAR) images have high reliability but typically require a great number of training samples. Consequently, it is very difficult and time‑consuming to manually collect enough training samples when extracting water body from multi‑temporal PolSAR images. To this problem, transfer learning is used to reduce the labor cost of querying new samples and improve the timeliness of water body extraction of multi‑temporal PolSAR images. Methods: Firstly, an optimal source domain image from multi‑temporal images is automatically selected according to the distribution difference between images, and the other images are taken as target domain images. Secondly, a group of training samples are queried in the source domain image as source sample set, and the same number of unlabeled samples are randomly sampled from each target domain image as their target domain sample set. And the knowledge of source domain samples is transferred to target domain samples via the transfer learning method. Finally, a random forest classifier‑based water body extraction model is trained using the target domain sample set, and is used for the water body extraction of target domain images. Results: We have conducted experiments using six PolSAR images and two kinds of transfer learning methods, the results show that:(1) The label transfer accuracy and the water body extraction accuracy are positively correlated. (2) Inductive transfer learning methods achieve higher label transfer accuracy and lower standard deviation. (3) A smaller distribution difference between source and target domain images indicate a greater transferability, and thus a better water body extraction accuracy. (4) The water body extraction results of inductive transfer learning methods have a higher rate of missing detection, while the results of transductive transfer learning methods have a higher rate of false detection. Conclusions: In the water body extraction of multi‑temporal PolSAR images, the use of transfer learning methods can significantly reduce the number of manually labeled samples needed to construct high‑performance classifiers, while maintaining the water body extraction accuracy at a high level. It has great application potentiality in the emergency response of flood disaster. [ABSTRACT FROM AUTHOR]

Objectives: Terrestrial laser scanning (TLS), ground‑based interferometric synthetic aperture radar (GB‑InSAR) and unmanned aerial vehicle(UAV) photography have been applied to emergency deformation monitoring of rockslide. The objectives are concerns about new problems of 3D visualization of deformation data. Methods: Firstly, we introduced iterative closest point method (ICP) to complete the registration between the TLS point cloud and UAV terrain model. Secondly, geometric mapping method was applied to achieve data fusion between GB‑InSAR deformation map and TLS point cloud. The mapping deviation is difficult to be corrected by conventional methods, because the lack of artificial control points on the rockslide. Visual interpretation and GB‑InSAR images simulation method were combined to extract the control points. The spatial coordinate transformation parameters were estimated with the least square method. The transformation model was established to correct the mismatch. Results: The proposed method was verified by the simulation data and the actual monitoring data of a landslide. The result shows that the accuracy reaches the pixel level which meet the needs of emergency monitoring. Conclusions: The proposed method is limited by the working conditions of emergency and the subjectivity of image interpretation. If the accuracy of control point selection was not narrowed, the final matching fusion position accuracy may reach pixel‑level accuracy. [ABSTRACT FROM AUTHOR]

Objectives: On January 8th 2022, an Ms 6.9 earthquake struck Menyuan, Qinghai Province, China. The epicenter is located in the middle section of Qilian-Haiyuan fault in the northeast margin of the Tibetan Plateau, which is a seismic gap. It is of great significance to study this earthquake for under-standing the tectonic deformation mechanism, strain release process and seismic hazards of the NE Tibetan Plateau.Methods: Differential synthetic aperture radar interferometry (D-InSAR) is used to obtain the coseismic deformation along the line of sight (LOS) of the Menyuan earthquake from Sentinel-1 satellite images. The co-registration of optically sensed images and correlation are used to obtain the EW coseismic deformation from the multi-temporal Sentinel-2 data. Based on the coseismic deformation field derived from InSAR interferograms, we invert the slip distribution. Results: The LOS deformation fields showed that the Menyuan earthquake caused a surface rupture of about 20 km long and the maximum LOS displacement was 0.75 m. The maximum EW displacement derived from optical images was 2.5 m. The slip distribution model showed that the Menyuan earthquake was dominated by a sinistral strike-slip. The coseismic slip is concentrated at a depth of 0 - 10 km. The maximum slip was 3.25 m at a depth of 4.89 km and the rake angle was 10.44°. The estimated total seismic moment is 1.07x1019 N.m, corresponding to a magnitude of Mw 6.6. Conclusions: According to the results from aftershock relocation, field investigation and tectonic background, we infer that the seismogenic fault is the Lenglongling fault and the eastern segment of Tuolaishan fault is also ruptured. The stress on the eastern segment of the Lenglongling fault and the western segment of the Tuolaishan fault is enhanced, implying that there is the possibility of strong earthquakes in this region in the future. [ABSTRACT FROM AUTHOR]

Objectives: Two earthquakes with a distance less than 40 km occurred in Menyuan, Qinghai Province, China, including Mw 5.9 earthquake on 26th January 2016 and Mw 6.7 earthquake on 8th January 2022. Methods: The interferometric synthetic aperture radar (InSAR) technique was utilized to process Sentinel-1A radar images from both ascending and descending tracks to retrieve high resolution coseismic surface displacement maps for these two Menyuan earthquakes. The coseismic surface displacements were then inverted to determine the fault geometry parameters and the non-uniform fault slip distributions with a dislocation model in an elastic half- space. Results: The results show that the 2016 Menyuan earthquake was a reverse event on a buried fault with a maximum surface displacement of 6.7 cm and 7.0 cm on the ascending and descending tracks, respectively. The maximum slip on the fault plane reached 0.53 m and concentrated at depths of 4 − 12 km. The 2022 Menyuan earthquake ruptured along the NWW - SEE strike with a maximum surface displacement of 78 cm in the satellite radar line of sight and a maximum fault slip of 3.5 m occurred at a depth of approximately 4 km. Conclusions: Considering the tectonic settings around the Lenglongling fault zone, the 2022 Menyuan earthquake is likely to occure in the western segment of the left‐lateral Lenglongling fault, with its western end extending westward to the Tuolaishan fault. By calculating the coulomb failure stress changes (CFS), it is inferred that the 2016 Menyuan earthquake generated considerable positive CFS which might have triggered the 2022 Menyuan earthquake. [ABSTRACT FROM AUTHOR]

Objectives: Integrated remote sensing application based on the space-air-ground investigation system has been promoted nationwide and has achieved certain results. However, as the work progresses, some urgent problems have been exposed, such as how to define hidden hazards and what types of them, what are the targeted identification technologies, and what is the focus of the next step, etc. Methods: By summarizing a large number of practical cases, the definition of hidden landslide hazards is further clarified. Potential landslides are divided into three categories: Deformation area, historical deformation failure area and potentially unstable slope, and we propose targeted identification techniques for different types. Results: The deformation zone is mainly identified by combining multi-temporal high resolution optical images and synthetic aperture radar interferometry technology (InSAR), high-resolution optical remote sensing images and light detection and ranging (LiDAR) can be used in the historical deformation and destruction region. Conclusions: Potential unstable slopes are difficult to identify by remote sensing and traditional artificial investigations. It is necessary to vigorously develop aerial and semi-aerial geophysical prospecting technology to identify the underground structure of the slope quickly, especially the spatial distribution of the base-cover interface and the groundwater status, and delineate dangerous areas through quantitative analysis of stability. [ABSTRACT FROM AUTHOR]

Objectives: In the research of time-series interferometric synthetic aperture radar (TS-InSAR) deformation monitoring, the conventional triangular reflectors are increasingly used as point targets. But the conventional triangular reflectors cannot meet the requirements of obtaining multi-angle deformation information on ascending and descending track scenes. They need to be tailored for specific datasets to meet monitoring needs. Methods: A set schemes for corner reflector design and accuracy verification are constructed. According to the characteristics of X-band radar satellite satellite, a two-way adjustable rectangular corner reflector is designed to meet the monitoring scene of the ascending and descending track. Due to the relationship between the radar cross section(RCS) and the size of the corner reflector, the side length of the bidirectional rectangular corner reflector is calculated to be 0.35 m. Ten bidirectional rectangular corner reflectors are placed on the side slope of the highway in the monitoring area. By artificially simulating the two-dimensional deformation, the detection accuracy and reliability of the TS-InSAR on the bidirectional rectangular corner reflector is verified. Results: Under the condition of effectively eliminating the sidelobe effect of the pixel, the designed bidirectional rectangular corner reflector maintains a high and stable backscatter ing intensity, which can be easily extracted. Through the analysis of the weather conditions at the time of imaging, heavy rainfall has a certain influence on the backscattering intensity of corner reflectors. Under the blind experimental conditions of artificially adjusting the height of the corner reflector, CR1 and CR2 are moved between May 13, 2018 and May 24, 2018. The root mean square error is 0.6 mm. The deformation detection accuracy reaches the sub-millimeter level. Conclusions: The TS-InSAR is highly sensitive to the deformation of the designed bidirectional rectangular corner reflector. The bidirectional rectangular corner reflector can meet the requirements of obtaining multi-angle deformation information on ascending and descending track scenes. [ABSTRACT FROM AUTHOR]

Objectives: The mountain excavation and city construction in Yan’an New District is one of the largest land‑creation projects in the Loess Plateau. Large‑scale mountain excavation has induced extensive land uplift in the region. At the same time, the excavation area is the main site for the urbanization of the new district, and the land uplift has seriously affected the sustainable development of urbanization in the new district. Methods: In this study, the small baseline subset synthetic aperture radar interferometry (SBAS‑InSAR) technology was utilized to process the ascending Sentinel‑1A data stacks covering the new district from May 2016 to October 2019 to retrieve long ‑ term ground deformation information after construction. The spatial pattern of ground uplift and its time deformation evolution characteristics was revealed, and the causes of land uplift were further quantitatively analyzed in detail. Results: The results show that significant land uplift of the new district after construction is mainly distributed in the excavation area that has the greatest impact on human engineering activities, with a maximum uplift rate of 24.8 mm/a. The ground lifting deformation develops rapidly in the early post‑construction period, and it tends to slow down and stabilize over time.Conclusions: The unloading rebound deformation caused by the rapid mountain excavation is the main internal factor causing land uplifting, and the distribution and magnitude of the land uplifting are the result of the combined effects of external factors such as excavation engineering, excavation thickness, and building load. [ABSTRACT FROM AUTHOR]

Objectives: Yan’an New District has attracted much attention because of its large‑scale mountain excavation and city construction project. Such mega‑scale land creation, rapid urban construction and complex geological conditions has induced wide land subsidence in the region. Methods: In this paper, the temporal and spatial evolution of land subsidence in Yan’ an New District was assessed using small baseline subset synthetic aperture radar interferometry (SBAS‑InSAR) to process the Sentinel‑1A images covering the period from May 2016 to October 2019. Then, the mechanism and evolution trend of land subsidence in the new district was analyzed. Results: Experimental results suggest that the spatial evolution pattern of land subsidence in Yan’an New District is highly relevant with the process of land creation project. The land subsidence area in the new district gradually expanded from the major urban area to the area of forest along with the land creation project, and its area decreased gradually with the passing of time. The temporal deformation evolution of land subsidence in the new district is related to the thickness of the fill, and all exhibit non‑linear subsidence with various velocities. The temporal evolution of land subsidence undergoes three processes of rapid‑slow‑steady. The greater the thickness of the fill, the larger the amount of sedimentation, the longer it takes for the sedimentation to stabilize. Conclusions: The main internal mechanism of land subsidence is the consolidation and compression deformation of the loess fill, which can be divided into three stages: Instantaneous settlement, consolidation settlement, and secondary consolidation settlement. Based on the characteristics of different stages combined with time series deformation, the land subsidence process and development trend can be analyzed. The research results provide a scientific reference for further monitoring and early warning, urban planning, and prevention and control of subsidence disasters in the new district. [ABSTRACT FROM AUTHOR]

Objectives: On 8th August 2017, a catastrophic earthquake of Ms 7.0 struck the County of Jiu-zhaigou in Sichuan Province, China. Because of its high magnitude and shallow epicenter, the earthquake caused grave casualties and property losses. Furthermore, the earthquake triggered numerous secondary mountain disasters such as rockfall and landslide. We use the synthetic aperture radar interferometry (In-SAR) technique to detect post-earthquake landslides in Jiuzhaigou area. Methods: Firstly, we carry out a quick detection across wide area using differential InSAR (DInSAR) technique with 6 ALOS-2 PALSAR-2 ascending images. DInSAR technique can observe the subtle deformation of landslide surface during the pe-riod between two SAR observations, so the active deforming landslides can be efficiently distinguished from stable areas. Nevertheless, the performance of DInSAR is often limited by temporal and spatial decor-relations induced by dense vegetation coverage, etc. To retrieve the temporal evolution of these landslides, detailed monitoring of specific landslides is carried out using time-series InSAR technique with 112 Sentinel-1 ascending and descending images. Time-series InSAR technique can overcome the limitation of tradi-tional DInSAR technique in spatial-temporal decorrelation, and improve the deformation measurement accuracy. Therefore, the combination of DInSAR and time-series InSAR is adopted to accurately monitor sin-gle landslide, and efficiently detect landslides in Jiuzhaigou area. Results: The results show that there are 13 landslides that can be detected by InSAR analyses in Jiuzhaigou area, including 7 landslides close to residential areas. The time series displacements retrieved indicate that several landslides have long-term linear deformations after the 2017 earthquake. These landslides severely threaten people's life and property and need to be monitored regularly. Finally, the reliability of detection results based on InSAR technology is verified by a field survey. Conclusions: The results can provide key supports for the geohazard mitigation and prevention in the Jiuzhaigou scenic area. We indicate that a strategy which combines DInSAR technique with time-series InSAR can efficiently and accurately detect and monitor landslides over a wide area. [ABSTRACT FROM AUTHOR]

In order to promote the wider use of ground-based synthetic aperture radar (GB-SAR) in the field of surface micro-deformation monitoring, this paper summarizes the main types of GB-SAR and introduces the working principles and important parameters of some representative systems. By taking the S-SAR (slope-synthetic aperture radar) ground-based radar independently developed by the China Academy of Safety and Production Science as an example，the key processing techniques utilized by the GB-SAR are introduced. Experiments verify the accuracy of GB-SAR. It also shows the application of GB-SAR in open pit mine slope monitoring，mine dumping site and landslide emergency rescue. Finally，the problems and suggestions of GB-SAR in monitoring are proposed. [ABSTRACT FROM AUTHOR]

The theory and method of mining subsidence monitoring based on D-InSAR(differential interferometric synthetic aperture radar) technology is a hot topic in the current field of mine deformation monitoring. At present, three-dimensional monitoring model of mining subsidence has been constructed successfully by some scholars through integrating single line of sight to D-InSAR technology and mining subsidence law. However, the model is only suitable for the researches on mining subsidence of horizontal coal seam rather than inclined coal seam. In view of the above defects，according to the relationship between deformation of line of sight monitored by D-InSAR and three-dimensional deformation, considering the surface subsidence law of inclined coal seam mining, the surface subsidence monitoring method for inclined coal seam based on single line of sight to D-InSAR technology is proposed. The simulation experiment shows that the absolute relative error of the subsidence monitoring method is less than 9. 53%, with an average of 1. 31%. The absolute relative error of the horizontal movement monitoring along North-South (SN) and East-West (EW) is less than 9. 78%，with an ave-rage of 3. 71%, which meets the mining subsidence monitoring accuracy requirements. The surface movement and deformation caused by the mining of the 9310 working face of Nantun coal mine in yanzhou mining area, Shandong province China during 2012-01-27 - 2012-02-07 are successfully monitored by this method in this paper [ABSTRACT FROM AUTHOR]