Search

Understanding the characteristics of the raindrop size distribution (DSD) is crucial to improve our knowledge of the microphysical processes of precipitation and to improve the accuracy of radar quantitative precipitation estimation (QPE). In this study, the spatial variability of DSD in different regions of Beijing and its influence on radar QPE are analyzed using 11 disdrometers. The DSD data are categorized into three regions: Urban, suburban, and mountainous according to their locations. The DSD exhibits evidently different characteristics in the urban, suburban, and mountain regions of Beijing. The average raindrop diameter is smaller in the urban region compared to the suburban region. The average rain rate and raindrop number concentration are lower in the mountainous region compared to both urban and suburban regions. The difference in DSD between urban and suburban regions is due to the difference in DSD for the same precipitation types, while the difference in DSD between mountain and plains (i.e., urban and suburban regions) is the combined effect of the convection/stratiform ratio and the difference of DSD for the same precipitation types. Three DSD-based polarimetric radar QPE estimators were retrieved and estimated. Among these three QPE estimators, R(ZH), R(Kdp), and R(Kdp, ZDR), R(Kdp, ZDR) performs best, followed by R(Kdp), and R(ZH) performs worst. R(Kdp) is more sensitive to the representative parameters, while R(ZH) and R(Kdp, ZDR) are more sensitive to observational error and systematic bias (i.e., calibration). [ABSTRACT FROM AUTHOR]

In order to solve the current problems that conventional video inspection can only detect, as an internal pipeline defect and drainage pipeline radar inspection device detects in a single direction and at radar frequency in water pipeline defect detection, a three-channel drainage pipeline ground penetrating radar (GPR) inspection device was designed and developed, the assembly and commissioning of the device prototype were completed, and an actual engineering test application was carried out. Focusing on the problem that the detection direction and depth of the single-channel detection device are limited, a three-channel drainage pipeline GPR inspection device is designed to realize the synchronous detection of the inside of the pipeline, the pipeline body, and the external environment of the pipeline, improving the detection depth and efficiency. According to the design scheme of the three-channel drainage pipeline GPR inspection device, the assembly of the device prototype was completed. The device contains three radar channels, the top of the main frequency of the antenna is 1.4 GHz, the two sides are 750 MHz, the video camera has a pixel count of 4 million, and the positioning accuracy is less than 1 mm, the waterproof grade is IP68, the detection accuracy of pipe deformation (slope) is 0.1°, the detection depth outside the pipe is 1.2 m, and the detection accuracy of corrosion thickness is 15 mm. In a practical application of the device, the Jianguomenqiao sewage pipeline in Beijing, China, was tested, resulting in the discovery of 87 defects, including 39 loose soil areas at the bottom of the pipe exterior, 40 void areas, and 8 cavities. [ABSTRACT FROM AUTHOR]

Severe hailstorms frequently occurred in Beijing between May and August 2021, leading to extensive hail damage. These hailstorms were observed by radar and satellite data, and reported by surface observers. In this study, the spectral and cloud microphysical characteristics of typical Beijing events in 2021 were analyzed using Himawari-8 satellite products and ground-based S-band weather radar data obtained from the Beijing Meteorological Bureau. The relationship between Himawari-8 brightness temperature differences (BTD) and radar reflectivity was also investigated. The results revealed that the significant spectral depression of brightness temperatures (BTs) in hail clouds was observed by a satellite. Furthermore, the stronger the radar reflectivity was, the more rapidly BTD decreased, with a nonlinear relationship between them. The results of cloud physical characteristics show that, for cloud-top heights above 10 km, the cloud effective radius was about 25 μm, with a cloud-top temperature of 225 K during these hail events. By means of Gaussian fitting, the BT threshold value (11.2 μm) was determined by satellite at 230 K, with a BTD focused on 1.9 K when hailstorms occurred. These results will help us better understand the characteristics of hailstorms, while also providing information for future hail suppression in Beijing. [ABSTRACT FROM AUTHOR]

The performance of general circulation models (GCMs) in simulating horizontal winds is important because the distribution and variation in horizontal winds are central to investigating atmospheric dynamic characteristics and processes. Also, horizontal wind data can be used to extract some of the required information on gravity waves, tides, and planetary waves. In this context, the present paper evaluates the capability of the Institute of Atmospheric Physics atmospheric general circulation model high-top version (IAP-HAGCM) in simulating the horizontal winds and tides of the troposphere and lower stratosphere by presenting a climatological and statistical comparison against observations of the powerful Beijing mesosphere–stratosphere–troposphere (MST) radar (39.78°N, 116.95°E) during 2012–2014. The results illustrated that the IAP-HAGCM can successfully reproduce the time–altitude distribution of the monthly mean zonal wind and diurnal tide amplitude, albeit with some underestimation. The mean correlation coefficients and root-mean-square error for the zonal (meridional) winds were 0.94 (0.73) and 6.60 m s−1 (2.90 m s–1), respectively. Additionally, the IAP-HAGCM can capture the temporal variation in both the zonal and meridional winds. It is worth noting that, compared with the seven coupled model intercomparison project phase 6 (CMIP6) models, the IAP-HAGCM performs better in meridional wind simulations below 15 km. However, there are discrepancies in altitudinal ranges with large wind velocities, such as the westerly jet, in the transition region of the troposphere and stratosphere, and in February, April, July, and September. It is suggested that model users should take advantage of the model's simulation ability by combining this information regarding when and where it is optimal with their own research purposes. Moreover, the evaluation results in this paper can also serve as a reference for guiding improvements of the IAP-HAGCM. [ABSTRACT FROM AUTHOR]

We studied the characteristics of quasi-two-day wave (QTDW) using the meridional wind in the mesosphere and lower thermosphere (MLT) obtained from a meteor radar over Kototabang (KB, 0.2 ° S, 100.3 ° E) from 2003 to 2012. Atmospheric oscillations have a crucial impact on atmospheric dynamics, which contributes to more accurate space weather forecasting, thus providing a more secure space environment for human space exploration activities such as remote sensing and satellite navigation. QTDWs are typical atmospheric oscillations in the upper stratosphere, mesosphere and lower thermosphere. The occurrence time, amplitudes, periods and vertical wavelengths of QTDW events are analyzed statistically. Data obtained from the TIMED Doppler Interferometer (TIDI), which can measure wind and temperature and is onboard the Thermosphere, Ionosphere, Mesosphere Energetics and Dynamics (TIMED) satellite, are used to analyze the global distribution and spatial structure of QTDWs with different zonal wavenumbers. The characteristics of the QTDWs over KB are compared with the QTDWs at the middle latitudes using the meridional wind data from a meteor radar over Wuhan (114.4 ° E, 30.6 ° N), Beijing (116.5 ° E, 39.9 ° N) and Mohe (121.1 ° E, 50.1 ° N). The amplitudes of the QTDW and spectral analysis are calculated by the least squares fitting method. Our results demonstrate that QTDWs are present almost all year around over KB. The occurrence time, amplitudes, periods and vertical wavelengths of QTDW events with different zonal wavenumbers are determined in this study. We also find that the statistical characteristics of the QTDWs in KB are different from those at middle latitudes. The westward zonal wavenumber −4 (W4) events gradually increase with increasing latitude, whereas westward zonal wavenumbers −1, −2, and −3 (W1, W2 and W3, respectively) events all decrease with increasing latitude. [ABSTRACT FROM AUTHOR]

InSAR (Interferometric Synthetic Aperture Radar) is widely recognized as a crucial remote sensing tool for monitoring various geological disasters because it provides all-day and all-weather monitoring. Nevertheless, the current interpretation methods for InSAR heavily depend on the interpreter's experience, which hinders efficiency and fails to meet the requirements for the timely detection of geologic hazards. Furthermore, the results obtained through current InSAR processing carry inherent noise interference, further complicating the interpretation process. To address those issues, this paper proposes an approach that enables automatic and rapid identification of deformation zones. The proposed method leverages IPTA (Interferometric Point Target Analysis) technology for SAR data processing. It combines the power of HNSW (Hierarchical Navigable Small Word) and DBSCAN (Density-Based Spatial Clustering of Applications with Noise) algorithms to cluster deformation results. Compared with traditional methods, the computational efficiency of the proposed method is improved by 11.26 times, and spatial noise is suppressed. Additionally, the clustering results are fused with slope units determined using DEM (Digital Elevation Model), which facilitates the automatic identification of slopes experiencing deformation. The experimental verification in the western mountainous area of Beijing has identified 716 hidden danger areas, and this method is superior to the traditional technology in speed and automation. [ABSTRACT FROM AUTHOR]

Accurately estimating building heights is crucial for various applications, including urban planning, climate studies, population estimation, and environmental assessment. However, this remains a challenging task, particularly for large areas. Satellite-based Light Detection and Ranging (LiDAR) has shown promise, but it often faces difficulties in distinguishing building photons from other ground objects. To address this challenge, we propose a novel method that incorporates building footprints, relative positions of building and ground photons, and a self-adaptive buffer for building photon selection. We employ the Ice, Cloud, and Land Elevation Satellite 2 (ICESat-2) photon-counting LiDAR, specifically the ICESat-2/ATL03 data, along with building footprints obtained from the New York City (NYC) Open Data platform. The proposed approach was applied to estimate the heights of 17,399 buildings in NYC, and the results showed strong consistency with the reference building heights. The root mean square error (RMSE) was 8.1 m, and for 71% of the buildings, the mean absolute error (MAE) was less than 3 m. Furthermore, we conducted an extensive evaluation of the proposed approach and thoroughly investigated the influence of terrain, region, building height, building density, and parameter selection. We also verified the effectiveness of our approach in an experimental area in Beijing and compared it with other existing methods. By leveraging ICESat-2 LiDAR data, building footprints, and advanced selection techniques, the proposed approach demonstrates the potential to accurately estimate building heights over broad areas. [ABSTRACT FROM AUTHOR]

Revealing the land subsidence in Beijing, China, induced by the massive groundwater extraction in the past three decades, is important to mitigate the hazards and protect the residences and infrastructure. Many SAR (Synthetic Aperture Radar) datasets have been successfully applied to reveal the land subsidence over Beijing in previous research, while few works were achieved on land subsidence revealed by time-series InSAR (Interferometric Synthetic Aperture Radar) with Gaofen-3 SAR images. In this study, we successfully perform the time-series InSAR analysis with Gaofen-3 SAR images to extract the land subsidence in Beijing from 2020 to 2021. The Sentinel-1 SAR images were used to assess the accuracy of Gaofen-3 images. The subsidence scale and extent are consistent in detected major subsidence bowls between the two datasets. The spatial–temporal evolution and the deceleration of Beijing land subsidence were revealed by comparing with the Sentinel-1 results from 2017 to 2020. Moreover, we evaluated the interferometric performance of Gaofen-3 satellite SAR imagery and analyzed the main factors that mostly influence the coherence and quality of interferograms. Our results proved that the long perpendicular baselines decrease the coherence seriously over the study area, and the artifacts induced by inaccurate orbit information reduce the quality of the Gaofen-3 interferograms. Refining and removing the two main artifacts could improve the quality of interferograms formed by Gaofen-3 SAR images. [ABSTRACT FROM AUTHOR]

Accurate, rapid, and automatic local climate zone (LCZ) mapping is essential for urban climatology and studies in terms of urban heat islands. Remotely sensed imageries incorporated with machine learning algorithms are widely utilized in LCZ labeling. Nevertheless, large-scale LCZ mapping is still challenging due to the complex vertical structure of underlying urban surfaces. This study proposed a new method of LCZ labeling that uses a random forest classifier and multi-source remotely sensed data, including Sentinel 1A Synthetic Aperture Radar (SAR), Sentinel 2 Multispectral Instrument, and Luojia1-01 night-time light data. In particular, leaf-on and -off imageries and surface thermal dynamics were utilized to enhance LCZ labeling. Additionally, we systematically evaluated how daytime and night-time features influence the performance of the classification procedure. Upon examination, the results for Beijing, China, were confirmed to be robust and refined; the Overall Accuracy (OA) value of the proposed method was 88.86%. The accuracy of LCZs 1–9 was considerably increased when using the land surface temperature feature. Among these, the Producer Accuracy (PA) value of LCZ 3 (compact low-rise) significantly increased by 16.10%. Notably, it was found that NTL largely contributed to the classification concerning LCZ 3 (compact low-rise) and LCZ A/B (dense trees). The performance of integrating leaf-on and -off imageries for LCZ labeling was better than merely uses of leaf-on or -off imageries (the OA value increased by 4.75% compared with the single use of leaf-on imagery and by 3.62% with that of leaf-off imagery). Future studies that use social media big data and Very-High-Resolution imageries are required for LCZ mapping. This study shows that combining multispectral, SAR, and night-time light data can improve the performance of the random forest classifier in general, as these data sources capture significant information about surface roughness, surface thermal feature, and night-time features. Moreover, it is found that incorporating both leaf-on and leaf-off remotely sensed imageries can improve LCZ mapping. [ABSTRACT FROM AUTHOR]

It is commonly believed that clear-air echoes detected by weather radars are caused by atmobios migration. However, clear-air echoes are sometimes inconsistently related to the activity of living creatures. In some cases, the characteristics of radar products seem to conform to biological scattering, but the movement of echoes cannot be observed. For these reasons, we sought to expand the cause of clear-air echoes from a Chinese Doppler S-band Weather Radar (CINRAD/SA) in Beijing. Some contradictions were discovered in a case which diverged from previous conclusions. It was found that the progression and movement of clear-air echoes do not conform to the rules of biological activities. The frequency distribution of dual-wavelength ratio peaks is 21.5 dB, which is in accordance with Villars–Weisskopf's turbulence theory. From 1 May to 20 May, the 58% dual-wavelength ratio between the S-band and the X-band was distributed between 18 dB and 24 dB. These results show that more than half of the clear-air echoes of CINRAD/SA at night were caused by turbulence in Beijing. A new model of troposcatter propagation, the reflecting-layers model, was then introduced to explain the radar observations. According to the reflecting-layers model, the echoes' diurnal variation and reflectivity characteristics are attributed to the effects of turbulent mixing. Excessive turbulent mixing affects the generation of the reflective layer, thereby weakening the echo signal. It is necessary to re-examine the position of turbulence in clear-air echoes. [ABSTRACT FROM AUTHOR]

Surface subsidence is a serious threat to human life, buildings and traffic in Beijing. Surface subsidence is closely related to human activities, and human activities in Beijing area showed a decreasing trend during the Corona Virus Disease 2019 (COVID-19). To study surface subsidence in Beijing before and after the COVID-19 outbreak and its causes, a total of 51 Sentinel-1A SAR images covering Beijing from January 2018 to April 2022 were selected to derive subsidence information by Time Series Interferometry Synthetic Aperture Radar (TS-InSAR). The results of surface subsidence in Beijing demonstrate that Changping, Chaoyang, Tongzhou and Daxing Districts exhibited the most serious subsidence phenomenon before the COVID-19 outbreak. The four main subsidence areas form an anti-Beijing Bay that surrounds other important urban areas. The maximum subsidence rate reached −57.0 mm/year. After the COVID-19 outbreak, the main subsidence area was separated into three giant subsidence funnels and several small subsidence funnels. During this period, the maximum subsidence rate was reduced to −43.0 mm/year. Human activity decrease with the COVID-19 outbreak. This study effectively analysed the influence of natural factors on surface subsidence after excluding most of the human factors. The following conclusions are obtained from the analysis: (1) Groundwater level changes, Beijing's geological structure and infrastructure construction are the main reasons for surface subsidence in Beijing. (2) Seasonal changes in rainfall and temperature indirectly affect groundwater level changes, thereby affecting surface subsidence in the area. (3) The COVID-19 outbreak in early 2020 reduced the payload of Beijing's transportation facilities. It also slowed down the progress of various infrastructure construction projects in Beijing. These scenarios affected the pressure on the soft land base in Beijing and reduced the surface subsidence trend to some extent. [ABSTRACT FROM AUTHOR]

Foreign Object Debris (FOD) refers to any foreign material on the airfield that may injure and threaten the aircraft and airport system. Due to the complex background on the airfield pavement and weak target echoes in long-distance monitoring, it is not easy to detect objects of various types and sizes. The existing FOD radar system's detection method has a short effective range, and the detectable objects' radar cross-section intensity is no less than −20 dBsm. In this paper, we propose an integrated FOD automatic target detection algorithm for millimeter-wave (MMW) surveillance radar to improve small target detection under long-range conditions of over 660 m. The signal form of FOD and a clutter model of ground clutter received by millimeter-wave radar are primarily utilized and established theoretically. The runway edge detection means that it is employed based on the in-continuity features as the runway region of interest during the automatic extraction step. Following the clutter map constant false alarm detection algorithm, we utilize a time-domain algorithm that functions as the vital detection processor. Moreover, an explicit definition of the FOD detection performance is developed in a characteristic quantitative way. This criterion involves an absolute reference value for all FOD radar systems. The well-designed FOD frequency-modulated continuous-wave MMW surveillance radar is utilized, and actual experiments are carried out in a real airport in Beijing, China. The results validate the proposed method's effectiveness and the superior performance of FOD target detection in long-range situations. [ABSTRACT FROM AUTHOR]

This paper compares the wind fields measured by the meteor radar at Mohe, Beijing, Wuhan, and Sanya stations and horizontal wind model (HWM14) predictions. HWM14 appears to successfully reproduce the height-time distribution of the monthly mean zonal winds, although large discrepancies occur in wind speed between the model and measurement, especially in the summer and winter months. For meridional wind, the consistency between model prediction and radar observation is worse than that of zonal wind. The consistency between radar measurements and model prediction at Sanya station is worse than other sites located at higher latitudes. Harmonic analysis reveals large discrepancies in diurnal, semidiurnal, and terdiurnal tides extracted from meteor radar observations and HWM14 predictions. [ABSTRACT FROM AUTHOR]

The values of ratio a of the linear relationship between specific attenuation and specific differential phase vary significantly in convective storms as a result of resonance scattering. The best-linear-fit ratio a at X band is determined using the modified attenuation correction algorithm based on differential phase and attenuation, as well as the premise that reflectivity is unattenuated in S band radar detection. Meanwhile, the systemic reflectivity bias between the X band radar and S band radar and water layer attenuation (ZW) on the wet antenna cover of the X band radar are also considered. The good performance of the modified correction algorithm is demonstrated in a moderate rainfall event. The data were collected by four X band dual-polarization (X-POL) radar sites, namely, BJXCP, BJXFS, BJXSY, and BJXTZ, and a China's New Generation Weather Radar (CINRAD/SA radar) site, BJSDX, in Beijing on 20 July 2016. Ratio a is calculated for each volume scan of the X band radar, with a mean value of 0.26 dB deg−1 varying from 0.20 to 0.31 dB deg−1. The average values of systemic reflectivity bias between the X band radar (at BJXCP, BJXFS, BJXSY, and BJXTZ) and S band radar (at BJSDX) are 0, −3, 2, and 0 dB, respectively. The experimentally determined ZW is in substantial agreement with the theoretically calculated ones, and their values are an order of magnitude smaller than rain attenuation. The comparison of the modified attenuation correction algorithm and the empirical-fixed-ratio correction algorithm is further evaluated at the X-POL radar. It is shown that the modified attenuation correction algorithm in the present paper provides higher correction accuracy for rain attenuation than the empirical-fixed-ratio correction algorithm. [ABSTRACT FROM AUTHOR]

A Chinese Ka-band solid-state transmitter cloud radar (CR) was employed to investigate clouds and precipitation over the Tibetan Plateau (TP) in summer 2014. After performing quality control and merging the output from the three work modes of CR, the Doppler spectral density data were analyzed to estimate the vertical air motion (Vair) and particle fall speed (Vfall) in clouds. The reproduced radar reflectivity was employed to retrieve the ice water content (IWC) and ice effective radius (Re). The cloud type classification algorithm was developed to classify cloud conditions into eight different categories. The vertical and daily variations of reflectivity, Vair, Vfall, IWC, and Re were then analyzed, and as a result, four conclusions were formulated. First, the clouds with reflectivity less than −10 dBZ were mainly located above 7 km, and the precipitable particles were formed below this layer (rain echo top) due to the abundant supercooled water therein. Second, the weak updraft in the range of 2–4 km caused the high occurrence of weak clouds during 04:00–12:00 Beijing local time (BT), and the rain echo top did not exceed 7 km due to the existing downdraft in 4–6 km. After 14:00 BT, convective clouds quickly developed, and the averaged updraft, reflectivity, and the echo top notably increased. Third, particular cloud types commonly exhibited weak reflectivity, low IWC, and less vertical variations, while the others more obvious vertical variations with larger IWC and Re. Last, compared with the radar sensitivity and range from radar, Vair biases that were introduced by the Doppler spectra broadening due to turbulence, wind shear, and radar beam width, could be neglected. [ABSTRACT FROM AUTHOR]

Geological hazards often occur in mountainous areas and are sudden and hidden, so it is important to identify and assess geological hazards. In this paper, the western mountainous area of Beijing was selected as the study area. We conducted research on landslides, collapses, and unstable slopes in the study area. The surface deformation of the study area was monitored by multi-temporal interferometric synthetic aperture radar (MT-InSAR), using a combination of multi-looking point selection and permanent scatterer (PS) point selection methods. Random forest (RF), support vector machine (SVM), convolutional neural network (CNN), and recurrent neural network (RNN) models were selected for the assessment of geological hazard susceptibility. Sixteen geological hazard-influencing factors were collected, and their information values were calculated using their features. Multicollinearity analysis with the relief-F method was used to calculate the correlation and importance of the factors for factor selection. The results show that the deformation rate along the line-of-sight (LOS) direction is between −44 mm/year and 28 mm/year. A total of 60 geological hazards were identified by combining surface deformation with optical imagery and other data, including 7 collapses, 25 unstable slopes, and 28 landslides. Forty-eight of the identified geological hazards are not recorded in the Beijing geological hazards list. The most effective model in the study area was RF. The percentage of geological hazard susceptibility zoning in the study area is as follows: very low susceptibility 27.40%, low susceptibility 28.06%, moderate susceptibility 21.19%, high susceptibility 13.80%, very high susceptibility 9.57%. [ABSTRACT FROM AUTHOR]

Based on three years of data observed by Beijing MST radar between 2012 and 2014, this study gives the distribution characteristics of the turbulence energy dissipation rate ε and the vertical turbulence diffusion coefficient K z in the troposphere–lower stratosphere over Beijing based on the spectral width method. Additionally, from the perspective of atmospheric stability expressed by the Brunt–Väisälä frequency N and the gradient Richardson number Ri, the seasonal variations of turbulence parameters at different heights are analyzed. The results show that ε ( K z ) is distributed in the range of 10−4 to 10−1.5 m2 s−3 (100 to 101 m2 s−1), with median values of 10−3 to 10−2.5 m2 s−3 (100.4 to 100.75 m2 s−1). In the height range of −1 km to about +1 km from the tropopause, the profiles of ε have no noticeable seasonal difference from October to next May, and ε increases sharply with the slope of linear regression ∇log(ε)/∇h = 0.25 (m2 s−3) km−1. The seasonal variation of turbulence parameters has noticeable differences at different atmospheric layers. Furthermore, the atmospheric static/dynamic stability and turbulence intensity are the influencing factors of turbulence parameters. From the tropopause to 14 km, the values of ε and N both have a peak from October to next May, and this peak layer disappears from July to August. Within the range of 2–4 km from the ground, K z and the frequency of Ri < 0.25 have a significant positive correlation, with a Pearson correlation coefficient greater than 0.5, and both have large value layers in March to May and September to October. The turbulence intensity above 15 km is likely the main factor affecting the turbulence dissipation rate and vertical turbulence diffusion coefficient. [ABSTRACT FROM AUTHOR]

The millimeter-wave cloud radar, ceilometer, and disdrometer have been widely used to observe clouds and precipitation. However, there are some drawbacks when those three instruments are solely employed due to their own limitations, such as the fact that radars usually suffer from signal attenuation and ceilometers/disdrometers cannot provide measurements of the hydrometeors of aloft clouds and precipitation. Thus, in this paper, we developed an integrated technology by combining and utilizing the advantages of three instruments together to investigate the vertical structure and diurnal variation of warm clouds and precipitation, and the raindrop size distribution. Specifically, the technology consists of appropriate data processing, quality control, and retrieval methods. It was implemented to study the warm clouds and precipitation in South China during the pre-flood season of 2016. The results showed that the hydrometeors of warm clouds and precipitation were mainly distributed below 2.5 km and most of the rainfall events were very light with a rain rate less than 1 mm h−1, however, the stronger precipitation primarily contributed the accumulated rain amount. Furthermore, a rising trend of cloud base height from 1000 to 1900 BJT was found. The cloud top height and cloud thickness gradually increased from 1200 BJT to reach a maximum at 1600 BJT (Beijing Standard Time, UTC+8), and then decreased until 2000 BJT. Also, three periods of the apparent rainfall on the ground of the day, namely, 0400–0700 BJT, 1400–1800 BJT, and 2300–2400 BJT were observed. During three periods, the raindrops had wider size spectra, higher number concentrations, larger rain rates, and higher water contents than at other times. The hydrometeor type, size, and concentration were gradually changed in the vertical orientation. The raindrop size distributions of warm precipitation in the air and on the ground were different, which can be expressed by γ distributions N(D) = 1.49 × 104D−0.9484exp(−6.79D) in the air and N(D) = 1.875 × 103D0.862exp(−2.444D) on the ground, where D and N(D) denote the diameter and number concentration of the raindrops, respectively. [ABSTRACT FROM AUTHOR]

The long-term over-exploitation of groundwater has not only caused the compaction of aquifer thickness and surface deformation but has also further aggravated the loss of groundwater storage (GWS) in Beijing plain. The South-to-North Water Diversion Project (SNWDP) furnishes a new source of water for Beijing. By reviewing related studies, it was found that there are few studies on the realization of GWS estimation based on InSAR technology considering the lag effect. Therefore, in this study, firstly, the long-time series deformation characteristics of Beijing plain were obtained from 46 RadarSat-2 images using small baseline subset interferometric synthetic aperture radar (SBAS-InSAR) technology. Secondly, the seasonal components of surface deformation and hydraulic head change were extracted by means of multichannel singular spectrum analysis (MSSA), verifying the separation accuracy by means of Monto Carlo-SSA (MC-SSA). Finally, for the hydrodynamic delay (aquifer water supply/drainage) of the complex aquifer system, we introduced the time lag cross-correlation (TLCC) approach to correct the hysteresis response of seasonal deformation relative to the variation of the aquifer system head, so as to realize the estimation of aquifer storage properties and GWS loss, even unrecoverable GWS (UGWS). The results showed that the average annual variation of total GWS (TGWS) in Beijing plain was −6.702 × 107 m3, of which the depletion volume of UGWS was −6.168 × 107 m3, accounting for 92.03% of the TGWS. On a temporal scale, the depletion of UGWS lagged behind the total head change, with about one year of lag time. On a spatial scale, in contrast to the north of Beijing plain, the depletion of UGWS in the south only recovered briefly after 2015 and then continued to decline. This further indicated that the process of the decline of middle-deep confined head and long-term GWS loss caused by over-exploitation of groundwater was irreversible. These findings are of great significance to optimize the allocation of groundwater resources, reduce the harm of land subsidence and protect groundwater resources. [ABSTRACT FROM AUTHOR]

Cloud top height (CTH) is an essential pareter for the general circulation model in understanding the impact of clouds on the Earth's radiation budget and global climate change. This paper compares the CTH products, derived from the Moderate Resolution Imaging Spectroradiometer (MODIS), onboard the Aqua and Terra satellites with ground-based Ka band radar data in Beijing from 2014 to 2017. The aim was to investigate the data accuracy and the difference in CTH measurements between passive satellite data and active ground-based radar data. The results show that MODIS, on average, underestimates CTH relative to radar by −1.08 ± 2.48 km, but with a median difference of −0.65 km and about 48% of differences are within 1 km. Statistically, MODIS CTHs which are greater than 6 km show lower discrepancy to radar CTH than those of MODIS CTHs less than 4 km. The CTH difference is independent of cloud fraction and cloud layer. It shows strong dependence on cloud depth, decreasing as cloud depth increases. There is a tendency for MODIS to underestimate high thin clouds but overestimate low thin clouds relative to radar. Total ozone, SO2, CO, NO2, aerosol PM10, total water vapor and temperature inversion show unobvious influences in the CTH discrepancy. It is shown that the MODIS CO2-slicing technique performs much better than IRW (infrared window) technique when cloud layer is higher than 2 km. The average difference calculated from all comparisons by CO2-slicing technique and IRW technique is 0.09 ± 1.58 km, and −2.20 ± 2.73 km, respectively. [ABSTRACT FROM AUTHOR]

Vertical wind shear (VWS) is one of the key meteorological factors in modulating ground-level particulate matter with an aerodynamic diameter of 2.5 µm or less (PM2.5). Due to the lack of high-resolution vertical wind measurements, how the VWS affects ground-level PM2.5 remains highly debated. Here we employed the wind profiling observations from the fine-time-resolution radar wind profiler (RWP), together with hourly ground-level PM2.5 measurements, to explore the wind features in the planetary boundary layer (PBL) and their association with aerosols in Beijing for the period from December 1, 2018, to February 28, 2019. Overall, southerly wind anomalies almost dominated throughout the whole PBL or even beyond the PBL under polluted conditions during the course of a day, as totally opposed to the northerly wind anomalies in the PBL under clean conditions. Besides, the ground-level PM2.5 pollution exhibited a strong dependence on the VWS. A much weaker VWS was observed in the lower part of the PBL under polluted conditions, compared with that under clean conditions, which could be due to the strong ground-level PM2.5 accumulation induced by weak vertical mixing in the PBL. Notably, weak northbound transboundary PM2.5 pollution mainly appeared within the PBL, where relatively small VWS dominated. Above the PBL, strong northerlies winds also favored the long-range transport of aerosols, which in turn deteriorated the air quality in Beijing as well. This was well corroborated by the synoptic-scale circulation and backward trajectory analysis. Therefore, we argued here that not only the wind speed in the vertical but the VWS were important for the investigation of aerosol pollution formation mechanism in Beijing. Also, our findings offer wider insights into the role of VWS from RWP in modulating the variation of PM2.5, which deserves explicit consideration in the forecast of air quality in the future. [ABSTRACT FROM AUTHOR]

The coordinated development of the Beijing-Tianjin-Hebei has become a national strategy with Beijing and Tianjin as twin engines driving the regional development. However, the Beijing-Tianjin-Hebei region has suffered dramatic ground subsidence during last two to three decades, mainly due to long-term groundwater withdrawal. Although, annual spirit leveling has been conducted routinely in some parts of Beijing and Tianjin, and InSAR technique has also been used to monitor ground subsidence in some local areas of the region, there is a lack of a complete survey of ground subsidence over the whole region. In this paper, we report a research on mapping ground subsidence in the Beijing-Tianjin-Hebei region over a long time span from 1992 to 2014. Three SAR datasets from four satellites are used: ERS-1/2 SAR images from 1992 to 2000, ENVISAT ASAR images from 2003 to 2010, and RADARSAT-2 images from 2012 to 2014. An improved multi-temporal InSAR method, namely "Multiple-master Coherent Target Small-Baseline InSAR" (MCTSB-InSAR), has been developed to process the datasets. A unique feature of MCTSB-InSAR is the adjustment process useful for wide area monitoring which provides an integrated solution for both calibration of InSAR-derived deformation and the harmonization of the deformation estimates from overlapping SAR frames. Three maps of the subsidence rate corresponding to the three periods over the wide Beijing-Tianjin-Hebei region are generated, with respective accuracy of 8.7 mm/year (1992-2000), 4.7 mm/year (2003-2010), and 5.4 mm/year (2012-2014) validated by more than 120 leveling measurements. The spatial-temporal characteristics of the development of ground subsidence in Beijing and Tianjin are analyzed. This research represents a first-ever effort on mapping ground subsidence over very large area and over long time span in China. The result is of significance to serve the decision-making on ground subsidence mitigation in the Beijing-Tianjin-Hebei region. [ABSTRACT FROM AUTHOR]

Beijing-Tianjin-Hebei (BTH) has been suffering from severe groundwater storage (GWS) consumption and land subsidence (LS) for a long period. The overexploitation of groundwater brings about severe land subsidence, which affects the safety and development of BTH. In this paper, we utilized multi-frame synthetic aperture radar datasets obtained by the Rardarsat-2 satellite to monitor land subsidence's temporal and spatial distribution in the BTH from 2012 to 2016 based on multi-temporal interferometric synthetic aperture radar (MT-InSAR). In addition, we also employed the Gravity Recovery and Climate Experiment (GRACE) mascon datasets acquired by the Center for Space Research (CSR) and Jet Propulsion Laboratory (JPL) to obtain the GWS anomalies (GWSA) of BTH from 2003 to 2016. Then we evaluate the accuracy of the results obtained. Furthermore, we explored the relationship between the regional GWSA and the average cumulative subsidence in the BTH. The total volume change of subsidence is 59.46% of the total volume change of groundwater storage. Moreover, the long-term decreasing trend of the GWSA (14.221 mm/year) and average cumulative subsidence (17.382 mm/year) show a relatively high consistency. Finally, we analyze the heterogeneity of GWS change (GWSC) and LS change (LSC) in the four typical areas by the Lorenz curve model. The implementation of the South-to-North Water Diversion Project (MSWDP) affects the heterogeneity of GWSC and LSC. It can be seen that the largest heterogeneity of LSC lags behind the GWSC in the Tianjin-Langfang-Hengshui-Baoding area. The largest uneven subsidence in Beijing and Tianjin occurred in 2015, and the largest uneven subsidence in Hengshui-Baoding occurred in 2014. After that, the heterogeneity of subsidence gradually tends to stable. [ABSTRACT FROM AUTHOR]

Based on long-term lidar (light detection and ranging) observations at Yanqing (40.5°N, 116°E) and Pingquan (41°N, 118.7°E), cloud events occurred in the vicinity of the stratopause above Beijing were reported for the first time. These events occurred with tenuous and sparse layers within the altitude range of 33–65 km, and the maximum VBSC value ranged from 1 × 10 − 10 m − 1 sr − 1 to 5.5 × 10 − 9 m − 1 sr − 1 . Considering temperature and water vapor measurements from SABER/TIMED, the occurrence mechanism of these lidar-observed cloud events was examined. It was found that some cloud layers resulted from the nucleation of water vapor due to the local meteorological changes in the middle atmosphere, while other lidar-observed clouds could comprise floating clusters of cosmic dust, hydrate droplets, volcanic ash, space traffic exhaust, etc. These cloud events are rare cloud-like phenomena in the middle atmosphere observed by lidars at midlatitudes in China; they differ from NLCs and PSCs in terms of altitude distribution and seasonal variation, and the relevant microphysics processes behind their occurrence are likely meaningful to meteorology at midlatitudes. [ABSTRACT FROM AUTHOR]

Beijing's One Million Mu Plain Afforestation Project involves planting large areas with the exotic North American tree species Fraxinus pennsylvanica Marsh (ash). As an exotic tree species, ash is very vulnerable to infestations by the emerald ash borer (EAB), a native Chinese wood borer pest. In the early stage of an EAB infestation, attacked trees show no obvious sign. Once the stand has reached the late damage stage, death occurs rapidly. Therefore, there is a need for efficient early detection methods of EAB stress over large areas. The combination of unmanned aerial vehicle (UAV)-based hyperspectral imaging (HI) with light detection and ranging (LiDAR) is a promising practical approach for monitoring insect disturbance. In this study, we identified the most useful narrow-band spectral HI data and 3D LiDAR data for the early detection of EAB stress in ash. UAV-HI data of different infested stages (healthy, light, moderate and severe) of EAB in the 400–1000 nm range were collected from ash canopies and were processed by Partial Least Squares–Variable Importance in Projection (PLS-VIP) to identify the maximally sensitive bands. Band R678 nm had the highest PLS-VIP scores and the most robust classification ability. We combined this band with band R776 nm to develop an innovative normalized difference vegetation index (NDVI(776,678)) to estimate EAB stress. LiDAR data were used to segment individual trees and supplement the HI data. The new NDVI(776,678) identified different stages of EAB stress, with a producer's accuracy of 90% for healthy trees, 76.25% for light infestation, 58.33% for moderate infestation, and 100% for severe infestation, with an overall accuracy of 82.90% when combined with UAV-HI and LiDAR. [ABSTRACT FROM AUTHOR]

In urban areas, deformation of transportation infrastructures may lead to serious safety accidents. Timely and accurate monitoring of the structural deformation is critical for prevention of transportation accidents and assurance of construction quality, particularly in areas with regional land subsidence, such as the city of Beijing. In this study, we proposed a method for the detection of seasonal deformation of highway overpasses using the integration of persistent scatterers Interferometric Synthetic Aperture Radar (PS-InSAR) techniques and seasonal indices, i.e., deformation concentration degree (DCD) and deformation concentration period (DCP) indices. Taking eastern Beijing urban area as a case study area, we first used the PS-InSAR technique to derive time series surface deformation based on 55 TerraSAR-X images during 2010–2016. Then, we proposed DCD and DCP indices to characterize seasonal deformation of 25 highway overpasses in the study area, with DCD representing to what degree the annual deformation is distributed in a year, and DCP representing the period on which deformation concentrates in the year. Our results showed that the maximum annual deformation rate reached −141.3 mm/year in Beijing urban area, and the PS-InSAR measurements agreed well with levelling measurements (R2 > 0.97). For PS pixels with DCD ≥ 0.3, the monthly deformation showed obvious seasonal patterns with deformation values during some months greater than those during the other months. DCP revealed that the settlement during autumn and winter was more serious than that in spring and summer. The seasonal patterns seemed to be related to the location, structure, and construction age of the overpasses. The upper-level overpasses, the newly constructed overpasses, and those located in the subsidence area (rate < −40 mm/year) tended to show a greater seasonal pattern. The seasonal deformation variations were also affected by groundwater-level fluctuation, temperature, and compressible layer. [ABSTRACT FROM AUTHOR]

Beijing is a major city suffering from land subsidence due to long-term over-exploitation of groundwater. The South-to-North Water Diversion Project (SNWDP), however, has had a significant impact on the structure of water consumption since the end of 2014, and it is changing the status of land subsidence in Beijing. In this study, we employed Persistent Scatterers Synthetic Aperture Radar Interferometric (PS-InSAR) to investigate the decadal evolution of land subsidence in Beijing with 100 TerraSAR-X stripmap images collected from April 2010 to December 2019. The water resources, historic climate and urban construction data were compiled for the years of 2010 to 2019 to analyze changes in groundwater level, human activity, surface geology, active faults and land subsidence patterns. The results show that the changes in the water supply structure are correlated to a rise in groundwater level after 2015. These changes include an increase in the water supply from the SNWDP, a reduction in groundwater exploitation, the optimization of water consumption, replacing recycled water for environmental water and a reduction in the use of water for agriculture. Land subsidence in the study area was concentrated in the eastern regions, trending towards a decreasing velocity starting about two years after the commencement of SNWDP in 2015. Uneven subsidence in the land subsidence area was related to excavations of underground soil, and the construction of Line 6 and Line 7 led to rapid nonlinear subsidence. Our results have scientific significance for reducing subsidence hazards in the context of SNWDP and urban expansion. [ABSTRACT FROM AUTHOR]

Since the 1930s, due to the rapid development of the city and the increase of population, the demand from Beijing residents for water resources has gradually increased. Land deformation in the Beijing Plain is a serious issue. In order to warn of, and mitigate, disasters, it is urgently necessary to obtain the latest rate, extent, and temporal evolution of land subsidence in Beijing. Firstly, the temporal and spatial distribution characteristics of land deformation in Beijing during 2003–2020 were unveiled using the time-series interferometric synthetic aperture radar (InSAR) technique and two different satellite datasets, sentinel-1a/1b and ENVISAT ASAR. By means of combining calibration of InSAR results with the global positioning system (GPS), we studied the evolutionary process of long-term land subsidence in Beijing. The precision of our InSAR annual subsidence results is less than 10 mm. Land subsidence in Beijing is unevenly distributed, and so five main land subsidence zones were monitored. The time-series results showed that the rate of land subsidence rate continued to increase from 2003 to 2015, but has gradually shown a slowing trend from 2015 to 2020. Further, we used the quadratic polynomial fitting method to interpolate the time-series deformation results from 2010 to 2015, and compared these with GPS. The results demonstrated that although the InSAR observation method is not strictly registered with GPS in time, its deformation trend is consistent. In addition, the calibrated long time-series was consistent with the three deformation stages of land subsidence evolution in Beijing. Finally, we analyzed the deformation information obtained by InSAR technology in combination with land use type data, precipitation and groundwater data. The results demonstrated that the central area is mostly stable, and land deformation in the northeast is obvious and uneven. In addition, land use type and precipitation have little influence on land subsidence. Changes in land subsidence were closely related to changes in groundwater level, and seasonal variations in deformation correlated with precipitation. [ABSTRACT FROM AUTHOR]

The Lightning Mapping Imager (LMI) onboard the Fengyun-4A (FY-4A) satellite is the first independently developed satellite-borne lightning imager in China. It enables continuous lightning detection in China and surrounding areas, regardless of weather conditions. The FY-4A LMI uses a Charge-Coupled Device (CCD) array for lightning detection, and the accuracy of lightning positioning is influenced by cloud top height (CTH). In this study, we proposed an ellipsoid CTH parallax correction (ECPC) model for lightning positioning applicable to FY-4A LMI. The model utilizes CTH data from the Advanced Geosynchronous Radiation Imager (AGRI) on FY-4A to correct the lightning positioning data. According to the model, when the CTH is 12 km, the maximum deviation in lightning positioning caused by CTH in Beijing is approximately 0.1177° in the east–west direction and 0.0530° in the north–south direction, corresponding to a horizontal deviation of 13.1558 km, which exceeds the size of a single ground detection unit of the geostationary satellite lightning imager. Therefore, it is necessary to be corrected. A comparison with data from the Beijing Broadband Lightning Network (BLNET) and radar data shows that the corrected LMI data exhibit spatial distribution that is closer to the simultaneous BLNET lightning positioning data. The coordinate differences between the two datasets are significantly reduced, indicating higher consistency with radar data. The correction algorithm decreases the LMI lightning location deviation caused by CTH, thereby improving the accuracy and reliability of satellite lightning positioning data. The proposed ECPC model can be used for the real-time correction of lightning data when CTH is obtained at the same time, and it can be also used for the post-correction of space-based lightning detection with other cloud top height data. [ABSTRACT FROM AUTHOR]

Using ground-based remote sensing equipment (wind profile radar and microwave radiometer) data and ground automatic station data, a ground-based remote sensing sounding profile system (FAS) is constructed, which aims to make use of its advantages of high resolution, high accuracy, and low cost to make up for the lack of space–time density in existing conventional sounding layer information. The retrieval results of remote sensing sounding profiles in Beijing from May 2021 to September 2022 were tested and evaluated. The results show that the correlation coefficient between FAS and conventional sounding specific humidity is 0.89, the root–mean–square deviation is 1.53 g/kg, and the evolution trends of different data sources of convective available potential energy (CAPE) and vertical wind shear are synchronous. A case study was conducted to evaluate the effectiveness of 40 severe convective processes in the Beijing Plain area. The results show that, due to the minute-level time resolution of FAS, the retrieved convective parameters could track the evolution trend of the atmospheric state with high timeliness, dynamically describe the configuration of thermodynamic parameters, and indicate the time-varying local convective potential and instability level. Therefore, it has certain short-term forecasting significance for the occurrence time, intensity, and convection type. [ABSTRACT FROM AUTHOR]

Rapid simulation of land subsidence can provide an effective means of facilitating its management and control. This paper innovatively introduced a back-propagation (BP) neural network subsidence simulation model which was optimized by the sparrow search algorithm (SSA), hereinafter referred to as the SSA-BP model, to simulate land subsidence in Tongzhou District, Beijing. First, based on the time series interferometric synthetic aperture radar (InSAR) monitoring, different technologies such as spatial analysis, Google Earth Engine (GEE), and machine learning were utilized together with a variety of multi-source spatial data, including groundwater level, compressible layer thickness, Visible Infrared Imager Radiometer Suite (VIIRS) nighttime light images, and the OpenStreetMap (OSM) road distribution. Furthermore, we determined the optimal grid scale for land subsidence research and worked out a multifactor-driven SSA-BP land subsidence simulation model for which sensitivity analysis was performed with great care. Main conclusions: (1) From November 2010 to January 2020, the average annual surface displacement rate in Beijing's subcentre, Tongzhou, ranged from −133.9 to +3.9 mm/year. (2) Our SSA-BP land subsidence simulation model allows for a relatively high overall accuracy. The average root mean square error (RMSE) is 5.00 mm/year, the average mean absolute error (MAE) is 3.80 mm/year, and the average coefficient of determination (R2) is 0.83. (3) Sensitivity analysis shows that the SSA-BP model driven by multiple factors has a far better simulation performance than the model driven by any single weighting factor, and displays the highest sensitivity to the groundwater level factor among all the weighting factors. In terms of subdividing the study area, our SSA-BP land subsidence model runs with multifunctional zones displayed a higher accuracy than that without. This paper made use of a machine learning model to simulate land subsidence in Beijing's Tongzhou District and explored the sensitivity of different factors to land subsidence, which is helpful for its scientific prevention. [ABSTRACT FROM AUTHOR]

Excessive exploitation of the groundwater has resulted in obvious three-dimensional (3D) deformation features on the surface of the Beijing Plain. This paper, by combining Interferometric Synthetic Aperture Radar (InSAR) and Global Positioning System (GPS) technologies, has obtained time-series information of the 3D surface deformation in the Beijing Plain, analyzing its spatial distribution characteristics. On this basis, the relationship between different controlling factors with the 3D deformation of the surface has been analyzed as well. The following results are obtained: (1) From 2013 to 2018, the land subsidence, which generally showed the trend of slowing down, was mainly concentrated in the eastern, northern, and southern regions of Beijing Plain, with multiple subsidence centers. (2) Under the International Terrestrial Reference Frame 2005 (ITRF2005), the horizontal direction of all GPS points in the plain is basically the same, with the dominant movement direction being NE112.5°~NE113.8°. Under the Eurasian reference frame, the horizontal movement rate of GPS points significantly decreases. The movement rate and direction of each point are not characteristic of overall trend activity. (3) The distribution and extent of the 3D surface deformation in the Beijing Plain are controlled by the basement structure. Part of the subsided area corresponds to a Quaternary depression formed at the junction of active faults disrupting the area. Similarly, the distribution of horizontal deformation in the E-W and N-S directions of the plain is controlled by the regional basement structure comprising major faults bounding horizontal deformation. (4) Groundwater exploitation is the main cause of the 3D surface deformation in the Beijing Plain. The groundwater funnels of the second and third confined aquifer are in suitable agreement with the land subsidence. The horizontal movement in the Beijing Plain is either directed toward the center of the groundwater or the land subsidence funnel, and the deformation is directed from areas with higher to areas with lower groundwater levels. [ABSTRACT FROM AUTHOR]

The seasonal variations of raindrop size distribution (DSD) and rainfall are investigated using three-year (2016–2018) observations from a two-dimensional video disdrometer (2DVD) located at a suburban station (40.13° N, 116.62° E, ~30 m AMSL) in Beijing, China. The annual distribution of rainfall presents a unimodal distribution with a peak in summer with total rainfall of 966.6 mm, followed by fall. Rain rate (R), mass-weighted mean diameter (Dm), and raindrop concentration (Nt) are stratified into six regimes to study their seasonal variation and relative rainfall contribution to the total seasonal rainfall. Heavy drizzle/light rain (R2: 0.2~2.5 mm h−1) has the maximum occurrence frequency throughout the year, while the total rainfall in summer is primarily from heavy rain (R4: 10~50 mm h−1). The rainfall for all seasons is contributed primarily from small raindrops (Dm2: 1.0~2.0 mm). The distribution of occurrence frequency of Nt and the relative rainfall contribution exhibit similar behavior during four seasons with Nt of 10~1000 m−3 registering the maximum occurrence and rainfall contributions. Rainfall in Beijing is dominated by stratiform rain (SR) throughout the year. There is no convective rainfall (CR) in winter, i.e., it occurs most often during summer. DSD of SR has minor seasonal differences, but varies significantly in CR. The mean values of log10Nw (Nw: mm−1m−3, the generalized intercept parameter) and Dm of CR indicate that the CR during spring and fall in Beijing is neither continental nor maritime, at the same time, the CR in summer is close to the maritime-like cluster. The radar reflectivity (Z) and rain rate () relationship (Z = R) showed seasonal differences, but were close to the standard NEXRAD Z-R relationship in summer. The shape of raindrops observed from 2DVD was more spherical than the shape obtained from previous experiments, and the effect of different axis ratio relations on polarimetric radar measurements was investigated through T-matrix-based scattering simulations. [ABSTRACT FROM AUTHOR]

This study investigates the characteristics of space-borne Lightning Mapping Imager (LMI) lightning products and their relationships with cloud properties using ground-based total lightning observations from the Beijing Broadband Lightning Network (BLNET) and cloud information from S-band Doppler radar data. LMI showed generally consistent lightning spatial distributions with those of BLNET, and yielded a considerable lightning detection capability over regions with complex terrain. The ratios between the LMI events, groups and flashes were approximately 9:3:1, and the number of LMI-detected flashes was roughly one order of magnitude smaller than the number of BLNET-detected flashes. However, in different convective episodes, the LMI detection capability was likely to be affected by cloud properties, especially in strongly electrified convective episodes associated with frequent lightning discharging and thick cloud depth. As a result, LMI tended to detect lightning flashes located in weaker and shallower cloud portions associated with fewer cloud shielding effects. With reference to the BLNET total lightning data as the ground truth of observation (both intra-cloud lightning and cloud-to-ground lightning flashes), the LMI event-based detection efficiency (DE) was estimated to reach 28% under rational spatiotemporal matching criteria (1.5 s and 65 km) over Beijing. In terms of LMI flash-based DE, it was much reduced compared with event-based DE. The LMI flash-based ranged between 1.5% and 3.5% with 1.5 s and 35–65 km matching scales. For 330 ms and 35 km, the spatiotemporal matching criteria used to evaluate Geostationary Lightning Mapper (GLM), the LMI flash-based DE was smaller (<1%). [ABSTRACT FROM AUTHOR]

Based on the resonance fluorescence scattering mechanism, a narrowband sodium (Na) lidar can measure temperature and wind in the mesosphere and lower thermosphere (MLT) region. By using a narrowband spectral filter, background light noise during the day can be suppressed, allowing for continuous observations. To obtain full-diurnal-cycle temperature and wind measurement results, a complex and precise retrieval process is required, along with necessary corrections to minimize measurement errors. This paper introduces the design of a data acquisition unit for three frequencies in three directions of the Na lidar system in the Chinese Meridian Project (Phase II) and investigates the calibration and retrieval methods for obtaining diurnal temperature and horizontal wind in the MLT region, using a Na Doppler lidar with Faraday anomalous dispersion optical filter (FADOF). Furthermore, these methods are applied to observations conducted by a Na lidar in Beijing, China. The wind and temperature results over full diurnal cycles obtained from the all-solid-state Na Doppler lidar are reported for the first time and compared with temperature measurements from satellite, as well as wind observations from a meteor radar. The comparison demonstrates a reasonable agreement between the results, indicating the rationality of the lidar-retrieved results and the feasibility and effectiveness of the data correction and retrieval method. [ABSTRACT FROM AUTHOR]

Groundwater overexploitation is a critical issue in the North China Plain (NCP), resulting in groundwater level decline and surface subsidence for the last half-century. This problem, however, has been greatly alleviated by the South-to-North Water Diversion (SNWD) Project since 2015. Monitoring of this process has been steadily improved in recent years using water level and geodetic observations. Here, we characterize the water storage change at the Huairou groundwater reserve site (HGRS) in Beijing due to the SNWD by combining Interferometric Synthetic Aperture Radar (InSAR) data of the Sentinel-1 satellites, continuous Global Positioning System (GPS) data, and well water level data observed during the same time. InSAR observations revealed subsidence up to ~400 mm in the Beijing plain but uplift at ~40 mm in the HGRS during 2015–2019, and more than 70% of the uplift occurred from October 2018 to January 2019. By integrating the most significant uplift deformation during October 2018 to January 2019 with water level observations at the same time, we estimated the storativity of the confined aquifer system at HGRS as 1. 68 – 7. 82 × 10 − 3 , weighing in the correction for effective stress and surface deformation for various situations. Based on the estimated aquifer storativity and the observed water level change in the unconfined and confined aquifer, the recharged water storage for the confined and unconfined aquifers was estimated as 1. 20 – 1. 39 × 10 7   m 3 and ~ 2. 86 × 10 8   m 3 from 6 October 2018 to 22 January 2019, respectively, which is about 4% and 91% of the surface water recharge through river channels in the same period due to the SNWD Project. Our study demonstrates that integration of geodetic and hydrological data can provide crucial information for the assessment of groundwater circulation and assistance of groundwater management. [ABSTRACT FROM AUTHOR]

Land subsidence is one of the most important geological hazards in Beijing, China, and its scope and magnitude have been growing rapidly over the past few decades, mainly due to long-term groundwater withdrawal. Interferometric Synthetic Aperture Radar (InSAR) has been used to monitor the deformation in Beijing, but there is a lack of analysis of the long-term spatiotemporal evolution of land subsidence. This study focused on detecting and characterizing spatiotemporal changes in subsidence in the Beijing Plain by using Persistent Scatterer Interferometry (PSI) and geographic spatial analysis. Land subsidence during 2003–2015 was monitored by using ENVISAT ASAR (2003–2010), RADARSAT-2 (2011–2015) and TerraSAR-X (2010–2015) images, with results that are consistent with independent leveling measurements. The radar-based deformation velocity ranged from −136.9 to +15.2 mm/year during 2003–2010, and −149.4 to +8.9 mm/year during 2011–2015 relative to the reference point. The main subsidence areas include Chaoyang, Tongzhou, Shunyi and Changping districts, where seven subsidence bowls were observed between 2003 and 2015. Equal Fan Analysis Method (EFAM) shows that the maximum extensive direction was eastward, with a growing speed of 11.30 km2/year. Areas of differential subsidence were mostly located at the boundaries of the seven subsidence bowls, as indicated by the subsidence rate slope. Notably, the area of greatest subsidence was generally consistent with the patterns of groundwater decline in the Beijing Plain. [ABSTRACT FROM AUTHOR]

Beijing is severely affected by land subsidence, and rapid urbanisation and building construction might accelerate the land subsidence process. Based on 39 Envisat Advanced Synthetic Aperture Radar (ASAR) images acquired between 2003–2010, 55 TerraSAR-X images acquired between 2010–2016, and urban building information, we analysed the relationship between land subsidence and buildings at the regional, block, and building scales. The results show that the surface displacement rate in the Beijing urban area ranged from −109 mm/year to +13 mm/year between 2003–2010, and from −151 mm/year to +19 mm/year between 2010–2016; two subsidence bowls were mainly distributed in the eastern part of the Chaoyang District. The displacement rate agreed well with the levelling measurements, with an average bias of less than six mm/year. At the regional scale, the spatial pattern of land subsidence was mainly controlled by groundwater extraction, compressible layer thickness, and geological faults. Subsidence centres were located in the area around ground water funnels with a compressible layer depth of 50–70 m. The block-scale analysis demonstrated a clear correlation between the block construction age and the spatial unevenness of subsidence. The blocks constructed between 1998–2005 and after 2005 showed considerably more subsidence unevenness and temporal instability than the blocks constructed before 1998 during both time periods. The examination of the new blocks showed that the spatial unevenness increased with building volume variability. For the 16 blocks with a high building volume, variability, and subsidence unevenness, the building-scale analysis showed a positive relationship between building volume and settlement in most blocks, although the R2 was lower than 0.5. The results indicate that intense building construction in urban areas could cause differential settlement at the block scale in Beijing, while the settlement of single buildings could be influenced by the integrated effects of building volume, foundation structures, and the hydrogeological background. [ABSTRACT FROM AUTHOR]

This study demonstrates the capability of a cloud model in simulating a real-world tornado using observed radiosonde data that define a homogeneous background. A reasonable simulation of a tornado event in Beijing, China, on 21 July 2012 is obtained. The simulation reveals the evolution of a descending reflectivity core (DRC) that has commonalities with radar observations, which retracts upward right before tornadogenesis. Tornadogenesis can be divided into three steps: the downward development of mesocyclone vortex, the upward development of tornado vortex, and the eventual downward development of condensation funnel cloud. This bottom-up development provides a numerical evidence for the growing support for a bottom-up, rapid tornadogenesis process as revealed by the state-of-the-art mobile X-band phase-array radar observations. The evolution of the simulated tornado features two replacement processes of three near-surface vortices coupled with the same midlevel updraft. The first replacement occurs during the intensification of the tornado before its maturity. The second replacement occurs during the tornado's demise, when the connection between the midlevel mesocyclone and the near-surface vortex is cut off by a strong downdraft. This work shows the potential of idealized tornado simulations and three-dimensional illustrations in investigating the spiral nature and evolution of tornadoes. [ABSTRACT FROM AUTHOR]

Identifying Persistent Scatterers (PSs) is one of the key processing steps of the Persistent Scatterer Interferometric Synthetic Aperture Radar (PS-InSAR) technique. The number, density, and reliability of identified PSs directly affect the monitoring accuracy of land subsidence, especially in higher density urban environments. As a result of the side-looking viewing geometry of SAR, the layover effect poses a major challenge to the PS identification. This research proposes joint modeling of the PS-InSAR technique and RELAX algorithm for SAR tomography (PS-InSAR+RELAX) to detect single and double scatterers and to improve the identification and reliability of PSs. It has been demonstrated that RELAX improves separation of the scatterers when compared to two other spectral analysis methods for SAR tomography, Beam-Forming (BF) and Singular Value Decomposition (SVD). RELAX exhibits the least noise when the number of baseline changes from 15 to 30, and it can separate the scatterers at a lower Normal-Slant-Range (NSR) height than the two other methods. As RELAX can better identify, separate, and then filter out layover scatterers, the number and density of PSs identified by PS-InSAR+RELAX is reduced and visually simplified, suggesting that the method can effectively reduce the influence of the layover effect on the PS identification. Also, the PSs identified by PS-InSAR+RELAX are more coherent than those identified by the traditional PS-InSAR technique. The proposed technique has been applied to Sentinel-1A data acquired from 2014 to 2016, to monitor land subsidence in the city of Beijing, China. When evaluated against the leveling measurements, PS-InSAR+RELAX performs better than the traditional PS-InSAR technique, with the correlation coefficients (r) of r = 0.98 and r = 0.95, respectively. [ABSTRACT FROM AUTHOR]

In this study, the temperature and relative humidity profiles retrieved from five ground-based microwave radiometers in Beijing were assimilated into the rapid-refresh multi-scale analysis and prediction system-short term (RMAPS-ST). The precipitation bifurcation prediction that occurred in Beijing on 4 May 2019 was selected as a case to evaluate the impact of their assimilation. For this purpose, two experiments were set. The Control experiment only assimilated conventional observations and radar data, while the microwave radiometers profilers (MWRPS) experiment assimilated conventional observations, the ground-based microwave radiometer profiles and radar data into the RMAPS-ST model. The results show that in comparison with the Control test, the MWRPS test made reasonable adjustments for the thermal conditions in time, better reproducing the weak heat island phenomenon in the observation prior to the rainfall. Thus, assimilating MWRPS improved the skills of the precipitation forecast in both the distribution and the intensity of rainfall precipitation, capable of predicting the process of belt-shaped radar echo splitting and the precipitation bifurcation in the urban area of Beijing. The assimilation of the ground-based microwave radiometer profiles improved the skills of the quantitative precipitation forecast to a certain extent. Among multiple cycle experiments, the onset of 0600 UTC cycle closest to the beginning of rainfall performed best by assimilating the ground-based microwave radiometer profiles. [ABSTRACT FROM AUTHOR]

Ground-based microwave radiometers (MWRPS) can provide continuous atmospheric temperature and relative humidity profiles for a weather prediction model. We investigated the impact of assimilation of ground-based microwave radiometers based on the rapid-refresh multiscale analysis and prediction system-short term (RMAPS-ST). In this study, five MWRP-retrieved profiles were assimilated for the precipitation enhancement that occurred in Beijing on 21 May 2020. To evaluate the influence of their assimilation, two experiments with and without the MWRPS assimilation were set. Compared to the control experiment, which only assimilated conventional observations and radar data, the MWRPS experiment, which assimilated conventional observations, the ground-based microwave radiometer profiles and the radar data, had a positive impact on the forecasts of the RMAPS-ST. The results show that in comparison with the control test, the MWRPS experiment reproduced the heat island phenomenon in the observation better. The MWRPS assimilation reduced the bias and RMSE of two-meter temperature and two-meter specific humidity forecasting in the 0–12 h of the forecast range. Furthermore, assimilating the MWRPS improved both the distribution and the intensity of the hourly rainfall forecast, as compared with that of the control experiment, with observations that predicted the process of the precipitation enhancement in the urban area of Beijing. [ABSTRACT FROM AUTHOR]

We present the mapping at fine spatial scale of aerosol optical properties using a mobile laboratory equipped with LIDAR (Light Detection And Ranging), sun photometer and in situ instruments for performing on-road measurements. The mobile campaign was conducted from 9 May to 19 May 2017 and had the main objective of mapping the distribution of pollutants in the Beijing and North China Plain (NCP) region. The highest AOD (Aerosol Optical Depth) at 440 nm of 1.34 and 1.9 were recorded during two heavy pollution episodes on 18 May and 19 May 2017, respectively. The lowest Planetary Boundary Layer (PBL) heights (0.5–1.5 km) were recorded during the heavy pollution events, correlating with the highest AOD and southern winds. The transport of desert dust from the Gobi Desert was captured during the mobile measurements, impacting Beijing during 9–13 May 2017. Exploring the NCP outside Beijing provided datasets for regions with scarce ground measurements and allowed the mapping of high aerosol concentrations when passing polluted cities in the NCP (Baoding, Tianjin and Tangshan) and along the Binhai New Area. For the first time, we provide mass concentration profiles from the synergy of LIDAR, sun photometer and in situ measurements. The case study along the Binhai New Area revealed mean extinction coefficients of 0.14 ± 0.10 km−1 at 532 nm and a mass concentration of 80 ± 62 μg/m3 in the PBL (<2 km). The highest extinction (0.56 km−1) and mass concentrations (404 μg/m3) were found in the industrial Binhai New Area. The PM10 and PM2.5 fractions of the total mass concentration profiles were separated using the columnar size distribution, derived from the sun photometer measurements. This study offers unique mobile datasets of the aerosol optical properties in the NCP for future applications, such as satellite validation and air quality studies. [ABSTRACT FROM AUTHOR]

Land subsidence threatens the stable operation of urban rail transit, including subways. Obtaining deformation information during the entire life-cycle of a subway becomes a necessary means to guarantee urban safety. Restricted by sensor life and cost, the single-sensor Multi-temporal Interferometric Synthetic Aperture Radar (MTI) technology has been unable to meet the needs of long-term sequence, high-resolution deformation monitoring, especially of linear objects. The multi-sensor MTI time-series fusion (MMTI-TSF) techniques has been proposed to solve this problem, but rarely mentioned. In this paper, an improved MMTI-TSF is systematically explained and its limitations are discussed. Taking the Beijing Subway Network (BSN) as a case study, through cross-validation and timing verification, we find that the improved MMTI-TSF results have higher accuracy (R2 of 98% and, Root Mean Squared Error (RMSE) of 4mm), and compared with 38 leveling points, the fitting effect of the time series is good. We analyzed the characteristics of deformation along the BSN over a 15-year periods. The results suggest that there is a higher risk of instability in the eastern section of Beijing Subway Line 6 (L6). The land subsidence characteristics along the subway lines are related to its position from the subsidence center, and the edge of the subsidence center presents a segmented feature. [ABSTRACT FROM AUTHOR]

Clustering of multi-source geospatial big data provides opportunities to comprehensively describe urban structures. Most existing studies focus only on the clustering of a single type of geospatial big data, which leads to biased results. Although multi-view subspace clustering methods are advantageous for fusing multi-source geospatial big data, exploiting a robust shared subspace in high-dimensional, non-uniform, and noisy geospatial big data remains a challenge. Therefore, we developed a method with adaptive graphs to constrain multi-view subspace clustering of multi-source geospatial big data (agc2msc). First, for each type of data, high-dimensional and noisy original features were projected into a low-dimensional latent representation using autoencoder networks. Then, adaptive graph constraints were used to fuse the latent representations of multi-source data into a shared subspace representation, which preserved the neighboring relationships of data points. Finally, the shared subspace representation was used to obtain the clustering results by employing a spectral clustering algorithm. Experiments on four benchmark datasets showed that agc2msc outperformed nine state-of-the-art methods. agc2msc was applied to infer urban land use types in Beijing using the taxi GPS trajectory, bus smart card transaction, and points of interest datasets. The clustering results may provide useful calibration and reference for urban planning. [ABSTRACT FROM AUTHOR]

Fourteen-month precipitation measurements from a second-generation PARSIVEL disdrometer deployed in Beijing, northern China, were analyzed to investigate the microphysical structure of raindrop size distribution and its implications on polarimetric radar applications. Rainfall types are classified and analyzed in the domain of median volume diameter D 0 and the normalized intercept parameter N w . The separation line between convective and stratiform rain is almost equivalent to rain rate at 8.6 mm h−1 and radar reflectivity at 36.8 dBZ. Convective rain in Beijing shows distinct seasonal variations in log 10 N w – D 0 domain. X-band dual-polarization variables are simulated using the T-matrix method to derive radar-based quantitative precipitation estimation (QPE) estimators, and rainfall products at hourly scale are evaluated for four radar QPE estimators using collocated but independent rain gauge observations. This study also combines the advantages of individual estimators based on the thresholds on polarimetric variables. Results show that the blended QPE estimator has better performance than others. The rainfall microphysical analysis presented in this study is expected to facilitate the development of a high-resolution X-band radar network for urban QPE applications. [ABSTRACT FROM AUTHOR]

In this study, Terra-MODIS (Moderate Resolution Imaging Spectroradiometer) Collections 6 and 6.1 (C6 & C6.1) aerosol optical depth (AOD) retrievals with the recommended high-quality flag (QF = 3) were retrieved from Dark-Target (DT), Deep-Blue (DB) and merged DT and DB (DTB) level–2 AOD products for verification against Aerosol Robotic Network (AERONET) Version 3 Level 2.0 AOD data obtained from 2004–2014 for three sites located in the Beijing-Tianjin-Hebei (BTH) region. These are: Beijing, located over mixed bright urban surfaces, XiangHe located over suburban surfaces, and Xinglong located over hilly and vegetated surfaces. The AOD retrievals were also validated over different land-cover types defined by static monthly NDVI (Normalized Difference Vegetation Index) values obtained from the Terra-MODIS level-3 product (MOD13A3). These include non-vegetated surfaces (NVS, NDVI < 0.2), partially vegetated surfaces (PVS, 0.2 ≤ NDVI ≤ 0.3), moderately vegetated surfaces (MVS, 0.3 < NDVI < 0.5) and densely vegetated surfaces (DVS, NDVI ≥ 0.5). Results show that the DT, DB, and DTB-collocated retrievals achieve a high correlation coefficient of ~ 0.90–0.97, 0.89–0.95, and 0.86–0.95, respectively, with AERONET AOD. The DT C6 and C6.1 collocated retrievals were comparable at XiangHe and Xinglong, whereas at Beijing, the percentage of collocated retrievals within the expected error (↔EE) increased from 21.4% to 35.5%, the root mean square error (RMSE) decreased from 0.37 to 0.24, and the relative percent mean error (RPME) decreased from 49% to 27%. These results suggest significant relative improvement in the DT C6.1 product. The percentage of DB-collocated AOD retrievals ↔EE was greater than 70% at Beijing and Xinglong, whereas less than 66% was observed at XiangHe. Similar to DT AOD, DTB AOD retrievals performed well at XiangHe and Xinglong compared with Beijing. Regionally, DB C6 and C6.1-collocated retrievals performed better than DT and DTB in terms of good quality retrievals and relatively small errors. For diverse vegetated surfaces, DT-collocated retrievals reported small errors and good quality retrievals only for NVS and DVS, whereas larger errors were reported for PVS. MVS. DB contains good quality AOD retrievals over PVS, MVS, and DVS compared with NVS. DTB C6.1 collocated retrievals were better than C6 over NVS, PVS, and DVS. C6.1 is substantially improved overall, compared with C6 at local and regional scales, and over diverse vegetated surfaces. [ABSTRACT FROM AUTHOR]

In metropolitan areas, the static load of high-rise buildings may result in uneven settlement, which seriously threatens residents' living safety. Studying the response relationship between the additional stress of high-rise buildings and foundation settlement plays an important role in ensuring the safe development of metropolitan cities. Firstly, based on Persistent Scatterers Interferometric Aperture Radar (PS-InSAR) technology, we used 68 descending TerraSAR-X images to obtain the surface settlement in the study area from April 2010 to October 2018, which were validated with leveling benchmark monitoring results. Secondly, we calculated the additional stress of the building loads to quantify its effect on the uneven settlement in the Central Business District (CBD) of Beijing. Finally, two sets of characteristic points were selected to analyze the response relationships between foundation settlement and additional stress generated by building loads. The findings show: (1) The surface settlement rate varied from −145.2 to 24 mm/year in the Beijing Plain. The InSAR results agree well with the monitoring results derived from the leveling benchmark; the Pearson correlation coefficients were 0.98 and 0.95 in 2011–2013 and 2015–2016, respectively. (2) The stress results show that the depth of the influence of the static load of high-rise buildings was 74.9 m underground in the CBD. (3) The spatial distribution pattern of the additional stress is consistent with the foundation settlement. A characteristic point with greater additional stress in the same group has a higher foundation settlement rate. This relationship has also been found between the uneven foundation settlement and additional stress gradients. These findings provide scientific support for mitigating economic losses due to foundation settlement caused by additional stresses derived from building loads. [ABSTRACT FROM AUTHOR]

In this study, based on the National Centers for Environmental Prediction (NCEP) Final Analysis (FNL) data, the reliability and performances of their application on clean days and polluted days (based on the PM2.5 mass concentrations) in Beijing were assessed. Conventional meteorological factors and diagnostic physical quantities from the NCEP/FNL data were compared with the L-band radar observations in Beijing in the autumns and winters of 2017–2019. The results indicate that the prediction reliability of the temperature was the best compared with those of the relative humidity and wind speed. It is worth noting that the relative humidity was lower and the near-surface wind speed was higher on polluted days from the NCEP/FNL data than from the observations. As far as diagnostic physical quantity is concerned, it was revealed that the temperature inversion intensity depicted by the NCEP/FNL data was significantly lower than that from the observations, especially on polluted days. For example, the difference in the temperature inversion intensity between the NCEP/FNL data and the observation ranged from −0.56 to −0.77 °C on polluted days. In addition, the difference in the wind shears between the NCEP/FNL reanalysis data and the observations increased to 0.40 m/s in the lower boundary layer on polluted days compared with that on clean days. Therefore, it is suggested that the underestimation of the relative humidity and temperature inversion intensity, and the overestimation of the near-surface wind speed should be seriously considered in simulating the air quality in the model, particularly on polluted days, which should be focused on more in future model developments. [ABSTRACT FROM AUTHOR]

NOX (NOX = NO + NO2) emissions measurements in Beijing are of great significance because they can aid in understanding how NOX pollution develops in mega-cities throughout China. However, NOX emissions in mega-cities are difficult to measure due to changes in wind patterns and moving sources on roads during measurement. To obtain good spatial coverage on different ring roads in Beijing over a short amount of time, two mobile differential optical absorption spectroscopy (DOAS) instruments were used to measure NOX emission flux from April 18th to 26th, 2018. In addition, a wind profile radar provided simultaneous wind field measurements for altitudes between 50 m and 1 km for each ring road measurement. We first determined NOX emission flux of different ring roads using wind field averages from measured wind data. The results showed that the NOX emission flux of Beijing's fifth ring road, which represented the urban part, varied from (19.29 ± 5.26) × 1024 molec./s to (36.46 ± 12.86) × 1024 molec./s. On April 20th, NOX emission flux for the third ring was slightly higher than the fourth ring because the two ring roads were measured at different time periods. We then analyzed the NOX emission flux error budget and error sensitivity. The main error source was the wind field uncertainty. For some measurements, the main emission flux error source was either wind speed uncertainty or wind direction uncertainty, but not both. As Beijing's NOX emissions came from road vehicle exhaust, we found that emission flux error had a more diverse sensitivity to wind direction uncertainty, which improved our knowledge on this topic. The NOX emission flux error sensitivity study indicated that more accurate measurements of the wind field are crucial for effective NOX emission flux measurements in Chinese mega-cities. Obtaining actual time and high resolved wind measurements is an advantage for mega-cities' NOX emission flux measurements. The emission flux errors caused by wind direction and wind speed uncertainties were clearly distinguished. Other sensitivity studies indicated that NOX/NO2 ratio uncertainty dominated flux errors when the NOX/NO2 ratio uncertainty was >0.4. Using two mobile-DOAS and wind profile radars to measure NOx emission flux improved the quality of the emission flux measuring results. This approach could be applied to many other mega-cities in China and in others countries. [ABSTRACT FROM AUTHOR]