Search

Ozone plays an important role in the Earth's atmosphere. It is mainly formed in the tropical stratosphere and is transported by the Brewer–Dobson Circulation to higher latitudes. In the stratosphere, ozone can filter the incoming solar ultraviolet radiation, thus protecting life at the surface. Although tropospheric ozone accounts for only ~10%, it is a powerful GHG and pollutant, harmful to the health of the environment and living beings. Several studies have highlighted biomass burning as a major contributor to the tropospheric ozone budget. Our study focuses on the Natal site (5.40°S, 35.40°W, Brazil), one of the oldest ozone-observing stations in Brazil, which is expected to be influenced by fire plumes in Africa and Brazil. Many studies that examined ozone trends used the total atmospheric columns of ozone, but it is important to assess ozone separately in the troposphere and the stratosphere. In this study, we have used radiosonde ozone profiles and daily TCO measurements to evaluate the variability and changes of both tropospheric and stratospheric ozone separately. The dataset in this study comprises daily total columns of colocalized ozone and weekly ozone profiles collected between 1998 and 2019. The tropospheric columns were estimated by integrating ozone profiles measured by ozone sondes up to the tropopause height. The amount of ozone in the stratosphere was then deduced by subtracting the tropospheric ozone amount from the total amount of ozone measured by the Dobson spectrometer. It was assumed that the amount of ozone in the mesosphere is negligible. This produced three distinct time series of ozone: tropospheric and stratospheric columns as well as total columns. The present study aims to apply a new decomposition method named Empirical Adaptive Wavelet Decomposition (EAWD) that is used to identify the different modes of variability present in the analyzed signal. This is achieved by summing up the most significant Intrinsic Mode Functions (IMF). The Fourier spectrum of the original signal is broken down into spectral bands that frame each IMF obtained by the Empirical Modal Decomposition (EMD). Then, the Empirical Wavelet Transform (EWT) is applied to each interval. Unlike other methods like EMD and multi-linear regression (MLR), the EAWD technique has an advantage in providing better frequency resolution and thus overcoming the phenomenon of mode-mixing, as well as detecting possible breakpoints in the trend mode. The obtained ozone datasets were analyzed using three methods: MLR, EMD, and EAWD. The EAWD algorithm exhibited the advantage of retrieving ~90% to 95% of ozone variability and detecting possible breakpoints in its trend component. Overall, the MRL and EAWD methods showed almost similar trends, a decrease in the stratosphere ozone (−1.3 ± 0.8%) and an increase in the tropospheric ozone (+4.9 ± 1.3%). This study shows the relevance of combining data to separately analyze tropospheric and stratospheric ozone variability and trends. It highlights the advantage of the EAWD algorithm in detecting modes of variability in a geophysical signal without prior knowledge of the underlying forcings. [ABSTRACT FROM AUTHOR]

High-speed railways in the Beijing–Tianjin–Hebei (BTH) Plain are gradually becoming more widespread, covering a greater area. The operational safety of high-speed railways is influenced by the continuous development of land subsidence. It is necessary to predict the subsidence along the high-speed railways; thus, this work is of critical importance to the safety of high-speed railway operation. In this study, we processed Sentinel-1A data using the Persistent Scatterer Interferometric Synthetic Aperture Radar (PS-InSAR) technique to acquire the land subsidence in the typical BTH area. Then, we combined the Empirical Mode Decomposition (EMD) and Gradient Boosting Decision Tree (GBDT) methods (EMD-GBDT) to forecast land subsidence along high-speed railways. The results revealed that some parts of the high-speed railways in the BTH plain had passed through or approached the land subsidence area; the maximum cumulative subsidence of the Beijing–Shanghai, Tianjin–Baoding and Shijiazhuang–Jinan high-speed railways reached 326 mm, 384 mm and 350 mm, respectively. The forecasting accuracy for land subsidence along high-speed railways was enhanced by the EMD-GBDT model. The Root Mean Square Error (RMSE) and Mean Absolute Error (MAE) were 0.38 mm to 0.56 mm and 0.23 mm to 0.38 mm, respectively. [ABSTRACT FROM AUTHOR]

Ozone plays an important role in the Earth’s atmosphere. It is mainly formed in the tropical stratosphere and is transported by the Brewer–Dobson Circulation to higher latitudes. In the stratosphere, ozone can filter the incoming solar ultraviolet radiation, thus protecting life at the surface. Although tropospheric ozone accounts for only ~10%, it is a powerful GHG and pollutant, harmful to the health of the environment and living beings. Several studies have highlighted biomass burning as a major contributor to the tropospheric ozone budget. Our study focuses on the Natal site (5.40°S, 35.40°W, Brazil), one of the oldest ozone-observing stations in Brazil, which is expected to be influenced by fire plumes in Africa and Brazil. Many studies that examined ozone trends used the total atmospheric columns of ozone, but it is important to assess ozone separately in the troposphere and the stratosphere. In this study, we have used radiosonde ozone profiles and daily TCO measurements to evaluate the variability and changes of both tropospheric and stratospheric ozone separately. The dataset in this study comprises daily total columns of colocalized ozone and weekly ozone profiles collected between 1998 and 2019. The tropospheric columns were estimated by integrating ozone profiles measured by ozone sondes up to the tropopause height. The amount of ozone in the stratosphere was then deduced by subtracting the tropospheric ozone amount from the total amount of ozone measured by the Dobson spectrometer. It was assumed that the amount of ozone in the mesosphere is negligible. This produced three distinct time series of ozone: tropospheric and stratospheric columns as well as total columns. The present study aims to apply a new decomposition method named Empirical Adaptive Wavelet Decomposition (EAWD) that is used to identify the different modes of variability present in the analyzed signal. This is achieved by summing up the most significant Intrinsic Mode Functions (IMF). The Fourier spectrum of the original signal is broken down into spectral bands that frame each IMF obtained by the Empirical Modal Decomposition (EMD). Then, the Empirical Wavelet Transform (EWT) is applied to each interval. Unlike other methods like EMD and multi-linear regression (MLR), the EAWD technique has an advantage in providing better frequency resolution and thus overcoming the phenomenon of mode-mixing, as well as detecting possible breakpoints in the trend mode. The obtained ozone datasets were analyzed using three methods: MLR, EMD, and EAWD. The EAWD algorithm exhibited the advantage of retrieving ~90% to 95% of ozone variability and detecting possible breakpoints in its trend component. Overall, the MRL and EAWD methods showed almost similar trends, a decrease in the stratosphere ozone (−1.3 ± 0.8%) and an increase in the tropospheric ozone (+4.9 ± 1.3%). This study shows the relevance of combining data to separately analyze tropospheric and stratospheric ozone variability and trends. It highlights the advantage of the EAWD algorithm in detecting modes of variability in a geophysical signal without prior knowledge of the underlying forcings.

High-speed railways in the Beijing–Tianjin–Hebei (BTH) Plain are gradually becoming more widespread, covering a greater area. The operational safety of high-speed railways is influenced by the continuous development of land subsidence. It is necessary to predict the subsidence along the high-speed railways; thus, this work is of critical importance to the safety of high-speed railway operation. In this study, we processed Sentinel-1A data using the Persistent Scatterer Interferometric Synthetic Aperture Radar (PS-InSAR) technique to acquire the land subsidence in the typical BTH area. Then, we combined the Empirical Mode Decomposition (EMD) and Gradient Boosting Decision Tree (GBDT) methods (EMD-GBDT) to forecast land subsidence along high-speed railways. The results revealed that some parts of the high-speed railways in the BTH plain had passed through or approached the land subsidence area; the maximum cumulative subsidence of the Beijing–Shanghai, Tianjin–Baoding and Shijiazhuang–Jinan high-speed railways reached 326 mm, 384 mm and 350 mm, respectively. The forecasting accuracy for land subsidence along high-speed railways was enhanced by the EMD-GBDT model. The Root Mean Square Error (RMSE) and Mean Absolute Error (MAE) were 0.38 mm to 0.56 mm and 0.23 mm to 0.38 mm, respectively.

Micro-Doppler is regarded as a unique signature of a target with micro-motions. The sophisticated recognition of the cone-shaped targets can be realized through the micro-Doppler effect. However, it is difficult to extract the micro-motion features perpendicular to the radar line of sight (LOS) effectively. In this paper, a micro-Doppler parameters extraction method of the cone-shaped targets is put forward based on the rotating antenna. First, a new radar configuration is proposed, in which an antenna rotates uniformly on a fixed circle, thus producing Doppler frequency shift. Second, the expression of the micro-Doppler frequency shift induced by the precession cone-shaped target is derived. Then, the micro-Doppler curves of point scatterers at the cone top and bottom are separated by the smoothness of the curves, and the empirical mode decomposition (EMD) method is utilized for the detection and estimation of the coning frequency. Finally, the micro-motion components perpendicular to the radar LOS are inverted by the peak of micro-Doppler frequency curve. Simulation results prove the effectiveness and robustness of the proposed method. [ABSTRACT FROM AUTHOR]

The development of technology has strongly affected regional urbanization. With development of mobile communication technology, intelligent devices have become increasingly widely used in people's lives. The application of big data in urban computing is multidimensional; it has been involved in different fields, such as urban planning, network optimization, intelligent transportation, energy consumption and so on. Data analysis becomes particularly important for wireless networks. In this paper, a method for analyzing cellular traffic data was proposed. Firstly, a method to extract trend components, periodic components and essential components from complex traffic time series was proposed. Secondly, we introduced causality data mining. Different from traditional time series causality analysis, the depth of causal mining was increased. We conducted causality verification on different components of time series and the results showed that the causal relationship between base stations is different in trend component, periodic component and essential component in urban wireless network. This is crucial for urban planning and network management. Thirdly, DIC-ST: a spatial temporal time series prediction based on decomposition and integration system with causal structure learning was proposed by combining GCN. Final results showed that the proposed method significantly improves the accuracy of cellular traffic prediction. At the same time, this method can play a crucial role for urban computing in network management, intelligent transportation, base station siting and energy consumption when combined with remote sensing map information. [ABSTRACT FROM AUTHOR]

Micro-Doppler is regarded as a unique signature of a target with micro-motions. The sophisticated recognition of the cone-shaped targets can be realized through the micro-Doppler effect. However, it is difficult to extract the micro-motion features perpendicular to the radar line of sight (LOS) effectively. In this paper, a micro-Doppler parameters extraction method of the cone-shaped targets is put forward based on the rotating antenna. First, a new radar configuration is proposed, in which an antenna rotates uniformly on a fixed circle, thus producing Doppler frequency shift. Second, the expression of the micro-Doppler frequency shift induced by the precession cone-shaped target is derived. Then, the micro-Doppler curves of point scatterers at the cone top and bottom are separated by the smoothness of the curves, and the empirical mode decomposition (EMD) method is utilized for the detection and estimation of the coning frequency. Finally, the micro-motion components perpendicular to the radar LOS are inverted by the peak of micro-Doppler frequency curve. Simulation results prove the effectiveness and robustness of the proposed method.

The development of technology has strongly affected regional urbanization. With development of mobile communication technology, intelligent devices have become increasingly widely used in people’s lives. The application of big data in urban computing is multidimensional; it has been involved in different fields, such as urban planning, network optimization, intelligent transportation, energy consumption and so on. Data analysis becomes particularly important for wireless networks. In this paper, a method for analyzing cellular traffic data was proposed. Firstly, a method to extract trend components, periodic components and essential components from complex traffic time series was proposed. Secondly, we introduced causality data mining. Different from traditional time series causality analysis, the depth of causal mining was increased. We conducted causality verification on different components of time series and the results showed that the causal relationship between base stations is different in trend component, periodic component and essential component in urban wireless network. This is crucial for urban planning and network management. Thirdly, DIC-ST: a spatial temporal time series prediction based on decomposition and integration system with causal structure learning was proposed by combining GCN. Final results showed that the proposed method significantly improves the accuracy of cellular traffic prediction. At the same time, this method can play a crucial role for urban computing in network management, intelligent transportation, base station siting and energy consumption when combined with remote sensing map information.

To improve the accuracy of pulse matching and the mapping quality of lightning discharges, the Pearson correlation method combined with empirical mode decomposition (EMD) is introduced for discharge electric field pulse matching. This paper uses the new method to locate the lightning channels of an intra-cloud (IC) lightning flash and a cloud-to-ground (CG) lightning flash and analyzes the location results for the two lightning flashes. The results show that this method has a good performance in lightning location. Compared with the pulse-peak feature matching method, the positioning results of the new method are significantly improved, which is mainly due to the much larger number of positioning points (matched pulses). The number of located radiation sources has increased by nearly a factor of seven, which can significantly improve the continuity of the lightning channel and clearly distinguish the developmental characteristics. In the CG flash, there were three negative recoil streamers in the positive leader channel. After the three negative recoil streamers were finished, taking approximately 1 ms, 12 ms, and 2 ms, respectively, the negative leader channel underwent a K-process. The three negative recoil streamers are not connected to the K-processes in the negative leader channel. We think that the three negative recoil streamers may have triggered the three K-processes, respectively.

The global navigation satellite system (GNSS) has seen tremendous advances in measurement precision and accuracy, and it allows researchers to perform geodynamics and geophysics studies through the analysis of GNSS time series. Moreover, GNSS time series not only contain geophysical signals, but also unmodeled errors and other nuisance parameters, which affect the performance in the estimation of site coordinates and related parameters. As the number of globally distributed GNSS reference stations increases, GNSS time series analysis software should be developed with more flexible format support, better human–machine interaction, and with powerful noise reduction analysis. To meet this requirement, a new software named GNSS time series noise reduction software (GNSS-TS-NRS) was written in MATLAB and was developed. GNSS-TS-NRS allows users to perform noise reduction analysis and spatial filtering on common mode errors and to visualize GNSS position time series. The functions’ related theoretical background of GNSS-TS-NRS were introduced. Firstly, we showed the theoretical background algorithms of the noise reduction analysis (empirical mode decomposition (EMD), ensemble empirical mode decomposition (EEMD)). We also developed three improved algorithms based on EMD for noise reduction, and the results of the test example showed our proposed methods with better effect. Secondly, the spatial filtering model supported five algorithms on a separate common model error: The stacking filter method, weighted stacking filter method, correlation weighted superposition filtering method, distance weighted filtering method, and principal component analysis, as well as with batch processing. Finally, the developed software also enabled other functions, including outlier detection, correlation coefficient calculation, spectrum analysis, and distribution estimation. The main goal of the manuscript is to share the software with the scientific community to introduce new users to the GNSS time series noise reduction and application.

Rice crop monitoring is an important activity for crop management. This study aimed to develop a phenology-based classification approach for the assessment of rice cropping systems in Mekong Delta, Vietnam, using Moderate Resolution Imaging Spectroradiometer (MODIS) data. The data were processed from December 2000, to December 2012, using empirical mode decomposition (EMD) in three main steps: (1) data pre-processing to construct the smooth MODIS enhanced vegetation index (EVI) time-series data; (2) rice crop classification; and (3) accuracy assessment. The comparisons between the classification maps and the ground reference data indicated overall accuracies and Kappa coefficients, respectively, of 81.4% and 0.75 for 2002, 80.6% and 0.74 for 2006 and 85.5% and 0.81 for 2012. The results by comparisons between MODIS-derived rice area and rice area statistics were slightly overestimated, with a relative error in area (REA) from 0.9–15.9%. There was, however, a close correlation between the two datasets (R2 ≥ 0.89). From 2001 to 2012, the areas of triple-cropped rice increased approximately 31.6%, while those of the single-cropped rain-fed rice, double-cropped irrigated rice and double-cropped rain-fed rice decreased roughly −5.0%, −19.2% and −7.4%, respectively. This study demonstrates the validity of such an approach for rice-crop monitoring with MODIS data and could be transferable to other regions.

The Ice, Cloud and Land Elevation Satellite-2 (ICESat-2) mission employs a micro-pulse photon-counting LiDAR system for mapping and monitoring the biomass and carbon of terrestrial ecosystems over large areas. In preparation for ICESat-2 data processing and applications, this paper aimed to develop and validate an effective algorithm for better estimating ground elevation and vegetation height from photon-counting LiDAR data. Our new proposed algorithm consists of three key steps. Firstly, the noise photons were filtered out using a noise removal algorithm based on localized statistical analysis. Secondly, we classified the signal photons into canopy photons and ground photons by conducting a series of operations, including elevation frequency histogram building, empirical mode decomposition (EMD), and progressive densification. At the same time, we also identified the top of canopy (TOC) photons from canopy photons by percentile statistics method. Thereafter, the ground and TOC surfaces were generated from ground photons and TOC photons by cubic spline interpolation, respectively. Finally, the ground elevation and vegetation height were estimated by retrieved ground and TOC surfaces. The results indicate that the noise removal algorithm is effective in identifying background noise and preserving signal photons. The retrieved ground elevation is more accurate than the retrieved vegetation height, and the results of nighttime data are better than those of the corresponding daytime data. Specifically, the root-mean-square error (RMSE) values of ground elevation estimates range from 2.25 to 6.45 m for daytime data and 2.03 to 6.03 m for nighttime data. The RMSE values of vegetation height estimates range from 4.63 to 8.92 m for daytime data and 4.55 to 8.65 m for nighttime data. Our algorithm performs better than the previous algorithms in estimating ground elevation and vegetation height due to lower RMSE values. Additionally, the results also illuminate that the photon classification algorithm effectively reduces the negative effects of slope and vegetation coverage. Overall, our paper provides an effective solution for estimating ground elevation and vegetation height from micro-pulse photon-counting LiDAR data.

Rice crop monitoring is an important activity for crop management. This study aimed to develop a phenology-based classification approach for the assessment of rice cropping systems in Mekong Delta, Vietnam, using Moderate Resolution Imaging Spectroradiometer (MODIS) data. The data were processed from December 2000, to December 2012, using empirical mode decomposition (EMD) in three main steps: (1) data pre-processing to construct the smooth MODIS enhanced vegetation index (EVI) time-series data; (2) rice crop classification; and (3) accuracy assessment. The comparisons between the classification maps and the ground reference data indicated overall accuracies and Kappa coefficients, respectively, of 81.4% and 0.75 for 2002, 80.6% and 0.74 for 2006 and 85.5% and 0.81 for 2012. The results by comparisons between MODIS-derived rice area and rice area statistics were slightly overestimated, with a relative error in area (REA) from 0.9-15.9%. There was, however, a close correlation between the two datasets (R2 ⩾ 0.89). From 2001 to 2012, the areas of triple-cropped rice increased approximately 31.6%, while those of the single-cropped rain-fed rice, double-cropped irrigated rice and double-cropped rain-fed rice decreased roughly -5.0%, -19.2% and -7.4%, respectively. This study demonstrates the validity of such an approach for rice-crop monitoring with MODIS data and could be transferable to other regions. [ABSTRACT FROM AUTHOR]

To improve the accuracy of pulse matching and the mapping quality of lightning discharges, the Pearson correlation method combined with empirical mode decomposition (EMD) is introduced for discharge electric field pulse matching. This paper uses the new method to locate the lightning channels of an intra-cloud (IC) lightning flash and a cloud-to-ground (CG) lightning flash and analyzes the location results for the two lightning flashes. The results show that this method has a good performance in lightning location. Compared with the pulse-peak feature matching method, the positioning results of the new method are significantly improved, which is mainly due to the much larger number of positioning points (matched pulses). The number of located radiation sources has increased by nearly a factor of seven, which can significantly improve the continuity of the lightning channel and clearly distinguish the developmental characteristics. In the CG flash, there were three negative recoil streamers in the positive leader channel. After the three negative recoil streamers were finished, taking approximately 1 ms, 12 ms, and 2 ms, respectively, the negative leader channel underwent a K-process. The three negative recoil streamers are not connected to the K-processes in the negative leader channel. We think that the three negative recoil streamers may have triggered the three K-processes, respectively. [ABSTRACT FROM AUTHOR]

The global navigation satellite system (GNSS) has seen tremendous advances in measurement precision and accuracy, and it allows researchers to perform geodynamics and geophysics studies through the analysis of GNSS time series. Moreover, GNSS time series not only contain geophysical signals, but also unmodeled errors and other nuisance parameters, which affect the performance in the estimation of site coordinates and related parameters. As the number of globally distributed GNSS reference stations increases, GNSS time series analysis software should be developed with more flexible format support, better human–machine interaction, and with powerful noise reduction analysis. To meet this requirement, a new software named GNSS time series noise reduction software (GNSS-TS-NRS) was written in MATLAB and was developed. GNSS-TS-NRS allows users to perform noise reduction analysis and spatial filtering on common mode errors and to visualize GNSS position time series. The functions' related theoretical background of GNSS-TS-NRS were introduced. Firstly, we showed the theoretical background algorithms of the noise reduction analysis (empirical mode decomposition (EMD), ensemble empirical mode decomposition (EEMD)). We also developed three improved algorithms based on EMD for noise reduction, and the results of the test example showed our proposed methods with better effect. Secondly, the spatial filtering model supported five algorithms on a separate common model error: The stacking filter method, weighted stacking filter method, correlation weighted superposition filtering method, distance weighted filtering method, and principal component analysis, as well as with batch processing. Finally, the developed software also enabled other functions, including outlier detection, correlation coefficient calculation, spectrum analysis, and distribution estimation. The main goal of the manuscript is to share the software with the scientific community to introduce new users to the GNSS time series noise reduction and application. [ABSTRACT FROM AUTHOR]

The Ice, Cloud and Land Elevation Satellite-2 (ICESat-2) mission employs a micro-pulse photon-counting LiDAR system for mapping and monitoring the biomass and carbon of terrestrial ecosystems over large areas. In preparation for ICESat-2 data processing and applications, this paper aimed to develop and validate an effective algorithm for better estimating ground elevation and vegetation height from photon-counting LiDAR data. Our new proposed algorithm consists of three key steps. Firstly, the noise photons were filtered out using a noise removal algorithm based on localized statistical analysis. Secondly, we classified the signal photons into canopy photons and ground photons by conducting a series of operations, including elevation frequency histogram building, empirical mode decomposition (EMD), and progressive densification. At the same time, we also identified the top of canopy (TOC) photons from canopy photons by percentile statistics method. Thereafter, the ground and TOC surfaces were generated from ground photons and TOC photons by cubic spline interpolation, respectively. Finally, the ground elevation and vegetation height were estimated by retrieved ground and TOC surfaces. The results indicate that the noise removal algorithm is effective in identifying background noise and preserving signal photons. The retrieved ground elevation is more accurate than the retrieved vegetation height, and the results of nighttime data are better than those of the corresponding daytime data. Specifically, the root-mean-square error (RMSE) values of ground elevation estimates range from 2.25 to 6.45 m for daytime data and 2.03 to 6.03 m for nighttime data. The RMSE values of vegetation height estimates range from 4.63 to 8.92 m for daytime data and 4.55 to 8.65 m for nighttime data. Our algorithm performs better than the previous algorithms in estimating ground elevation and vegetation height due to lower RMSE values. Additionally, the results also illuminate that the photon classification algorithm effectively reduces the negative effects of slope and vegetation coverage. Overall, our paper provides an effective solution for estimating ground elevation and vegetation height from micro-pulse photon-counting LiDAR data. [ABSTRACT FROM AUTHOR]